数学建模社区-数学中国
标题:
模糊C均值聚类
[打印本页]
作者:
509
时间:
2013-8-20 16:42
标题:
模糊C均值聚类
[m,n]=size(data(:,:,1));
q! E) A6 X9 O2 [- Y- U( T
x=zeros(m*n,3);
( r7 L* h9 I0 d* {) j
for i=1:m
9 c- z% k6 t4 r8 X) U3 r
for j=1:n
+ j5 w' i' \% |3 A* X2 D
x(i,1)=j;
0 z8 Z4 \7 {/ [% H* C9 B
x(i,2)=i;
& c0 c! v$ P0 b( A% @1 ]# K0 w+ {
x(i,3)=data(i,j,1);
% e1 h& _+ }$ v4 N% G
end
) r6 L+ k+ h! S1 t) T
end
8 e0 [. e% g( j0 j3 _- k9 T- O
X = zscore(xdata);
p* a; C: l; ?2 Z& m, x2 j
8 n7 B& B4 ]0 Q7 I% M
- R/ P/ w# h6 m$ V% x
%****************************进行模糊C均值聚类******************************
- x5 d: f& L2 s. {
% 设置幂指数为3，最大迭代次数为200，目标函数的终止容限为1e-6，不显示中间迭代过程
) C6 k* H- U( D: E7 g1 K$ t* o" ?
options = [3, 200, 1e-6, 0];
" w" F' A1 V* Q T1 P: r
% 调用fcm函数进行模糊C均值聚类，返回类中心坐标矩阵center，隶属度矩阵U，目标函数值obj_fcn
6 Z# D, o8 `0 b+ q! H
[center,U,obj_fcn] = fcm(X,3,options);
+ |2 ?: b& z) ~5 l# D3 e
8 P) ?" r/ u- x8 S4 a: G
; I' z6 P& W& G" D6 ]: D( }
%****************查看聚类结果*****************************
( ~$ O* y3 ^ s
id1 = find(U(1,:) == max(U));
W# H1 p* K& o( W" J+ m* Y5 Z
id2 = find(U(2,:) == max(U));
. J* a6 B1 T' a
id3 = find(U(3,:) == max(U));
( L6 j7 \9 z& [; k7 I9 g' Y
city(id1)
1 I6 M. l6 R5 t
city(id2)
- y& W! `: `( q- A8 @
city(id3)
2 r9 f) q9 m! i4 Z: Z2 B4 u4 Q
& j4 f5 K) e, _! Y9 k8 R7 k0 I
为什么一查看聚类结果
4 _4 e# n7 I- a* b P! Q
就出现错误呢
- \; E$ J. i& Y
请大神指点一下为什么
" g( Q5 j; O$ R* z
总提示指数超过矩阵维数
! z! C6 _3 C, P$ ?8 ~" @
. d/ e& o, V6 n) N' p
作者:
我爱照镜子
时间:
2013-8-20 16:54
那三个表情是什么啊
作者:
净心、精心
时间:
2013-8-20 18:13
我爱照镜子 发表于 2013-8-20 16:54
4 W z! `0 M: \& v5 U
那三个表情是什么啊
) v- }, s# ~) g u' u( n5 v
：：：：：：：：：
欢迎光临 数学建模社区-数学中国 (http://www.madio.net/)
Powered by Discuz! X2.5
4 W z! `0 M: \& v5 U