常微分方程的数值解法-Read.DOCVIP

常微分方程的数值解法 在自然科学的许多领域中，都会遇到常微分方程的求解问题。然而，我们知道，只有少数十分简单的微分方程能够用初等方法求得它们的解，多数情形只能利用近似方法求解。在常微分方程课中已经讲过的级数解法，逐步逼近法等就是近似解法。这些方法可以给出解的近似表达式，通常称为近似解析方法。还有一类近似方法称为数值方法，它可以给出解在一些离散点上的近似值。利用计算机解微分方程主要使用数值方法。 我们考虑一阶常微分方程初值问题 在区间[a, b]上的解，其中f (x, y)为x, y的已知函数，y0为给定的初始值，将上述问题的精确解记为y(x)。数值方法的基本思想是：在解的存在区间上取n + 1个节点 这里差，i = 0,1, …, n称为由xi到xi+1的步长。这些hi可以不相等，但一般取成相等的，这时。在这些节点上采用离散化方法，（通常用数值积分、微分。泰勒展开等）将上述初值问题化成关于离散变量的相应问题。把这个相应问题的解yn作为y(xn)的近似值。这样求得的yn就是上述初值问题在节点xn上的数值解。一般说来，不同的离散化导致不同的方法。§1 欧拉法与改进欧拉法1.欧拉法 1．对常微分方程初始问题 用数值方法求解时，我们总是认为(9.1)、(9.2)的解存在且唯一。 欧拉法是解初值问题的最简单的数值方法。从（9.2）式由于y (x0) = y0已给定，因而可以算出