模糊C均值聚类

509

[m,n]=size(data(:,:,1));
0 c1 G. F& f+ X# k6 ex=zeros(m*n,3);
, Q3 l; L& P+ x4 U5 K' A( gfor i=1:m
    for j=1:n
1 X9 [- o4 g& D, b( z2 }7 ~       x(i,1)=j;
0 h6 B5 D7 R" n) }# ~       x(i,2)=i;
       x(i,3)=data(i,j,1);
    end
end
X = zscore(xdata);


%****************************进行模糊C均值聚类******************************
% 设置幂指数为3，最大迭代次数为200，目标函数的终止容限为1e-6，不显示中间迭代过程
$ _% L2 C% X$ `# L. p) Z+ A: b4 Boptions = [3, 200, 1e-6, 0];
% 调用fcm函数进行模糊C均值聚类，返回类中心坐标矩阵center，隶属度矩阵U，目标函数值obj_fcn
[center,U,obj_fcn] = fcm(X,3,options);
+ N$ B" W0 V( i. `

* ?! Q: V4 G$ H! j* \6 O4 m$ w%****************查看聚类结果*****************************
2 ^6 c! d: ^+ y2 P4 Mid1 = find(U(1,:) == max(U));  
' j7 L/ v" }9 b5 @3 V. Aid2 = find(U(2,:) == max(U));  
id3 = find(U(3,:) == max(U));   
city(id1)  
  j) Z/ I+ R! s6 kcity(id2)  
city(id3)  

# G; M) V  h: n+ f: s" h为什么一查看聚类结果
4 S$ w  U' {+ O4 ^1 d; D就出现错误呢
请大神指点一下为什么
总提示指数超过矩阵维数

; L1 c: [9 z2 q

我爱照镜子

那三个表情是什么啊

净心、精心

我爱照镜子 发表于 2013-8-20 16:54
+ x2 n1 K4 a+ J6 L; b那三个表情是什么啊

：：：：：：：：：