模糊C均值聚类

509

[m,n]=size(data(:,:,1));
x=zeros(m*n,3);
3 c" k: P- G! X0 ~8 Tfor i=1:m
    for j=1:n
1 C1 g' m# l9 L. d+ w       x(i,1)=j;
       x(i,2)=i;
% z- S) i( ?6 y+ r" ^2 i# C7 }       x(i,3)=data(i,j,1);
    end
: G' d, a* d, S; F% k7 Lend
X = zscore(xdata);


. b) z$ X- L8 f( d9 S5 H' v%****************************进行模糊C均值聚类******************************
% 设置幂指数为3，最大迭代次数为200，目标函数的终止容限为1e-6，不显示中间迭代过程
options = [3, 200, 1e-6, 0];
% 调用fcm函数进行模糊C均值聚类，返回类中心坐标矩阵center，隶属度矩阵U，目标函数值obj_fcn
[center,U,obj_fcn] = fcm(X,3,options);

! i: ]: c: v8 B( x  ]0 v
- n8 N# N& N' k%****************查看聚类结果*****************************
id1 = find(U(1,:) == max(U));  
5 X0 ^# x; N; D2 q5 [7 y: r8 p% zid2 = find(U(2,:) == max(U));  
/ f" y, n: \. c- g: Z) vid3 = find(U(3,:) == max(U));   
city(id1)  
city(id2)  
city(id3)  
5 [' \0 W2 v' A4 @& h
为什么一查看聚类结果
就出现错误呢
$ E/ n# p8 d3 U$ h+ U+ Q2 G0 L请大神指点一下为什么
5 u  t$ ^  d! I- e: n" m总提示指数超过矩阵维数
( K* H& b* G5 \9 B  u* {- V

我爱照镜子

那三个表情是什么啊

净心、精心

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

：：：：：：：：：