常微分方程的数值解法

《常微分方程的数值解法》由会员上传分享，免费在线阅读，更多相关内容在应用文档-天天文库。

1、第九章常微分方程的数值解法在自然科学的许多领域中，都会遇到常微分方程的求解问题。然而，我们知道，只有少数十分简单的微分方程能够用初等方法求得它们的解，多数情形只能利用近似方法求解。在常微分方程课中已经讲过的级数解法，逐步逼近法等就是近似解法。这些方法可以给出解的近似表达式，通常称为近似解析方法。还有一类近似方法称为数值方法，它可以给出解在一些离散点上的近似值。利用计算机解微分方程主要使用数值方法。我们考虑一阶常微分方程初值问题在区间[a,b]上的解，其中f(x,y)为x,y的已知函数，y0为给定的初始值，将上述问题的精确解记为y(x)。数值方法的基本思想是：在解的存在区间上取n+1个节点这里

2、差，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已给定，因而可以算出设x1=h充分小，则近似地有：(9.3)记从而

3、我们可以取作为y(x1)的近似值。利用y1及f(x1,y1)又可以算出y(x2)的近似值：一般地，在任意点xn+1=(n+1)h处y(x)的近似值由下式给出（9.4）这就是欧拉法的计算公式，h称为步长。不难看出，近似解的误差首先是由差商近似代替微商（见(9.3)）引起的，这种近似代替所产生的误差称为截断误差。还有一种误差称为舍入误差，这种误差是由于利用（9.4）进行计算时数值舍入引起的。例9.1用欧拉法求初值问题当h=0.02时在区间[0,0.10]上的数值解。解把代入欧拉法计算公式。就得具体计算结果如下表：nxnyny(xn)en=y(xn)-yn001.00001.0000010.020

4、.98200.98250.000520.040.96500.96600.000530.060.94890.95030.001440.080.93360.93540.001850.100.91920.9230.0021在上表中y(xn)列，乃是初值问题（9.5）、（9.6）的真解在xn上的值。为近似值yn的误差。从表中可以看出，随着n的增大，误差也在增大，所以说，欧拉法计算简便，对一些问题有较大的使用价值，但是，它的误差较大，所得的数值解精确度不高。2．改进欧拉法为了构造比较精确的数值方法，我们从另一角度重新分析一下初值问题。一般说来，一阶方程的初值问题与积分方程（9.7）是等价的，当x=x1

5、时，（9.8）要得到y(x1)的值，就必须计算出（9.8）式右端的积分。但积分式中含有未知函数，无法直接计算,只好借助于数值积分。假如用矩形法进行数值积分则因此有可见，用矩形法计算右端的积分与用欧拉法计出的结果完全相同。因此也可以说欧拉法的精度之所以很低是由于采用矩形法计算右端积分的结果。可以想象，用梯形公式计算（9.8）式右端的积分，可期望得到较高的精度。这时将这个结果代入（9.3）并将其中的y(x1)用y1近似代替，则得这里得到了一个含有y1的方程式，如果能从中解出y1，用它作为y(x1)的近似值，可以认为比用欧拉法得出的结果要好些。仿照求y1的方法，可以逐个地求出y2,y3，…。一般地

6、当求出yn以后，要求yn+1，则可归结为解方程：这个方法称为梯形法则。用梯形法则求解，需要解含有yn+1的方程式，这常常很不容易。为此，在实际计算时，可将欧拉法与梯形法则相结合，计算公式为（9.9）这就是先用欧拉法由(xn,yn)得出y(xn+1)的初始近似值，然后用（9.9）中第二式进行迭代，反复改进这个近似值，直到（e为所允许的误差）为止，并把取作y(xn+1)的近似值yn+1。这个方法就叫改进欧拉方法。显然，应用改进欧拉法，如果序列收敛，它的极限便满足方程即序列的极限可取作yn+1。可以证明，如果有界，则只要h取得适当小，上述序列必定收敛。这样当h取得充分小，就可保证上述迭代过程收敛到

7、一个解。当步长h取得适当时，欧拉方法算出的值已是较好的近似，因此改进欧拉法收敛很快，通常只需二、三次迭代即可。如果迭代很多步仍不收敛，这表明表长h选的过大，应缩小步长后再计算。通常把（9.9）叫做预报校正公式，其中第一式叫预报公式，第二式叫校正公式。这个公式还可写为(9.9)’3．公式的截断误差现在来考察两个公式的截断误差：有多大？这里假定前一步得的结果yn=y(xn)是准确的。写出y(xn+1)的泰勒展开式