数学建模社区-数学中国
标题:
模糊C均值聚类
[打印本页]
作者:
509
时间:
2013-8-20 16:42
标题:
模糊C均值聚类
[m,n]=size(data(:,:,1));
/ ^* w% w& C" g P; o, F" r
x=zeros(m*n,3);
& b4 P1 J1 z. J0 q* q; ?( ^; e5 ~+ G
for i=1:m
! k) ^$ f; e7 |" F
for j=1:n
6 ^( o4 Z0 J1 m4 e
x(i,1)=j;
3 w/ @ x- x5 F, t. C2 F
x(i,2)=i;
7 b1 d. v. {) m" E. n/ J
x(i,3)=data(i,j,1);
, R6 o" r2 k& e) ]) l5 h# ^
end
5 w# @& Z J7 Z$ T
end
4 D9 a# a$ Y% t2 c
X = zscore(xdata);
" @3 M1 L$ ]2 o- K
" N2 \0 U4 T" V; ]4 M
) l% N7 G' {" \
%****************************进行模糊C均值聚类******************************
1 b5 Z* k2 m( _: V/ n4 u
% 设置幂指数为3，最大迭代次数为200，目标函数的终止容限为1e-6，不显示中间迭代过程
3 Y% G( N! p6 O4 w
options = [3, 200, 1e-6, 0];
! \: Q0 I& m1 d
% 调用fcm函数进行模糊C均值聚类，返回类中心坐标矩阵center，隶属度矩阵U，目标函数值obj_fcn
8 ?. @! q: o, v& G& L, R' T+ \( B
[center,U,obj_fcn] = fcm(X,3,options);
$ Y: s. ^' A9 o! U
! O& X6 X0 Z( J2 b# J( ^! r
7 B; k- ~5 z! G9 r' ^' z
%****************查看聚类结果*****************************
" M; p* @5 C- H9 B
id1 = find(U(1,:) == max(U));
. z/ R7 m0 k/ I9 F! V0 N3 b
id2 = find(U(2,:) == max(U));
- q% \8 t5 m+ X! X6 ^ L
id3 = find(U(3,:) == max(U));
' N/ J7 U7 O9 I; p" w# f
city(id1)
5 L: u5 p! y3 r! a2 [
city(id2)
: h) i5 K) B/ ?7 \3 d/ N( h. y
city(id3)
9 f( |( \8 I# o; |% ?
/ h X# B- C1 r2 W) f
为什么一查看聚类结果
- C/ s3 \7 E+ m; S0 V! M" E K% M; n
就出现错误呢
& I3 o: E8 z3 H# f6 {
请大神指点一下为什么
% G5 i0 p6 i! U7 A) }$ j9 q
总提示指数超过矩阵维数
" z( | ]; h! Q7 V- M, N! g
2 b! p6 R; F( n+ S
作者:
我爱照镜子
时间:
2013-8-20 16:54
那三个表情是什么啊
作者:
净心、精心
时间:
2013-8-20 18:13
我爱照镜子 发表于 2013-8-20 16:54
; J7 E; X3 \' N' j$ _
那三个表情是什么啊
5 h: G. }# i, d7 g% M, {
：：：：：：：：：
欢迎光临 数学建模社区-数学中国 (http://www.madio.net/)
Powered by Discuz! X2.5
; J7 E; X3 \' N' j$ _