- 在线时间
- 25 小时
- 最后登录
- 2016-11-8
- 注册时间
- 2013-1-12
- 听众数
- 8
- 收听数
- 0
- 能力
- 0 分
- 体力
- 348 点
- 威望
- 0 点
- 阅读权限
- 30
- 积分
- 145
- 相册
- 0
- 日志
- 0
- 记录
- 0
- 帖子
- 80
- 主题
- 3
- 精华
- 0
- 分享
- 0
- 好友
- 5
升级   22.5% TA的每日心情 | 怒 2016-9-15 09:07 |
|---|
签到天数: 50 天 [LV.5]常住居民I
- 自我介绍
- 东北农大本科生
 |
本帖最后由 袁海亮 于 2013-8-20 16:15 编辑 - _) S5 G' U5 a8 k9 p
5 o8 ]3 y3 _4 C2 O( m
[m,n]=size(data(:,:,1));
6 k" B1 }# ]. Q# @# sx=zeros(m*n,3);
0 A8 T6 r" s! Z9 U& w$ W" E9 E$ D, [" gfor i=1:m* }" _3 [) D* K- Z: [
for j=1:n, ?5 a0 I; q/ B3 K
x(i,1)=j;+ H1 {0 n% j( W# a# Q
x(i,2)=i;: K4 R4 `2 k% V/ h
x(i,3)=data(i,j,1);
: {! c3 T& @& |2 h7 l- z& X end
# O$ g) u' `, O- hend- V) ]5 H& x7 [9 y
X = zscore(xdata);
8 x0 j: r3 V ]1 A4 n! k; }% n6 C* K- Q) C2 ?1 f; d1 ^
. {, L* F6 L$ s! t
%****************************进行模糊C均值聚类******************************* p9 r+ K7 [+ K# e. L1 a# V( t
% 设置幂指数为3,最大迭代次数为200,目标函数的终止容限为1e-6,不显示中间迭代过程
# N- w& X: o" Aoptions = [3, 200, 1e-6, 0];: }) [- T. z* m5 {& L
% 调用fcm函数进行模糊C均值聚类,返回类中心坐标矩阵center,隶属度矩阵U,目标函数值obj_fcn4 ^+ `2 |# a3 B
[center,U,obj_fcn] = fcm(X,3,options); `$ [9 h4 A6 m
+ [! L. X6 Q$ ]6 ]
" ^8 E/ C% r6 [9 F4 e%****************查看聚类结果*****************************
% d* V, M$ G4 h# h/ ?id1 = find(U(1,:) == max(U)); , g! ?6 U4 `: W* F! \+ H# ^( F
id2 = find(U(2,:) == max(U));
5 x8 g" Y2 r ?1 p6 I! R) Iid3 = find(U(3,:) == max(U)); ( O9 V4 K5 B" q" c
city(id1) - s: D4 Z# T2 {$ |1 I6 u. A
city(id2)
2 l" H. ^& z; bcity(id3)
8 H/ C. l+ y" A. Y: W! @) S; m/ C- J/ t: c: n: t
为什么一查看聚类结果
* E& q: {2 I9 V+ g8 ?+ B& A9 t就出现错误呢: S3 @, H" q7 V6 B
请大神指点一下为什么
3 i9 h* ~- x& Z% b" L' r+ \总提示指数超过矩阵维数% k* v1 _8 e9 j1 y/ G: d
' z: b: U0 b2 t" h2 k# D. z
/ w% E; {. o+ h- c6 b$ ?0 k
|
zan
|