数电-第六章 常微分方程初值问题的数值解法.pptVIP

* 第六章 常微分方程初值问题的数值解法 6.1 欧拉方法　 　 　　　 　　 6.2 计算公式的误差分析 　　 6.3 龙格－库塔方法 　　　　 6.4 向一阶方程组与高阶方程的推广 问题的提出 数值求解方法 6.1 欧拉方法 6.1.1 欧拉公式与改进欧拉公式 算法： 选择不同的数值积分公式来求 近似值就得到初值问题的各种数值解法 1.欧拉公式 这称为欧拉公式 2.后退欧拉公式 这称为后退欧拉公式 后退欧拉公式是一个隐式公式，通常采用迭代法求解。 例6.1 以 h=0.1为步长，用欧拉法求常微分方程初值问题 …. …. ….. …… 9.362x10-3 0.909362 0.900000 0.1 0 1 1 0 y(xi)-yi y(xi) yi xi 立表 (见表6-1（书121页）) 6.1.2 梯形公式与改进欧拉公式 3.梯形公式 ---梯形公式也是隐式单步法公式 图1 梯形公式 用梯形公式计算时，通常取欧拉公式的解作为迭代初值进行迭代计算，即采用下式 (1) (2) 4.改进欧拉公式 这称为改进欧拉公式 例6.2 仍取步长h = 0.1，采用改进欧拉法重新计算例 6.1 的 常微分方程初值问题。 (见表6-2（书125页）) 这时改进欧拉公式为 解 …. …. ….. …… 1.363x10-4