在线等回复

509

[m,n]=size(data(:,:,1));
x=zeros(m*n,3);
/ a* V/ A; m* R1 V) nfor i=1:m
    for j=1:n
+ T0 X: Q: q" A  L% w& U; q       x(i,1)=j;
) Z+ k4 j2 |- T       x(i,2)=i;
       x(i,3)=data(i,j,1);
; o9 g& b  r5 _    end
end
# M* `$ t7 K9 k" x  Y* x& x) DX = zscore(xdata);
0 t: S$ n: \) Q) _# ~* m
1 z; C8 b% Z3 D) u
%****************************进行模糊C均值聚类******************************
$ ]6 ]+ C) j  V( n# ~% 设置幂指数为3，最大迭代次数为200，目标函数的终止容限为1e-6，不显示中间迭代过程
options = [3, 200, 1e-6, 0];
% 调用fcm函数进行模糊C均值聚类，返回类中心坐标矩阵center，隶属度矩阵U，目标函数值obj_fcn
[center,U,obj_fcn] = fcm(X,3,options);
' o4 x0 P( U+ u+ C! C, u- J* h

%****************查看聚类结果*****************************
id1 = find(U(1,:) == max(U));  
id2 = find(U(2,:) == max(U));  
1 t9 V3 s( t8 H( j) C! w8 I* mid3 = find(U(3,:) == max(U));   
city(id1)  
- W/ ^) W$ V. c. {  \city(id2)  
/ e7 W# u% V) e% w, @city(id3)  
6 R# t3 K$ m8 j' J+ H
为什么一查看聚类结果
  L' A  f( K; S4 e2 i7 c就出现错误呢
请大神指点一下为什么
总提示指数超过矩阵维数
2 i4 P7 |/ K7 X
5 Y1 A. a( H! A2 Z

净心、精心

http://www.madio.net/thread-198267-1-1.html
! h. V% X7 z8 s  y. e, B( t! ~亲，在有问题时，记得标题要写相关的困惑和哪一方面的知识，这样才能能快的得到帮助，要不一般看到你这标题，没人愿意点开看的