计算方法课件-第六章.pptVIP

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