Ch常微分方程的数值方法剖析.ppt

* 第5章 常微分方程的数值方法 近似解析法、数值方法 初值问题 本章主要研究求解一阶常微分方程初值问题的几个常用的数值方法。 且f(x,y)满足李普希茨（Lipschitz)条件，即存在常数L，使 则由常微分的理论知道，初值问题（5.1）在区间[a,b]上存在唯一解。 1.求解区间[a,b]的离散化 5.1 建立常微分数值方法的基本思想与途径 微分方程数值解法的基本思想： 把求解区间和方程离散化，求出方程的解y(x)在一系列离散点上的近似值。 因此，不同的离散方式就产生不同的数值解法。 2.将微分方程离散化 （1）差商逼近法：用差商代替导数。 （2）数值积分法： （3）Taylor展开法： 用 处的向前、向后差商分别代替（5.1）左边的微商，实现微分算子离散化。即 5.2 欧拉（Euler)方法及其截断误差和阶 5.2.1 Euler公式 Euler方法是一种最 简单的显式单步法。 5.2.1 Euler公式 显式Euler公式 隐式Euler公式 梯形公式 （隐式公式） 差分公式 单步法 5.2.2 梯形公式的计算 下面以 梯形公式为例，介绍隐式公式的迭代算法。 对于隐式方法，如果f(x,y)是y的线性函数，则隐式公式可显式计算。 但当f(x,y)是y的非线性函数时，如 5.2.2 梯形公式的计算