第五章 常微分方程的数值解法

1、第五章 常微分方程的数值解法,许多工程实际问题的数学模型可以用常微分方程来描述，但是，除了常系数线性微分方程和一些特殊的微分方程可以用解析方法求解以外，在生产实践和科学研究中遇到的绝大多数微分方程往往比较复杂，不能给出解的解析表达式；有时候即使能用解析表达式来表示，又因为计算量太大而不实用，有些是已经有了求解的基本方法的典型方程，但实际使用时也是有困难的。,很多实际问题只需要求某些点上的函数值，不必求出函数表达式。,因此，我们有必要研究常微分方程的数值近似解法。本章主要研究一阶常微分方程初值问题的几个常用解法。,例：某伞兵与降落伞质量为90kg，当伞张开时以12m/s速度垂直下降。假设空气阻力

2、与伞的下降速度成正比，在速度为6m/s时，测得的空气阻力为353N。试求伞兵开伞后，第1秒末、第2秒末、直到第5秒末时各时刻的速度。,解：,设阻力为 ，则,由力学知识可得，,本章主要内容,一阶常微分方程初值问题,欧拉方法、梯形格式、改进的欧拉方法,龙格库塔（Runge-Kutta）方法,第一节 引言,一阶方程的初值问题：,由常微分方程理论，在 满足一定条件时存在唯一解函数 。,要计算出解函数 y(x) 在一系列节点 a = x0 x1 xn= b 处的近似值 。,常微分方程的数值解,在这些节点上采用离散化方法，（通常用数值积分、微分、 泰勒展开等）将上述初值问题化成关于离散变量的相应问题。 把

3、这个相应问题的解 作为 的近似值。这样求得的 就是上述初值问题在节点 上的数值解。一般说来，不同的 离散化导致不同的方法。,第二节 欧拉方法,1、欧拉方法的导出,对于上面的初值问题，若将函数 在 处的导数用差商来表示，则有，,解：根据欧拉公式可以得到：,此外，方程的真解为：,2、举例,例1 用欧拉法求初值问题 取步长h = 0.2。,求解结果如下：,根据已知条件：曲线y(x)上的点（x0,y0）及该点处曲线的导数f(x0,y0)，则可以得到过该点的直线：,该直线与xx1的交点P1，则P1的纵坐标y1为：,就用y1作为y(x1)的近似值,逐次进行后可以得到一条折线P0 P1Pn，该折线看作是初值

4、问题的积分曲线的近似，因此欧拉方法也称为欧拉折线法。,3、欧拉方法的几何意义,欧拉法虽然形式简单，计算方便。 但从上述几何意义上得知，由Euler法所得的折线明 显偏离了积分曲线，可见此方法非常粗糙,精度也低， 即误差太大。特别当的曲线曲率较大时，欧拉法的效果 更差。,二、梯形格式,对微分方程两边积分得，,欧拉方法,左矩形公式,后退的欧拉方法,右矩形公式,梯形方法,梯形公式,和欧拉公式相比较，梯形公式精度有所提高；但梯形公式是隐式的，在计算yi+1时候也只用到前一步的值yi，若yi已知，将yi代入公式求解时，一般不能直接得到yi+1，而需要通过其他方法（比如迭代法）求解，计算量较大。,实际中，

5、将欧拉公式和梯形公式联合使用，从而得到改进的欧拉方法。,三、改进的欧拉方法,将欧拉公式和梯形公式联合使用，先用欧拉公式得出一个y(xi+1)的粗糙近似值，称为预估值，然后对预估值使用梯形公式对它进行精确化，得到较为精确的近似值yi+1，称之为校正值，计算公式为：,这样的预估校正系统称为改进的欧拉方法,为了便于编写程序，常将上面的公式改写为如下式：,改进的欧拉方法与梯形方法具有同样的精度级别，但改进的欧拉方法为显式格式，使用更方便。,解：根据改进的欧拉公式得：,此外，方程的真解为：,举例,例1 用改进的欧拉法求初值问题 取步长h = 0.2。,求解结果如下：,例：某伞兵与降落伞质量为90kg，当

6、伞张开时以12m/s速度垂直下降。假设空气阻力与伞的下降速度成正比，在速度为6m/s时，测得的空气阻力为353N。试求伞兵开伞后，第1秒末、第2秒末、直到第5秒末时各时刻的速度。,0 12.0000 1.0000 13.9600 2.0000 14.6395 3.0000 14.8750 4.0000 14.9567 5.0000 14.9850,0 12.0000 1.0000 13.3197 2.0000 14.0589 3.0000 14.4729 4.0000 14.7048 5.0000 14.8346,欧拉法,改进的欧拉法,第三节 龙格库塔方法,一、龙格库塔法的基本思想,1、预备知

7、识,2、推导,称为 在 上的平均斜率。,问题转化为如何对 进行数值计算。,特例：,1）取一个点 上的斜率值 作为平均斜率 的近似值。,欧拉公式,精度较低,2）用两个点 和 上的斜率 和 的算术平均值作为平均斜率 的近似值。,改进的欧拉公式,精度稍高,根据前面的分析得到： 如果在区间 上多预测几个点的斜率值，然后取其加权平均作为平均斜率的近似值，有可能构造出具有更高精度的计算公式，此即为龙格库塔算法的基本思想。,二、龙格库塔公式,1、二阶龙格库塔公式,选择参数 使得算法具有2阶精度,满足该式的参数不止一组，而是一簇，所有满 足条件的公式通称为二阶龙格库塔公式。,改进的欧拉公式是二阶龙格库塔公式的其中一种。,2、四阶龙格库塔公式,二阶的龙格库塔公式是由两个点的斜率加权平均（或线性组合）得到的，同样，用四个点上的斜率线性组合就可得到四阶的龙格库塔公式。,经典的四阶的龙格库塔公式：,举例,例1 用经典龙格库塔公式求初值问题 取步长h = 0.2。,求解结果如下：,四阶龙格库塔法的精确度高于欧拉法和改进的欧拉法。,必须说明，由于四阶龙格库塔法有较高的精确度，大家自然会想到推导更高阶的龙格库塔公式。但实践证明，高于四阶的龙格库塔公式，不但计算量大，而且精