- 在线时间
- 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
- 自我介绍
- 东北农大本科生
 |
[m,n]=size(data(:,:,1));
) G3 u" E9 w, L; Y8 x$ i4 Ax=zeros(m*n,3);$ v3 e, t% A) \1 z) g9 R
for i=1:m4 U+ h8 X% g7 Q% ~ Q
for j=1:n6 s! o3 u1 J7 }+ G3 h
x(i,1)=j;
, w4 p, T7 F3 m+ U5 { x(i,2)=i;0 U$ b* F" p- Y( p
x(i,3)=data(i,j,1); . P( |! y4 @ @7 T1 ]/ u! q( ~
end
/ z+ d2 K( v& c+ oend" o/ p* a! ]( z% N
X = zscore(xdata);0 o( J W" w+ Q: Q' l
+ f2 r! I! z9 P7 {3 l) c
! C, L T7 y, s8 }%****************************进行模糊C均值聚类******************************
8 z3 i2 `6 t8 c0 R# P7 r0 ^5 G% 设置幂指数为3,最大迭代次数为200,目标函数的终止容限为1e-6,不显示中间迭代过程
B" Q7 y! w% r- Moptions = [3, 200, 1e-6, 0];
; t! e) [2 _+ k$ {% 调用fcm函数进行模糊C均值聚类,返回类中心坐标矩阵center,隶属度矩阵U,目标函数值obj_fcn
7 }0 I; G$ Z& l1 r3 x$ q[center,U,obj_fcn] = fcm(X,3,options);9 O: T( G8 S N5 b; O. m
8 C+ p* G/ l/ Q. c4 K5 x' C0 Y2 f/ F4 {. C
%****************查看聚类结果*****************************
$ `, o, h1 ]) r1 A2 Lid1 = find(U(1,:) == max(U)); , k4 v9 }! q1 d* P# U8 t
id2 = find(U(2,:) == max(U));
9 b. G: h0 a/ k& W* }+ Vid3 = find(U(3,:) == max(U));
1 X [& D9 x0 @2 r# t; acity(id1)
' x' l7 S) V, M! k! p: [% x0 |city(id2) + N* ]7 A* }1 Z/ k; {( g
city(id3)
7 O; {) i& R1 x3 X& d4 |- ~3 C+ _% d x u9 Z# L6 D1 p) u
为什么一查看聚类结果
3 z% l% e- n" G9 |就出现错误呢
7 X7 R- S5 Y$ B请大神指点一下为什么
7 [9 h" Y6 r4 [: f; N总提示指数超过矩阵维数5 k5 K. n: |+ T+ j: M: H
) h6 R# `& @' p" m" } |
zan
|