用python写一段k-means聚类算法代码
时间: 2023-09-29 11:08:23 浏览: 107
好的,以下是一段基于Python的K-Means聚类算法代码:
```python
import random
import numpy as np
import matplotlib.pyplot as plt
# 生成测试数据
def generate_data(n, k):
data = []
for i in range(k):
center = np.random.uniform(-10, 10, size=2)
for j in range(n):
x = np.random.normal(center[0], 1)
y = np.random.normal(center[1], 1)
data.append([x, y, i])
return np.array(data)
# 计算两点之间的距离
def distance(point1, point2):
return np.sqrt(np.sum((point1 - point2) ** 2))
# K-Means算法
def k_means(data, k, max_iter=100):
# 随机选择k个中心点
centers = data[random.sample(range(len(data)), k)][:, :-1]
for i in range(max_iter):
# 初始化每个簇的点
clusters = [[] for _ in range(k)]
# 分配每个点到最近的簇
for point in data:
distances = [distance(point[:-1], center) for center in centers]
cluster_index = np.argmin(distances)
clusters[cluster_index].append(point)
# 更新每个簇的中心点
for j in range(k):
centers[j] = np.mean(clusters[j], axis=0)[:-1]
return clusters, centers
# 绘制聚类结果
def plot_clusters(clusters, centers):
plt.figure(figsize=(8, 6))
colors = ['r', 'g', 'b', 'c', 'm', 'y', 'k']
for i in range(len(clusters)):
cluster = np.array(clusters[i])
plt.scatter(cluster[:, 0], cluster[:, 1], c=colors[i], alpha=0.5)
plt.scatter(centers[:, 0], centers[:, 1], s=100, marker='*', c='black')
plt.show()
# 测试
data = generate_data(100, 3)
clusters, centers = k_means(data, 3)
plot_clusters(clusters, centers)
```
在此代码中,我们首先生成了一些测试数据,然后使用K-Means算法将其聚类成k个簇,并绘制聚类结果。在该算法中,我们首先随机选择k个中心点,然后迭代地将每个点分配到最近的簇中,并更新每个簇的中心点。最终,我们返回每个簇的点和中心点,并绘制聚类结果。
阅读全文