- 在线时间
- 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
- 自我介绍
- 东北农大本科生
 |
本帖最后由 袁海亮 于 2013-8-20 16:15 编辑 2 K/ L3 ?6 `6 p- C0 Q) S1 y
1 Q9 G* B: [- `3 k6 |' t0 q
[m,n]=size(data(:,:,1));& q) n% c3 k) s# X* d
x=zeros(m*n,3);
$ r1 V" n$ g* P, f' mfor i=1:m7 C/ G8 J. D/ v4 V' C
for j=1:n
" e6 p/ a' C4 G) b4 H& N x(i,1)=j;+ V: O( |2 M/ ?1 w& g: X
x(i,2)=i;7 o& t) v9 Z9 e5 ?4 o
x(i,3)=data(i,j,1);
0 H2 C1 x: e, y end! h' c' R* Q/ s+ V6 S
end+ w- @, ]# p9 t& k
X = zscore(xdata);
/ X* j9 R9 s" {+ @ y7 ?7 o8 W- g8 ^4 S, _% j' N$ }
% ^) P% N' u' O5 t& Y%****************************进行模糊C均值聚类******************************
; x6 w- a9 o* B' ]9 _1 }7 u% 设置幂指数为3,最大迭代次数为200,目标函数的终止容限为1e-6,不显示中间迭代过程3 A/ x4 e1 \" d+ C
options = [3, 200, 1e-6, 0];
$ W# f# X; a; c% 调用fcm函数进行模糊C均值聚类,返回类中心坐标矩阵center,隶属度矩阵U,目标函数值obj_fcn0 H) f0 c- O1 O
[center,U,obj_fcn] = fcm(X,3,options);* i5 Q( m8 e! w
# N. W! P0 ?; n; u# k2 K+ s, o
* d, \0 u& R( _- B( J8 ]%****************查看聚类结果*****************************
- T% _- w( I$ t, ~id1 = find(U(1,:) == max(U));
- Y5 N1 Q8 U* L4 T3 w* |id2 = find(U(2,:) == max(U)); . g9 a8 V- T: ~9 B" C9 |0 e* ]
id3 = find(U(3,:) == max(U)); # G {+ f* t8 e! Z* b, N( `
city(id1) : c0 I. Q' P- K+ B5 p3 w9 |2 u8 J
city(id2) ) H: ^; L0 W- E5 \) T* T6 U
city(id3) 3 }# n& o/ v/ ~* f7 C
2 U% }" P+ V" y9 r& I; L! E$ h; V
为什么一查看聚类结果+ @. M L; Y: h' h7 E0 G' B0 }
就出现错误呢8 Q$ d5 l5 b$ c. x% b' `# a4 _
请大神指点一下为什么* P" e: o. {% @* z; ]# A
总提示指数超过矩阵维数
" t2 f; I) I. r0 {' z, H5 V2 b+ m1 W' ^" k# U: E; O8 n6 z N5 M* B
6 _& |' L, z' D/ R( o2 @$ q$ u |
zan
|