- 在线时间
- 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
- 自我介绍
- 东北农大本科生
 |
[m,n]=size(data(:,:,1));$ w+ f" a. [+ Q8 \' s. t+ [
x=zeros(m*n,3);2 L! c7 |* d3 |! U
for i=1:m0 w3 z4 K- s, B; p) y6 t
for j=1:n
1 j0 |6 Z7 h& H) y x(i,1)=j;
1 q5 y/ f* _, o' J0 v2 i) s% X x(i,2)=i;
5 M8 N6 m5 `/ Q0 O* z+ b z x(i,3)=data(i,j,1); " e, M2 \* ]& i/ S5 f) ]
end
* R7 J3 u. a$ D5 Z* k3 [6 ?end
( l9 A9 g/ D$ t( y3 H, hX = zscore(xdata);' ?: e+ j7 Y2 E! n+ ~. B7 I
% H) {& _ D" `8 H
5 L# o5 Z2 g S* z, f%****************************进行模糊C均值聚类******************************3 U [9 o$ D# i' W
% 设置幂指数为3,最大迭代次数为200,目标函数的终止容限为1e-6,不显示中间迭代过程) r+ v }: c* P
options = [3, 200, 1e-6, 0];3 P8 G e" V) g% j* h, c1 W( \
% 调用fcm函数进行模糊C均值聚类,返回类中心坐标矩阵center,隶属度矩阵U,目标函数值obj_fcn
, ~5 d/ T/ o" r; X[center,U,obj_fcn] = fcm(X,3,options);
6 f( B0 a9 r" M
# r6 x6 x' E; e7 K7 u O/ l$ w# M6 `/ s" |" X; X2 t q( s/ {
%****************查看聚类结果*****************************
. D$ t8 f E* _id1 = find(U(1,:) == max(U));
, J4 @- U/ f$ b% i: p, I# P @3 yid2 = find(U(2,:) == max(U)); 1 W7 ^7 Z- k' Y" W
id3 = find(U(3,:) == max(U)); , }8 t" U1 a9 @1 I! x4 I6 y
city(id1) 7 c* B. e. X4 l
city(id2)
! H, I2 p% V7 hcity(id3)
; d# K" E, l) g, E8 M* N T4 S1 p9 R1 D2 Q g6 b' @
为什么一查看聚类结果. I6 Z2 F. n0 _( o& S! M/ l
就出现错误呢- t% k: K% ]* T& [* g, @
请大神指点一下为什么1 H V* T2 q6 K1 O9 }4 O- [
总提示指数超过矩阵维数
' y+ g' u+ k* E
* @- a7 E9 ]$ x' P- ?7 X. |" E4 \# X |
zan
|