第九章常微分方程数值解法初次修改稿1.ppt

第九章 常微分方程数值解法 再离散化，即可得如下二阶Adams显式公式 ? 截断误差 ? 二阶Adams隐式公式: 梯形公式 再离散化，即可得如下计算公式 F(x)用过节点 xn , xn+1的线性插值函数代替 ? 三阶Adams显式公式: F(x)用过节点 xn , xn?1 , xn ?2的2次插值函数 L2(x)代替 再离散化，得三阶Adams显式公式 ? 三阶Adams显式公式截断误差 ? 三阶Adams隐式公式: F(x)用过节点 xn+1 , xn , xn ?1的2次插值函数 L2(x)代替 再离散化，得三阶Adams隐式公式 ? 三阶Adams隐式公式截断误差 ? 四阶Adams显式公式 截断误差 ? 四阶Adams隐式公式 截断误差 ? Milne公式 F(x)用过节点 xn , xn?1 , xn ?2的2次插值函数 L2(x)代替 再离散化，得Milne公式 Milne公式是四阶四步显式方法, 截断误差为 ? 基于Taylor展开方法能导出所有的常用线性多步法, 下面以线性二步法为例来说明这种方法 ? 线性二步法的一般格式 ? 适当选取参数?0, ?1, ??1, ?0, ?1, 使其局部截断误差达到最高阶(使公式的阶尽可能高) ? 局部截断误差 : 假设第n步以前的计算是精确的, yi=y(xi) ( i ? n), 用多步法计算一步所产生的误差, 即 y(xn+1) 与 yn+1 的差. 为符号简单记, 把 简记为 由Taylor展开式得 ? 局部截断误差 要使局部截断误差的阶高, Rn+1的前面几项应尽可能多为0, 因只有5个未知数, 可令前5项为0. 其中 fi=f (xi , yi), i=n ? 1, n, n+1. Simpson公式 两步四阶隐式方法 局部截断误差 ? 答案: 该公式为(r+1)阶Adams显式公式 ? 适当选取?0, ?1, …, ?r 使下列公式的阶尽可能高 适当选取?0, ?1使公式的阶尽可能高. 例 r=1的情形 局部截断误差 得二阶Adams公式 要使阶尽可能高, 须令 局部截断误差 ? 答案: 该公式为(r+1)阶Adams隐式公式. ? 适当选取 ??1, ?0, ?1, …, ?r?1 使下列公式的阶尽可能高 Hamming公式 四阶三步隐式公式 ? 适当选取?0, ?2, ??1, ?0, ?1, 使下列公式的阶尽可能高 ? 答案: 所求公式为 截断误差 例 分别用四阶Adams显式和隐式公式求解初值问题 的数值解, 取 h=0.1. ? 利用线性多步法求解初值问题, 必须先用其他方法算出开始几个点的近似值. 一般可用同阶的单步法. ? 从以上例子看到同阶的Adams方法, 隐式方法要比显式方法误差小, 这可以从两种方法的局部截断误差主项的系数大小来加以解释. ? 先用Euler法求得一个初步的近似值，记为 ，称之为预测值，然后用它替代梯形法右端的yn+1 再直接计算 fn+1, 得到校正值 yn+1, 这样建立的预测－校正系统称为改进的Euler法： 预测 校正 ? 改进的Euler法 ? 实践表明, 改进Euler法的精度明显高于Euler法. 利用Taylor展开可以证明, 改进Euler的局部截断误差为O ( h3 ), 故它是二阶方法. 它有下列平均化形式： 例 用改进Euler法求解 解 取步长h=0.1, 计算结果见下表 改进Euler法公式为 y (xn) 1.0954 1.1832 1.2649 1.3416 1.4142 1.4832 1.5492 1.6125 1.6733 1.7321 xn 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 yn 1.0959 1.1841 1.2662 1.3434 1.4164 1.4860 1.5525 1.6165 1.6782 1.7379 ? 同Euler法的计算结果比较, 改进Euler法明显改善了精度. ? 改进Euler法的几何意义 ( K1 + K2)/2 K2 K1 xn xn+1 考察差商 根据微分中值定理，存在点?，利用所给方程 y?=f (x, y) 得 称 K*=f (?, y(? ))为区间[xn , xn+1]上的平均斜率, 这样只要对平均斜率 K*提供一种算法，相应地便导出一种计算格式. §3 Runge-Kutta法 ? Runge-Kutta法的基本思想 梯形法： 改进Euler法： 其中 Euler法：