我有一些矩阵值的复杂数据Z(f)和f∈{f0,f1,…}是一个 3x3 矩阵。我的工作流程中需要逆。在我的结果遇到一些问题后,我将问题缩小到反转,这似乎由 numpy.linalg.inv 执行得不太准确,而 scipy.linalg.inv 执行的范围也较小。Z(fi)Z−1(f)
的图中所见,在中间值应用 numpy 反转两次不会产生原始数据。在图中,scipy 反演似乎没有问题,但使用 scipy 计算也不会在的这些值下完全得到统一矩阵。Z12103<f<105ZZ−1Zf

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