- 在线时间
- 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 编辑 ! L- X, V" w) Q) \! \
+ ]0 f# {0 u5 ^- r B7 V2 l5 }[m,n]=size(data(:,:,1));
" L+ {. J( ^) e. @- j+ z6 i7 c* gx=zeros(m*n,3); S, D% w2 D$ p" b9 ` S3 e
for i=1:m& {* b0 t3 u$ t' M
for j=1:n
' A( `1 t9 |( K V x(i,1)=j;: i+ q9 \ j4 U( H
x(i,2)=i; p- a6 f k' c4 f* q8 ^
x(i,3)=data(i,j,1);
/ [, K! z3 v2 J3 w9 C w end
) K3 q) ^9 @5 f: ?# X/ t3 eend- [% |6 A ?7 F1 g4 @
X = zscore(xdata);
7 m* R' p; e/ n7 P
, D; z$ |7 X, @3 [) l) O3 p$ q( L6 T+ I
%****************************进行模糊C均值聚类******************************
: K* ?9 b; G# ]+ f8 P: s4 @% 设置幂指数为3,最大迭代次数为200,目标函数的终止容限为1e-6,不显示中间迭代过程% a6 ]% b5 s# |' r w3 ?
options = [3, 200, 1e-6, 0]; b" t4 s6 |" Y4 T5 P$ I9 [
% 调用fcm函数进行模糊C均值聚类,返回类中心坐标矩阵center,隶属度矩阵U,目标函数值obj_fcn5 H/ l$ n" U* ] A
[center,U,obj_fcn] = fcm(X,3,options);
2 U. [- `. S* T" Z7 c& ?
' P& |3 x' j9 X: J: w e+ ~- z8 S' m# t
%****************查看聚类结果*****************************4 ?, |4 M: k( W
id1 = find(U(1,:) == max(U)); , ^7 D8 q' ?9 ?! s9 n; a
id2 = find(U(2,:) == max(U));
! n: |9 f1 O' l( s0 C5 ?id3 = find(U(3,:) == max(U)); 0 P8 a1 e% h* D! s8 t1 \3 b
city(id1) - O M; ?7 ?7 [6 P( |# E o$ k
city(id2) ! ~9 v* R. R \
city(id3) 6 t: [$ ^- p) f3 o: c; V
9 z2 Y! w9 H( d& Y' i为什么一查看聚类结果
1 K4 w" c% U. @. E就出现错误呢3 h* ]& X2 a k2 y9 N! v
请大神指点一下为什么
- S) v2 g. X( Y! S+ M总提示指数超过矩阵维数( U9 b4 e; S: K$ F$ _: t2 N& l7 [
0 ^8 `, a) X2 K
7 T$ u# I+ |* F+ g8 I |
zan
|