- 在线时间
- 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 编辑
U0 m9 X/ ?" g' y, P; Q; b) n( `1 {" }( p6 M
[m,n]=size(data(:,:,1));
$ w' z. Q& a9 U4 o/ m" Nx=zeros(m*n,3);' @ n* E! P r7 Z* \2 c
for i=1:m
" Q8 v/ V1 k" h1 V* J7 ~ for j=1:n$ Y5 G" Y6 F4 _% I ^8 O4 E
x(i,1)=j;( l' p1 ^, p" s7 [2 ?$ K
x(i,2)=i;4 _4 ?3 j M7 T' ^+ r& x9 z
x(i,3)=data(i,j,1);
* r( S8 ^1 A; x0 N! n8 h) Y end
7 s% k$ e1 V8 V' H% I) [" g% [! _end7 i" c0 U4 P/ Y% x* Y$ S4 I( ], a
X = zscore(xdata);
* @0 b% u+ {& R& H) R/ c) C3 [( w
* ]) V. r# @1 Z0 B" B. y' g+ ^8 B5 I
%****************************进行模糊C均值聚类******************************( o% F, S( z* a( b1 D
% 设置幂指数为3,最大迭代次数为200,目标函数的终止容限为1e-6,不显示中间迭代过程, n) Q6 J: p m& ]' e7 f6 d
options = [3, 200, 1e-6, 0];! O- _/ W w1 ~5 Q- a0 Z
% 调用fcm函数进行模糊C均值聚类,返回类中心坐标矩阵center,隶属度矩阵U,目标函数值obj_fcn
( F2 c/ t4 e" |2 M' v* O[center,U,obj_fcn] = fcm(X,3,options);
: n- i5 V4 x3 o+ E7 |5 \& X8 }2 B6 W+ Q/ |4 F) q+ I2 I+ P
' \1 X* ?3 P0 `' W- [%****************查看聚类结果*****************************4 k/ N4 `+ ]4 `" B% e
id1 = find(U(1,:) == max(U));
% `, f2 y' w, z# Y* cid2 = find(U(2,:) == max(U));
, i2 u1 M/ U6 ?' C5 c9 gid3 = find(U(3,:) == max(U));
/ ?5 B4 O3 [% n( s! R' }# Kcity(id1)
0 y; ^, i; _" H' F" scity(id2)
4 E7 ]) R# V' _5 V# s- h! ucity(id3) / v: r5 j4 Z9 y% b. b! d8 _3 i
( B A' z' \: \+ m7 d为什么一查看聚类结果0 O F1 o! E0 Y% K/ |, L) n
就出现错误呢5 j" i- }8 _/ y0 E& I7 W* C+ s
请大神指点一下为什么0 e: v5 ^6 Z) P+ W
总提示指数超过矩阵维数$ w$ h/ m w3 ^# w! {& G9 ^
$ Q* t8 B* H. [- j
, y2 K6 v1 M" k |
zan
|