《常微分方程课程计》指导书3.docVIP

第3章 数值算法之一：单步法 接下来的章节，我们将详细讲解常微分方程的数值求解方法，特别是编程实现方程的初值问题，这也将成为解决实际问题的必要基础和课程设计的主要内容。 本章重点讲解一阶常微分方程的数值算法中最简单的一类方法：单步法。首先介绍Euler法、后退Euler法与梯形法，并分析单步法的局部截断误差，然后给出了改进的Euler法。 3.1　简单的单步法及基本概念 3. 求初值问题(1.2.1)的一种最简单方法是将节点的导数 代替，于是(1.2 　　　　 　　　(3.1 从出发，由(3.1.1)求得再将代入(3.1.1)右端，得到的近似，一般写成 　　　　　　　　 　　(3.1 该数值解法称为解微分方程初值问题的Euler法. Euler法的几何意义如图3-1所示。初值问题(1.2.1)的解曲线y=y(x)过点，从出发，以为斜率作一段直线，与直线交点于，显然有，再从出发，以为斜率作直线推进到上一点，其余类推，这样得到解曲线的一条近似曲线，它就是折线. 图3-1 Euler法的几何意义显示 Euler法也可利用的Taylor展开式得到，由 　　　　　　(3.1 略去余项，以，就得到近似计算公式(3.1.2).另外，还可对(1.2.1 　　　　　　　　　 　　(3.1.4) 若右端积分用左矩形公式，用，，则得(3.1.2).如果在(3.1.4)的积分中用右矩形公式，则得