- 在线时间
- 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 编辑
q' b. i( Z7 h' @4 f7 g! [1 Z7 o5 B/ z* \9 p
[m,n]=size(data(:,:,1));
1 T1 z h5 L& ^; v) V/ U4 @x=zeros(m*n,3);' U! {: e8 t1 y1 R1 y6 W% {
for i=1:m
/ d! i3 }, j' u. ` for j=1:n3 M( t7 S9 S" G1 P
x(i,1)=j;
- ]6 ?+ {8 }! I1 D x(i,2)=i;8 ~0 n' e" [% a P/ q! v* I
x(i,3)=data(i,j,1); ( G4 r/ a) d7 y0 S
end: O1 Y9 M$ ^9 n4 b2 d; ?
end
% C0 b% `) k! ]( [7 o4 oX = zscore(xdata);% h" m$ |8 E) P# A
* W; D8 O+ }8 `. y
, f/ L+ M% `3 a7 U( i: b9 C3 P7 {9 S%****************************进行模糊C均值聚类******************************
3 T* o" J' [. I2 R& G' }% 设置幂指数为3,最大迭代次数为200,目标函数的终止容限为1e-6,不显示中间迭代过程
5 u4 X3 g$ ?; p- Toptions = [3, 200, 1e-6, 0];
& D; Y) _5 o& t, [% G% 调用fcm函数进行模糊C均值聚类,返回类中心坐标矩阵center,隶属度矩阵U,目标函数值obj_fcn2 }5 j4 b! M5 M, P+ n
[center,U,obj_fcn] = fcm(X,3,options);. L5 c- ^+ Q: x( n/ L6 t* Z7 O
' h: x5 I- p, n3 C6 R1 V/ o
2 M# \2 f3 ]; |% A$ r%****************查看聚类结果*****************************5 w! b8 a. h5 \
id1 = find(U(1,:) == max(U)); 2 w) ~4 @ a- U& F
id2 = find(U(2,:) == max(U)); 5 G& p* J1 U' }9 T
id3 = find(U(3,:) == max(U)); ) V2 L/ M0 W) C5 K2 ]
city(id1) : Y3 Z7 t9 Z M1 V5 d7 c
city(id2) 2 \8 `1 n! `# q
city(id3) 6 p$ e; \7 _# K* W
2 r+ f( \! }3 P6 H' Y- K$ o为什么一查看聚类结果
" `4 L& `) c8 d% z7 f: U3 M/ `就出现错误呢" ~" `. s0 Z0 e5 ]
请大神指点一下为什么
" e. b' O. i" e总提示指数超过矩阵维数
" e8 z% \" C6 K+ b! @
' f! ]( x2 q0 u
3 E/ p, C3 C7 L, q |
zan
|