在python中反转矩阵

Question

我想要反转矩阵Q + 1e-5 * np.eye（d）（大小d X d）并使用以下代码来获得近似结果。

Q = X.dot(X.T) # X is a large sparse matrix, Q is singular
P = np.linalg.inv(Q+1e-5*np.eye(d))

但我得到了这个：

P=[[ nan  nan  nan ...,  nan  nan  nan]
  [ nan  nan  nan ...,  nan  nan  nan]
  [ nan  nan  nan ...,  nan  nan  nan]
   ...,
  [ nan  nan  nan ...,  nan  nan  nan]
  [ nan  nan  nan ...,  nan  nan  nan]
  [ nan  nan  nan ...,  nan  nan  nan]]

有人有任何想法吗？

Answer 1

根据定义，奇异矩阵没有逆矩阵。根据这个例子，当你试图反转一个不能被反转的矩阵时，看起来像numpy的编写者选择给你一个未定义的结果而不是抛出错误。