- 在线时间
- 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));2 s$ m" V5 I$ ~% f! O" B& F
x=zeros(m*n,3);% C) k6 e& X" I; l
for i=1:m
9 y: V) f* o& ^! l' e for j=1:n
|3 G; E7 J, N x(i,1)=j;. N/ ~3 c) d/ B" {7 \
x(i,2)=i;
. p. o: F9 o- J- C2 N/ e x(i,3)=data(i,j,1); 4 `* I" |2 Q% i5 Z+ }2 G7 C5 q2 A
end
% L- f O: z) z+ G% x0 _; nend
8 q2 f+ ~0 \, I( jX = zscore(xdata);7 q! R Z( \ W% ]% X9 K
" \0 D; Z" ]# o8 y0 ^9 `9 S
1 k: V. z" Q: d* Q. F%****************************进行模糊C均值聚类******************************( a5 ^( ?1 O Z8 ^( i
% 设置幂指数为3,最大迭代次数为200,目标函数的终止容限为1e-6,不显示中间迭代过程/ `- a4 \! b" h+ _2 x: ^; x
options = [3, 200, 1e-6, 0];
1 O m2 L4 B; Z+ L% 调用fcm函数进行模糊C均值聚类,返回类中心坐标矩阵center,隶属度矩阵U,目标函数值obj_fcn
[8 n: z6 F& \2 C1 [[center,U,obj_fcn] = fcm(X,3,options);. P/ _! V8 K* {* g4 g1 |8 Z4 v( W
" `" {) {+ ^8 [3 e2 Q- ?$ \
$ w" L/ k4 H0 r/ N
%****************查看聚类结果*****************************
. L% h; `* n6 z. @0 Wid1 = find(U(1,:) == max(U)); ! T0 a0 K% a5 k( M( A0 d* k
id2 = find(U(2,:) == max(U));
3 P. ~7 A$ S! \# Eid3 = find(U(3,:) == max(U));
# Y1 {' y9 t! I2 q+ Zcity(id1)
2 X( C! R7 h5 b! C) Y' [% Zcity(id2)
1 K2 W; e* h( S* f. _) hcity(id3) 2 ^. I+ V- }/ M+ ~
N' ?" q) I, e; H# e
为什么一查看聚类结果0 b# q) ^* |1 J# E+ u
就出现错误呢0 u$ E* ]/ E1 ^$ u
请大神指点一下为什么
6 F# _4 S% T/ ?2 Q- t( K3 W& c7 @+ g总提示指数超过矩阵维数
( m! t6 r$ y5 a# `0 ~9 O/ ]: }, S2 _' A; I9 ]% }
|
zan
|