数学建模社区-数学中国
标题:
在线等回复
[打印本页]
作者:
509
时间:
2013-8-20 15:55
标题:
在线等回复
本帖最后由 袁海亮 于 2013-8-20 16:15 编辑
: y6 l- [7 ~, r( w4 i7 x
- W6 F' J9 t$ I2 ^9 i, A. C! j
[m,n]=size(data(:,:,1));
0 h5 U" C, y- H& s
x=zeros(m*n,3);
3 o" W: F7 Q, e! h0 m) {% I5 @7 U
for i=1:m
) [8 E; r. P5 n, V
for j=1:n
0 x" q. Z; K" o7 _' V
x(i,1)=j;
* `+ H9 M: u( F# I
x(i,2)=i;
7 P8 B& x$ r( g( J; |# ]3 M" E9 s
x(i,3)=data(i,j,1);
! U& K' e* }$ ^! _( Y7 V& b: H
end
8 S u! c: P! e( x
end
) U u( W8 Q0 g0 f4 t5 v# L
X = zscore(xdata);
$ q4 b; P* B5 ?4 D! B; t$ B
- a% J9 n0 a! F( i- T" X
: [5 x' o; k4 v$ E7 b& W
%****************************进行模糊C均值聚类******************************
7 f C7 e% |4 E, p& y" k* A
% 设置幂指数为3，最大迭代次数为200，目标函数的终止容限为1e-6，不显示中间迭代过程
6 w8 F, u: P3 r; h
options = [3, 200, 1e-6, 0];
9 u; X! p8 k, Z. Y) N
% 调用fcm函数进行模糊C均值聚类，返回类中心坐标矩阵center，隶属度矩阵U，目标函数值obj_fcn
5 B. L0 b# h2 q3 c; b
[center,U,obj_fcn] = fcm(X,3,options);
6 S. O4 T6 _$ V
% n7 u7 r! B+ l8 z( K6 V
+ ?- a& ]1 e4 \0 q" x; B8 F) V; H
%****************查看聚类结果*****************************
$ O" Y) s+ E3 Z# _6 q% s
id1 = find(U(1,:) == max(U));
$ d: r% y" ?0 r( d/ B# j
id2 = find(U(2,:) == max(U));
) s# S: Q* H# d3 \+ I" e5 x
id3 = find(U(3,:) == max(U));
2 u/ O" X% K1 ]0 i. x
city(id1)
# m1 K# `. H0 ^# y1 b4 F
city(id2)
8 {5 ^* y. @! A, i0 N; R! @
city(id3)
3 I& j' x, R5 y; s- W# Y
7 _$ k+ T) l4 ?: W( s/ F+ H( E
为什么一查看聚类结果
: K) M" g+ V3 N( {, H
就出现错误呢
. q3 k9 ^4 R* s9 M
请大神指点一下为什么
; J) ]/ Q. Y* b% t' \
总提示指数超过矩阵维数
) K/ M) a1 W; L, y/ o: Q
" r' o, i7 ]# b5 g2 }/ R
, m; R. @/ `0 C ]
作者:
净心、精心
时间:
2013-8-20 16:40
http://www.madio.net/thread-198267-1-1.html
4 L( G, N: \7 y, R5 E' Q
亲，在有问题时，记得标题要写相关的困惑和哪一方面的知识，这样才能能快的得到帮助，要不一般看到你这标题，没人愿意点开看的
欢迎光临 数学建模社区-数学中国 (http://www.madio.net/)
Powered by Discuz! X2.5