计算差分方程的收敛点_PDE有限差分方法(4)——稳定性

今天讲差分格式的稳定性，主要讲【初值稳定性】. 关于它的一些基本介绍就不说了, 因为这是学习笔记而不是用来讲课的notes，只是用来回顾一下主要内容.

参考书：
(1) J.W. Thomas - Numerical Partial Differential Equations_ Finite Difference Methods (1995, Springer)
(2) 张强《偏微分方程的有限差分方法》科学出版社，2019年1月版。
(3) K. W. Morton, D. F. Mayers - Numerical solution of partial differential equations (2005, Cambridge University Press)

本次的主要内容：

定义 [初值稳定性] 考虑齐次线性差分格式

(

可以和时间层数与网格函数有关; 若与网格函数有关, 则PDE是非线性的)

给定离散范数

当

趋于零时, 若数值解满足有界性

其中界定常数

与

和

均无关, 则称差分格式按

具有初值稳定性.

注：若K与T有关, 则称差分格式具有短时间的数值稳定性; 若定界常数与T无关, 则称差分格式有长时间的初值稳定性.

例1 当且仅当

时, 模型问题(HP)的全显格式具有最大模初值稳定性, 此时, 离散最大模原理成立.

证明：全显格式是

它的系数和等于1, 且

时系数非负. 那么

于是

另一方面, 当

时, 右端出现负系数, 凸组合的系数结构不复存在. 反例: 考虑如下问题：

那么

进一步可推出

由于

则数值解趋于无穷, 全显格式不稳定. QED

例2 模型问题(HP)的全隐格式无条件具有最大模的初值稳定性.

证明：全隐格式的差分方程是

均具有隐式的凸组合系数结构, 即等号右端系数都是正的, 左端系数不低于右端系数之和. 这样类型的数值格式满足离散最大模原理. 这是因为:

<