- 在线时间
- 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 编辑
5 b, J+ ?! x7 A$ N3 ^2 c
$ G! S. \3 @* x! n- A& p9 P[m,n]=size(data(:,:,1));
# [7 ] y7 v+ Z2 f6 D8 _x=zeros(m*n,3);, U- P1 P! |/ F3 ?+ z- P
for i=1:m
5 c7 k- l, Y% n# e: V for j=1:n
6 Q$ f3 ?8 T! L- F- |2 }1 N5 D x(i,1)=j;
3 j4 F- x, y' f; p+ g* O1 N x(i,2)=i;
) }( T) P) A4 I0 s+ L x(i,3)=data(i,j,1); . \+ Y5 N. B% ]
end6 e6 K/ T6 y$ ]- X8 |" }. P
end) `! u) g2 Q1 O4 ^5 j$ V
X = zscore(xdata);/ j' T' j3 m0 C
7 _' W! C* r) k" w* \* o- u+ O
2 m4 `# ]" l6 o2 O/ u%****************************进行模糊C均值聚类******************************5 I8 A& {: l# Q4 q+ k
% 设置幂指数为3,最大迭代次数为200,目标函数的终止容限为1e-6,不显示中间迭代过程! } C B1 A* S0 J; _! {- x8 d
options = [3, 200, 1e-6, 0];
5 ] g$ t5 W& [$ L' A( p% 调用fcm函数进行模糊C均值聚类,返回类中心坐标矩阵center,隶属度矩阵U,目标函数值obj_fcn
6 Q3 i: p7 U& E4 ]2 [0 y' T9 F[center,U,obj_fcn] = fcm(X,3,options);+ W d x y7 v3 z' |
' p: a% o( \' i [. [ Y
0 ?/ [6 p( P T. Q5 d4 B%****************查看聚类结果*****************************$ O+ @2 i$ A6 j& ]! V0 o
id1 = find(U(1,:) == max(U));
( }6 C- { Z" F# U! i' Hid2 = find(U(2,:) == max(U));
' }( |& B* l; X8 ^' N: d0 n! yid3 = find(U(3,:) == max(U));
8 B. u2 L/ Y+ {4 Y' mcity(id1) , G3 K ^8 I5 w" I4 @ }5 x
city(id2)
* L! V% l! w9 f9 C. N0 Z# Fcity(id3) # \' z2 \! B0 C w; b1 R5 s V
- [% n; q$ H$ v1 J为什么一查看聚类结果9 q4 o- j& p; x% f; {
就出现错误呢
R! q% |! E, P' }$ [/ c请大神指点一下为什么
( w: d* J5 Q" I3 K1 C0 [总提示指数超过矩阵维数
( p- N/ `: b1 }' x
9 b2 G4 c: B" D$ a* W5 c1 _5 g4 j8 `) l4 C5 s
|
zan
|