第6讲 常微分方程数值解法.pptVIP

第6章 常微分方程数值解法 绪论 在工程和科学计算中，所建立的各种常微分方程的初值或边值问题，除很少几类的特殊方程能给出解析解，绝大多数的方程是很难甚至不可能给出解析解的，其主要原因在于积分工具的局限性。因此，人们转向用数值方法去解常微分方程，并获得相当大的成功，讨论和研究常微分方程的数值解法是有重要意义的。 6.1 初值问题的Euler方法 初值问题的Euler方法 初值问题的Euler方法 初值问题的Euler方法 初值问题的Euler方法 初值问题的Euler方法 初值问题的Euler方法 初值问题的Euler方法 初值问题的Euler方法 初值问题的Euler方法 6.1.2 误差概述 误差概述 误差概述 误差概述 6.1.3 数值稳定性分析 数值稳定性分析 定义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方法 四阶Runge-Kutta方法 6.2.3 R-K法的稳定性 R-K法的稳定性 R-K法的稳定性 6.2.5 隐式R-K法 隐式