- 在线时间
- 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));# O: ?1 d- ?2 ]+ H( y- y
x=zeros(m*n,3);
o* g, O$ p% S* F% hfor i=1:m
3 Z, K. ^/ L$ w; [: _ for j=1:n
- u9 q# @: y% @( z. Q x(i,1)=j;% y1 t% s7 I7 c1 F! A. N- [7 `; _
x(i,2)=i;
a0 T9 n' o0 `& m x(i,3)=data(i,j,1);
' w6 S1 e3 c) u3 D( V3 E end- I0 |: o9 x& O; Z f, l- w
end
) w( j6 E: |8 k8 ~+ I |& \% dX = zscore(xdata);
1 A( H: {( o+ m" d2 [2 [1 D" a u" I0 _+ g7 H0 l8 k
& b" C9 o6 p' F0 V( E%****************************进行模糊C均值聚类******************************
% E" b5 a% s- s" K# G" U% 设置幂指数为3,最大迭代次数为200,目标函数的终止容限为1e-6,不显示中间迭代过程
4 |. J, r$ v7 h! V# H/ r5 Ioptions = [3, 200, 1e-6, 0];/ Y4 M4 g I, l* }
% 调用fcm函数进行模糊C均值聚类,返回类中心坐标矩阵center,隶属度矩阵U,目标函数值obj_fcn) v) l8 l. [/ H2 y3 g3 b; J T
[center,U,obj_fcn] = fcm(X,3,options); C3 z) \* p3 p, p/ A" @
5 W1 E6 g" W% H- H! I) T
9 g. Y9 u- G7 P; F7 ~0 @%****************查看聚类结果*****************************
0 S) `( H& `2 n* d$ c6 Eid1 = find(U(1,:) == max(U));
' [# A% v( {) u# w: Y2 `4 n% bid2 = find(U(2,:) == max(U));
* Y6 m) d" R( R$ J% C" j3 ?id3 = find(U(3,:) == max(U));
0 e: T0 n) W! }$ `: n* j0 I* `city(id1) 4 w$ v3 I( G1 o
city(id2) $ s$ R3 {. F9 B" S( y. B# X
city(id3)
" Y, g; u9 }9 o$ S, z$ ~2 f1 A& W1 q5 o, Z4 J9 L4 a
为什么一查看聚类结果
/ D. p$ z& y) W: P1 k就出现错误呢9 u" Y$ Q" q3 |! ?
请大神指点一下为什么
& A1 T+ l+ S! k0 ?9 g总提示指数超过矩阵维数
5 W: o8 ~% L& y% d2 X8 O. ?. w' v, ^+ S! j8 L' G$ H
|
zan
|