- 在线时间
- 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 编辑 8 G* V4 r* S' n' d
% o, v2 V+ N* m2 ^' T1 r9 S
[m,n]=size(data(:,:,1));
. U; p# h0 G8 G' X/ Ex=zeros(m*n,3);
; R& ` l; p- L* ]for i=1:m2 D1 \' r7 C Q: W0 U4 [9 t0 d
for j=1:n
( v% X; a) K2 X: E$ v( _+ P x(i,1)=j; m2 u6 A* e# i2 C( I2 ^* Z# S' X7 d
x(i,2)=i;2 o4 \. K* o' {. ]! n" K% x6 L
x(i,3)=data(i,j,1); C" a! |5 h$ |4 a
end
1 q9 ~! I3 C2 v% u9 A5 [7 Uend
( o4 z9 N2 h4 u9 ~2 mX = zscore(xdata);
. k5 s" X) Q! s' r: r1 n
. ~7 g: ^' x7 O9 d
3 p& V6 n: u3 V7 q/ C%****************************进行模糊C均值聚类******************************
/ E w7 Z' u- ?' q# ]: H! C% 设置幂指数为3,最大迭代次数为200,目标函数的终止容限为1e-6,不显示中间迭代过程
# e6 D% F! ?! b' [; i m; poptions = [3, 200, 1e-6, 0];- {' X* C! o1 P8 |1 d+ s: G
% 调用fcm函数进行模糊C均值聚类,返回类中心坐标矩阵center,隶属度矩阵U,目标函数值obj_fcn
4 }; `3 j2 A3 L6 X$ W[center,U,obj_fcn] = fcm(X,3,options);- B" P4 F4 i' d: Y0 J! @
6 @6 l+ x. q% c2 U) t' l
/ a. P+ M0 C4 w5 |: r. _. a
%****************查看聚类结果*****************************3 D a' n, |0 [3 g* S0 {
id1 = find(U(1,:) == max(U)); r. M5 A6 t; F
id2 = find(U(2,:) == max(U)); + _8 [' t8 U7 }
id3 = find(U(3,:) == max(U));
& l. r) x, {/ I9 Y8 n6 lcity(id1)
# Y9 |. `4 Y5 F9 Z9 ncity(id2) 8 a8 v0 B7 o, K) v! q
city(id3) - c* g. l; p1 o% ^! e9 N7 x
x M3 Z+ F9 l8 W为什么一查看聚类结果
, i* O* ]. F2 H' |/ [* }/ O# k就出现错误呢1 k! z+ A# v/ C5 V- O( i5 \$ T3 ?4 G
请大神指点一下为什么% W7 z8 R, ?( O* p
总提示指数超过矩阵维数" I/ F: }8 F Y q( b: V
6 n/ g/ i* a1 k% W% [# L
. ~' A) {) C5 n |
zan
|