- 在线时间
- 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 编辑
9 z, ~* e3 P2 d2 o+ C9 |( e% n/ t% u( G# ^' ^
[m,n]=size(data(:,:,1));
G: r* s2 ]- i/ y% wx=zeros(m*n,3);$ K, b* N6 B3 B/ q
for i=1:m
) v3 z: z. @- A for j=1:n9 H k" M3 @ `8 a) V
x(i,1)=j;# e5 A$ ~ V5 I
x(i,2)=i;+ X* ^/ V/ M$ n% u
x(i,3)=data(i,j,1);
- e& e# w6 j3 S' \( c2 ~1 ~ end' G+ c' a4 J& s6 p
end
. x' t* H7 V5 z FX = zscore(xdata);
7 B" W8 O( F9 u1 R, D7 U/ m+ {5 C; _7 x/ t8 I2 F% C" P
! {- }* e6 I* m7 u
%****************************进行模糊C均值聚类******************************' f/ D( P+ T/ D9 C) {
% 设置幂指数为3,最大迭代次数为200,目标函数的终止容限为1e-6,不显示中间迭代过程
+ o6 O. a- P( xoptions = [3, 200, 1e-6, 0];: ?1 M( U' n/ t
% 调用fcm函数进行模糊C均值聚类,返回类中心坐标矩阵center,隶属度矩阵U,目标函数值obj_fcn
* P0 m: q+ T, `, F. }[center,U,obj_fcn] = fcm(X,3,options);
& j) E. s. ~. I7 x1 K3 O8 E& g) V# Z# R6 f9 a( o0 q0 Y
4 T% {/ Y2 d' C. L w' V+ u%****************查看聚类结果*****************************
/ H, S5 k4 j# v, D# T L. R' T6 m. zid1 = find(U(1,:) == max(U)); + Z" f8 y+ U7 Y8 e
id2 = find(U(2,:) == max(U)); ) k! `! q4 y/ e v. F( K
id3 = find(U(3,:) == max(U)); ) D3 a' B+ i- o+ i+ I& ?
city(id1) # E5 T6 \, U( d( r( G/ j
city(id2) / I3 X) m# R3 q0 H$ {# r
city(id3) * L D) S) l2 v' ]" j
! ^- e( k+ ]0 i2 z8 E
为什么一查看聚类结果
" Z( N) m2 k* E$ s就出现错误呢 p" u" I8 y' @6 Z/ e( v8 z/ b' K
请大神指点一下为什么% F; w8 K, G# E# p. m3 l3 \
总提示指数超过矩阵维数
, }% m/ C: Y5 s# W; S5 _! D2 [1 s' \4 N* t. \
% O1 [( u; n" P |
zan
|