如何使用python获得方阵的伪行列式

Question

我有一个矩阵，它在我计算朴素贝叶斯分类器的奇异测试中失败。我正在处理ln(det(sigma))等式的一部分。

if np.linalg.cond(covarianceMatrix) < 1/sys.float_info.epsilon:
    return np.log(np.linalg.det(covarianceMatrix))
else:
    return a pseudo determinant

当协方差矩阵是奇异矩阵时，我必须找到伪行列式。我怎样才能做到这一点？

Answer 1

你能使用numpy的pinv来计算伪逆，然后用它来计算行列式吗？ http://www.sosmath.com/matrix/inverse/inverse.html