当前位置: 首页 > news >正文

K-means

K均值算法(K-means)聚类

【关键词】K个种子,均值

一、K-means算法原理

聚类的概念:一种无监督的学习,事先不知道类别,自动将相似的对象归到同一个簇中。

K-Means算法是一种聚类分析(cluster analysis)的算法,其主要是来计算数据聚集的算法,主要通过不断地取离种子点最近均值的算法。

K-Means算法主要解决的问题如下图所示。我们可以看到,在图的左边有一些点,我们用肉眼可以看出来有四个点群,但是我们怎么通过计算机程序找出这几个点群来呢?于是就出现了我们的K-Means算法

这个算法其实很简单,如下图所示:

从上图中,我们可以看到,A,B,C,D,E是五个在图中点。而灰色的点是我们的种子点,也就是我们用来找点群的点。有两个种子点,所以 K = 2 K=2 K=2

然后,K-Means的算法如下:

1.随机在图中取K(这里 K = 2 K=2 K=2 )个种子点。

2.然后对图中的所有点求到这K个种子点的距离,假如点Pi离种子点Si最近,那么Pi属于Si点群。(上图中,我们可以看到A,B属于上面的种子点,C,D,E属于下面中部的种子点)

3.接下来,我们要移动种子点到属于他的“点群”的中心。(见图上的第三步)

4.然后重复第2)和第3)步,直到,种子点没有移动(我们可以看到图中的第四步上面的种子点聚合了A,B,C,下面的种子点聚合了D,E)。

这个算法很简单,重点说一下“求点群中心的算法”:欧氏距离(Euclidean Distance):差的平方和的平方根

d i j = ∑ k = 1 n ( x i k − x j k ) 2 d_{ij}=\sqrt {\sum _{k=1}^{n}(x_{ik}-x_{jk})^{2}} dij=k=1n(xikxjk)2

K-Means主要最重大的缺陷-都和初始值有关:

K是事先给定的,这个K值的选定是非常难以估计的。很多时候,事先并不知道给定的数据集应该分成多少个类别才最合适。(ISODATA算法通过类的自动合并和分裂,得到较为合理的类型数目K)

K-Means算法需要用初始随机种子点来搞,这个随机种子点太重要,不同的随机种子点会有得到完全不同的结果。(K-Means + + ++ ++ 算法可以用来解决这个问题,其可以有效地选择初始点)

总结:K-Means算法步骤:

1.从数据中选择k个对象作为初始聚类中心;

2.计算每个聚类对象到聚类中心的距离来划分;

3.再次计算每个聚类中心

4.计算标准测度函数,直到达到最大迭代次数,则停止,否则,继续操作。

5.确定最优的聚类中心

K-Means算法应用

看到这里,你会说,K-Means算法看来很简单,而且好像就是在玩坐标点,没什么真实用处。而且,这个算法缺陷很多,还不如人工呢。是的,前面的例子只是玩二维坐标点,的确没什么意思。但是你想一下下面的几个问题:

1)如果不是二维的,是多维的,如5维的,那么,就只能用计算机来计算了。

2)二维坐标点的X,Y坐标,其实是一种向量,是一种数学抽象。现实世界中很多属性是可以抽象成向量的,比如,我们的年龄,我们的喜好,我们的商品,等等,能抽象成向量的目的就是可以让计算机知道某两个属性间的距离。如:我们认为,18岁的人离24岁的人的距离要比离12岁的距离要近,鞋子这个商品离衣服这个商品的距离要比电脑要近,等等。

二、实战

重要参数:

·n_clusters:聚类的个数

重要属性:

·cluster_centers_:[n_clusters,n_features]的数组,表示聚类中心点的坐标

·labels_:每个样本点的标签

1、聚类实例

导包,使用make_blobs生成随机点cluster_std

In [ ]:

1

建立模型,训练数据,并进行数据预测,使用相同数据

In [ ]:1

In [ ]: 1

无监督的情况下进行计算,预测 现在机器学习没有目标

In [ ]:

1

绘制图形,显示聚类结果kmeans.cluster_centers

In [ ]:1

In [ ]: 1

2、实战,三问中国足球几多愁?

