- 在线时间
- 25 小时
- 最后登录
- 2016-11-8
- 注册时间
- 2013-1-12
- 听众数
- 8
- 收听数
- 0
- 能力
- 0 分
- 体力
- 348 点
- 威望
- 0 点
- 阅读权限
- 30
- 积分
- 145
- 相册
- 0
- 日志
- 0
- 记录
- 0
- 帖子
- 80
- 主题
- 3
- 精华
- 0
- 分享
- 0
- 好友
- 5
升级   22.5% TA的每日心情 | 怒 2016-9-15 09:07 |
|---|
签到天数: 50 天 [LV.5]常住居民I
- 自我介绍
- 东北农大本科生
 |
[m,n]=size(data(:,:,1));& z8 I( D& S! d) t
x=zeros(m*n,3);
T/ l( T0 j- ] |$ N# v: Jfor i=1:m
" m1 v+ A, F7 r for j=1:n9 D$ z' n$ `9 i& c$ E* E6 U) E8 G
x(i,1)=j;
" f* g( \4 E/ N! ^/ k x(i,2)=i;
; }4 e( w1 d- S) A) e: | x(i,3)=data(i,j,1); & ^# `2 l1 i- Y
end* R3 R$ [2 o3 V+ W1 g h+ U3 M
end
: C+ T! ]3 v( [3 XX = zscore(xdata);' e; }3 I2 m% U# O9 D
1 X% z& b5 o2 F& U0 G
1 c, g/ P3 q- z2 |%****************************进行模糊C均值聚类******************************
4 O4 Q0 Y# ^" h" V5 U; Q4 U% 设置幂指数为3,最大迭代次数为200,目标函数的终止容限为1e-6,不显示中间迭代过程
% w8 C% G3 O. Z( [& _% Woptions = [3, 200, 1e-6, 0];
- M, f: _1 b6 @; O% 调用fcm函数进行模糊C均值聚类,返回类中心坐标矩阵center,隶属度矩阵U,目标函数值obj_fcn. f2 j* I5 T" W0 o2 h( }
[center,U,obj_fcn] = fcm(X,3,options);- Y9 K# F( t% O, U3 A
3 F& z7 j: ~% `# a# ^
3 r! G; O- z8 U2 a3 s% [%****************查看聚类结果*****************************9 j. N5 g S2 I
id1 = find(U(1,:) == max(U)); ) r; y% R# w9 A) y6 w% j* E1 z' C
id2 = find(U(2,:) == max(U));
# k9 f( D6 j: ^% a# ^( ^9 Uid3 = find(U(3,:) == max(U));
0 `; Y% D- @$ q4 _- ]city(id1) ! h' T/ t9 x5 v5 b% }) H
city(id2)
$ d7 ?! F& p. p y6 ucity(id3)
, C4 H" R& H# @9 A c, j3 {% d% ?9 Z9 h% p w. w, | G. I/ B Y
为什么一查看聚类结果
9 x) C) n- t+ \- S7 b% d就出现错误呢0 _$ E* B8 A0 j& d( s
请大神指点一下为什么5 C7 n5 n4 F% j: f3 V. I- }
总提示指数超过矩阵维数
4 v: |8 u8 X9 K& s
5 B) D5 F6 y7 u& p- g3 E |
zan
|