第9章5-7节多步方程组与高阶方程.pptVIP

§5 线性多步法 1 线性多步法的一般公式 2 阿当姆斯显式与隐式公式 3 米尔尼方法与辛普森方法 4 汉明方法 5 构造多步法公式的注记和例  §6 线性多步法的相容性，收敛性与稳定性 1.相容性 与单步法一样，我们定义若局部截断误差Tn+k的价p≥1,则称k步法与微分方程（1.1）相容.   定理5 线性多步法(5.1)式与微分方程（1.1） 相容的充分必要条件是   2.收敛性  定义9  定理6 若线性多步法是收敛的，则它一定 是相容的.  例9 用线性二步法  但是方法（6.6）的解并不收敛. 在方法（6.6）中取 y0=0,y1=h 此时方法（6.6）化为  用归纳法可以证明  定义10 若线性多步法（5.1）式的第一特征多项式的根都在单位圆内，或单位圆上，且在单位圆上的根为单根，则称线性多步法（5.1）满足根条件.  稳定性与绝对稳定性 1）稳定性 设线性多步法（6.5）有扰动 经过扰动后的解为Zn,n=0,…,N,它满足  定义11 对初值问题（1.1）（1.2),由方法（6.5）得到的解为yn,n=0,…,N由于有扰动 使得方程（6.9）的解为zn,n=0,…,N,若存在常数 C及h0,对所有h∈(0,h0),当 则称线性多步法（5.1）是稳定的。  定理8 线性多步法（5.1）稳定的充分必要条件是它满足根条件。  2）绝对稳定性 定义12 对于给定的 ，若 的零点满足 则称线性多步法（5.1）关于此 值是绝对稳定的。 若在 的复平面的某个区域R中对所有 值线性多步法（5.1）都是绝对稳定的，在R外，方法是不稳定的，则称R为 线性多步法（5.1）的绝对稳定域。 R与实轴的交集称为线性多步法（5.1）的绝对稳定区间。 §7 方程组和高阶方程 2 化高阶方程为一阶方程组 3 刚性方程组 （6.1） 对多步法 （5.1） 引入多项式 分别称为多步法（5.1）的第一特征多项式 和第二特征多项式. （6.2） （6.3） （6.4） 由于用k步法(5.1)需要k个初值，而微分方程（1.1） 只给出一个初值y(x0)=y0,因此还要给出k-1个初值才能用 k步法(5.1)，既需要由相应的单步法求出y1,…,yk-1. （6.5） 注意：相容的线性多步法不一定收敛. （6.6） 方法（6.6）的解并不收敛. 定理7 线性多步法（5.1）是相容的， 则线性多步法（6.5）收敛的充分必要条件是 线性多步法（5.1）满足根条件. 在例9中， 不满足根条件，因此方法（6.6）的解并不收敛. （6.9） 1 一阶方程组 前面研究了单个方程 的数值解法，只要把 和 理解为向量，那么，所提供的各种计算公式即可应用 到一阶方程组的情形. 考察一阶方程组 的初值问题，初始条件给为 则上述方程组的初值问题可表示为 （7.1） 求解这一初值问题的四阶龙格-库塔公式为 式中 若采用向量的记号，记 或表示为 其中 这里 是第 个因变量 在节点 的近似值. 考察两个方程的特殊情形： 这时四阶龙格-库塔公式具有形式 其中 （7.2） （7.3） 这是一步法，利用节点 上的值 ，由（7.3） 式顺序计算 ，然后代入（7.2） 式即可求得节点 上的 . （7.3） 高阶微分方程（或方程组）的初值问题，原则上总可 以归结为一阶方程组来求解. 例如，考察下列 阶微分方程 （7.4） 初始条件为 （7.5） 只要引进新的变量 即可将 阶方程（7.4）化为如下的一阶方程组： （7.6） 初始条件（7.5）则相应地化为 （7.7） 初值问题（7.4）,（7.5）和（7.6）,（7.7）是彼此等价的. （7.4） （7.5） 在逐步推进的求解过程中，计算 之前事实上已经 求出了一系列的近似值 如果充分利用前面多步的信息来预测 ，则可以期望获得较高的精度. 这就是构造所谓线性多步法的基本思想. 本节主要介绍基于泰勒展开的构造方法. 构造多步法的主要途径是基于数值积分方法和基于泰 勒展开方法，前者可直接