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

how to calculate inverse square root of a square matrix

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

给定一个邻接方阵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)))