常微分方程初值问题的数值解法课件

常微分方程初值问题的数值解法§9.2

Euler方法

§9.3

Runge–Kutta公式§9.4

单步法的进一步讨论§9.5线性多步法§9.1引言数值算例8/1/20231常微分方程初值问题的数值解法§9.2Euler方法§9.定义：初值问题的单步显式方法，若对于任意固定的有近似解yn

满足极限

，则称该单步法收敛。

1收敛性定义：Remark：从定义可知，若格式收敛，整体截断误差en=y(xn)-yn必然趋于零。§9.4单步法的进一步讨论8/1/20232定义：初值问题的单步显式方法，若对于任意固定的2整体截断误差与局部截断误差的关系定理：若初值问题的单步方法之局部截断误差为且单步法中函数关于y满足lipschitz条件，则有8/1/202332整体截断误差与局部截断误差的关系定理：若初值问题的单步方二.

3相容性单步法局部截断误差：8/1/20234二.

3相容性单步法局部截断误差：7/28/20234由于假设为连续函数，因而上式可以表示为定义：如果当时，近似方程能逼近微分方程，则称数值公式与原微分方程相容。相容性定义8/1/20235由于假设为连续函数，因而上定义：如结论：若显式单步法的阶大于或等于1，则该单步法与微分方程相容；反之，如果单步法与微分方程相容，且关于h满足Lipschitz条件，则单步法至少为一阶方法。相容性与收敛阶的关系8/1/20236结论：若显式单步法的阶大于或等于1，则该单步法与微分方程相容定理：设增量函数在区域

中连续，并对变量y满足利普希茨条件，则单步法收敛的充要条件为相容性条件成立。Remark：在满足定理的条件下，Euler方法，Euler预估-校正格式，Runge-Kutta方法等都与原微分方程相容。4单步法收敛的条件8/1/20237定理：设增量函数在区域Remark：在满足定定义1用一个数值方法求解微分方程初值问题时，对给定步长h>0，若在计算时引入误差（也称扰动），由此引起计算后面的时误差绝对值均不增加，则称这个数值方法是绝对稳定的。

单步法收敛性概念以及定理都是在计算过程中无舍入误差的前提条件下建立的。5稳定性注：由稳定性定义可以看出方法是否稳定依赖于方程的右端函数，即方法是否稳定是指对于某个问题该方法是否稳定。8/1/20238定义1用一个数值方法求解微分方程初值问题时，对给定步长h针对模型方程研究稳定性设f(x,y)关于y满足Lipschitz条件，这样就可以针对如下模型方程研究方法的稳定性：其中λ为复常数，为使微分方程自身稳定，假定

定义1

设步长h>0的单步法用于求解模型方程时，中由引起的误差满足，则称单步法对于所用步长h和复数是绝对稳定的。若在计算时有误差,但在计算后面的8/1/20239针对模型方程研究稳定性设f(x,y)关于y满足LipschiRemark1：在上面的定义中，可以取小于或等于关系符。取小于号是为了和线性多步法相一致。Remark2：单步法是否稳定，与模型方程中的复数以及所用步长h有关。若对复平面上的某个区域G，当时，单步法绝对稳定，则称G为单步法的绝对稳定区域，G与实轴的交集为绝对稳定区间。关于模型稳定的说明8/1/202310Remark1：在上面的定义中，可以取小于或等于关系符。取小①Euler显式公式是保证绝对稳定性对步长h所加限制当为实数时，得到用h表示的绝对稳定的区间(-2,0)6常用公式的稳定性8/1/202311①Euler显式公式是保证绝对稳定性对步长h所加限制当为得绝对稳定区域②隐式Euler公式8/1/202312得绝对稳定区域②隐式Euler公式7/28/202312③梯形公式Back8/1/202313③梯形公式Back7/28/202313线性多步法的基本思想：如果充分利用前面多步的信息预测yn+k，则可期望获得较高精度。1线性多步法有关概念K步线性多步法一般形式为

其中为常数，不全为零。§9.5线性多步法8/1/202314线性多步法的基本思想：如果充分利用前面多步的信息预测yn+k若则为隐式方法，若则为显式方法Remark：R－K方法是增加一些非节点处的函数值提高单步法的精度，这样使计算量增加了许多。线性多步法每步只需要计算一个函数值。1线性多步法有关概念8/1/202315若则为隐式方法，若则为显式方法Re

对于隐式公式（），f(x,y)一般是非线性函数，故难以求解到yn+k的显示表达式，故常用迭代法求解：其中任意给出，s＝0,1,2,…,迭代到满足给定精度要求。可以证明，当f(x,y)满足Lipschitz条件或时，只要，迭代关系式就是收敛的。线性多步法有关概念（续）8/1/202316对于隐式公式（），f(x,y)一般是定义处的局部截断误差为

线性多步法2线性多步法局部截断误差若,则称线性多步法为p阶方法。8/1/202317定义处的局部截断误差为线性多步法2线性多步法局若线性多步法为p阶方法，则称为主局部截断误差系数。即主局部截断误差为8/1/202318若线性多步法为p阶方法，则关于局部截断误差定义的说明8/1/202319关于局部截断误差定义的说明7/28/202319利用微分中值定理

其中介于与之间。说明2故在的假定下，若

（显示公式），则：8/1/202320利用微分中值定理其中介于与之间。说若,且即为p阶方法，

则当y(x)充分可微时，说明38/1/202321若,且即为p阶方法，则当y即的首项与的首项相同，因此两种局部截断误差的定义相同。

Remark2：可以证明，显示线性多步法的整体截断误差比局部截断误差低一阶。Remark1：可以利用此处的截断误差定义分析前面的单步隐式方法。对于Euler方法，其主局部截断误差为，而对于梯形方法，其主局部截断误差为。说明48/1/202322即的首项与的首项相同，因此两种局部将方程两端从积分得

构造p次Lagrange插值多项式：3用数值积分法构造8/1/202323将方程两端从积分得构造其中

公式建立8/1/202324其中公式建立7/28/202324系数计算8/1/202325系数计算7/28/202325取可得到Adams显式公式

具体公式1：Adams显式公式

8/1/202326取可得到Adams显式公式具体公式1具体公式2：Adams隐式公式取k=0,j=1可得到Adams隐式公式再用n＋1代替n，得到8/1/202327具体公式2：Adams隐式公式取k=0,j=1可得到Adam取k＝1，j=1，得到Nyström

显式公式：具体公式3：Nyström显式公式8/1/202328取k＝1，j=1，得到Nyström显式公式：具体公式3：线性多步法

的局部截断误差为

局部截断误差8/1/202329线性的局部截断误差为局部截断误差7/28/202329局部截断误差8/1/202330局部截断误差7/28/202330对于Adams显式公式与隐式公式，由于显式(j=0,k=1)在[0,1]恒正,隐式(j=1，k=0):在[-1,0]恒负为某中间点，E(Explicit），I（Implicit）。Adams公式局部截断误差8/1/202331对于Adams显式公式与隐式公式，由于显式(j=0,k=1)当p=3时,

局部截断误差表明，在y(x)具有p+2阶连续导数的条件下，p+1步Adams显式方法与p步Adams隐式方法的局部截断误差是O(hp+2)，即它们是p＋1阶方法。特别地，当p＝3时，Adams显、隐方法都是四阶的。8/1/202332当p=3时,局部截断误差表明，在y(x)具有p+2阶连续导Taylor展开法更具一般性。例：用Taylor展开法构造下述公式，使其为四阶方法，并求其局部截断误差的主项。4用Taylor展开构造8/1/202333Taylor