在线等回复

509

: T6 F7 S8 Z7 F: m5 K  [( m
[m,n]=size(data(:,:,1));
# D  Y# q/ G9 X6 }" d! |, gx=zeros(m*n,3);
for i=1:m
3 V. A0 v; E) a0 Q    for j=1:n
* W, b4 }6 p' Z% ^- Y$ B, c       x(i,1)=j;
* ?/ @* @7 N& z* V2 n! q  U) [       x(i,2)=i;
# K4 O) B9 M( R       x(i,3)=data(i,j,1);
    end
end
X = zscore(xdata);
) r# g" g+ M1 Z8 H0 J% E& e
2 j/ _3 [) @& i. v- p7 B
1 V, y( T( t" I- l7 w" X%****************************进行模糊C均值聚类******************************
! l- z7 s3 h9 y) `$ \9 p% 设置幂指数为3，最大迭代次数为200，目标函数的终止容限为1e-6，不显示中间迭代过程
% [& w# F+ w, koptions = [3, 200, 1e-6, 0];
9 W) l4 ^3 x9 b5 M/ i$ l% 调用fcm函数进行模糊C均值聚类，返回类中心坐标矩阵center，隶属度矩阵U，目标函数值obj_fcn
$ J% t) A$ s  C4 y4 I" [[center,U,obj_fcn] = fcm(X,3,options);
9 ^2 f$ b# b2 }2 i: L
" P9 q/ W6 {  {( V. V& Z
- X  `2 y8 H/ K# G2 j  Z7 R%****************查看聚类结果*****************************
  E; N: _: W: J* @6 _id1 = find(U(1,:) == max(U));  
id2 = find(U(2,:) == max(U));  
id3 = find(U(3,:) == max(U));   
/ e2 w; R: t5 |. }. \4 Fcity(id1)  
city(id2)  
city(id3)  
; w' k* }$ W9 E) ]3 G/ i
" W! l% F% v$ U( X. h  R( f$ D为什么一查看聚类结果
就出现错误呢
9 c, p# q. ?3 T: S  F% G请大神指点一下为什么
( H, V* O( M% m! o2 r总提示指数超过矩阵维数
% a, O" H' A$ S, t3 c

净心、精心

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