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

作者: ray0215 (OKOK It's good!)   2016-06-22 17:37:51
小弟我最近需要計算相關係數(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]])
每個都十五次方
我想應該不是正確解答吧@@
請問要怎麼解決
先謝謝版上的大大了
作者: bigpigbigpig (To littlepig with love)   2016-06-22 18:26:00
此矩陣的行列式 = 3.19867e-10,李組長眉頭一皺...
作者: alen332l (alen3321)   2016-06-22 19:40:00
先確定你是否真的需要反矩陣?如果是這樣:Ax = b 已知A,b求x,不一定要用反矩陣可否提供更宏觀的題目?
作者: s06yji3 (阿南)   2016-06-22 20:18:00
行列式為0沒有反矩陣。1e-10可視為0了吧@@
作者: starcloud (傷心嗎?)   2016-07-12 04:17:00
這是奇異矩陣 所以請用 pinv來解

Links booklink

Contact Us: admin [ a t ] ucptt.com