python 谱聚类　幂迭代

最新推荐文章于 2024-05-30 11:11:00 发布

ASD991936157

最新推荐文章于 2024-05-30 11:11:00 发布

阅读量1.1k

点赞数

本文链接：https://blog.csdn.net/ASD991936157/article/details/54317393

版权

本文延续前文，详细介绍了如何使用Python实现幂迭代谱聚类算法，通过实例展示了聚类过程，对于数据集的分类提供了有效的方法。

摘要由CSDN通过智能技术生成

紧接着上篇博文，实现了幂迭代聚类算法：

# encoding=utf-8
import numpy as np
import matplotlib.pyplot as plt
from numpy import linalg as LA
from sklearn.cluster import KMeans
from sklearn.metrics.pairwise import rbf_kernel
from sklearn.preprocessing import normalize
from sklearn.datasets import make_blobs


def similarity(points):
    """
    亲和矩阵
    :param points:
    :return:
    """
    res = rbf_kernel(points)
    for i in range(len(res)):
        res[i, i] = 0
    return res


def spectral(points, k):
    """
    谱聚类
    :param points:
    :param k:
    :return:
    """
    W = similarity(points)
    Dn = np.diag(np.power(np.sum(W, axis=1), -0.5))
    L = np.eye(len(points)) - np.dot(np.dot(Dn, W), Dn)
    eigvals, eigvecs = LA.eig(L)
    indices = np.argsort(eigvals)[:k]
    subvecs = normalize(eigvecs[:, indices])
    return KMeans(n_clusters=