- 在线时间
- 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 编辑
3 [3 ^6 U3 K9 @+ F5 d; l- J" U: I
2 m1 w2 W% w# `8 ^[m,n]=size(data(:,:,1));/ X) O6 L' @" P* `+ u# c
x=zeros(m*n,3);6 S1 D: j! J' T% f7 R# i6 N8 W2 c
for i=1:m
( V/ Q: i2 ]1 U3 `8 X: E, e& G for j=1:n" R. G4 D; ?7 h' s$ I. _
x(i,1)=j;
9 A# y" D9 p' s6 ^1 {# _+ U x(i,2)=i;
, r' U3 T% J' f9 B9 ? x(i,3)=data(i,j,1);
2 T& ]' _7 D _; U9 @ end
$ z! k( b7 p* Qend5 B, l% ?1 J2 q$ L# b
X = zscore(xdata);( d; K; u2 {; p# y" q
$ W" M+ J$ k( p- G9 a0 X& a. j
' I" Y, `3 U, d0 l%****************************进行模糊C均值聚类******************************* n. z, J; I+ ? `9 i3 u# c% v
% 设置幂指数为3,最大迭代次数为200,目标函数的终止容限为1e-6,不显示中间迭代过程* Z r9 W3 |! R' P, h0 i- C) e
options = [3, 200, 1e-6, 0];
1 V. i( _" E3 ~( b8 t/ k% V1 Z% 调用fcm函数进行模糊C均值聚类,返回类中心坐标矩阵center,隶属度矩阵U,目标函数值obj_fcn& s2 v' _- w7 w
[center,U,obj_fcn] = fcm(X,3,options);% O5 A$ T5 ?' u; l) ]
' a2 o) M& M: w9 l, f- p+ l2 R! l a# ^$ D, _, v2 I! p
%****************查看聚类结果*****************************1 z* W! x5 @/ e4 e
id1 = find(U(1,:) == max(U)); 5 f# ?5 _: n. H: E+ a
id2 = find(U(2,:) == max(U));
$ ~. |' p' w a0 Q, c! y& B! ]( h( Hid3 = find(U(3,:) == max(U));
4 I0 P+ y$ c F8 n* U0 mcity(id1) 3 y) t2 U/ k$ v9 G% M
city(id2) / d8 B" Z" X8 b
city(id3)
. w3 B$ _1 H! A8 q! v
- E j4 e: C9 [% n6 m3 S, I为什么一查看聚类结果/ S; j9 M/ S$ y& y7 G* B
就出现错误呢4 Y& `. U6 {& Y) ?
请大神指点一下为什么& S w' W' c/ z; C5 \
总提示指数超过矩阵维数, H! m9 M9 e# J* R( @, J
* Y$ r/ V& @! V! p' Y# i+ `. n
( g5 o9 d$ S3 w. W2 v X
|
zan
|