常微分方程数值解法

1、i.常微分方程初值问题数值解法本章讨论常微分方程初值问题数值解法，主要是差分法。解微分方程的所谓差分法的要点如下：首先是区域的离散，即将连续的求解区域离散化成有限个网格点。其次是方程的离散，例如用差商代替微商，或者对微分方程积分使之变成积分方程，然后数值积分，或者。最后得到网格点上的近似解所满足的一个差分方程，解之即得差分解。i.1 常微分方程差分法考虑常微分方程初值问题：求函数满足 （i.1a） (i.1b)其中是定义在区域: , 上的函数，和是给定的常数。我们假设对满足Lipschitz条件，即存在常数使得 （i.2）这一条件保证了（i.1）的解是适定的，即存在，唯一，而且连续依赖于初值。

2、通常情况下，（i.1）的精确解不可能用简单的解析表达式给出，只能求近似解。本章讨论常微分方程最常用的近似数值解法-差分方法。先来讨论最简单的Euler法。为此，首先将求解区域离散化为若干个离散点： （i.3）其中，称为步长。在微积分课程中我们熟知，微商（即导数）是差商的极限。反过来，差商就是微商的近似。在处，在（i.1a）中用向前差商代替微商，便得 如果忽略误差项，再换个记号，用代替便得到 一般地，我们有 （i.4）从(i.1b) 给出的初始值出发，由上式可以依次算出上的差分解。下面我们用数值积分法重新导出 Euler法以及其它几种方法。为此，在区间上积分常微分方程（i.1a），得 （i.5）

3、用各种数值积分公式计算（i.5）中的积分，便导致各种不同的差分法。例如，若用左矩形公式就得到 Euler 法（i.4）。如果用右矩形公式，便得到下面的： （i.6）类似地，如果用梯形公式，就得到 (i.7)当关于是非线性函数的时候，不能由（i.6）和(i.7)从直接算出，称这一类方法为隐式，通常采用某种迭代法求解。例如，将一般的隐式方法写成 (i.8)则可以利用如下的迭代法由算出： (i.9)关于的迭代通常只需进行很少几步就可以满足精度要求了。为了避免对隐式方法进行迭代的麻烦，比如说对于改进的Euler方法（i .7），可以采用某种预估法近似算出，然后再用（i .7）作校正，这就导致所谓预估校

4、正法。下面给出一个例子： （i.10）这是一个多步法，即计算节点上的近似值时，除了用到前一点的近似值之外，还要用到，甚至可能用到。而用前面的各种Euler法计算节点上的近似值时，只用到，因此称之为单步法。下面给出另一个多步法的例子。在区间上积分（i.1a），得 用Simpson公式(即把被积函数看作二次函数)近似计算积分，便得到 （i.11）用多步法（i.10）或（i.11）计算时，必须先用某种单步法由计算出，称为造表头。然后再逐次算出。一般说来，多步法比Euler法等简单的单步法精度要高一些。下面我们讨论一类所谓Ronge-Kutta法。他们是单步法，但是其精度可以与多步法比美。最常用的是下

5、面的标准Ronge-Kutta法和Gear法： （i.12） (i.13)从几何上， Ronge-Kutta法可以粗略地解释为：在区间中选取若干个点(可以重复)，仅仅利用在区间内可以得到的所有信息，依次给出函数在这些点上尽可能精确的的近似值，然后把它们组合起来，尽可能精确地近似计算（i.5）中的积分。Ronge-Kutta法一般的构造方式如下。选定常数，令 （i.14）选取这些待定常数的原则是：将（i.14）在作Taylor展开，然后按照的幂重新整理，使得 （i.15）与微分方程（i.1a）的解在处的Taylor展式 （i.16）有尽可能多的项重合，即要求 这里等等。按照（i.14）构造出的都

6、是显式Runge-Kutta方法, 在每一个的表达式中只出现。如果允许在某一个的表达式中出现，则可以导出隐式Runge-Kutta方法。 i.2 常微分方程组与高阶常微分方程 先来考虑下面的常微分方程组初值问题 (i.17)利用向量记号,上式可以改写为 (i.18)上节中各方法都可以直接应用到常微分方程组(i.18)。例如，Euler方法成为 再来考虑高阶常微分方程 (i.19)这时,可以令 (i.20) (i.21) (i.22)于是可以把高阶常微分方程(i.19)化成一阶常微分方程组(i.18)。i.3 收敛性与稳定性截断误差 粗略地说，截断误差可以定义为将微分方程解带入到差分方程后得到的

