模糊C均值聚类

509

[m,n]=size(data(:,:,1));
& ~3 ?, d( P2 P7 S- w( e  R6 ~x=zeros(m*n,3);
4 o' C, [8 s- ^; a- e+ R) q: N- j6 dfor i=1:m
    for j=1:n
       x(i,1)=j;
  T: [3 o; c) C; P3 S       x(i,2)=i;
0 i4 V- v' i' i5 P! P: @1 R% ^       x(i,3)=data(i,j,1);
    end
- {8 ]) o( N6 f0 d/ Jend
% G4 T+ e; D+ SX = zscore(xdata);

# }8 c2 k: d; N+ q$ \  z6 M
+ S1 x: a' ~/ }# ]%****************************进行模糊C均值聚类******************************
% 设置幂指数为3，最大迭代次数为200，目标函数的终止容限为1e-6，不显示中间迭代过程
options = [3, 200, 1e-6, 0];
, n7 T8 V% I8 d# B) s% 调用fcm函数进行模糊C均值聚类，返回类中心坐标矩阵center，隶属度矩阵U，目标函数值obj_fcn
7 u; m; o# `7 i+ D  `[center,U,obj_fcn] = fcm(X,3,options);
  D- e' k2 N' N% V3 h7 s
7 \+ P5 h- y) i
% G" \) d# |+ a! d%****************查看聚类结果*****************************
6 F; a- H. q7 m. s1 c$ pid1 = find(U(1,:) == max(U));  
$ Y6 M! Y( ?! bid2 = find(U(2,:) == max(U));  
id3 = find(U(3,:) == max(U));   
1 H  H" t; j" _& `' i( F+ Ecity(id1)  
city(id2)  
city(id3)  
3 {0 {1 K1 Y/ Z$ A3 `9 g2 R0 M+ H0 K
4 y/ S4 c0 D1 s" Q# ^为什么一查看聚类结果
就出现错误呢
' s- s( w5 C! [9 [: |请大神指点一下为什么
总提示指数超过矩阵维数

- c& a* q8 U1 g) i. p) C

我爱照镜子

那三个表情是什么啊

净心、精心

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

9 H: t/ j/ }2 w# m( y+ g：：：：：：：：：