- 在线时间
- 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));3 n% I& Q7 ^( k* m$ |# [
x=zeros(m*n,3); Z4 v0 j! \4 B/ Q, r
for i=1:m! |6 r6 [! F7 y) {+ v
for j=1:n$ W+ Y' p: O7 A: A5 d0 X g
x(i,1)=j;
, _1 \) P1 B. i# m( ?$ O/ V x(i,2)=i;
) v9 ^/ u4 E8 Z$ P. X k- p x(i,3)=data(i,j,1); 1 l) I R' B4 p4 o$ P& [
end
) [% S* M& x2 e% r5 Y o. S6 Fend# w' P( G( F7 V+ X, o K# g
X = zscore(xdata);
/ S( w, ]1 a' Q
0 i# {, ~8 l. o: q$ a: L' d$ ?+ R/ c$ V+ @9 G2 N( V# @: c- A
%****************************进行模糊C均值聚类******************************
4 Q* q5 U2 e* E" ^' U$ O, B2 ?: z% 设置幂指数为3,最大迭代次数为200,目标函数的终止容限为1e-6,不显示中间迭代过程
" e/ E8 K. e! F+ ioptions = [3, 200, 1e-6, 0];, \3 K, q0 |5 @/ A9 x% f1 r/ g
% 调用fcm函数进行模糊C均值聚类,返回类中心坐标矩阵center,隶属度矩阵U,目标函数值obj_fcn
4 e: N. Z' {. l7 _& o[center,U,obj_fcn] = fcm(X,3,options);* Y7 G, h& @* K8 O& d, P( C/ A: S
; B* Y4 }2 M" B' K
/ {9 D1 A5 J( f- H! }
%****************查看聚类结果***************************** g: n4 y0 ^& [: m3 g u: l
id1 = find(U(1,:) == max(U)); * S7 t7 l; O' a1 L7 r+ i: h, a$ s' Z
id2 = find(U(2,:) == max(U)); ' c4 ^6 d9 H2 C, [* ^- I
id3 = find(U(3,:) == max(U)); 8 ?# [# K' N0 B
city(id1) $ r1 w0 t+ m1 ]
city(id2) # L' F: T/ {" ` ]' m
city(id3) ! i7 T1 n4 `1 d- `, D
3 V/ W9 y* t9 u, o% J+ V3 W' w
为什么一查看聚类结果
5 h7 ?5 h- G' d; h8 l$ K) ~就出现错误呢6 `& x" L2 N; B0 j, ?0 R
请大神指点一下为什么
4 F9 N4 H# {, A$ I* W$ x& J) h' A总提示指数超过矩阵维数
* Y" o j7 N7 o5 k7 D5 a' _2 E2 X5 U2 D! Y
|
zan
|