模糊C均值聚类

509

[m,n]=size(data(:,:,1));
# t6 E! y+ W% j0 ]6 ~: Dx=zeros(m*n,3);
for i=1:m
& n: M% D/ W0 ?5 m) I! G5 `" W1 y    for j=1:n
       x(i,1)=j;
( c" @; y. X# `0 Y       x(i,2)=i;
       x(i,3)=data(i,j,1);
    end
end
! Y$ c2 l2 n2 b& Y9 |0 mX = zscore(xdata);
# m. q( H2 T4 V' a) G; O
  _, c& q; X6 k
%****************************进行模糊C均值聚类******************************
: j0 ]2 G$ v% [" t( |% 设置幂指数为3，最大迭代次数为200，目标函数的终止容限为1e-6，不显示中间迭代过程
options = [3, 200, 1e-6, 0];
% 调用fcm函数进行模糊C均值聚类，返回类中心坐标矩阵center，隶属度矩阵U，目标函数值obj_fcn
$ K$ }" f# t( J0 m9 G" `' S1 i[center,U,obj_fcn] = fcm(X,3,options);
" B" q5 B5 P0 v. j/ K( P

%****************查看聚类结果*****************************
id1 = find(U(1,:) == max(U));  
id2 = find(U(2,:) == max(U));  
0 Z9 T! H& t* ?5 Rid3 = find(U(3,:) == max(U));   
$ ^& Q( a8 q! R+ f, ccity(id1)  
city(id2)  
city(id3)  
9 z% ~$ d) E" Q& {  L2 `* A5 y6 R
为什么一查看聚类结果
: t; h1 C% G% h4 H就出现错误呢
4 T; B" K, U3 u9 W/ M0 C0 j# Q请大神指点一下为什么
/ Z6 H( w) E  x1 [总提示指数超过矩阵维数

我爱照镜子

那三个表情是什么啊

净心、精心

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

：：：：：：：：：