模糊C均值聚类

509

[m,n]=size(data(:,:,1));
x=zeros(m*n,3);
6 J2 O9 S& |0 R7 n/ [  M7 ^for i=1:m
    for j=1:n
4 @: C% ~8 N% W2 e2 ~/ i# u       x(i,1)=j;
       x(i,2)=i;
5 p  r: v/ j* l- h       x(i,3)=data(i,j,1);
8 n; Y- [% j+ A3 q6 z    end
end
3 y3 i1 P# {+ r# G! {4 k5 n1 rX = zscore(xdata);
( ~1 V- i# Y+ p" l; o

%****************************进行模糊C均值聚类******************************
% 设置幂指数为3，最大迭代次数为200，目标函数的终止容限为1e-6，不显示中间迭代过程
options = [3, 200, 1e-6, 0];
. @3 p5 u( d: V& P* Q- w; ]$ j% 调用fcm函数进行模糊C均值聚类，返回类中心坐标矩阵center，隶属度矩阵U，目标函数值obj_fcn
[center,U,obj_fcn] = fcm(X,3,options);
/ e0 G1 T) H  ~0 ~
# V$ _8 Q0 l: G) M: h$ ]$ S
%****************查看聚类结果*****************************
id1 = find(U(1,:) == max(U));  
id2 = find(U(2,:) == max(U));  
6 m$ g% R3 R, H! F/ sid3 = find(U(3,:) == max(U));   
9 H5 j# N* L- O' z* w9 Kcity(id1)  
& p3 Y7 @- P6 o% |$ ucity(id2)  
! g/ ~# l8 y7 E' U0 X9 Z/ e% ?city(id3)  
  v9 Z* I/ z4 ?
为什么一查看聚类结果
就出现错误呢
请大神指点一下为什么
总提示指数超过矩阵维数
' h8 s6 ~) v, H2 ]$ ~9 }( u

我爱照镜子

那三个表情是什么啊

净心、精心

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

：：：：：：：：：