- 在线时间
- 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 编辑
% H3 G1 @. J+ g2 t) k) ]* Z( F9 ^) b- |9 p: q0 n8 Q8 K# L
[m,n]=size(data(:,:,1)); i+ U- u* m; X0 Z
x=zeros(m*n,3); L# u% x: O/ I$ Y6 I5 c$ I
for i=1:m
+ G0 p5 R5 z9 c# }+ @ for j=1:n% v C l# U4 C, i$ g Y# y5 G- {
x(i,1)=j;
* S' i) |- _( q5 W6 T0 n& ` x(i,2)=i;* r9 K4 @+ u7 C' R2 H- S
x(i,3)=data(i,j,1); , m9 R3 u7 a- r% D5 _
end
$ T) }9 K' |; n6 Wend
+ o5 q7 y2 r& S+ tX = zscore(xdata);
~2 G4 ?% i* Q6 |; J& o1 X1 X9 k( Y' c9 T% L- s; l* m/ s
- S4 ~0 d# } c5 V; \4 T%****************************进行模糊C均值聚类******************************
& u0 N' r) V% M7 \% 设置幂指数为3,最大迭代次数为200,目标函数的终止容限为1e-6,不显示中间迭代过程
& a/ z+ A, z- H: Z& Noptions = [3, 200, 1e-6, 0];) V# m( q. F% R y
% 调用fcm函数进行模糊C均值聚类,返回类中心坐标矩阵center,隶属度矩阵U,目标函数值obj_fcn4 T- X" Y( n% ]2 `9 o. j+ G
[center,U,obj_fcn] = fcm(X,3,options);2 @( F4 R# E9 O8 A/ V& p
! K8 w3 w7 T5 E
: {5 T7 O" t" R%****************查看聚类结果*****************************
7 O5 A) o& f7 pid1 = find(U(1,:) == max(U));
# ~$ e* q7 y$ a( V xid2 = find(U(2,:) == max(U));
+ S0 j1 t: a. p/ jid3 = find(U(3,:) == max(U)); 9 m) n* C* l# P
city(id1)
0 E: M4 H9 i9 ^4 W5 Gcity(id2) % a9 p2 R+ E* k2 B* U
city(id3)
+ Q- z$ y- m9 {* r2 _, j3 V& U8 w2 M2 d! f7 u
为什么一查看聚类结果
! I0 S/ h. s! k, e7 f( E0 E4 ?( J就出现错误呢
6 K9 |# A# H: Y6 {/ l4 c( D! u请大神指点一下为什么
: ^1 x4 ^% C |! _总提示指数超过矩阵维数1 Z0 l' M. g, L4 r% p
/ C3 A! h+ W- Z, B! Q( T i
9 S& y7 d$ K& m' W- f5 x: q# w* a
|
zan
|