- 在线时间
- 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 编辑 2 r6 k0 s( ?; N+ j2 w/ v
. M! l6 F* |6 M5 _' `7 q6 v[m,n]=size(data(:,:,1));4 b3 L4 \: b8 {$ v; E
x=zeros(m*n,3);
; u- q2 r" S' ^8 x' efor i=1:m5 v5 _: u5 J/ R! x; l
for j=1:n8 ]( S5 l/ H* `. f2 w) E
x(i,1)=j;) g3 o" K! m* k
x(i,2)=i;) [9 p" `6 m& F9 O/ h) o% [/ Y
x(i,3)=data(i,j,1);
7 V( e- T4 \8 w, b+ w) @3 O/ ~$ F end/ ]. w3 i" x) W" H: z6 }
end
: @$ n7 x' r" LX = zscore(xdata);$ i6 S" N% Z' ?% ]% K, b
& D( I! [8 m& U' g
% q, h% p$ y9 P2 q%****************************进行模糊C均值聚类******************************
, \7 }/ ^3 R8 x' I% 设置幂指数为3,最大迭代次数为200,目标函数的终止容限为1e-6,不显示中间迭代过程
0 b$ J$ r, f4 S9 ]1 e) q7 O; ~options = [3, 200, 1e-6, 0];
# @) ?' K0 V* q3 h Z0 V' Z/ P% 调用fcm函数进行模糊C均值聚类,返回类中心坐标矩阵center,隶属度矩阵U,目标函数值obj_fcn2 V U# s3 T. |! T+ p3 k
[center,U,obj_fcn] = fcm(X,3,options);
1 W/ v$ t7 I* u7 y# E C4 M, j( g8 s! q8 P4 |# n, O1 n) ?
& v- k2 C( {5 S( x3 H. N5 E- _% _%****************查看聚类结果****************************** s/ H3 l4 ]3 w2 ^7 G8 A( H
id1 = find(U(1,:) == max(U));
3 L4 m' I: h0 G9 G: {" t# u4 |3 [id2 = find(U(2,:) == max(U));
1 K3 `! L; K1 S3 E) c1 z4 `id3 = find(U(3,:) == max(U));
) n5 W0 b- z2 t$ b8 x7 G# u+ Q4 Gcity(id1) 0 r0 E* I) s9 U6 g
city(id2)
& N+ S1 \' ^2 c; Dcity(id3) 8 w( v+ p- c8 X7 W/ E: l
r! n' J8 H; c" J! X为什么一查看聚类结果; k" ?7 O4 v" H
就出现错误呢; [* ?! W8 ` A1 \
请大神指点一下为什么4 ]+ Z" [& {4 k# [' W/ G# n
总提示指数超过矩阵维数' x0 z& ^9 z* [- e( T `+ l, ~& s
3 X* [$ D" o7 J; \9 T' X+ ^; P! h- v5 i) ?; Q9 m
|
zan
|