第8章 常微分方程数值解法 本章主要内容： 1．欧拉法、改进欧拉法 2

《第8章 常微分方程数值解法 本章主要内容： 1．欧拉法、改进欧拉法 2》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、第8章常微分方程数值解法本章主要内容：1．欧拉法、改进欧拉法.2．龙格-库塔法。3．单步法的收敛性与稳定性。重点、难点一、微分方程的数值解法在工程技术或自然科学中，我们会遇到的许多微分方程的问题，而我们只能对其中具有较简单形式的微分方程才能够求出它们的精确解。对于大量的微分方程问题我们需要考虑求它们的满足一定精度要求的近似解的方法，称为微分方程的数值解法。本章我们主要讨论常微分方程初值问题的数值解法。数值解法的基本思想是：在常微分方程初值问题解的存在区间[a,b]内，取n+1个节点a=x0＜x1＜…＜xN=b（其中差hn=xn–xn-1称为步长，一般取h为常数，即

2、等步长），在这些节点上把常微分方程的初值问题离散化为差分方程的相应问题，再求出这些点的上的差分方程值作为相应的微分方程的近似值（满足精度要求）。二、欧拉法与改进欧拉法欧拉法与改进欧拉法是用数值积分方法对微分方程进行离散化的一种方法。将常微分方程变为　1．欧拉法（欧拉折线法）欧拉法是求解常微分方程初值问题的一种最简单的数值解法。欧拉法的基本思想：用左矩阵公式计算（＊）式右端积分，则得欧拉法的计算公式为：欧拉法局部截断误差或简记为O（h2）。8我们在计算时应注意欧拉法是一阶方法，计算误差较大。欧拉法的几何意义：过点A0（x0，y0），A1（x1，y1），…，An（xn

3、，yn），斜率分别为f（x0，y0），f（x1，y1），…，f（xn，yn）所连接的一条折线，所以欧拉法亦称为欧拉折线法。例1用欧拉法解初值问题 在x＝0(0.2)1处的近似解。（计算过程保留4位小数）。【思路】用欧拉法求解常微分方程的初值问题时，首先熟练掌握欧拉公式的一般形式，根据具体题目写出找出欧拉公式的迭代式，并根据初始条件和所给步长进行迭代求解。解∵f(x，y)＝－2xy，h＝0.2，欧拉公式为：列表计算如下：nxnyny(xn)y(xn)-yn0011010.210.9608-0.039220.40.920.8521-0.067930.60.77280.

4、6977-0.075140.80.58730.5273-0.06510.39940.3679-0.031582．改进欧拉法改进欧拉法比欧拉法的计算准确，是对欧拉法的改进。改进欧拉法的基本思想：用梯形公式计算（＊）式右端积分，则得改进欧拉法的计算公式为：利用改进欧拉法计算常微分方程初值问题时，我们应注意此公式为隐式表达式，需要对它进行迭代求解。计算时可以采用一次迭代和多次迭代，因此，就有改进欧拉法预估-校正法公式和反复迭代的改进欧拉法预估-校正法公式。改进欧拉法预估-校正法公式：反复迭代的改进欧拉法预估-校正法公式：改进欧拉法的局部截断误差或简记为O（h3）。从局部

5、截断误差的形式看，改进欧拉法是二阶方法，因此，它比欧拉法更精确。例2用预估-校正法求初值问题在x=0（0.2）1的解。【思路】掌握预估-校正法的计算公式，根据已知条件迭代求解。解步长h=0.2，将代入预估-校正公式，整理得列表计算如下：8nxnyny(xn)001110.20.807200.8046320.40.636900.6314530.60.490480.4891840.80.377800.37720510.291030.29100例3用改进欧拉法求解例1的初值问题,要求。【思路】掌握改进欧拉法的计算公式，根据已知条件迭代求解，并检验迭代解是否满足精度要求，

6、若满足则确定此解为常微分方程在某点的近似解。解将代入改进欧拉法的计算公式得：列表计算如下：ny(xn)y(xn)-yn0011111010.210.960.96160.96150.9608-0.000720.40.88460.85230.85490.85490.8521-0.002830.60.71810.70030.70250.70220.6977-0.004540.80.53370.53250.53270.53270.5273-0.0054510.36220.37500.37250.37300.3679-0.00518三、龙格-库塔法1．龙格-库塔法龙格-库塔

7、法具有精度高、收敛、稳定，不需要计算高阶导数等优点，是求解微分方程初值问题的一组著名的显示单步方法，广泛应用于求解常微分方程的初值问题。本章我们介绍了二、三、四阶龙格-库塔法。龙格-库塔法的基本思想：在计算初值问题的数值解时，考虑均差，则由微分中值定理可得，由初值问题可得公式为：上式中称为区间上的平均斜率。如果给平均斜率一种计算方法，就可得到计算y（xn+1）的近似值yn+1的公式。如果仅取处的斜率值作为平均斜率的近似值，则得到的的公式为欧拉公式；如果取处的斜率值，的平均值作为平均斜率的近似值，则得到的的公式改进欧拉公式。㈠二阶龙格-库塔法的公式和局部截断误差：在

8、上式中选择