病态矩阵向量问题的正则化

计算科学 线性代数 线性求解器 条件数
2021-12-18 09:55:41

我有一个由积分微分系统引起的线性*问题,并写道:

(I+λA)x=b
在哪里A是一个真正的完整矩阵,大小n×n, 但不是对称的并且具有一些零特征值。我知道这些特征值对应什么。

增加(但仍远小于 1)时,的条件数(我们称之为)恶化得非常快(到和更多)。即使使用直接求解器(幸运的是对我来说就足够了),结果也变得毫无意义。λI+λAK1010n=100

我不确定,但我相信这与的零特征值有关。有没有办法可以利用我对相应特征向量的知识来规范数值分辨率?A

(*注意,这实际上是我正在解决的定点线性问题,原问题有λ=λ(x)

编辑: 本身确实有一个非常大的条件数(由 numpy 用 L2-norm 进行数值评估), (当然应该解析为 0,因为如上所述,的特征值为 0)。它的最大特征值的模数约为 600,这里是A1019σmin=1.8e07σmax=3.2e11σminAn=102

在此处输入图像描述

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