常微分方程数值解法试卷.ppt

第7章 常微分方程数值解法 绪论 在工程和科学计算中，所建立的各种常微分方程的初值或边值问题，除很少几类的特殊方程能给出解析解，绝大多数的方程是很难甚至不可能给出解析解的，其主要原因在于积分工具的局限性。因此，人们转向用数值方法去解常微分方程，并获得相当大的成功，讨论和研究常微分方程的数值解法是有重要意义的。 初值问题的Euler方法 Question : 上述方法得到的解是一个什么形式？ 与高等数学中常微分方程的解的形式一样么？有什么区别 初值问题的Euler方法 初值问题的Euler方法 初值问题的Euler方法 初值问题的Euler方法 初值问题的Euler方法 初值问题的Euler方法 初值问题的Euler方法 误差 改进的Euler方法局部截断误差为O(h3) 整体截断误差 Euler方法：O(h) 改进的Euler方法：O(h^2) 误差概述 数值稳定性分析 数值稳定性分析 定义6.1.3 若某数值算法的绝对稳定性区域包含hλ平面上的左半平面Re(hλ)＜0，则称该方法是A稳定的。 隐式Euler法是A稳定的。 6.2 Runge-Kutta方法 Runge-Kutta方法 Runge-Kutta方法 Runge-Kutta方法 6.2.2 四阶Runge-Kutta方法 四阶