如何计算方阵的平方根倒数

问题描述

任何人都可以确认这是否是计算谱聚类的度矩阵的平方根平方根的正确方法

给定一个邻接方阵A,邻接矩阵

D 是度矩阵,L 是拉普拉斯矩阵

D = np.diag(1/np.sqrt(np.sum(A,axis=1)))

L = D @ A @ D

或者我应该使用

D = A**-.5

我在两种方法的结果矩阵中都看到了 inf 值。

解决方法

这对我有用.. 避免任何除以零的情况。

D = np.diag(1. / np.sqrt(sum(A,1)))