模糊C均值聚类

509

[m,n]=size(data(:,:,1));
x=zeros(m*n,3);
for i=1:m
- }# h  s( P+ j    for j=1:n
8 a* y( v) c6 [       x(i,1)=j;
+ a/ [/ r* Z* h  @5 X, e       x(i,2)=i;
       x(i,3)=data(i,j,1);
7 \3 _. S9 e6 f/ P" a; z    end
end
X = zscore(xdata);


7 z$ D# ^7 C' f9 ~( w( `' T%****************************进行模糊C均值聚类******************************
% 设置幂指数为3，最大迭代次数为200，目标函数的终止容限为1e-6，不显示中间迭代过程
options = [3, 200, 1e-6, 0];
3 J- x' ^+ \3 ?! H8 f) s; M% 调用fcm函数进行模糊C均值聚类，返回类中心坐标矩阵center，隶属度矩阵U，目标函数值obj_fcn
; i( y8 Y" z! G- h[center,U,obj_fcn] = fcm(X,3,options);


%****************查看聚类结果*****************************
id1 = find(U(1,:) == max(U));  
- I' t) Y* ?: K6 jid2 = find(U(2,:) == max(U));  
id3 = find(U(3,:) == max(U));   
$ y" X& [5 i7 L# J" }' y5 icity(id1)  
( K5 _( k. v% c2 ], xcity(id2)  
city(id3)  
- I+ e! R# D' c# A7 a
为什么一查看聚类结果
就出现错误呢
请大神指点一下为什么
总提示指数超过矩阵维数

% n/ a2 B% L4 f2 _

净心、精心

我爱照镜子 发表于 2013-8-20 16:54
那三个表情是什么啊

：：：：：：：：：

我爱照镜子

那三个表情是什么啊