数学建模社区-数学中国
标题:
在线等回复
[打印本页]
作者:
509
时间:
2013-8-20 15:55
标题:
在线等回复
本帖最后由 袁海亮 于 2013-8-20 16:15 编辑
2 u0 }* U# s5 h# D. {
: D) }8 W7 d3 \3 t0 A: J
[m,n]=size(data(:,:,1));
5 j# l w0 c! m+ v
x=zeros(m*n,3);
7 h9 ]) f8 H j( R( B
for i=1:m
* G% l" [1 U1 g, c0 H# M3 B2 V
for j=1:n
9 ]+ [; f3 Y4 S! P, O4 l0 Z
x(i,1)=j;
8 d! L7 u6 O1 _, R. h0 a9 w* J
x(i,2)=i;
1 p1 Y, w% p! F
x(i,3)=data(i,j,1);
0 q* L; Q; |! e# Z+ ~2 y
end
# n6 U' y0 P* a: p
end
' r3 b& M2 h8 U4 D( U6 l1 _
X = zscore(xdata);
' ^5 b' A$ }5 @6 k* \2 z: B3 E
, T- i$ Q" E% |$ S3 Z: y- ~
8 t& `# N% \* Q% T
%****************************进行模糊C均值聚类******************************
% M7 l" z9 Z2 O, P- c- @4 B
% 设置幂指数为3，最大迭代次数为200，目标函数的终止容限为1e-6，不显示中间迭代过程
* z+ P9 F3 O* N4 q5 P6 D
options = [3, 200, 1e-6, 0];
# `7 r6 L! z7 Z# U
% 调用fcm函数进行模糊C均值聚类，返回类中心坐标矩阵center，隶属度矩阵U，目标函数值obj_fcn
$ @. y: L- v- X/ J: F$ J1 i
[center,U,obj_fcn] = fcm(X,3,options);
4 p1 F; Q2 r8 P% r: ?
5 J# e: [+ w) K- k/ n7 r2 q' [
( r( r0 W# t7 W" i* b
%****************查看聚类结果*****************************
5 f1 ?" Y4 M5 T2 _3 S) J7 y
id1 = find(U(1,:) == max(U));
* B: h+ j% H4 T& C; z, F! \
id2 = find(U(2,:) == max(U));
. {+ o! v0 F1 ~& u2 H! p
id3 = find(U(3,:) == max(U));
( H. [! v" b) k/ I& K& q2 E
city(id1)
7 i# t6 M7 W( I' ]
city(id2)
, G) g4 }$ X) ~# C# \0 z. [) X/ C/ m
city(id3)
9 l9 ?9 r- V: [! |. X5 S
/ ?4 s8 _, {9 {& d5 v. H$ a
为什么一查看聚类结果
1 `: m6 P) c5 w$ K2 V) Q: s
就出现错误呢
$ X8 y$ Y6 l4 w3 t$ \
请大神指点一下为什么
+ `! E! K+ P& ^/ \ r/ m y
总提示指数超过矩阵维数
: \4 t6 m/ v( q& ]8 b
. m6 @* \4 |* s
6 T/ b7 u D/ N3 C
作者:
净心、精心
时间:
2013-8-20 16:40
http://www.madio.net/thread-198267-1-1.html
b7 p( x' i2 D) P' K" s
亲，在有问题时，记得标题要写相关的困惑和哪一方面的知识，这样才能能快的得到帮助，要不一般看到你这标题，没人愿意点开看的
欢迎光临 数学建模社区-数学中国 (http://www.madio.net/)
Powered by Discuz! X2.5