【推荐】数值计算方法：常微分方程初值问题的数值解法.pptVIP

第六章 常微分方程 初值问题的数值解法 6.1 欧拉方法 6.2 龙格－库塔方法 1 问题的提出 数值求解方法 6.1 欧拉方法 6.1.1 引言 2 6.1.1 欧拉公式与后退欧拉公式与梯形公式 算法： 3 选择不同的数值积分公式来求 近似值就得到初值问题的各种数值解法 1 . 欧拉公式 这称为 欧拉公式 4 2. 后退欧拉公式 这称为 后退欧拉公式 后退欧拉公式是一个隐式公式，通常采用迭代法求解。 5 3. 梯形公式 --- 梯形公式 也是隐式单步法公式 梯形公式 6 用梯形公式计算时，通常取欧拉公式的解作为迭代初值进行迭 代计算，即采用下式 (1) (2) 7 例 6.1 以 h= 0.2 为步长，用欧拉公式、后退欧拉公式与梯形公式 求常微分方程初值问题 8 计算结果见表 6.1 9 6.1.3 改进欧拉公式 这称为 改进欧拉公式 10 例 6.2 仍取步长 h = 0.2 ，采用改进欧拉法重新计算例 6.1 的 常微分方程初值问题。 ( 计算结果见表 6-2) 这时改进欧拉公式为 解 11 6.1.4 计算公式的误差分析 定义 6.1 为该方法的 整体截断误差 . 如果 则称该方法是收敛的 . 12 设单步显式公式的一般形式为 一般地 , 微分方程初值问题精确解不满足 (6.1.15), 即 定义 6.2 称 13 截断误差的估计 设 y ( x ) ? C 3 [ x 0 , b ] , 则 （ 1 ） 对欧拉公式 ，有局部截断误差 因此，欧拉公式的局部截断误差为 O ( h 2 ) 14 对 (6.1.16) 中的积分 , 用习题五第 2 题中的右矩形公式 , 得 其中 定义 6.3 称 （ 2 ）对后退欧拉公式，有如下分析： 15 因此， 后退欧拉公式的局部截断误差为 O ( h 2 ) 其中 （ 3 ） 对梯形公式 ，对 (6.1.16) 中的积分使用数值积分得 16 其局部截断误差为 因此， 梯形公式的局部截断误差为 O ( h 3 ) （ 4 ）改进欧拉公式可以改写为 其局部截断误差为 17 因为 因此， 改进欧拉公式的局部截断误差为 O ( h 3 ) 而 从而改进欧拉公式的局部截断误差为 18