- 在线时间
- 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));
" x- A& j( x3 Yx=zeros(m*n,3);
8 U# `, c' s, Q9 Dfor i=1:m
% }" L6 u& q& \4 @& U% o8 D) s9 F for j=1:n
' k: D2 K+ x9 k j# M0 } x(i,1)=j;
. S6 E- Y/ s) E% u# X x(i,2)=i;
$ \ V+ u. y+ M$ \* P x(i,3)=data(i,j,1);
5 Z+ U" D( | g1 u% U end
3 `0 j- d2 ~0 ]3 ?" j# gend
% t! A# M9 f! R1 xX = zscore(xdata);$ e( }% A0 e8 x' H* c
1 y0 t3 i7 [1 z# ?. m: G* L( u4 y3 H0 Q& s6 q; t, I
%****************************进行模糊C均值聚类******************************
, x+ F& b+ b( t% 设置幂指数为3,最大迭代次数为200,目标函数的终止容限为1e-6,不显示中间迭代过程0 V# k! N& b1 @4 j7 V
options = [3, 200, 1e-6, 0];
) R7 E* Y' p# @% 调用fcm函数进行模糊C均值聚类,返回类中心坐标矩阵center,隶属度矩阵U,目标函数值obj_fcn" R6 U7 y$ t9 ]) F9 N. N
[center,U,obj_fcn] = fcm(X,3,options);: m( R0 N8 f. `& K! ~
5 ~, d7 V! _ k& H! d1 V6 J3 A" c
* d# U; M' I" @. ~, u ]
%****************查看聚类结果****************************** B9 b. D4 Q+ b- i
id1 = find(U(1,:) == max(U)); " ?! u" o- K+ A/ k& d
id2 = find(U(2,:) == max(U)); % D& r6 r* C: u- [: I
id3 = find(U(3,:) == max(U));
% w& C; d: W s- Z: c6 W4 J+ zcity(id1) ; n H" j9 i v' K$ v
city(id2)
+ }" C% T3 E o0 }9 Hcity(id3) 9 f/ w- y( m: S' g3 p5 b; R7 N
S2 K1 J u/ D为什么一查看聚类结果
( H9 J9 _2 F3 Y4 h6 Y: j就出现错误呢
3 k0 s! _; {2 r, b4 c2 k- A请大神指点一下为什么2 o/ q: Z, ?! H3 T9 {8 G+ G4 U+ \
总提示指数超过矩阵维数
" j7 |9 i |* V5 p
% F# x" g4 g# |* w& o; _9 ^5 s, ] |
zan
|