作者ray0215 (❺❺❻❻)
看板Python
标题[问题] 关於反矩阵inv计算爆值是正常的吗?
时间Wed Jun 22 17:37:51 2016
小弟我最近需要计算相关系数(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://webptt.com/cn.aspx?n=bbs/Python/M.1466588274.A.BAD.html
1F：→ bigpigbigpig: 此矩阵的行列式 = 3.19867e-10，李组长眉头一皱... 06/22 18:26
好像真的可能这样欸
2F：推 alen332l: 先确定你是否真的需要反矩阵? 06/22 19:40
3F：→ alen332l: 如果是这样：Ax = b 已知A,b求x，不一定要用反矩阵 06/22 19:41
4F：→ alen332l: 可否提供更宏观的题目? 06/22 19:41
我要计算R^(-1)-I 的值
R = correlation matrix
I = idendity matrix
照理说R^(-1)和I两边的值不会差太大@@
就因为反矩阵的值和单位矩阵差太大 後面就不用计算了
※ 编辑: ray0215 (140.112.63.249), 06/22/2016 20:01:18
5F：推 s06yji3: 行列式为0没有反矩阵。1e-10可视为0了吧@@ 06/22 20:18
6F：推 starcloud: 这是奇异矩阵 所以请用 pinv来解 07/12 04:17