抛物型方程的差分方法.pptxVIP

微分方程数值解法

第2章抛物型方程的差分方法2.1差分格式建立的基础2.2显式差分格式2.3隐式差分格式2.4解三对角形方程的追赶法2.5差分格式的稳定性和收敛性2.6非线性抛物型方程的差分解法举例2.7二维抛物型方程的差分格式2.8交替方向的隐式差分格式(ADI格式)

本章，我们研究线性抛物型方程的差分解法，主要讨论差分方程的构造方法和有关的理论问题以及研究方法等，重点在于一维线性抛物型方程的差分方法，对于非线性以及多维抛物型方程的差分解法也进行了研究。其中,为平面上某一区域。(2.1)众所周知，一维线性抛物型方程的一般形式为

(2)初边值问题(或称混合问题)通常考虑的定解问题有：(1)初值问题(或称Cauchy问题)在区域上求函数，使满足(2.2)为给定的初始函数。(2.3)(2.4)在区域上求函数，使满足边值条件初值条件

2.1差分格式建立的基础为了构造微分方程(2.1)的有限差分逼近，首先将求解区域用二组平行于轴和轴的直线构成的网格覆盖，网格边长在方向为，在方向为(如图2.1所示)。分别称为空间方向和时间方向的步长，网格线的交点称为网格的结点。对初值问题来说，网格是在上的结点称为边界结点，属于内的结点称为内部结点。

对于初边值问题，设，则网格是研究导数的差商近似表达式。为此对二元函数定义，且假定具有我们需要的有界偏导数。在上的结点称为边界结点，属于内的结点称为内部结点。差分方程就是在网格点上求出微分方程解的近似值的一种方法，因此又称为网格法。构造逼近微分方程的差分方程的方法。

向后差商(2.7)(2.5)中心差商向前差商则在处对的一阶偏导数有三个可能的近似:(2.6)由Taylor展开，有

显然，用差商近似导数存在误差，令(2.8)则关于导数的近似差商表达式，也可以通过线性算子作为推导工具得到，定义：截断误差,阶为用向后差商近似导数的截断误差阶也为而中心差商近似导数的截断误差阶为为方向偏导数算子

为方向位移算子，为方向平均算子，其中:方向的差分算子:(2.9)前差算子:,(2.10)后差算子:,中心差算子:(2.11),,

建立差分算子和导数算子之间的关系，由Talyor展开，有由得(2.12)或者(2.13)同理有

因为故(2.14)同理(2.15)因为(2.16)则(2.17)式(2.14)，(2.15)，(2.17)分别给出了偏导数算子关于前差、后差、中心差的级数表达式双曲正弦3246

(2.18.1)(2.18.2)(2.18.3)利用这些关系式就可给出偏导数的差分表达式返回

又由可得二阶偏导数的差分表达式(2.19.1)(2.19.2)(2.19.3)返回返回

对于三阶、四阶偏导数的差分表达式为(2.20.1)(2.20.2)(2.20.3)(2.21.1)(2.21.2)(2.21.3)

从以上这些偏导数的差分表达式，我们可以得到偏导数的各种精度的近似表达式。且又由二阶导数的前差表达式(2.19.1)，得因此在的前差表达式中取第一项，则有即截断误差阶为。

现在研究构造微分方程(2.1)的差分方程的方法，为此记微分方程(2.1)为(2.22)L是关于的线性算子，。包括二个相邻时间层的网格结点的差分方程可以从Talor展开式推出返回

设，于是(2.23)如果算子L不依赖于t，即