关于极限和正交积分规则的两个摄动矩阵的条件数

计算科学 条件数 调理
2021-12-08 03:18:38

我有一个关于两个不同扰动矩阵的条件数的问题。从让开始A是一个 spd 矩阵,其元素定义为ai,j=ΩRdφiφjdΩ出现在有限元方法中φi近似空间的基函数。此外让cond2(A)是条件数A.

  1. 我的第一个问题是关于矩阵元素的数值积分以及它们如何影响条件数。我的想法是写A~=A+E其中的元素E说明数值积分和精确积分之间的误差。有什么办法可以给出不等式关系cond2(A~)cond2(A)?

  2. 我的第二个问题是关于条件数内的限制。假设一个矩阵E有元素ei,j=Ω~Rdφ~iφ~jdΩΩ~Ω=. 条件数如何对矩阵的限制做出反应,即它是否成立cond2(A+limα0αE)=cond2(A)或者这两者之间是否存在任何不平等关系?我首先尝试使用

    cond2(A+limα0αE)cond2(A)+cond2(limα0αE)
    但我意识到零矩阵的条件数定义为. 因此,不平等是没有用的。

1个回答

1.对于矩阵A具有不同的特征值,添加扰动δA结果 [1] 改变大小的特征值(到一阶)

δλi=(X1δAX)ii,
所以|δλi|κ(X)δA. 条件数的变化由 其中是绝对值的最小和最大特征值。成正比 ,你能得到的唯一问题是如果你的矩阵几乎是奇异的并且扰动变成奇异的,那么会很大。κ(A)
δ(κ(A))=κ(A)(δ(logλn)δ(logλ1))κ(A)δA(|λn|1+|λ1|1),
λ1,λnδAδ(κ(A))

2.是不正确的考虑,那么两者都是良条件的,,但κ(A+B)κ(A)+κ(B)A=diag(1,1+ϵ)B=diag(1,1)κ(A)=1+ϵκ(B)=1κ(A+B)=ϵ1

  • [1]用于正向和反向模式自动微分的矩阵导数结果的扩展集合,M. Giles,(pdf)