lapack 奇异矩阵

计算科学 拉帕克 矩阵
2021-12-10 04:49:38

我想找到一个条件来确定矩阵是否可逆。天真地,我计算了行列式,看看它是否为零。但后来我意识到即使对于非常小的值(1e-100),lapack 仍然能够找到一个反转矩阵。

我在这里读到分解可用于计算矩阵的秩。QR

分解是检查矩阵是否可逆的最佳方法吗?它会考虑数字错误吗?有没有更有效的方法?QR

我还发现了这个链接,它给出了一个倒数条件数这个数字是什么意思?我知道如果那么矩阵是奇异的!但是截止点是什么,对于计算机κκ=

1个回答

LU 分解仅在矩阵可逆(不是 LUP)时存在。在 Lapack 中执行此操作的常用方法是尝试获得 LU 分解。许多 Lapack 库会给你一个信息代码,说明失败的原因。如果您的矩阵不可逆 LU 分解输出“信息”将不是 0,您可以检查此标志来确定您的矩阵是否可逆。