理论教育 差分方程收敛性-偏微分方程数值解法:土建类

差分方程收敛性-偏微分方程数值解法:土建类

时间:2023-10-26 理论教育 版权反馈
【摘要】:差分方程解的收敛性问题,简单地说,就是当步长h无限缩小时,差分方程的解是否逼近微分方程的解。为讨论差分方程解的收敛性问题先引入两个概念。定义一 称由Rj=Ly-Lhy 所确定的Rj为差分方程的截断误差。下面我们对微分方程式和差分方程式进行讨论。

差分方程收敛性-偏微分方程数值解法:土建类

差分方程解的收敛性问题,简单地说,就是当步长h无限缩小时,差分方程的解是否逼近微分方程的解。这个问题是有限差分法中的一个重要问题,若不能保证当h→0时,εj=yxj)-yj→0,则所得的差分解没有实用价值。为讨论差分方程解的收敛性问题先引入两个概念。

定义一 称由

Rjy=Lyxj)-Lhyxj) (1.2.22)

所确定的Rjy)为差分方程的截断误差。

定义二 若存在一个与步长h无关的常数K,使不等式

|Rjy)|≤Khp

成立,即

Rjy=Ohp

则称Lh具有p阶精度。

下面我们对微分方程式(1.2.17)和差分方程式(1.2.19)进行讨论。

y′的差分表达式(1.2.8)和y″的差分表达式(1.2.10)代入,得

其中,xj-1ξxj+1xj-1ηxj+1。显然,当h→0时,Rjy)→0。事实上,只要yx)在[ab]上具有直到二阶的连续导数,则由微分中值定理容易推得

于是,对于[ab]上具有二阶连续导数的任何函数,差分算子Lh逼近微分算子L,且由式(1.2.23)可知,在yx)具有四阶连续导数时,Lh就有二阶精度。考虑到

其中,ajbjcj由式(1.2.20)给出。

若令

根据式(1.2.20),考虑到

为了推导方便,把条件qx)≥0改成(www.daowen.com)

qx)≥k>0,x∈[ab]

bj=2+h2qj,得

于是

Rjy)→0,得ε→0。由此可以得出,对于方程式(1.2.17),若该方程的解yx)在[ab]上具有二阶连续导数,且px),qx)和rx)都是在[ab]上的连续函数,qx)≥k>0,则由差分方程(1.2.19)所确定的差分解必收敛于微分方程的解,即

进一步,若方程式(1.2.17)的解yx)在[ab]上具有四阶连续导数,则对于式(1.2.19),有

yxj)-yj=Oh2

此即对差分解的误差估计。

事实上,若令

则由式(1.2.23),可得

由式(1.2.24)得

yxj)-yj=Oh2

上面考虑的边值问题为第一类边界条件,对于第二类和第三类的边界条件,因其中包含导数,必须用适当的有限差分格式来表示导数。例如,若要求误差项是Oh),则有

若要求误差项是Oh2),则有

不难看出,在以上两种情况中,对于前一种情况,差分方程的系数矩阵是三对角的,而后一种情况,作简单的变形,也可将差分方程的系数矩阵化为三对角矩阵。

免责声明:以上内容源自网络,版权归原作者所有,如有侵犯您的原创版权请告知,我们将尽快删除相关内容。

我要反馈