- 在线时间
- 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 编辑 9 T& C O5 |* @2 N' W- |
7 F0 N* L- G' n/ Q: e* f2 C
[m,n]=size(data(:,:,1));
0 y: g) Y; n! y- c# ox=zeros(m*n,3);
5 K5 b# y% p0 o. h& p4 v8 pfor i=1:m
& R, [8 L, A" x9 |) L- v7 H t for j=1:n- O; l7 _9 m6 A
x(i,1)=j;; C. H# b# w" d4 \+ I
x(i,2)=i;
$ j( D2 X9 P, Z1 N. a7 J x(i,3)=data(i,j,1); / q4 g4 V+ l/ V" ^3 X* [2 `) l. q
end2 p- `( a) a& D& |+ c* a- r: @
end. c, y2 z* s+ j" _$ K; q
X = zscore(xdata);
% ~4 A- U0 Z" Z' E0 D7 ~ ~
$ s( k. j- n3 { m# {1 s. }9 p7 F: M W& a- p D0 ~
%****************************进行模糊C均值聚类******************************# Z" `2 T0 U! J) o9 ^5 _) y
% 设置幂指数为3,最大迭代次数为200,目标函数的终止容限为1e-6,不显示中间迭代过程
) t, ?: q( m( x0 j# P. N# z% {options = [3, 200, 1e-6, 0];* Q" k; k4 z) ^- n" O0 Q; y$ r+ `
% 调用fcm函数进行模糊C均值聚类,返回类中心坐标矩阵center,隶属度矩阵U,目标函数值obj_fcn
* q- a% g+ [ S6 Z7 u4 U[center,U,obj_fcn] = fcm(X,3,options);( u8 c2 ^: d; w( d
, O ~/ Y, z6 X
- Q- J3 x G0 H7 N- t3 O9 E$ ^%****************查看聚类结果*****************************
; k# ]. L8 q+ J- c& Did1 = find(U(1,:) == max(U)); # U( d+ y& \+ t. ^1 g) z
id2 = find(U(2,:) == max(U));
4 S1 i7 Q5 T; }id3 = find(U(3,:) == max(U)); " u: k" Z' b( u6 m! a# z3 o! ~
city(id1) . m. R& n0 z* ~, d1 m- [% F P
city(id2) 3 T# c0 E# V9 [$ Q
city(id3) & s5 [, n7 D1 i8 C+ k: {' `9 o
+ ]0 h+ j$ e5 v6 K- u1 W
为什么一查看聚类结果" O3 t: B3 Y. X! }6 c
就出现错误呢
2 `5 j. k7 H( u/ K6 X' Y6 V请大神指点一下为什么
q7 b1 D; \* B总提示指数超过矩阵维数
: A! c" W) }0 C+ w8 L! L g p$ h. h! k- S0 @4 g- \: I7 O
, X- y) D" l$ w( p; c) i/ F' q: Z |
zan
|