- 在线时间
- 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 编辑 4 @; L5 W, g( L$ Z, ?+ r! `3 w" A5 k
3 w6 q' h/ {- z. n[m,n]=size(data(:,:,1));$ W! A0 l1 d) G; Z
x=zeros(m*n,3);
; N8 L* S! [5 Q! f+ I: Hfor i=1:m2 v0 M& f- n! ^) u5 D/ M8 d, R
for j=1:n$ d4 N) b$ Y6 u( l; B/ W
x(i,1)=j;
9 B; Z, K9 ^, T* b x(i,2)=i;
2 p/ s" Z& S& `# P- G/ X4 P x(i,3)=data(i,j,1); 2 G/ _- j" j8 d7 ~0 K
end
9 t- T# A$ f l; y6 Z: yend
$ I ^# Z+ P9 c* w4 O" X" D( {. a: PX = zscore(xdata);
- X1 n# s6 @) z% u$ d
$ a* h+ N# u* G' E* @- l _3 F0 X3 |, g. [
%****************************进行模糊C均值聚类******************************
5 G& Y' {+ M6 R4 x6 i$ X% 设置幂指数为3,最大迭代次数为200,目标函数的终止容限为1e-6,不显示中间迭代过程 d9 k% o3 f S, I
options = [3, 200, 1e-6, 0];
8 o9 y# V4 G# f4 d' k6 ^% 调用fcm函数进行模糊C均值聚类,返回类中心坐标矩阵center,隶属度矩阵U,目标函数值obj_fcn
7 l& @; p* y& D* T* _[center,U,obj_fcn] = fcm(X,3,options);
+ [$ n% O! b7 k( J8 y
) z( D8 r3 I9 j& ~4 Z" n/ @; j0 P/ w* Z7 i. f
%****************查看聚类结果*****************************5 m8 p2 z% D2 v2 C
id1 = find(U(1,:) == max(U)); 2 t9 D. O$ Z( r* m& d2 o& S) C1 B# ?
id2 = find(U(2,:) == max(U)); ' M; J" D. _1 r5 ?1 W7 L2 c$ y
id3 = find(U(3,:) == max(U));
3 r( G+ i8 e% U1 W) dcity(id1) , O! T9 U' Q [! r
city(id2) 5 X, r3 f& T& c
city(id3) * q5 H* X0 a7 z! H
6 c" n+ N Q5 g" d# `% x" k5 }: {2 G
为什么一查看聚类结果: U" u0 A% y& |' @/ r0 j7 S
就出现错误呢+ E$ ^3 ]# o* }! [8 ~
请大神指点一下为什么
" [$ g! }+ m" ~7 {总提示指数超过矩阵维数
8 E6 Q5 `6 R2 h4 o4 s- l8 K7 [& S7 M( s" c) A
2 Z- y9 \. F/ _! N+ n( V
|
zan
|