第九章+常微分方程初值问题数值解.docVIP

常微分方程初值问题 的数值解法 在自然科学的工程技术的许多领域中，常会遇到常微分方程初值问题，但这种问题大多数情况下不存在初等形式的解析解，只能用近似方法来求解。近似解法主要有两类：一类叫做近似解析方法，它能给出解得近似表达式，例如熟知的级数解法和逐次逼近法等；另一类近似解法称为数值解法，它可以给出解在一些离散点上的近似值。数值方法便于电子计算机求解微分方程。本章讨论初值问题常用的数值解法，并介绍有关的基本理论。 假设给定一阶常微分方程的初值问题： （9.1）（9.2） 其中f为x, y的已知函数，α为给定的初始值。 由微分方程的理论可知，如果f (x,y)在区域，内连续，且关于y满足李普希兹(Lipschitz)条件，即存在常数L，使 对所有的及任何，均成立，则初值问题（9.1），（9.2）有连续可微的解y(x)存在且唯一。所谓初值问题的数值解，则是问题的解y(x)在一系列点 处的值的近似值(n=0,1,…,N)。这里相邻两个节点之间的距离通常称为步长，通常将步长取为常数h。 §9.1 欧拉方法与改进的欧拉方法 欧拉(Euler)方法是最简单的数值方法，由于它的精确程度较差，已不常用于实际计算。但构造这个方法的基本原理，对于构造一般的数值方法具有普遍意义，因此首先对它进行讨论。 欧拉方法