- 在线时间
- 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 编辑
5 J! l' H2 [& x. m* U; m- ]: j2 K' k5 `
[m,n]=size(data(:,:,1));
0 A( p9 \" T: c/ J2 ax=zeros(m*n,3);
4 {3 q! }! ^3 z' U7 M: ~. L4 O# Kfor i=1:m1 K' i; {/ j4 A0 i" v
for j=1:n9 ]! m: W8 [ P/ ^' s9 U3 d
x(i,1)=j;. Z, m+ j% J1 B4 v* c2 P
x(i,2)=i;& F# B& F# U# x' |% J0 @- Z
x(i,3)=data(i,j,1);
3 r& T8 p0 G* g% T end. g+ a# W# H7 @6 }( L) d$ R
end
% q3 h* e! m2 TX = zscore(xdata);3 A& Y( d6 U: S& }
. S6 Q2 z% z2 Q9 d7 }% r% U. {( U. s6 B' m# I: X
%****************************进行模糊C均值聚类******************************, l- b7 j) s/ Z8 M) W) v" l) W
% 设置幂指数为3,最大迭代次数为200,目标函数的终止容限为1e-6,不显示中间迭代过程
3 A, E9 F2 ?- Z0 Y! M! N% boptions = [3, 200, 1e-6, 0]; r" {0 N6 E4 \! z
% 调用fcm函数进行模糊C均值聚类,返回类中心坐标矩阵center,隶属度矩阵U,目标函数值obj_fcn5 S, Q0 O1 W8 F4 v
[center,U,obj_fcn] = fcm(X,3,options);7 W1 m3 w6 L4 G6 s
4 \& B1 Z7 Q! ~
& {+ ~/ G5 Y* F: _! R
%****************查看聚类结果*****************************3 I3 T8 j. x% j6 D3 Q0 `, T( T+ v
id1 = find(U(1,:) == max(U)); , F+ n" s- [: C# S0 h/ i
id2 = find(U(2,:) == max(U)); / H5 Z' j' r. b% N" v9 A. u. |; i
id3 = find(U(3,:) == max(U)); $ S' b5 |' X$ x
city(id1)
! A( r% u6 h/ h0 R1 acity(id2)
& j/ T2 B/ {$ z# T3 vcity(id3)
8 r" R/ K8 D5 U4 k6 t. p7 c# @' |
为什么一查看聚类结果
- h$ A) z2 N( c1 n% G$ P就出现错误呢
{( t }+ y5 @2 c& h: u! a$ O5 ] L请大神指点一下为什么) ]8 l6 h+ V+ {
总提示指数超过矩阵维数
F/ l" T) {. I
8 u" W e* b% W$ C' ]. S! O, l9 Z4 z+ S6 W1 {
|
zan
|