- 在线时间
- 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 编辑
% V7 D% y2 R5 `' U1 r3 J) ? a$ V8 [- C+ L9 i* {! M1 k5 h
[m,n]=size(data(:,:,1));
( S: Y' `! Z% Ix=zeros(m*n,3);2 A) ?5 ]' x) C; A8 l2 ~; g
for i=1:m. l& B+ J5 H9 l" K. e I" i
for j=1:n
. E4 {4 Z l/ s x(i,1)=j;* D& |. K* r# _0 T' k
x(i,2)=i;
: n" Y8 z( A# e: A2 {0 V x(i,3)=data(i,j,1); 0 N( M6 i0 x" U6 i* D4 ~/ D0 l! Q
end
( _/ w& B, V7 ~# _1 c+ tend
+ V+ O) Y$ H0 j" [X = zscore(xdata);/ I" I W. t+ Y
1 b/ O7 t7 k5 E+ l: O' ], `, V/ C3 t! T% |1 Q
%****************************进行模糊C均值聚类******************************/ P5 v. |0 e0 G" I# U% M2 x% L
% 设置幂指数为3,最大迭代次数为200,目标函数的终止容限为1e-6,不显示中间迭代过程
% j. T( J+ g; p6 Joptions = [3, 200, 1e-6, 0];
/ f$ @7 R4 e: o6 }2 L+ R5 Q# S2 ]5 D% 调用fcm函数进行模糊C均值聚类,返回类中心坐标矩阵center,隶属度矩阵U,目标函数值obj_fcn
+ \6 k* L& q. C6 F: Z- B1 L[center,U,obj_fcn] = fcm(X,3,options);6 ^0 S/ p) K) w! o
" k4 ]/ `* G) C& I
$ \ x$ b9 c: f4 b%****************查看聚类结果*****************************4 t) z0 G* {( z4 Z: S
id1 = find(U(1,:) == max(U));
1 n! X! y; n( B. D. rid2 = find(U(2,:) == max(U));
( ?/ K8 Z' D# Fid3 = find(U(3,:) == max(U)); / s7 q' ^- A9 d5 E# w! ?' Z
city(id1)
5 P# { e1 N9 ]! m- [7 ?city(id2) 0 L1 h' L0 S" t! d! K( m
city(id3)
; m8 s( M* w4 e* o3 b0 R. g: R& ?/ {5 o' e! n
为什么一查看聚类结果
. s e* n$ N$ C0 ]1 G就出现错误呢" ? W& c1 B4 {% l2 Q
请大神指点一下为什么
/ E; o/ ~' R- X6 X1 I总提示指数超过矩阵维数: [1 n: y) r& S; q% m2 c
* G8 M; R* P' @& g
5 G! Q6 P5 @1 v
|
zan
|