- 在线时间
- 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 编辑
! j, F7 o2 W( W* W$ z5 h, G) O6 x! b* E8 ^ e2 ^* L
[m,n]=size(data(:,:,1));7 N' J+ ?! c; l
x=zeros(m*n,3);! s+ v/ h" W: W0 b
for i=1:m
5 I8 b) @' q4 M7 u- { for j=1:n) n: C5 N. M: T$ Q" f
x(i,1)=j;
: C" g' Z3 b5 c( n x(i,2)=i;
# R/ S+ a: D% `1 w$ _ x(i,3)=data(i,j,1); 1 W- T% t+ a8 u! }5 e
end
& ]7 p# [5 ]; v2 Q; b8 j- s8 |end1 D" ?2 [( J8 L" W6 }% O7 A
X = zscore(xdata);6 j4 d6 T6 H- O1 M3 u
8 }$ g! D1 y7 w* H' a7 R2 F! Y8 U$ @/ S, Q
%****************************进行模糊C均值聚类******************************8 }! r7 Y" {6 T: p
% 设置幂指数为3,最大迭代次数为200,目标函数的终止容限为1e-6,不显示中间迭代过程6 u, f/ _* _. c F
options = [3, 200, 1e-6, 0];
; P% k* y' A3 V. P0 U. m$ ]% 调用fcm函数进行模糊C均值聚类,返回类中心坐标矩阵center,隶属度矩阵U,目标函数值obj_fcn. Z& m; |$ M* k0 _) G
[center,U,obj_fcn] = fcm(X,3,options);
; {$ W, ?. A& O7 _- Z8 W) @4 M; C u; {% t$ e
/ o( S3 W4 k/ u9 X2 q4 c" {
%****************查看聚类结果*****************************- {6 V- z" V: }6 f- ~ s
id1 = find(U(1,:) == max(U)); ( \$ `' k. W4 l4 g6 [# s% v
id2 = find(U(2,:) == max(U)); 7 ~: b/ a- M o+ R: W7 `
id3 = find(U(3,:) == max(U));
# U2 i9 }$ X' A" T) ~9 ncity(id1) / E& e: j( K8 W/ p$ L
city(id2)
0 ~* g% R9 q& c* Q' i7 W5 Ycity(id3)
1 k0 I& T$ y' h) B i% b' i, Z. J- D* ^/ u& U; H8 ~8 \- \
为什么一查看聚类结果* x6 b+ f" [2 S2 S! E3 h; ~& C) {
就出现错误呢
: i% H- C8 b- Q2 }7 X0 J& e请大神指点一下为什么
# F- \9 R. T; ~总提示指数超过矩阵维数
4 G: }8 a7 f+ K' H& S! _ }# X
- |0 o0 ]' l `# i8 F8 X) h# G; ]2 h/ K; \' W
|
zan
|