麻烦用 numpy 和 scipy 反转复杂矩阵

计算科学 线性代数 矩阵 scipy 麻木的
2021-12-28 13:48:06

我有一些矩阵值的复杂数据Z(f)f{f0,f1,}一个 3x3 矩阵。我的工作流程中需要逆在我的结果遇到一些问题后,我将问题缩小到反转,这似乎由 numpy.linalg.inv 执行得不太准确,而 scipy.linalg.inv 执行的范围也较小。Z(fi)Z1(f)

的图中所见,在中间值应用 numpy 反转两次不会产生原始数据在图中,scipy 反演似乎没有问题,但使用 scipy 计算也不会在的这些值下完全得到统一矩阵。Z12103<f<105ZZ1Zf

Z12

条件数不同,如下面的第二张图所示,但我想了解为什么我的问题会出现在中间值处,以及如何选择一种能够最小化此类问题的反演算法。Z(f)f0f

条件编号

0个回答
没有发现任何回复~