- 在线时间
- 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 编辑
/ D9 n1 u" p: i# n, q( y4 [' K' g' l
[m,n]=size(data(:,:,1));
; d- P7 j3 G. [x=zeros(m*n,3);
! {' H n* Y- f# k I' nfor i=1:m6 e0 m( Z; e3 D/ \) R4 u
for j=1:n
4 O* o/ v0 Z! v7 E" `% S3 P8 a& W x(i,1)=j;+ X! W, i( t! r7 U
x(i,2)=i;: F# j6 F/ @. y2 h0 l
x(i,3)=data(i,j,1); $ |1 |/ |1 }/ a% L
end
' ]6 i' _/ E$ oend
: h8 J4 q1 `& G4 o, W+ wX = zscore(xdata);' C: a0 g; u* y! l n' V' h
( t( y: x8 n/ Q
* t% ?/ n9 D7 Y9 Q6 ?' a; K
%****************************进行模糊C均值聚类******************************$ y& Q+ K+ m Z; j5 f/ m
% 设置幂指数为3,最大迭代次数为200,目标函数的终止容限为1e-6,不显示中间迭代过程! H$ g e% _. C+ l# Z
options = [3, 200, 1e-6, 0];, a# `6 }, f7 L0 @. u
% 调用fcm函数进行模糊C均值聚类,返回类中心坐标矩阵center,隶属度矩阵U,目标函数值obj_fcn
3 A# ~" X0 @. v5 b6 s# f2 \[center,U,obj_fcn] = fcm(X,3,options);
- E6 S6 z5 Z3 g$ a4 I8 U
/ G" B% U1 ^! g& P7 x" }; ]4 O! j
6 E; Q8 i+ f4 I# u- H/ y%****************查看聚类结果*****************************
2 j$ f5 d/ D6 k) F# E' Jid1 = find(U(1,:) == max(U));
" {* U1 Q$ _# j% H. s3 b' h, Zid2 = find(U(2,:) == max(U)); 7 q0 k: |8 N2 i( N7 w# e" Q7 D
id3 = find(U(3,:) == max(U));
, a/ n3 u3 ~0 Y( Qcity(id1)
* V9 y [" K4 p( `4 |& h$ Gcity(id2) " O: W* Q9 a9 s" ` b9 A3 X4 N
city(id3)
: n i/ c0 ~# ^2 j. ~. E0 l$ f" c4 E& q) j" ~/ C: A+ \
为什么一查看聚类结果
* d6 K, f9 O) s0 R就出现错误呢( s8 T9 m) ~1 S, R% G5 }8 K
请大神指点一下为什么+ E1 \, [& l' {! {' y L
总提示指数超过矩阵维数
: M) p6 {9 w% h! h, G. V- Y2 ~" M/ u f" U1 N% _4 l' D
. V9 z- i! S6 x9 o; @) V! R& O: C
|
zan
|