模糊C均值聚类

509

[m,n]=size(data(:,:,1));
x=zeros(m*n,3);
for i=1:m
    for j=1:n
" D) n7 Q+ t7 ]; q  v       x(i,1)=j;
       x(i,2)=i;
       x(i,3)=data(i,j,1);
    end
4 j; S; }; Q4 z0 z. V7 i/ F5 Yend
. p: e# u: E: _, oX = zscore(xdata);


%****************************进行模糊C均值聚类******************************
# Z# g, f  B- A8 q% O8 V3 ^4 A% 设置幂指数为3，最大迭代次数为200，目标函数的终止容限为1e-6，不显示中间迭代过程
* E6 P$ _0 W; m5 {options = [3, 200, 1e-6, 0];
- o7 V2 I& a$ A" d- ^. X: i& Y% {% 调用fcm函数进行模糊C均值聚类，返回类中心坐标矩阵center，隶属度矩阵U，目标函数值obj_fcn
. c9 f- E9 v% _6 u8 y1 c[center,U,obj_fcn] = fcm(X,3,options);
- A  e% \5 z( F& u2 ~; M
6 T  C9 h- G9 E& P
%****************查看聚类结果*****************************
$ a) g! [8 o$ p. S* R' Q4 ?' C7 P2 Gid1 = find(U(1,:) == max(U));  
  e) M8 G" ~3 ^' w" l2 D/ pid2 = find(U(2,:) == max(U));  
- {/ v. B1 T/ }3 W9 aid3 = find(U(3,:) == max(U));   
# r( f! t; o8 f, F0 mcity(id1)  
) b1 c0 ^% |2 Y- t9 mcity(id2)  
. F: T  N( ^5 r* }8 Qcity(id3)  
% n* M3 L; E5 M. N% i; P( {
为什么一查看聚类结果
就出现错误呢
请大神指点一下为什么
8 s) x" E+ u) `& O. I1 f  w总提示指数超过矩阵维数
/ V" n6 W9 A& D- B0 \. G+ A6 {5 M

净心、精心

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

：：：：：：：：：

我爱照镜子

那三个表情是什么啊