- 在线时间
- 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 编辑 0 o8 g* f& L2 u5 r& {8 z
* C) w+ q1 s( B# ~) ^$ p8 ?, ?[m,n]=size(data(:,:,1));
( Z% P+ f# d. I7 m( fx=zeros(m*n,3);
$ e$ D9 x6 ~, |+ e# R2 pfor i=1:m
9 U) {5 { [1 ~" F, n for j=1:n
+ u0 L" n; ^; i J$ o# n- i; T x(i,1)=j;
2 r3 ?$ c! q, ?/ h8 _' @/ F x(i,2)=i;5 L' n$ E* d P1 J5 [% l
x(i,3)=data(i,j,1); 8 w; |- v, Q+ j2 c1 v' R
end9 c: c8 `0 y) y- ?. {2 U
end
4 [$ h7 j& p1 c F- D/ y: wX = zscore(xdata);
) Y6 G- m: \' z
. A- B( p2 a; L+ F1 D
1 L' z. k" c- K, W( x%****************************进行模糊C均值聚类******************************
7 c( g2 v5 e T% 设置幂指数为3,最大迭代次数为200,目标函数的终止容限为1e-6,不显示中间迭代过程3 U- `8 r6 }+ A0 F
options = [3, 200, 1e-6, 0];
) [ i3 b" V1 ]# y4 F2 h v7 y1 r% 调用fcm函数进行模糊C均值聚类,返回类中心坐标矩阵center,隶属度矩阵U,目标函数值obj_fcn8 Q8 L8 J7 A/ A9 F+ T6 M
[center,U,obj_fcn] = fcm(X,3,options);
! Q/ d* M5 d# S& z- H( R, Z1 X0 u( G7 h4 Y {0 @
8 u. S: l& M: } ?& k" R%****************查看聚类结果*****************************
- M! V9 @ e) _' k' o& did1 = find(U(1,:) == max(U)); 8 b, K9 y+ ~* s- l8 d& ]4 t* J- F! r
id2 = find(U(2,:) == max(U)); . t7 |3 q) c2 |; M
id3 = find(U(3,:) == max(U)); - Q% g1 ^; X/ b
city(id1) # c& [5 c) O6 J& {9 T" Z- g
city(id2)
- N4 V- O5 x8 H" v: f1 bcity(id3) ; _: p- m7 J; [3 Z' `* ^
& E5 V! |& G" O8 v# R- E
为什么一查看聚类结果' {. j) f, S- N2 y! [ |, q% L
就出现错误呢
- [9 T/ [# d9 F y! V m请大神指点一下为什么6 M( O0 x! c, x+ n+ S" f
总提示指数超过矩阵维数
+ w1 a r0 x5 s# s, I: |# ~) ~9 j8 M5 k* k$ z1 Z
5 ?3 t# V- S" d0 q6 m
|
zan
|