特征值分析与傅立叶分析的稳定性及其等价性

计算科学 pde 有限差分 抛物线pde
2021-11-29 22:04:39

我有一个关于常数系数 pde 稳定性分析的问题。假设我正在查看 pde

ut=auxx
所以第一种方法是计算放大因子,这就是所谓的傅里叶方法。称它为Q(ζ)以及用于插入测试功能的显式方法eijhw我有ζ=hw
Q(ζ)=14aτ/h2sin2(ζ/2)
另一方面,矩阵的特征值分析I+τA, 在哪里A是的离散近似auxx给出矩阵的特征值I+τA与放大系数完全相同Q(ζ).

它们是否总是相同的,这些方法是否等效?所以,我可以计算任何一个,并确保它是有界的1?

第二个问题。自从A是对称的,这意味着maxλi=||I+τA||2. 但是,如果我使用非均匀网格,那么max|λi|||I+τA||2因此两者都不能保证稳定性,因为我只能证明max|λi|1. 那么在非均匀网格的情况下,稳定性的充分条件是什么?

1个回答

对问题 1 的回答:

如果A循环矩阵,则其特征向量的形式为vj,w=eijhw; 即,它们是离散傅里叶模式因此,当且仅当您处理循环矩阵时,傅立叶分析(或通常称为冯诺依曼分析)和特征值分析是等价的。

回答问题 2:

最大特征值模和欧几里得范数相等,如果A是一个正规矩阵如果A是非正态的,那么欧几里得范数可以更大,如你所说。稳定性的充分条件当然是I+τA21. 还有其他充分条件,具体取决于您所指的稳定性。PDE 离散化最重要的稳定性是Lax-Richtmeyer 稳定性,在您的情况下,这意味着 对于所有,其中是一个常数,可能取决于但不取决于有关讨论,请参见Leveque 文本的第 9 章。

(I+τA)n<C
nτ<TCTnτ

证明非均匀网格离散化的这种稳定性可能很困难。在实践中,通常考虑基于均匀网格的稳定性标准就足够了,例如并应用它,但将替换为的最小值整个网格。Δt(Δx)2/2ΔxΔx