模糊C均值聚类

509

[m,n]=size(data(:,:,1));
# D$ E( I" N' H8 s" _8 J( B; Zx=zeros(m*n,3);
for i=1:m
    for j=1:n
$ ]9 O" D& x( ]  i7 z& {% k       x(i,1)=j;
3 q1 E  H7 s7 {/ {0 D7 H       x(i,2)=i;
* {- w1 _3 H* X       x(i,3)=data(i,j,1);
7 a; D+ G( P/ K5 d& a/ k+ \    end
1 K( R) n0 l1 @% G; @end
3 U8 G& X) Z" z' }1 q( NX = zscore(xdata);
' H# O0 n# \$ E  t+ ^  p

%****************************进行模糊C均值聚类******************************
: _$ [6 {* ^6 p1 Q7 ]/ H$ @6 Q9 p% 设置幂指数为3，最大迭代次数为200，目标函数的终止容限为1e-6，不显示中间迭代过程
8 N; G& W  y2 |+ J6 ooptions = [3, 200, 1e-6, 0];
% 调用fcm函数进行模糊C均值聚类，返回类中心坐标矩阵center，隶属度矩阵U，目标函数值obj_fcn
+ v3 p0 J; [; \[center,U,obj_fcn] = fcm(X,3,options);
, `* _  f) w2 d3 G0 o
5 h( Y+ @' ]% J2 {1 b
; E1 n+ e! v/ s' m1 `4 s7 H%****************查看聚类结果*****************************
id1 = find(U(1,:) == max(U));  
id2 = find(U(2,:) == max(U));  
6 \5 j7 T9 ^4 b" x# aid3 = find(U(3,:) == max(U));   
2 V( [, C4 v5 Y. K1 N, Ucity(id1)  
) Z! J" e5 C; J" N6 f( _city(id2)  
city(id3)  

为什么一查看聚类结果
& h( d, ~6 J( H* S# I就出现错误呢
请大神指点一下为什么
总提示指数超过矩阵维数
" n8 C8 P# J! @7 P- @" p; G

我爱照镜子

那三个表情是什么啊

净心、精心

我爱照镜子 发表于 2013-8-20 16:54
0 Y& L, K- c1 O( w# c8 D+ P那三个表情是什么啊

4 R/ x; z5 J: [  G：：：：：：：：：