[問題] 關於反矩陣inv計算爆值是正常的嗎?

[ray0215 (❺❺❻❻)]

小弟我最近需要計算相關係數(correlation coefficient matri)矩陣的反矩陣

例如

     Martrix   array([[ 1.        ,  0.92531576,  0.81998247,  0.92586894],
                      [ 0.92531576,  1.        ,  0.70974912,  0.72358554],
                      [ 0.81998247,  0.70974912,  1.        ,  0.72669914],
                      [ 0.92586894,  0.72358554,  0.72669914,  1.        ]])

這樣用numpy.linalg.inv算出來的值

會計算出

array([[  6.51500296e+15,  -3.18743561e+15,  -7.89299760e+14,
         -3.15207310e+15],
       [ -3.18743561e+15,   1.55943840e+15,   3.86161323e+14,
          1.54213745e+15],
       [ -7.89299760e+14,   3.86161323e+14,   9.56245323e+13,
          3.81877116e+14],
       [ -3.15207310e+15,   1.54213745e+15,   3.81877116e+14,
          1.52502845e+15]])


每個都十五次方
我想應該不是正確解答吧@@

請問要怎麼解決

先謝謝版上的大大了


--
※ 發信站: 批踢踢實業坊(ptt.cc), 來自: 140.112.63.249
※ 文章網址: https://www.ptt.cc/bbs/Python/M.1466588274.A.BAD.html
→ bigpigbigpig: 此矩陣的行列式 = 3.19867e-10，李組長眉頭一皺...
好像真的可能這樣欸

推 alen332l: 先確定你是否真的需要反矩陣?
→ alen332l: 如果是這樣：Ax = b 已知A,b求x，不一定要用反矩陣
→ alen332l: 可否提供更宏觀的題目?
我要計算R^(-1)-I 的值
R = correlation matrix
I = idendity matrix
照理說R^(-1)和I兩邊的值不會差太大@@
就因為反矩陣的值和單位矩陣差太大 後面就不用計算了





※ 編輯: ray0215 (140.112.63.249), 06/22/2016 20:01:18
推 s06yji3: 行列式為0沒有反矩陣。1e-10可視為0了吧@@
推 starcloud: 這是奇異矩陣 所以請用 pinv來解