- 在线时间
- 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 编辑 1 r4 o: W: s$ G
, A. L: P3 s8 i4 `[m,n]=size(data(:,:,1));
) X5 T5 F$ ^3 ^7 Ex=zeros(m*n,3);
' a2 Z$ K! i9 G( l3 Hfor i=1:m
* D6 P! @# M; M0 K. P for j=1:n
3 ]6 s1 ^0 H' d4 c4 P# s: P' A x(i,1)=j;
x Z2 j2 d2 l! {% n" R x(i,2)=i;
3 o# c& I! |+ N' n S5 q: R x(i,3)=data(i,j,1);
2 Z/ i) U6 F8 \# s0 w8 O end
' t; o' u; ~& |, bend# D9 F r# j9 M0 g
X = zscore(xdata);
5 w1 {- z5 Z& i" E
$ e. K6 F% x, ^, c
! A! r0 W* Q9 B$ j% C* m" X%****************************进行模糊C均值聚类******************************: E9 S+ S5 G3 [ H0 e$ c
% 设置幂指数为3,最大迭代次数为200,目标函数的终止容限为1e-6,不显示中间迭代过程3 H, Y+ U, m L
options = [3, 200, 1e-6, 0];2 n9 u8 H3 x5 P3 F8 r; G* J
% 调用fcm函数进行模糊C均值聚类,返回类中心坐标矩阵center,隶属度矩阵U,目标函数值obj_fcn0 k% J6 P3 u: E, D( _4 W+ {- s
[center,U,obj_fcn] = fcm(X,3,options);; R) v |0 V/ O* t
4 w+ l# x# V3 L0 j6 i* R) F7 z( l3 R V" d" ]) d; T, _
%****************查看聚类结果*****************************# K* U2 B5 Q2 n9 V: q" ~ I
id1 = find(U(1,:) == max(U)); " l- y5 v" Q5 i: N9 T
id2 = find(U(2,:) == max(U)); & X% r5 B) i, E4 Q, `1 B
id3 = find(U(3,:) == max(U)); % m& Z; F! J# ]- Y) G( |
city(id1)
; g. Y j W2 T, V# \ b: tcity(id2) ' n! @. x! v9 _) o* s) ?9 ^; A9 Q
city(id3) ( d# z4 X) u3 n% ?+ @
& U5 O$ f: I8 y5 m% N
为什么一查看聚类结果: T' Q$ [) X( y5 s% P
就出现错误呢" `1 q4 u. T: ?. z
请大神指点一下为什么" H! Z# k0 C9 |3 [
总提示指数超过矩阵维数
s, u; V" I9 X, u( t' {
3 b$ Z2 C/ d; ^* h, b1 X( V) ]; }; A* ^: k, t
|
zan
|