第1章 常微分方程两点边值问题的差分方法
有限差分方法是用于求解微分方程定解问题的*广泛的数值方法,其基本思想是用离散的只含有有限个未知量的差分方程组去近似代替连续变量的微分方程和定解条件,并把差分方程组的解作为微分方程定解问题的近似解.常微分方程两点边值问题可以看成一维椭圆型方程定解问题,模型简单.本章研究此模型问题的差分解法,介绍微分方程数值解法中的一些基本概念、差分格式的极值原理分析方法和能量分析方法,以及提高数值解精度的Richardson外推法.
1.1 Dirichlet边值问题
考虑如下定解问题:
(1.1a)
(1.1b)
其中q(x).0,f(x)为已知函数,α和β为已知常数.
当q(x)≡0时,在方程(1.1a)中用s代替x,并在两边关于s从0到x积分一次,得到
在上式中用ξ代替x,再在两边关于ξ从0到x积分一次,并应用左边界条件u(0)=α,得到
再应用右边界条件u(L)=β,可得
因而(1.1)的解可表示为
要想求出某点处的值还需要借助于数值积分.当时,用同样的方法要想得到解的精确表达式是困难的,甚至是办不到的.读者可对q(x)≡1的情形试一试.
尽管难以求出精确解,但我们可以设法给出解的估计式.
1.1.1 基本微分不等式
本书中Cm[0,L]表示闭区间[0,L]上所有具有m阶连续导数的函数的集合.
设函数u∈C[0,L].记
如果函数u∈C1[0,L],则进一步记
引理1.1(I)设函数,则有
(1.2)
(II)设函数v∈C2[0,L],且v(0)=0,v(L)=0,则有
(1.3)
(III)设函数v∈C1[0,L],且v(0)=v(L)=0,则有
(IV)设函数v∈C1[0,L],且v(0)=v(L)=0,则对任意.>0有
(1.4)
(V)设函数v∈C1[0,L],则对任意.>0有
(1.5)
证明 (I)由分部积分直接可得(1.2).
(II)由(1.2)易得(1.3).
(III)对于任意的x∈(0,L),有
(1.6)
(1.7)
将(1.6)和(1.7)两端平方并应用Cauchy-Schwarz不等式,得到
(1.8)
(1.9)
将(1.8)乘以L.x,将(1.9)乘以x,并将结果相加,得
(1.10)
注意到当x∈(0,L)时,
由(1.10)易得
将上式两边开方,得
易知
对(1.10)式两端关于x积分,得
两边开方得
(IV)对任意.>0,有
将以上两式相加并除以2,得到
因而(1.4)成立.
(V)设x∈[0,L]使得
当y∈[x,L]时,
由以上两式得到
将上式两边关于y从0到L求积分,得到
易得
因而(1.5)成立.
引理证毕.
1.1.2 解的先验估计式
我们给出齐次边值问题解的先验估计式.
定理1.1设函数v∈C2[0,L]为两点边值问题
(1.11a)
(1.11b)
的解,其中q(x).0,则有
(1.12)
(1.13)
证明(I)将(1.11a)两端同乘以v(x),并关于x在(0,L)上积分,得
(1.14)
注意到(1.11b),由引理1.1有
由q(x).0,有
此外,应用Cauchy-Schwarz不等式,有
将以上三式代入(1.14),得
再次应用引理1.1,有
于是
(II)注意到
由(1.12)及引理1.1,得
定理证毕.
称(1.12)和(1.13)为两点边值问题(1.11)解的先验估计式.
展开