- 在线时间
- 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));
& Q& K, G4 }8 u" I/ i3 F8 Ox=zeros(m*n,3);
. U8 h7 W9 q( ~( xfor i=1:m! v) J# i% c9 l* W% B) b0 \
for j=1:n' N: w1 U! I0 t
x(i,1)=j;
9 m) j0 g# e" J x(i,2)=i;2 P$ D. c# q: q
x(i,3)=data(i,j,1); + k8 J1 ~3 Z7 O1 }
end
( c' I+ A3 a/ x8 k6 Dend
; Z- l- n; j- ~" M/ k% ZX = zscore(xdata);
4 I: |- B& F( u8 P+ u+ p
, @2 o& b; Y8 Y! l* h
+ n4 A) }, J% M' o1 b%****************************进行模糊C均值聚类******************************
* d* B* t9 j8 y% 设置幂指数为3,最大迭代次数为200,目标函数的终止容限为1e-6,不显示中间迭代过程8 |5 i8 P- R1 W* _- g& q0 N, U
options = [3, 200, 1e-6, 0];
. k' P2 @3 b5 Q `% 调用fcm函数进行模糊C均值聚类,返回类中心坐标矩阵center,隶属度矩阵U,目标函数值obj_fcn
- ?3 J2 d' i. _: j- U! L$ Y[center,U,obj_fcn] = fcm(X,3,options);
6 Z1 F$ \- y# c9 T
2 F( Y3 q) Q! s" Q
, p5 c+ i; h) ?! X! ?%****************查看聚类结果*****************************$ e: A3 e+ }& S3 _
id1 = find(U(1,:) == max(U));
' c y( E% i9 ~" a4 Gid2 = find(U(2,:) == max(U));
$ c6 ^4 D/ v+ V1 L0 r& Q5 sid3 = find(U(3,:) == max(U));
/ [$ ]+ u& z1 x; k- d2 }8 Ecity(id1)
1 R5 Z& s8 V/ Y( {9 }+ O) R6 rcity(id2) / n. h7 O+ `+ x* B
city(id3)
. }1 ? `6 F+ @4 x; |" c2 t- [- |2 o: u4 h( s$ x2 C
为什么一查看聚类结果+ y3 h" t* P1 |' J( F
就出现错误呢' J; [1 f( C8 `& y& b4 _! E
请大神指点一下为什么
: V/ ^, c$ f1 a2 Z总提示指数超过矩阵维数" K: C; V R* d8 o- i! i7 U9 o
* H$ M+ F, i3 G& ^. {2 e, m4 u u( o
|
zan
|