数学建模社区-数学中国
标题:
在线等回复
[打印本页]
作者:
509
时间:
2013-8-20 15:55
标题:
在线等回复
本帖最后由 袁海亮 于 2013-8-20 16:15 编辑
, G0 r8 R' M" x- f- H4 x
# j3 l+ S: v; o/ Q& H) z1 v/ Y. ]
[m,n]=size(data(:,:,1));
' E, v( y8 [) X" u
x=zeros(m*n,3);
$ f& f& y) S! E* p, e
for i=1:m
9 V" K* V9 G, m1 n# N
for j=1:n
1 k$ }/ K: ], K" ^
x(i,1)=j;
/ \9 c3 y+ Q& K5 E) M
x(i,2)=i;
0 E' d+ ?6 j* ]* r9 Z% f
x(i,3)=data(i,j,1);
# p+ C7 n/ K+ q) s2 u
end
* V6 f2 o3 I, A: f4 a% p
end
% Q) k' X' g- S
X = zscore(xdata);
, S- F$ v* m/ v8 v7 x2 s' l S7 R
# T& w" k: L! T0 G0 v$ _
( c2 C+ O: [" o! f
%****************************进行模糊C均值聚类******************************
2 e* O# w. P' K [
% 设置幂指数为3，最大迭代次数为200，目标函数的终止容限为1e-6，不显示中间迭代过程
/ s- Z" ^7 k; y( q+ n0 [$ ^1 i
options = [3, 200, 1e-6, 0];
0 c) n' ^& f5 W
% 调用fcm函数进行模糊C均值聚类，返回类中心坐标矩阵center，隶属度矩阵U，目标函数值obj_fcn
5 ], t3 d4 b/ q# Z
[center,U,obj_fcn] = fcm(X,3,options);
1 D; k3 v9 ]7 D% K m$ ^. H
8 j$ e$ C& S6 D3 Y* q
+ f: M# N- ]& F2 X% ^9 g
%****************查看聚类结果*****************************
4 m% Q0 c% U( g, F, `2 }
id1 = find(U(1,:) == max(U));
& I6 h2 Z2 ~7 b5 |: x
id2 = find(U(2,:) == max(U));
! V$ a, q% ^3 L: s. \# r. T
id3 = find(U(3,:) == max(U));
+ O8 _# T& X7 W6 m8 g) ^* `
city(id1)
7 o- i1 ?9 H: a$ _$ F, p
city(id2)
9 V2 r* b0 Q; `1 J" r
city(id3)
7 S \) l, ]! d4 Z, `% z
0 \0 ^6 a8 |9 |) L2 p& g
为什么一查看聚类结果
' W Y0 U9 b3 X/ q+ C" A
就出现错误呢
7 L- A; t) h& K' P
请大神指点一下为什么
- D. t. u4 q8 E
总提示指数超过矩阵维数
4 N4 A% N& F+ l( p. U. ?8 a
9 p: r$ h& ]! Q! b1 I
0 O. |! v5 H- l4 T7 f/ n8 Q: O- B& `
作者:
净心、精心
时间:
2013-8-20 16:40
http://www.madio.net/thread-198267-1-1.html
0 k& d9 E9 q9 h
亲，在有问题时，记得标题要写相关的困惑和哪一方面的知识，这样才能能快的得到帮助，要不一般看到你这标题，没人愿意点开看的
欢迎光临 数学建模社区-数学中国 (http://www.madio.net/)
Powered by Discuz! X2.5