加载中…
个人资料
  • 博客等级:
  • 博客积分:
  • 博客访问:
  • 关注人气:
  • 获赠金笔:0支
  • 赠出金笔:0支
  • 荣誉徽章:
正文 字体大小:

转载:MATLAB K-Means 聚类函数Kmeans 初值设置

(2011-10-11 20:49:24)
标签:

杂谈

K-means聚类算法采用是将N*P矩阵X划分为K个类,使得类内对象之间距离最大,而类之间距离最小。



使用方法:
Idx=Kmeans(X,K)
[Idx,C]=Kmeans(X,K)
[Idx,C,sumD]=Kmeans(X,K)
[Idx,C,sumD,D]=Kmeans(X,K)
[…]=Kmeans(…,’Param1’,Val1,’Param2’,Val2,…)


各输入输出参数介绍:


X N*P数据矩阵
K 表示将X划分为几类,为整数
Idx N*1向量,存储是每个点聚类标号
C K*P矩阵,存储是K个聚类质心位置
sumD 1*K和向量,存储是类间所有点与该类质心点距离之和
D N*K矩阵,存储是每个点与所有质心距离

[…]=Kmeans(…,'Param1',Val1,'Param2',Val2,…)
这其中参数Param1、Param2等,主要可以设置为如下:

1. ‘Distance’(距离测度)
‘sqEuclidean’ 欧式距离(默认时,采用此距离方式)
‘cityblock’ 绝度误差和,又称:L1
‘cosine’ 针对向量
‘correlation’   针对有时序关系
‘Hamming’ 只针对二进制数据

2. ‘Start’(初始质心位置选择方法)
‘sample’ 从X中随机选取K个质心点
‘uniform’ 根据X分布范围均匀随机生成K个质心
‘cluster’ 初始聚类阶段随机选择10%X子样本(此方法初始使用’sample’方法)
matrix 提供一K*P矩阵,作为初始质心位置集合

3. ‘Replicates’(聚类重复次数)   整数


              
使用案例:

data=
5.0 3.5 1.3 0.3 -1
5.5 2.6 4.4 1.2 0
6.7 3.1 5.6 2.4 1
5.0 3.3 1.4 0.2 -1
5.9 3.0 5.1 1.8 1
5.8 2.6 4.0 1.2 0

[Idx,C,sumD,D]=Kmeans(data,3,'dist','sqEuclidean','rep',4)



运行结果:
Idx =
    1
    2
    3
    1
    3
    2

C =
5.0000 3.4000 1.3500 0.2500 -1.0000
5.6500 2.6000 4.2000 1.2000       0
6.3000 3.0500 5.3500 2.1000 1.0000

sumD =
0.0300
0.1250
0.6300

D =
0.0150 11.4525 25.5350
12.0950 0.0625 3.5550
29.6650 5.7525 0.3150
0.0150 10.7525 24.9650
21.4350 2.3925 0.3150
10.2050 0.0625 4.0850

matlab-kmeans函数注释

X = [randn(100,2)+ones(100,2);...
     randn(100,2)-ones(100,2)]; 产生100个样本点,行指向每个样本,列是维变量值。
opts = statset('Display','final');

[idx,ctrs] = kmeans(X,2,'Distance','city','Replicates',5,'Options',opts);
%返回参数意义:[IDX,C,sumd,D]=kmeans()
IDX:每个样本点所在的类别
C:所聚类别的中心点坐标位置k*p,k是所聚类别
sumd:每个类内各点到中心点的距离之和
D:每个点到各类中心点的距离n*k

Matlab聚类分析中kmeans函数运行结果,请教为什么?

k=5;
[IDX,C,sumd,D] =kmeans(SCORE(:,1:3),k);
我想把主成份分析后的结果SCORE(:,1:3),大小为89*3,聚成5类,我的理解的运行结果应该是:
IDX是1~5的整数,表示归到了那一类;
C是每一类的质心位置,大小是5*3;
sumd是每一类中各点到质心的距离和,大小是1*5;
D是每个点到质心的位置,大小是89*1
 

0

阅读 收藏 喜欢 打印举报/Report
  

新浪BLOG意见反馈留言板 欢迎批评指正

新浪简介 | About Sina | 广告服务 | 联系我们 | 招聘信息 | 网站律师 | SINA English | 产品答疑

新浪公司 版权所有