模糊C均值聚类

509

[m,n]=size(data(:,:,1));
, O& ]7 ?& A/ H  _x=zeros(m*n,3);
for i=1:m
    for j=1:n
: k' H6 w+ J% e; h       x(i,1)=j;
       x(i,2)=i;
4 ^5 v9 V( I  B" I: h2 d9 K+ ]; ^1 `       x(i,3)=data(i,j,1);
    end
! g6 B$ P9 l  }+ ^/ Aend
X = zscore(xdata);
/ U" p7 z4 ]8 A8 Z

7 G" A5 h1 j! O2 ^0 w( ?2 }%****************************进行模糊C均值聚类******************************
  j6 v( x4 @4 [& ~# j2 a# e% 设置幂指数为3，最大迭代次数为200，目标函数的终止容限为1e-6，不显示中间迭代过程
options = [3, 200, 1e-6, 0];
% 调用fcm函数进行模糊C均值聚类，返回类中心坐标矩阵center，隶属度矩阵U，目标函数值obj_fcn
[center,U,obj_fcn] = fcm(X,3,options);


%****************查看聚类结果*****************************
id1 = find(U(1,:) == max(U));  
- }) T5 r# u- \; O3 E" {id2 = find(U(2,:) == max(U));  
) }' z1 h/ z, Y/ J* [" m8 `id3 = find(U(3,:) == max(U));   
city(id1)  
1 R+ C& P! ?- d2 Kcity(id2)  
city(id3)  

) C4 G  X$ p1 D% A为什么一查看聚类结果
  j/ K1 j2 ?! B- b9 B就出现错误呢
请大神指点一下为什么
总提示指数超过矩阵维数
2 O% A, W1 ^% p* w4 j, [- R' O7 {
* ?* A; X7 S/ X4 |/ x5 l0 R% h

我爱照镜子

那三个表情是什么啊

净心、精心

我爱照镜子 发表于 2013-8-20 16:54
' A1 v* N$ O! k  K那三个表情是什么啊

' Y2 a6 T# ~5 h4 R  r- {：：：：：：：：：