数值分析8.48.5(免费阅读).ppt

华长生制作 方法收敛的条件 单步法的稳定性 Over! * 第八章 常微分方程数值解 Solution of differential equation 在前面所讨论的方法中,在计算 时只用到前一步 的信息(单步法)，为提高截断误差的阶，每个时间步必须增加计算右端函数 的次数。当 的结构比较复杂时，计算量较大。现在指出另一个提高截断误差阶的办法，即构造这样的方法: 在计算 公式中,充分利用前几步得到的信息 及 ，但每进一步，只计算一次 的值。这样的方法称为多步方法，若函数值 以线性组合的形式出现于公式中，则称方法为线性多步方法。 线性多步方法 Linear multistep methods 初值问题: -----------(1) 称为Euler二步法. 一个特例 A special case : 1. (l 步) 线性多步法的一般格式为：A general form ---------(2) 当 ??1 ? 0 时，为隐式公式; implicit schemes ??1 = 0 则为显式公式。 Explicit schemes 系数 及 可根据方法的局部截断误差及阶确定，其定义如下. 线性多步方法的构造 The construction of linear multistep methods 构造多步法有多种途径，常用的有基于Taylor展开 的构造方法(待定系数法)和基于数值积分的构造方法。 根据公式要达到的精度(p 阶),即 ,确定 基于Taylor展开的构造方法: 具体做法: (参考 P224 例8.4.1) 选择系数 使得这格式的精度 达到p阶，则满足 3. 几个重要的线性多步法 Adams(阿当姆斯)方法, Milne(米纳)方法, Hamming(哈明)方法, Simpson(辛普生)方法. 基于数值积分的构造法 将 在 上积分，得到 只要近似地算出右边的积分 ，则可通过 近似y(xn+1) 。而选用不同近似式 Ik，可得到不同的计算公式。 多步方法的特点： (1) 因初始条件只有一个，多步方法的启动要借助 高阶的单步方法来开始. (2) 多步方法比较简单，只要在这几个点的函数 值的线性组合, 而且每步中所用函数值, 有些下一 步还可使用。 8.5 收敛性与稳定性分析 用上式的差分方程来逼近微分方程的初值问题是否 合理,就要看差分方程的解是否收敛到初值问题的 精确解. -----------(1) 定义8.5.1 对于欧拉方法，由于其增量函数 就是 ，故当 关于 满足利普希茨条件时它是收敛的. 再考察改进的欧拉方法，其增量函数由 假设 关于 满足利普希茨条件，记利普希茨常数为 , 设 为定数）， 上式表明 关于 的利普希茨常数 则由上式推得 因此改进的欧拉方法也是收敛的. 假设 关于 满足利普希茨条件，记利普希茨常数为 , 设 为定数）， 上式表明 关于 的利普希茨常数 因此改进的欧拉方法也是收敛的. 定义定义定义定义 若一种数值方法在节点值 上大小为 的扰动， 于以后各节点值 上产生的偏差均不超过 ，则 称该方法是稳定的. 定义 我们考虑一种简单情况，即仅初值有误差，而其他计算 步骤无误差。 设 是初值有误差后的计算值，则 则对于显式 Euler 公式 可以看出，当初始误差充分小，以后各步的误差也充分小。 有： 例： --------- (2) 试验方程： (其中 ) 设在节点值 上有一扰动值 ，它的传播使节点值 产生大小为 的扰动值， 假设用 按欧拉公式得出 的计算过程不再有新的误差， 则扰动值满足 可见扰动值满足原来的差分方程. 如果差分方程的解是不增长的，即有 则