- 在线时间
- 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));
9 M) B9 }/ n% N, P2 fx=zeros(m*n,3);
# @* G: _9 X: Lfor i=1:m% U9 h& n, A& Q$ F2 `5 R) t7 `
for j=1:n
% k* Y4 C0 ~% Y5 r4 u* G( E x(i,1)=j;
6 R3 m0 b7 `5 ]9 Q9 M x(i,2)=i;
* k( U& }# F. o9 f2 q# a: m x(i,3)=data(i,j,1);
9 u! [1 K8 N# v; e end/ m! ?. o0 g1 M! M0 ]
end- x1 Q" G3 O+ e& C" S4 v, ~
X = zscore(xdata);) `! m' {- k: i5 J G# l+ P& c
4 V5 c0 P$ x7 r& j+ V* i3 D* M+ Z5 M9 H# d0 f
%****************************进行模糊C均值聚类******************************
+ |; C5 Y) R# M, i% y% 设置幂指数为3,最大迭代次数为200,目标函数的终止容限为1e-6,不显示中间迭代过程( d/ X8 G2 E$ p! v6 w0 d5 ]0 r
options = [3, 200, 1e-6, 0];* l: m" e& l0 B7 d/ w! `
% 调用fcm函数进行模糊C均值聚类,返回类中心坐标矩阵center,隶属度矩阵U,目标函数值obj_fcn$ u. p. t R* H( q u3 w4 o% v0 k
[center,U,obj_fcn] = fcm(X,3,options);
; Y, b, C2 C& |$ z1 r+ c" y
3 S( n3 a& R8 \+ X- T% e7 K( N( {- E/ E4 e
%****************查看聚类结果*****************************
- f2 C. b# w1 ^% Nid1 = find(U(1,:) == max(U)); 6 T7 K* C) b% K( x/ }% S0 s2 M
id2 = find(U(2,:) == max(U)); / Q: K! Y" t0 a2 \0 v
id3 = find(U(3,:) == max(U)); # a' u4 }! y5 ]7 H4 a
city(id1) 1 j& _1 P7 h) E- k
city(id2) / ]5 Z5 o8 A: E3 f8 U( x( O1 o; h
city(id3)
6 _# q$ P5 u, {: ~! t
3 p# j8 ]* E0 S/ Y6 t6 M为什么一查看聚类结果' V; u" K0 [, s. `
就出现错误呢+ h I6 C2 Z( A$ p k
请大神指点一下为什么
; N) q& j6 e3 i- C3 M# b' i. P总提示指数超过矩阵维数7 |- K; \4 X4 o. M4 P0 N+ b! K
# l) `! X# Y# x6 Z% i |
zan
|