- 在线时间
- 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 编辑 - F- A) V/ b" I# L! G
8 @, M2 [* H4 m/ G. B[m,n]=size(data(:,:,1));
* t6 F/ _: r. kx=zeros(m*n,3);1 n' o$ o; k* e4 I; I0 u$ x8 _
for i=1:m. S' g, J: ~* a2 E8 h7 O a
for j=1:n! ]/ J. j/ i; {( @* b
x(i,1)=j;
/ g( E/ s$ V* _- O: i- {( y x(i,2)=i;
$ o; J6 n, H& ] x(i,3)=data(i,j,1); ' y$ ?' [2 U$ D% D4 a' w
end
5 t& y3 g* T4 _! T2 {+ lend+ j+ L. b& e! b1 v' ~
X = zscore(xdata);+ {4 C0 k" t, I& |% |& y3 m
" y1 c H: h; G4 Y- ?
; r) X& t( V: M% k; A%****************************进行模糊C均值聚类******************************
, k. D5 u% n4 D. B% 设置幂指数为3,最大迭代次数为200,目标函数的终止容限为1e-6,不显示中间迭代过程
c4 y% {; Q" h' J( L0 Ooptions = [3, 200, 1e-6, 0];
2 o1 F" N( f1 e8 @6 G% 调用fcm函数进行模糊C均值聚类,返回类中心坐标矩阵center,隶属度矩阵U,目标函数值obj_fcn8 E" H: W! C4 T& \7 \- r
[center,U,obj_fcn] = fcm(X,3,options);
. @. \* O# i7 L! E g3 ~0 I: T; N- g7 N5 b' R. ^/ n& {2 C
% b, g" J D3 `( I6 V: g# f%****************查看聚类结果*****************************
' y3 {5 d, H/ ]+ [7 E8 A5 K7 r0 ~id1 = find(U(1,:) == max(U));
9 A7 |! \ n( r* lid2 = find(U(2,:) == max(U)); * e% y+ t, \" _8 S# e: f
id3 = find(U(3,:) == max(U));
0 e" i) t% e! Q" ]city(id1) 6 l1 a( V1 Y8 e5 \8 ?7 ]$ e
city(id2) " h# J# ?3 W& ?. P
city(id3) * D5 `- m( s0 P& b1 o
+ t" P5 O$ A, B! M. l% D x为什么一查看聚类结果; K1 G; d! b& P) a7 D/ M9 G
就出现错误呢
) _8 g1 P" Z6 l- Y$ o. }5 U请大神指点一下为什么
* q9 K( S; e+ s$ M5 d$ U总提示指数超过矩阵维数 i8 s6 E5 y: i$ ]) E. V
4 ?6 j2 m% @3 ^. w& [$ i
: k+ f, d( x4 `/ `; e6 f |
zan
|