模糊C均值聚类

509

[m,n]=size(data(:,:,1));
x=zeros(m*n,3);
/ Y; z5 M# z- [for i=1:m
2 S3 P/ i+ H, S% M' h    for j=1:n
/ |1 q/ |7 X; \! q0 s2 p       x(i,1)=j;
! ?4 J" ^+ Q/ a( ^8 d' o  m. A3 o9 w       x(i,2)=i;
" n; H! s* ^6 {8 f9 _' f       x(i,3)=data(i,j,1);
    end
end
X = zscore(xdata);
) L& H  n' l+ G, |$ ?
; Y& @5 ~$ |; Y9 o1 m" D% n& F. P* `
: ~0 A8 m0 `3 Y%****************************进行模糊C均值聚类******************************
  S" b: K* [9 y$ A% 设置幂指数为3，最大迭代次数为200，目标函数的终止容限为1e-6，不显示中间迭代过程
options = [3, 200, 1e-6, 0];
- q4 `4 H% F2 |- l1 g: v% 调用fcm函数进行模糊C均值聚类，返回类中心坐标矩阵center，隶属度矩阵U，目标函数值obj_fcn
. k5 p1 `' k# Z5 K6 w8 C[center,U,obj_fcn] = fcm(X,3,options);
; s: c  }- l. u; ~+ r6 R% }3 W

* t  i$ t. @: l, b7 g( i%****************查看聚类结果*****************************
7 q+ Q1 v( ]4 g" ]* t* Aid1 = find(U(1,:) == max(U));  
id2 = find(U(2,:) == max(U));  
( X1 ^" J2 r, lid3 = find(U(3,:) == max(U));   
  Y/ t" k; D9 H8 pcity(id1)  
2 \1 }7 X; [$ q' K% V  L$ G6 |city(id2)  
city(id3)  

为什么一查看聚类结果
% ^. Y* H8 ?: ?3 {8 {就出现错误呢
; r% ?. n  i! H; O请大神指点一下为什么
  t0 m, g; j. ?+ v6 y总提示指数超过矩阵维数
( U2 @1 y- t% E

我爱照镜子

那三个表情是什么啊

净心、精心

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

2 Q9 g5 S/ ?( l1 h0 o9 t- Y：：：：：：：：：