在线等回复

509

+ T" a  c  t3 V1 h7 C/ I6 ?7 I
. a3 v2 ]' E8 n5 \! ]9 s[m,n]=size(data(:,:,1));
x=zeros(m*n,3);
for i=1:m
+ q8 T. k) ~3 k- k% x" E    for j=1:n
' u5 [2 |7 H5 S# S% M7 f1 ~5 u       x(i,1)=j;
       x(i,2)=i;
: g: z- F$ M& t* g+ e+ i       x(i,3)=data(i,j,1);
    end
end
X = zscore(xdata);

4 p) e) t: X2 O4 @; H
%****************************进行模糊C均值聚类******************************
% 设置幂指数为3，最大迭代次数为200，目标函数的终止容限为1e-6，不显示中间迭代过程
3 ^8 s0 j: \5 I1 \" l! Uoptions = [3, 200, 1e-6, 0];
% 调用fcm函数进行模糊C均值聚类，返回类中心坐标矩阵center，隶属度矩阵U，目标函数值obj_fcn
[center,U,obj_fcn] = fcm(X,3,options);

+ i4 r7 s/ X2 K6 V& {9 b8 d& ?
%****************查看聚类结果*****************************
  N5 ?. n) O/ X# w7 f0 Q- O% p- Gid1 = find(U(1,:) == max(U));  
2 U6 x8 b# Q8 pid2 = find(U(2,:) == max(U));  
id3 = find(U(3,:) == max(U));   
city(id1)  
city(id2)  
city(id3)  

6 {5 s# R& R6 c6 ]: F' \为什么一查看聚类结果
就出现错误呢
1 H. j( C! [1 O' X2 l请大神指点一下为什么
总提示指数超过矩阵维数

净心、精心

http://www.madio.net/thread-198267-1-1.html
亲，在有问题时，记得标题要写相关的困惑和哪一方面的知识，这样才能能快的得到帮助，要不一般看到你这标题，没人愿意点开看的