导包,3D图像需导包:from mpl_toolkits.mplot3d import Axes3D

读取数据

In [ ]:

1

列名修改为:“国家”,"2006世界杯","2010世界杯”,"2007亚洲杯”

In [ ]:

1

使用K-Means进行数据处理,对亚洲球队进行分组,分三组

In [ ]:

1

for循环打印输出分组后的球队,argwhere()

In [ ]:

1

绘制三维立体图形,,ax = = = plt.subplot(projectio on’3d’) \text {on'3d')} on’3d’)

ax.scatter3D()

In [ ]:

3、聚类实践与常见错误

导包,使用make_blobs创建样本点

In [ ]:

1

第一种错误,k值不合适,make_blobs默认中心点三个

In [ ]:

1

第二种错误,数据偏差

trans = [ [ 0.6 , − 0.6 ] , [ − 0.4 , 0.8 ] ] \text {trans}=[[0.6,-0.6],[-0.4,0.8]] trans=[[0.6,0.6],[0.4,0.8]]

X 2 = np.dot ( X , trans ) X2=\text {np.dot}(X,\text {trans}) X2=np.dot(X,trans)

In [ ]:

1

第三个错误:标准偏差不相同cluster_std

1

In [ ]:

In [ ]:

1

第四个错误:样本数量不同

1

In [ ]:

In [ ]:

1

4、K-Means图片颜色点分类

导包from sklearn.metrics import pairwise_distances_argmin

In [ ]:

1

加载图片/创建模型/训练数据/获取labels

In [ ]:

1

使用pairwise_distances_argmin算法生成随机labels

In [ ]:

1

创造方法,重新合成图片

In [ ]:

1

数据可视化

In [ ]:

1
1

In [ ]:

三、作业

1、分析ex7data2.mat文件

找出最佳聚类数目,并画出聚类的中心点

from scipy.io import loadmat

from sklearn.cluster import KMeans

data = 1 \text {data}=1 data=1 oadmat(‘…/data/ex7data2.mat’)

X = data [ ′ X ′ ] \mathrm {X}=\text {data}\left[{}^{\prime }\mathrm {X}^{\prime }\right] X=data[X]

X.shape

2、图片压缩

使用聚类压缩图片

img = = = plt.imread(‘…/data/bird_small.png’)

img_shape = = = img.shape

img_shape

相关文章:

  • 机器学习-简要与数据集加载
  • React Native【详解】搭建开发环境,创建项目,启动项目
  • 虚拟机连不上网只有lo地址
  • 高频PCB设计如何选择PCB层数?
  • 2025年LangChain(V0.3)开发与综合案例
  • Spring MVC 如何自动将请求参数映射到 Controller 方法的参数对象(POJO)上?
  • 把其他conda的env复制到自己电脑的conda上
  • 【ULMD】基于单峰标签生成和模态分解的多模态情感分析
  • Java大师成长计划之第13天:Java中的响应式编程
  • Python实现NOA星雀优化算法优化BP神经网络分类模型项目实战
  • P56-P60 统一委托,关联游戏UI,UI动画,延迟血条
  • CSS Border 三角形阴影与多重边框的制作
  • ROS2:自定义接口文件(无废话)
  • 第100+40步 ChatGPT学习:R语言实现多轮建模
  • 使用Homebrew下载配置git和连接GitHub(Mac版)
  • dubbo限流
  • VMware Fusion安装win11 arm;使用Mac远程连接到Win
  • 使用 OpenSSL 吊销 Kubernetes(k8s)的 kubeconfig 里的用户证书
  • 论文速读《Embodied-R: 基于强化学习激活预训练模型具身空间推理能力》
  • 如何添加或删除极狐GitLab 项目成员?
  • 新华时评:任凭风云变幻,中俄关系从容前行
  • 陕南多地供水形势严峻:有的已呼吁启用自备水井
  • 潘功胜:坚定支持汇金公司在必要时实施对股票市场指数基金的增持
  • 央行:5月8日起7天期逆回购操作利率由此前的1.50%调整为1.40%
  • 湖北十堰市委副秘书长管聪履新丹江口市代市长
  • 山东如意集团及实控人等被出具警示函,此前实控人已被罚十年禁止入市