- 在线时间
- 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 编辑
. a5 A; n' D6 C. m% Q4 u# b# ]2 J% @( u+ v0 E# x4 M8 D
[m,n]=size(data(:,:,1));
9 `4 Y3 r) I5 j Y4 [ B1 F# V5 }1 `x=zeros(m*n,3);% ]+ n! A) `# Y. X |3 x; ?& r' W
for i=1:m
* r U1 G$ {; R, F! H for j=1:n7 X# d, V, S/ E5 g1 q- l# `/ h7 |* s5 |4 O
x(i,1)=j;7 W9 r1 X+ N& q) Z
x(i,2)=i;
. s+ {+ ?3 E" Q9 z' C0 _# k6 S# m x(i,3)=data(i,j,1);
0 ?* `+ f5 Z! g4 N, a- B; U7 b" U5 p end
9 G: E Q( I) `end
. I) w# ?7 f! f9 A" WX = zscore(xdata);" x s" V0 o( Z! V
1 ]( f3 X7 t \5 O% A
, S# x) k% m0 o& R2 T1 @$ _8 A
%****************************进行模糊C均值聚类******************************9 p" ]6 R7 j7 K% ~
% 设置幂指数为3,最大迭代次数为200,目标函数的终止容限为1e-6,不显示中间迭代过程) ~; c# q9 s) s: v, S+ }
options = [3, 200, 1e-6, 0];7 B% u6 Q! i4 B% T4 l
% 调用fcm函数进行模糊C均值聚类,返回类中心坐标矩阵center,隶属度矩阵U,目标函数值obj_fcn) f% j, j' _! {- Y; L9 m
[center,U,obj_fcn] = fcm(X,3,options);
7 n3 b# ]4 c: I* }; f M2 l6 }3 Q$ j
% o! M; s: {' F! Q2 A; N5 }" E4 E# ^0 Y" {
%****************查看聚类结果*****************************
2 m ~! y9 f, q6 X2 H W+ W) xid1 = find(U(1,:) == max(U)); 6 z+ w, u6 u# `% \
id2 = find(U(2,:) == max(U));
9 h) J% m) u, }id3 = find(U(3,:) == max(U));
/ V% T9 H) X% x4 [1 O! q4 ccity(id1) : l2 I! F7 v. H9 P& A
city(id2) * N1 V: m/ Q/ _* b! M" v& k# Y. P% P6 x
city(id3)
6 m2 u4 i# t @6 V- r1 g! f
0 Y7 H$ U S+ ?4 h9 E4 Z V为什么一查看聚类结果7 M; f; \! K7 x* a' N6 t
就出现错误呢- M K- {! E, g& X1 {; T [: F
请大神指点一下为什么3 M$ e8 P0 b0 x
总提示指数超过矩阵维数
& T# Y5 [, E: A" b) l/ P, h6 S8 @$ a% c% S) h6 B
# V- G2 Q8 q3 r2 S
|
zan
|