- 在线时间
- 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 编辑 & a( C: L9 P2 s# W+ X! G) V
1 t: ^+ B9 I. Q0 {
[m,n]=size(data(:,:,1));
( w* G. E/ z$ N5 A) ^6 jx=zeros(m*n,3);
9 g0 Y, J. [8 s# U" Afor i=1:m9 g9 a7 j* d, G
for j=1:n
8 L- q( z$ Z/ n3 x x(i,1)=j;" j0 s. o- s2 ^4 E) I) |- d# N
x(i,2)=i;
. w' U7 {5 e/ i# h# p x(i,3)=data(i,j,1);
& U* q# \$ A: ~ end+ T# }' C8 I) d; W+ G. ^
end+ }* L4 J2 i6 v' I% ^
X = zscore(xdata);) x; r( z7 a, k- |- X
' t0 \* w) w+ o0 |/ ^8 Y4 Q5 u8 c, Y( M$ d7 {
%****************************进行模糊C均值聚类******************************( N. I Y' h6 I) p
% 设置幂指数为3,最大迭代次数为200,目标函数的终止容限为1e-6,不显示中间迭代过程6 u+ k- S0 F" ` F* g
options = [3, 200, 1e-6, 0];
: Z( _9 c$ K M2 d' T: m( M' H% 调用fcm函数进行模糊C均值聚类,返回类中心坐标矩阵center,隶属度矩阵U,目标函数值obj_fcn* p) S2 M. S7 F3 |
[center,U,obj_fcn] = fcm(X,3,options);7 H, z6 E) g( T" ^! a! T# ?% O v
4 ~) K3 v+ A0 |3 ^$ B2 o% a0 X8 g( o/ n% f
%****************查看聚类结果*****************************: `1 h8 s+ K, H$ l/ N
id1 = find(U(1,:) == max(U)); * h' j0 F1 D, ], F* h
id2 = find(U(2,:) == max(U));
: d8 a4 `3 p" M3 y& |id3 = find(U(3,:) == max(U)); " a- q7 V- `* u G( Z3 S1 ^0 p
city(id1) ; T3 k6 {% I/ Q$ Q2 _3 i- F
city(id2)
! |; `3 Z* q$ rcity(id3)
( i( A4 x; \( ?) F2 }/ M
- z8 a- w% A1 s为什么一查看聚类结果- i8 v- k/ V# T% x" L! X, r
就出现错误呢
' j, w+ a; f' I4 b3 {请大神指点一下为什么
. J/ a6 d; f Y总提示指数超过矩阵维数
1 q2 H1 ]2 z( j7 Z7 L3 g* j! u: p* A/ W- r$ ~% e# N
2 C, h8 F$ V4 l' N+ d
|
zan
|