数学建模社区-数学中国
标题:
在线等回复
[打印本页]
作者:
509
时间:
2013-8-20 15:55
标题:
在线等回复
本帖最后由 袁海亮 于 2013-8-20 16:15 编辑
~8 w! M0 L; ]7 L! F
$ V. R+ O* e! Z
[m,n]=size(data(:,:,1));
$ l" M8 J1 _* K& Q' E
x=zeros(m*n,3);
I3 j* k' F. n E& |1 i# l. p2 h
for i=1:m
: _. ^* ?4 w5 z0 i
for j=1:n
% A: Q& w7 O% t! h
x(i,1)=j;
* [0 V$ E0 P; ?0 k& b# h! L7 l
x(i,2)=i;
" g* Q M3 v8 J- J( x. r! M$ V, U
x(i,3)=data(i,j,1);
7 S" A6 g/ I" n. w; `9 B/ k# U
end
0 m- {7 z1 j# z- m5 J# N
end
; h. `, L! ?" A) w0 X# }+ D
X = zscore(xdata);
& b, o9 _1 O- G; a& ~
" ^& ~# y, M( g( b
7 i; F/ x" J, }/ l( A
%****************************进行模糊C均值聚类******************************
( q/ h3 x! M8 c- B% Q8 F
% 设置幂指数为3，最大迭代次数为200，目标函数的终止容限为1e-6，不显示中间迭代过程
2 ]# p/ k# v0 q$ m
options = [3, 200, 1e-6, 0];
2 M* F9 B* K# j
% 调用fcm函数进行模糊C均值聚类，返回类中心坐标矩阵center，隶属度矩阵U，目标函数值obj_fcn
. p- d2 u4 E4 p [0 S+ g! F
[center,U,obj_fcn] = fcm(X,3,options);
0 ~+ h/ e" X: X5 b2 `( o) o5 `
2 u" @! D: a0 V. `/ b$ X
: _: D. `0 z, J4 [
%****************查看聚类结果*****************************
: v& F& Y- L4 C5 @- h1 J& R+ N9 O% X
id1 = find(U(1,:) == max(U));
. L( O$ {5 }2 `# Z( N5 F |
id2 = find(U(2,:) == max(U));
! K2 w* d5 f/ {
id3 = find(U(3,:) == max(U));
5 J6 ]% B1 L+ l L' ~8 S
city(id1)
7 s; s! I) O, u1 U( ^
city(id2)
" s0 X; J6 }4 f+ e# r9 o
city(id3)
* Y7 t* y' Y5 D, b0 U# m3 W
( E! u3 N3 {4 o* u1 m- B. a
为什么一查看聚类结果
, w7 ?" f1 a2 {4 h: M3 ]
就出现错误呢
9 @4 A$ x; v9 e6 n
请大神指点一下为什么
. b) x' }7 `9 f6 w
总提示指数超过矩阵维数
- P( r" U+ s+ v3 B* E, D1 ~( g5 R
# Z5 A5 X6 O* _, R7 k
# N+ v6 p+ ~$ q' |
作者:
净心、精心
时间:
2013-8-20 16:40
http://www.madio.net/thread-198267-1-1.html
2 I9 |: d: }3 T. c4 e
亲，在有问题时，记得标题要写相关的困惑和哪一方面的知识，这样才能能快的得到帮助，要不一般看到你这标题，没人愿意点开看的
欢迎光临 数学建模社区-数学中国 (http://www.madio.net/)
Powered by Discuz! X2.5