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

*  §6 线性多步法的相容性，收敛性与稳定性 1.相容性 与单步法一样，我们定义若局部截断误差Tn+k的价p≥1,则称k步法与微分方程（1.1）相容. （6.1） *  对多步法 （5.1） 引入多项式 分别称为多步法（5.1）的第一特征多项式 和第二特征多项式. （6.2） （6.3） *  定理5 线性多步法(5.1)式与微分方程（1.1） 相容的充分必要条件是 （6.4） *  *  2.收敛性 由于用k步法(5.1)需要k个初值，而微分方程（1.1） 只给出一个初值y(x0)=y0,因此还要给出k-1个初值才能用 k步法(5.1)，既需要由相应的单步法求出y1,…,yk-1. （6.5） *  定义9 *  定理6 若线性多步法是收敛的，则它一定 是相容的. 注意：相容的线性多步法不一定收敛. *  例9 用线性二步法 （6.6） *  但是方法（6.6）的解并不收敛. 在方法（6.6）中取 y0=0,y1=h 此时方法（6.6）化为 *  用归纳法可以证明 方法（6.6）的解并不收敛. *  定义10 若线性多步法（5.1）式的第一特征多项式的根都在单位圆内，或单位圆上，且在单位圆上的根为单根，则称线性多步法（5.1）满足根条件. 定理7 线性多步法（5.1）是相容的， 则线性多步法（6.5）收敛的充分必要条件是 线性多步法（5.1）满足根条件. 在例9中， 不满足根条件，因此方法（6.6）的解并不收敛. *  稳定性与绝对稳定性 1）稳定性 设线性多步法（6.5）有扰动 经过扰动后的解为Zn,n=0,…,N,它满足 （6.9） *  定义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 方程组和高阶方程 1 一阶方程组 前面研究了单个方程 的数值解法，只要把 和 理解为向量，那么，所提供的各种计算公式即可应用 到一阶方程组的情形. 考察一阶方程组 的初值问题，初始条件给为 * 则上述方程组的初值问题可表示为 （7.1） 求解这一初值问题的四阶龙格-库塔公式为 式中 若采用向量的记号，记 * 或表示为 其中 * 这里 是第 个因变量 在节点 的近似值. 考察两个方程的特殊情形： * 这时四阶龙格-库塔公式具有形式 其中 （7.2） （7.3） * §5 线性多步法 在逐步推进的求解过程中，计算 之前事实上已经 求出了一系列的近似值 如果充分利用前面多步的信息来预测 ，则可以期望获得较高的精度. 这就是构造所谓线性多步法的基本思想. 本节主要介绍基于泰勒展开的构造方法. 构造多步法的主要途径是基于数值积分方法和基于泰 勒展开方法，前者可直接由方程 两端积分后 利用插值求积公式得到. * 1 线性多步法的一般公式 一般的线性多步法公式可表示为 （5.1） 其中 为 的近似， 计算时需先给出前面 个近似值 ,再由 （5.1）逐次求出 . 如果计算 时，除了使用 的值，还用到 的值，则称此方法为线性多步法. 为常数. 及 不全为零， 则称为线性 步法. 若 * 如果 ，称（5.1）为显式 步法，这时 可 直接由（5.1）算出. 如果 ，则（5.1）称为隐式 步法. 求解时与梯形法相同，要用迭代法方可算出 . （5.1） * 设 是初值问题(1.1)，(1.2)的准