数学建模社区-数学中国
标题:
模糊C均值聚类
[打印本页]
作者:
509
时间:
2013-8-20 16:42
标题:
模糊C均值聚类
[m,n]=size(data(:,:,1));
' Z+ b* g8 }! a1 E- [" D) x
x=zeros(m*n,3);
& O t3 ?4 P8 n- S6 o4 Y
for i=1:m
) |5 p1 ]9 ]% h" I' G
for j=1:n
2 X" {9 x# s) Q
x(i,1)=j;
$ l4 o+ n7 N; `- c
x(i,2)=i;
; `2 o$ X' L5 l: B4 W. c' `
x(i,3)=data(i,j,1);
/ h5 r' V) Q* f3 k: x
end
% j$ O3 W* A2 v1 B
end
+ x6 ]+ I8 J7 B1 j9 D
X = zscore(xdata);
/ D! ?2 n! N" q* r, g! m
~, L) C" \1 N- E) Y% V& d
' f: R5 ^/ w1 n7 @, h1 ?
%****************************进行模糊C均值聚类******************************
\. S. P* }2 f, n
% 设置幂指数为3，最大迭代次数为200，目标函数的终止容限为1e-6，不显示中间迭代过程
" N: o) H! `6 \6 a' d
options = [3, 200, 1e-6, 0];
' u, C$ o+ L9 u( z; P/ N2 Q* g! b
% 调用fcm函数进行模糊C均值聚类，返回类中心坐标矩阵center，隶属度矩阵U，目标函数值obj_fcn
; G ?4 F- X6 u5 r, a) Q9 k6 g
[center,U,obj_fcn] = fcm(X,3,options);
& H9 A' N/ o* D! @
% ]4 L5 f% E: f" p# B+ `- T
4 v; W2 g" F" [$ N. } V" e5 S0 M
%****************查看聚类结果*****************************
: M4 l' X. Z, w( j
id1 = find(U(1,:) == max(U));
# u+ P' \- w* r: m
id2 = find(U(2,:) == max(U));
6 d, l1 Q8 ]/ f; {0 H+ ]# x# V
id3 = find(U(3,:) == max(U));
+ X3 G: O3 E8 q2 {- k* D2 _0 I
city(id1)
# P: r# k8 y. z) K& x: E: A
city(id2)
0 t: F# P) a2 ]: J! \3 L
city(id3)
2 [/ j5 W$ i9 \8 ~+ @
# C0 ~0 ?, E" }3 q% T i' q
为什么一查看聚类结果
2 V7 ]8 J! T# u# J8 m
就出现错误呢
, y5 G- f) c9 ?% @
请大神指点一下为什么
n/ w4 y4 z {9 E
总提示指数超过矩阵维数
0 b3 Y1 G; X' q! z7 N
7 u3 Y/ S$ Z& t& C+ s" X+ @4 I6 ~
作者:
我爱照镜子
时间:
2013-8-20 16:54
那三个表情是什么啊
作者:
净心、精心
时间:
2013-8-20 18:13
我爱照镜子 发表于 2013-8-20 16:54
' S+ f* b( k# z
那三个表情是什么啊
* J. f7 {- B: k# J# s' O
：：：：：：：：：
欢迎光临 数学建模社区-数学中国 (http://www.madio.net/)
Powered by Discuz! X2.5
' S+ f* b( k# z