数学建模社区-数学中国
标题:
在线等回复
[打印本页]
作者:
509
时间:
2013-8-20 15:55
标题:
在线等回复
本帖最后由 袁海亮 于 2013-8-20 16:15 编辑
( h& w7 F" D- J2 a G3 \
/ n) u5 ^6 f" n7 t* R% ^
[m,n]=size(data(:,:,1));
6 ]! p$ I! v# I3 f% E2 W
x=zeros(m*n,3);
" `% A3 i+ n5 _$ s
for i=1:m
( l) k3 W- Q5 S# G& t6 `
for j=1:n
- P2 B2 l+ K6 F4 Z1 p3 C
x(i,1)=j;
- Z1 O1 n0 {# o. B8 o- V2 w
x(i,2)=i;
5 [/ b! b5 @+ u6 z9 U/ ?8 B
x(i,3)=data(i,j,1);
, ]- T; a& V' z6 H* I
end
' }& q" g* v/ p. z2 `5 ] [
end
2 n1 \4 O: D8 W$ g3 v% g. h% N; T% ]
X = zscore(xdata);
& O6 b9 D5 Y% X$ z# [( [6 C
) \0 H& I0 y/ D3 R4 w) R# ]9 {
& |1 y: m: {9 X
%****************************进行模糊C均值聚类******************************
( p1 \ S/ A0 ?( o" E1 S- r+ l' c6 V3 x
% 设置幂指数为3，最大迭代次数为200，目标函数的终止容限为1e-6，不显示中间迭代过程
7 o: i+ x' Q* S z* q0 z$ R
options = [3, 200, 1e-6, 0];
& t5 q6 I; [% k% s' U6 l1 M+ W, M# S
% 调用fcm函数进行模糊C均值聚类，返回类中心坐标矩阵center，隶属度矩阵U，目标函数值obj_fcn
8 C* t' q& z! z/ w; {0 D) R
[center,U,obj_fcn] = fcm(X,3,options);
& m+ j0 u8 }! n
- u9 |: `, w& u" }' C0 }
; ^3 D; T# i- V8 H3 ^
%****************查看聚类结果*****************************
' [9 I, Q) a9 O5 g% V, p
id1 = find(U(1,:) == max(U));
. \5 {4 E9 e2 x( d+ T/ j
id2 = find(U(2,:) == max(U));
6 i0 z) ^3 t4 M- S
id3 = find(U(3,:) == max(U));
* j [* D8 z; N! q6 i5 l" `
city(id1)
: B4 z# W' \3 d( j2 i, j
city(id2)
9 ^. g/ a* g7 z
city(id3)
# @* O9 E0 p1 W- _ H+ k
. V }% W1 @1 o1 o
为什么一查看聚类结果
7 |5 M9 ~6 ^5 r
就出现错误呢
A) @+ s: @' m E: J
请大神指点一下为什么
" v |% y' @( Z8 m1 O l2 b
总提示指数超过矩阵维数
2 {$ q0 \$ I$ [
$ x& Z0 [4 |: @) a- s" f% T; p
2 `. P1 n6 x0 y; w& P. o* o
作者:
净心、精心
时间:
2013-8-20 16:40
http://www.madio.net/thread-198267-1-1.html
; y) M9 y$ b$ M# x
亲，在有问题时，记得标题要写相关的困惑和哪一方面的知识，这样才能能快的得到帮助，要不一般看到你这标题，没人愿意点开看的
欢迎光临 数学建模社区-数学中国 (http://www.madio.net/)
Powered by Discuz! X2.5