在Python中反轉矩陣

Question

我想逆矩陣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]] 

不要任何人有任何想法？

Answer 1

根據定義，奇異矩陣沒有反。基於這個例子，看起來像numpy的作者選擇給你一個未定義的結果，而不是在你嘗試顛倒一個不能顛倒的矩陣時拋出錯誤。