模糊C均值聚类

509

[m,n]=size(data(:,:,1));
% }8 \  \+ P. P* g9 t9 Vx=zeros(m*n,3);
for i=1:m
    for j=1:n
       x(i,1)=j;
0 {8 U, k# v) u) _       x(i,2)=i;
       x(i,3)=data(i,j,1);
! X2 d( Z. X2 p% Y    end
end
X = zscore(xdata);
2 A. ^* s3 i) @! S. I( }

%****************************进行模糊C均值聚类******************************
3 j$ O4 x& P" F( _% 设置幂指数为3，最大迭代次数为200，目标函数的终止容限为1e-6，不显示中间迭代过程
) V0 r% j; O  ?. Xoptions = [3, 200, 1e-6, 0];
% 调用fcm函数进行模糊C均值聚类，返回类中心坐标矩阵center，隶属度矩阵U，目标函数值obj_fcn
[center,U,obj_fcn] = fcm(X,3,options);
/ r! h8 h4 V, Q+ ^; s3 z- v4 F: R

%****************查看聚类结果*****************************
  E" N% K: e8 ~( F8 Q  aid1 = find(U(1,:) == max(U));  
/ P( k: O  f8 rid2 = find(U(2,:) == max(U));  
2 f, x& A  O; @id3 = find(U(3,:) == max(U));   
city(id1)  
1 F/ \9 t8 u: k6 vcity(id2)  
city(id3)  
4 H% Y+ O7 \# R3 R( v/ H4 J
为什么一查看聚类结果
% g8 ?% U, W7 |& {' d, u就出现错误呢
请大神指点一下为什么
总提示指数超过矩阵维数

. E% e+ M$ E' x. A, \" u

我爱照镜子

那三个表情是什么啊

净心、精心

我爱照镜子 发表于 2013-8-20 16:54
% m$ Y6 e$ g  X/ {+ p那三个表情是什么啊

* |# r: N. g9 H0 I7 H- c: {：：：：：：：：：