- 在线时间
- 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 编辑 * p7 V- s, ^' Z; X& g) c/ W
& E5 L9 i( E u3 P) D% W$ L[m,n]=size(data(:,:,1));
/ W. j% }$ ^& d0 f2 Z) K8 _x=zeros(m*n,3);
& U# @$ f* n4 ~$ Jfor i=1:m0 f* }& s& T, K2 V/ h! d7 S A
for j=1:n
* p7 S( S3 B' g" {( U) H8 q: @ x(i,1)=j;: ?' p/ V" y, @8 s. }$ g# R2 p
x(i,2)=i;
3 }" y5 s5 x7 ~ x(i,3)=data(i,j,1); 0 F) o, t( J0 w: c3 j3 C; l3 n
end
7 ]7 M5 {# b w! y7 f. U0 jend5 _8 ?' I2 o4 N1 A
X = zscore(xdata);
" z/ x# |: f9 Q, f5 o, G! Q
7 N4 ?$ P" Z b
! N7 K& D6 R4 F* y8 j%****************************进行模糊C均值聚类******************************! F. u. V; F1 u' ^
% 设置幂指数为3,最大迭代次数为200,目标函数的终止容限为1e-6,不显示中间迭代过程. A: j! V5 I6 m. c5 }! r6 O
options = [3, 200, 1e-6, 0];
1 W: T+ w5 r- I. L. t/ ?% 调用fcm函数进行模糊C均值聚类,返回类中心坐标矩阵center,隶属度矩阵U,目标函数值obj_fcn
9 O) K3 o+ j" v[center,U,obj_fcn] = fcm(X,3,options);
" S8 m. M0 b j* d( S9 r* S
5 I' R- X4 f8 b. G( G4 A% @( ?2 m( X r+ ?
%****************查看聚类结果*****************************- N& X/ B* ?$ K1 ?
id1 = find(U(1,:) == max(U)); & y: n% L8 q) g9 q
id2 = find(U(2,:) == max(U)); / W3 Y3 R6 A+ ~& } U/ G1 T1 @
id3 = find(U(3,:) == max(U));
) y0 _: a7 ]8 s: N, U3 Pcity(id1) + x) f" I. N( X0 g/ |" A
city(id2)
5 ^( t- Y4 |- C; D( w$ Tcity(id3) , v/ ^# \, @; j3 d H0 \! A
6 ?. |. u- N" {/ `/ g3 p为什么一查看聚类结果+ u7 I2 _( N$ F& |* D
就出现错误呢2 n2 j: W+ D! x8 ^2 B% _4 A
请大神指点一下为什么 _, S+ H) L! x8 r
总提示指数超过矩阵维数
# C- e7 S9 i/ x* w" C. M/ P! n% @" H
: S2 C5 u# ~" \! v R% M
|
zan
|