- 在线时间
- 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 s0 H2 H8 R1 ~1 P
: x& X) H$ M! d e1 ? ^
[m,n]=size(data(:,:,1)); T2 k1 S1 \0 W* p; E/ ]9 R
x=zeros(m*n,3);2 |& h) S: b, D+ E; {& a3 r
for i=1:m4 R7 X' S+ X" L+ B
for j=1:n3 U- V7 {+ \9 h4 r. u$ L
x(i,1)=j;
7 p7 L( V/ n1 S- I8 f; H- n# z x(i,2)=i;& p# w8 l0 ^# `# B. L: b. Y
x(i,3)=data(i,j,1); 8 N& v8 _! f# p" z! N- n
end. f/ c5 Y0 T& z, z# _1 g
end. V: q& B4 P0 y1 K4 y' z: ]9 T' q- i0 i
X = zscore(xdata);8 o: K9 `( Y2 a, p- y' N
W1 \1 N8 C# V
, I3 W: h, U& \, A# b: l2 w%****************************进行模糊C均值聚类******************************
' K8 T8 M% w6 n' J9 |; E- V1 N% w% 设置幂指数为3,最大迭代次数为200,目标函数的终止容限为1e-6,不显示中间迭代过程& g, b* O8 }* i! [0 A
options = [3, 200, 1e-6, 0];, R7 j9 \& a- _* Q q2 z: p4 [
% 调用fcm函数进行模糊C均值聚类,返回类中心坐标矩阵center,隶属度矩阵U,目标函数值obj_fcn! @/ u" w. T* d7 b4 I8 _$ e
[center,U,obj_fcn] = fcm(X,3,options);. i! o6 M$ T3 i$ q( s% n2 L: O7 q9 p
$ T+ [8 h) T* j# o
0 ^, w2 B' o7 n9 D6 {8 u. F0 ^2 Z%****************查看聚类结果*****************************
5 M: V6 H% ]/ y+ l8 b* Uid1 = find(U(1,:) == max(U));
5 r, v9 t3 \# ]3 Rid2 = find(U(2,:) == max(U));
$ r3 j* A0 ~: |id3 = find(U(3,:) == max(U)); # ]" [- ?2 x* G/ g
city(id1) - d' v2 b5 v/ t& @+ \5 Q$ V
city(id2)
1 E/ ?- m; ~- w v) gcity(id3) ' x6 O1 }$ {" @" G7 }7 i1 E& h7 R
/ Y2 ^6 J7 n4 \9 G0 Y为什么一查看聚类结果
Y- F% K0 |, E3 e j: A P就出现错误呢2 j' ]# H- Y' C: r8 Q1 \
请大神指点一下为什么
# e$ O0 V6 I( B' T总提示指数超过矩阵维数
& N$ ^$ p6 p' d& F5 L
# K1 C. x8 @4 ^; b, Q; n7 m* S8 c
: c- U2 z4 F% O' p& V |
zan
|