7、误差，代表了微分方程与差分方程之间的误差。例如，由Taylor展式和微分方程(i.1a)得到 其中是区间上某个常数。与Euler法相比较， 定义余项为Euler法的截断误差，它关于是2阶的，记为。将上节中讨论的单步法写成一般形式，则可以定义截断误差为。对于每一个差分法在适当的点（例如，或者）作Taylor展开，就可以得到截断误差的阶。对多步法可以类似处理。上节中各方法的截断误差阶分别为： 表i.1 常微分方程差分法的截断误差阶 差分法 Euler隐式Euler改进Euler预估校正(i.10)Milne标准Ronge-Kutta Gear 截断误差 相容性 一个差分方法称为相容的，如果其截断误

8、差至少是一阶的。收敛性 实用中我们更关心的是近似解的收敛性，即时，是否有。在适当的条件下，例如步长足够小，右端函数和解足够光滑等等，可以证明以上讨论的各方法都是收敛的，并且有估计式 （i.23）其中常数与和有关，与和无关；阶数等于截断误差的阶数减去1。稳定性 收敛性考虑的是差分方程的精确解与微分方程的精确解之间的误差。但是在计算机上求差分解时，由于计算机字长的限制，不可避免地产生舍入误差。另外，方程中的系数和初值等等由于测量条件等限制，也会产生数据误差。而稳定性研究的就是舍入误差和数据误差对差分方程计算结果的影响，即计算过程中某一步的“差之毫厘”会不会导致后面结果的“失之千里”。根据误差来源和

9、误差衡量标准的不同，可以定义形形色色的稳定性。一般说来，隐式方法的稳定性好于同类型的显式方法（例如Euler方法与隐式Euler方法）。另外，为了保证稳定性，常常要求步长足够小。对于常微分方程数值解来说，最简单常用的是关于初值的稳定性。以单步法为例。称某差分法关于初值稳定，如果对同一个微分方程和所有足够小的，存在常数，使得从不同初值和出发的两个差分解和之间的误差满足 （i.24）其中。由于任意一对，都可以看作是初值，关于初值的稳定性其实是考察如下问题：假设某一步计算有误差，而其后的计算不再有误差，那么这一步的误差对以后结果的影响如何。本章所讨论的所有差分方法都是关于初值稳定的。不加区别地说差分

10、方法稳定时，通常指的是关于初值稳定。不难证明，如果一个差分方法是稳定且相容的，则一定收敛，即 当 （i.25）在差分解的实际计算中，每步都难以避免地舍入误差。因而，我们进一步考察差分方法的所谓绝对稳定性，即每一步的舍入误差积累起来，是否会对后面的运算结果产生太坏的影响。考察某一个差分法对非线性常微分方程的绝对稳定性是十分困难的事情，无法得到象对初值稳定性那样的漂亮结论。通常的做法是考虑模型方程 （i.26）可以认为这个模型方程的真解和差分解的性质代表了非线性常微分方程在某一个局部的真解和差分解的性质。（例如在处将线性化，将看作，忽略其他项。）记。关于某一个差分法的绝对稳定性的典型结果是：给出一

11、个区间，使得对于所有的，求解（i.25）的这个差分法都是绝对稳定的。这个区间称之为这个差分法的绝对稳定域。绝对稳定域越大，的取值范围就越大，的允许区值就越多，所代表的非线性常微分方程也就越多。 表i.2 常微分方程差分法的绝对稳定域 差分法 Euler隐式Euler改进EulerMilne标准R-K Gear 绝对稳定域（-2，0） (0,0) (-2.78,0)(-2.78,0)关于绝对稳定性的一个必要条件是：若一个多步法是相容并且稳定的，则它绝对稳定（即它的绝对稳定域非空）的一个必要条件是。这相当于仅当时，多步法才可能绝对稳定。而当时，在差分计算中误差增长很快。那么这时候是否就无法用差分法近似计算微分方程了呢？为此，注意到模型方程的解是。从而微分方程的解当（相应于）时是按指数增长的。因此，如果这时