- 在线时间
- 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));$ Z6 S' t8 ^2 U, U/ `% {6 p0 c
x=zeros(m*n,3);# B8 B& i9 L* i; ]* K
for i=1:m
z, I J' _$ B/ x! R for j=1:n
2 k2 J& Z5 i% p8 i) V4 F- ] x(i,1)=j;1 e! t, ?8 W4 K' L( y' j3 P
x(i,2)=i;5 M1 ^& e3 B6 a9 q
x(i,3)=data(i,j,1); 4 W3 _. }/ h3 Q2 P( }" B
end3 d# g# k1 x( A( b& p5 \- |6 x
end V4 X. K. v) h) D
X = zscore(xdata);
; F2 X0 w) e9 J$ T
3 Q7 n2 w" I. I+ _
8 |* H( O3 b3 e8 h5 a4 s7 r%****************************进行模糊C均值聚类******************************
0 @" j3 ^& l/ v% 设置幂指数为3,最大迭代次数为200,目标函数的终止容限为1e-6,不显示中间迭代过程
( z4 t( K3 b) K, [+ W ?options = [3, 200, 1e-6, 0]; v0 q1 Y7 ]- x( ?9 F2 d
% 调用fcm函数进行模糊C均值聚类,返回类中心坐标矩阵center,隶属度矩阵U,目标函数值obj_fcn
& o; s. t) E i: b7 m[center,U,obj_fcn] = fcm(X,3,options);3 u _3 Q$ v8 Z$ Y/ w
$ J2 ]' | P9 Z6 G
K0 C( s) z* b%****************查看聚类结果*****************************3 S$ J( s& l3 Q( ?1 `! }
id1 = find(U(1,:) == max(U)); ( C/ I2 |3 P6 r; y
id2 = find(U(2,:) == max(U)); - w6 ~3 J6 Z: z3 K2 S
id3 = find(U(3,:) == max(U)); " T$ m) Q/ a, M# \
city(id1)
! `& A- t2 H) c, u% acity(id2)
& B2 n& |5 Z; M& }- R. {city(id3) 3 G9 Y* a2 ]9 E6 N6 X
7 J9 y2 Y Z ~) C' G# u
为什么一查看聚类结果7 P( E( a) P! M9 D( W* @
就出现错误呢
; C: D2 R: O, b0 a) o请大神指点一下为什么
, _ a& A# e, J; b0 z( P2 q5 p总提示指数超过矩阵维数
* Y5 Q {$ ^& k6 e, Q+ ^4 O
- Y. `: K& J, s |
zan
|