- 在线时间
- 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 编辑
8 r5 k/ Y6 C/ H- f& f3 L3 V
# b) t4 k# Y/ I O7 k' h& x6 X" D* _[m,n]=size(data(:,:,1));/ k+ I) e$ }( X) N; ^# g ]" i
x=zeros(m*n,3);; n; [% p& G* v8 }( V. w
for i=1:m
2 t/ m( O( j8 ~ for j=1:n- H% D. Q* n- I- c
x(i,1)=j;5 K9 b: j* e& z( M% E$ l
x(i,2)=i;
4 O6 W1 x) |. A# `2 c x(i,3)=data(i,j,1); 0 p% X9 U9 }1 F( a
end
2 [5 [! H C) }5 ]/ u" S3 |! z+ mend
, d7 n. Y5 V: OX = zscore(xdata);
8 c5 e3 a8 ?. Q+ p" w. a% K- T0 L, r) y
% n7 n! l2 `/ }8 D* m- U
%****************************进行模糊C均值聚类******************************
- K- h4 R7 q' ?& `- G. V+ g/ ]% 设置幂指数为3,最大迭代次数为200,目标函数的终止容限为1e-6,不显示中间迭代过程* g3 f7 P( |7 {8 `8 F P9 a7 A* [
options = [3, 200, 1e-6, 0];, T R$ R, h4 I; B0 x3 P+ |
% 调用fcm函数进行模糊C均值聚类,返回类中心坐标矩阵center,隶属度矩阵U,目标函数值obj_fcn
5 ]" Q8 P. j1 W4 `# R[center,U,obj_fcn] = fcm(X,3,options);
D7 @( r- n1 I/ h0 z
* V5 I0 F) v/ Z, G
6 e3 ]0 w; ?) }$ i* J! F9 Q%****************查看聚类结果*****************************& I/ E- W7 G, m# J4 Z8 k; b% O
id1 = find(U(1,:) == max(U)); ; L! ]5 k) ?7 F$ i
id2 = find(U(2,:) == max(U)); 6 U; M: R# V9 [4 t8 o
id3 = find(U(3,:) == max(U)); . t' t3 s& }/ T* f8 V
city(id1) 0 l, m3 X8 X0 }! G& g' D/ w
city(id2)
- _" m# z/ u2 }" [- pcity(id3)
. R: q4 ?6 ^$ x: O% K
' x( A# I8 ]$ [7 @) G! p4 U为什么一查看聚类结果
* L9 H# I ]7 J) w( N9 U { s就出现错误呢- A+ e0 |( q" ^% I1 W
请大神指点一下为什么! C P5 r& k/ L7 U8 T8 `2 N
总提示指数超过矩阵维数
9 P% b! D# |: Y# P6 e4 L2 \1 K, u4 {0 e9 H |
4 W6 K" O3 |* Y9 p% c8 A. R+ ]
|
zan
|