- 在线时间
- 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));
# z. H% h+ j9 |- R9 [. h& ax=zeros(m*n,3);
3 g ?+ L+ C# cfor i=1:m1 n' k8 C& ^1 G. K8 ~5 W
for j=1:n
3 U4 s, o, M) @& V x(i,1)=j;
3 k, v, }6 z/ S0 w! o: s x(i,2)=i;
1 K' H. e7 M! X; k: b x(i,3)=data(i,j,1); + E( E1 @# ^5 i4 h/ z W3 M* P
end. [, w$ f7 k6 h- d
end
$ Y- ~: j* {) ?4 JX = zscore(xdata);& G2 c4 {4 O3 \/ Z
. \: |/ L$ x4 j4 G; Y) @: ]3 \
9 n% v" v( D# n. E%****************************进行模糊C均值聚类******************************6 y, }+ n% s4 O5 A" {* _& Z
% 设置幂指数为3,最大迭代次数为200,目标函数的终止容限为1e-6,不显示中间迭代过程
/ x- c- g/ E8 z+ _* D$ V9 voptions = [3, 200, 1e-6, 0];7 `+ O2 m/ F5 c4 j8 `
% 调用fcm函数进行模糊C均值聚类,返回类中心坐标矩阵center,隶属度矩阵U,目标函数值obj_fcn
# P3 r, o( \( P1 J0 ~% w$ V+ z[center,U,obj_fcn] = fcm(X,3,options);& G" `! V4 C! x) j' l9 }
" _/ p& }$ m) G2 t: W
( c/ N" U |# }- D7 C
%****************查看聚类结果*****************************$ |# l+ G+ k0 W$ }" W2 P
id1 = find(U(1,:) == max(U));
: n& ^, [! |; a& ^5 _" O k) F/ kid2 = find(U(2,:) == max(U)); $ {8 G, d* m+ l
id3 = find(U(3,:) == max(U));
9 E2 B% \* z0 L; {6 G% jcity(id1) ' N! h% Y9 s! q5 F, ^
city(id2)
8 b( }# |9 a+ t$ F5 @; pcity(id3)
) e: d9 W1 U* ^
9 C+ U: c7 [' X" z. E/ q. E为什么一查看聚类结果
, C9 c+ b7 J+ E就出现错误呢
7 u. Z% d6 p& O( v3 G2 Y' K, a/ ^请大神指点一下为什么4 {9 m8 K1 I6 F2 f
总提示指数超过矩阵维数
. I) E3 q( P0 R: w: S, a* n6 X
+ M6 E @# t& O7 r; T' N; _+ s |
zan
|