- 在线时间
- 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 编辑
' |% N) g" ?% N+ ^5 m/ _
, M$ g8 c0 ^( S" j) h/ d[m,n]=size(data(:,:,1));
# t9 O% _( K2 d* T! w/ |1 px=zeros(m*n,3);
2 y) U3 Z# P' s0 Sfor i=1:m
3 m# p+ A' x7 V3 N( x$ I for j=1:n# `' F: r, X$ a+ D% p5 w1 t1 z9 h* v
x(i,1)=j;
c; s% g; q- }- O2 l x(i,2)=i;! Z- u* O0 K# `
x(i,3)=data(i,j,1); : E4 |( v, x9 p5 B! ]8 w9 ^$ u
end
; T8 j5 ]8 m! Q. o. Rend+ R! T* W/ c7 r4 {8 ]0 l( @; D0 E
X = zscore(xdata);
( V/ }$ }4 B& o* _+ m: l) Q" `" B1 O; A* G* \' J. I3 h
* W5 `7 y' E8 M* L; R
%****************************进行模糊C均值聚类******************************
$ P7 v+ C9 @' z9 V$ w- p7 k% 设置幂指数为3,最大迭代次数为200,目标函数的终止容限为1e-6,不显示中间迭代过程
. p8 s: U2 v9 N# G5 s% ^0 v9 Hoptions = [3, 200, 1e-6, 0]; Y, e% P3 l! G$ X. t
% 调用fcm函数进行模糊C均值聚类,返回类中心坐标矩阵center,隶属度矩阵U,目标函数值obj_fcn
# K! C4 L) d+ c[center,U,obj_fcn] = fcm(X,3,options);, @0 i5 g/ d5 W6 y& j; w
" }6 v1 V0 q1 W$ R {
$ c' r* ~8 B# O B3 {& D
%****************查看聚类结果****************************** B( w. a3 }7 g0 l; E
id1 = find(U(1,:) == max(U));
4 O0 }: v% L: j8 T! Kid2 = find(U(2,:) == max(U));
9 U9 a/ Z+ G$ ~1 A# S' _id3 = find(U(3,:) == max(U));
% j" p3 H6 U" B) N/ scity(id1)
8 l- `8 Y5 I1 }1 B' N. ^" Dcity(id2)
8 l1 |4 z$ t" G% j& H/ r( P Z8 O- fcity(id3)
+ u5 @. f- A' f0 J/ }
7 O- ]: @/ B: U9 C& U5 k9 a9 c Q为什么一查看聚类结果
) v0 ?5 C- A7 C- a就出现错误呢0 ]; b# @6 C) V T
请大神指点一下为什么
( `& K( ]! Z" A2 c) A" `& k总提示指数超过矩阵维数
, Z( p' z3 X* z \ R
9 @2 Z. C5 L: ?9 l s
, w! q# y2 H8 V, x6 [3 t5 L |
zan
|