- 在线时间
- 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));, x' s& F* ^6 F% F* w
x=zeros(m*n,3);
$ p M6 a& q5 L% }) L1 Gfor i=1:m" e P& s5 ]3 c; D
for j=1:n
0 s8 C* ^1 [ f, s% `, { x(i,1)=j;- X' |- `# S6 O3 }! l# P4 q6 ^
x(i,2)=i;
. n4 ?% h3 C0 m9 e x(i,3)=data(i,j,1);
+ f4 b4 T% T ^) Y end
) y( v/ P! z5 ~end1 R; ~1 Y3 ~! f/ e: |
X = zscore(xdata);
, ]4 B6 U; q2 Q% E) n$ E2 D) ~9 E$ n# v1 u' f: y3 S
" y# t. C5 y6 m, N! ~
%****************************进行模糊C均值聚类******************************8 c/ O# u* c- q" `+ b$ h2 z* E
% 设置幂指数为3,最大迭代次数为200,目标函数的终止容限为1e-6,不显示中间迭代过程" o d: P" \. ?9 k6 I, U8 w
options = [3, 200, 1e-6, 0];; \& Y# A* M2 Y7 V
% 调用fcm函数进行模糊C均值聚类,返回类中心坐标矩阵center,隶属度矩阵U,目标函数值obj_fcn8 q6 a5 }0 t9 Y* D; F. a) m
[center,U,obj_fcn] = fcm(X,3,options);. Q1 o |" u- ]# l7 h& X
9 W8 _7 v5 h) M# b y9 w3 Y
3 H& M% S1 J5 C/ o%****************查看聚类结果*****************************0 h3 c" I( o: [; g
id1 = find(U(1,:) == max(U));
+ v: y8 ], o @1 Iid2 = find(U(2,:) == max(U));
0 q) |$ P, X0 n4 m) [id3 = find(U(3,:) == max(U)); , J& l8 w* m; K, v4 ?
city(id1)
2 n7 |( X! @5 Hcity(id2) 9 q7 Z# s; e+ B$ h4 e2 R$ G, k5 N
city(id3) ; N( C: E0 @! ?" t
! g. F/ _1 G( ]! c. |
为什么一查看聚类结果$ d9 V C- k! @! L2 q
就出现错误呢
: e8 E0 u5 e4 y8 i- l9 |7 e* l0 o请大神指点一下为什么 s2 E' o1 n q w5 }2 _6 X( L
总提示指数超过矩阵维数5 ~4 J) h! p6 J0 }0 }5 E
% K/ r& n" ]. t9 W |
zan
|