数学建模社区-数学中国
标题:
在线等回复
[打印本页]
作者:
509
时间:
2013-8-20 15:55
标题:
在线等回复
本帖最后由 袁海亮 于 2013-8-20 16:15 编辑
( h3 j; a- |0 z7 T9 _1 T# n
* ~ g$ A U0 g7 Z1 \' U
[m,n]=size(data(:,:,1));
" P% G. Z' s' d8 \) o0 b% u
x=zeros(m*n,3);
( F4 A7 j4 r1 R. D* f/ v
for i=1:m
. H4 \1 O% p0 x! B
for j=1:n
: o4 L7 h+ c# k6 r& `. Y8 a
x(i,1)=j;
4 v& f; }, H# m0 Q9 q) s& @2 G
x(i,2)=i;
* E3 z' I' G9 P: K# @- C, L+ f; b
x(i,3)=data(i,j,1);
/ V+ i% F6 u- E% q$ B9 Z
end
$ n/ g! U U M2 B' U8 [% H( ?7 }4 B
end
: T. m T& w) ^2 Y5 U" H
X = zscore(xdata);
7 S& [) [- O2 J9 }6 a# Z0 @
& A, S" n7 B, ~0 S
- O6 [- J+ P! I$ N% K
%****************************进行模糊C均值聚类******************************
) Q9 _* N- [# L$ w$ g1 ^& u4 W4 x/ P
% 设置幂指数为3，最大迭代次数为200，目标函数的终止容限为1e-6，不显示中间迭代过程
9 s$ l& c: l1 y
options = [3, 200, 1e-6, 0];
9 G9 q; j5 a* W1 `- o4 d7 a
% 调用fcm函数进行模糊C均值聚类，返回类中心坐标矩阵center，隶属度矩阵U，目标函数值obj_fcn
" N6 R1 ?8 j! H- @4 q E" M% T. L! \4 a% }
[center,U,obj_fcn] = fcm(X,3,options);
- x! K4 J/ x" q& X
M$ U# A# J9 c8 v9 V9 {
) b* J- H2 [: @5 k1 M- _* u+ N# A7 }
%****************查看聚类结果*****************************
; f" G( h7 v. R( G* X2 c0 w# ^
id1 = find(U(1,:) == max(U));
0 N5 M" I5 O5 X+ S( C; f" Z% R3 X
id2 = find(U(2,:) == max(U));
) z# o7 p+ `" C* n: t/ y6 M! d' @
id3 = find(U(3,:) == max(U));
2 `7 \( e/ x* `3 C8 W, N
city(id1)
9 {* `+ W$ Z+ Y$ J% P
city(id2)
* x w$ T7 t1 }; j2 N9 _' w
city(id3)
C w+ W0 |5 {0 w+ P# I3 _
5 k& p1 _. T- g) _, e$ m
为什么一查看聚类结果
( Y# T) i/ g4 E4 E
就出现错误呢
7 Z0 Q! ~- s# e0 _/ u
请大神指点一下为什么
1 ]+ l3 ]( }+ `5 v5 e
总提示指数超过矩阵维数
" s5 ?, i9 Z, N* @$ j
+ k3 r& z# w' w6 ~0 S. }
$ l t! A; \) |
作者:
净心、精心
时间:
2013-8-20 16:40
http://www.madio.net/thread-198267-1-1.html
# `7 r$ @7 }3 M H, D; O
亲，在有问题时，记得标题要写相关的困惑和哪一方面的知识，这样才能能快的得到帮助，要不一般看到你这标题，没人愿意点开看的
欢迎光临 数学建模社区-数学中国 (http://www.madio.net/)
Powered by Discuz! X2.5