在python中逆矩阵答案

【问题标题】：Inverse a matrix in python在python中逆矩阵
【发布时间】：2016-11-09 23:23:26
【问题描述】：

我想对一个矩阵 Q+1e-5*np.eye(d) (size 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]]

有人知道吗？

【问题讨论】：

你想反转一个奇异矩阵吗？
根据定义，你知道奇异矩阵没有逆矩阵，对吧？
（你的意思是你试图求解一个由奇异矩阵给出的联立线性方程组？）
我知道不可能逆奇异矩阵。但我不认为 Q+1e-5*np.eye(d) 也可以是单数。我只是想知道为什么结果是 NaN。

标签： python matrix-inverse

【解决方案1】：

根据定义，奇异矩阵没有逆矩阵。根据这个例子，numpy 的作者似乎选择给您一个未定义的结果，而不是在您尝试反转无法反转的矩阵时抛出错误。

【讨论】：