计算方法_07.pptVIP

第七章 常微分方程数值解法 西南交通大学软件学院 计算机基础教研室 戴克俭制作 第七章 常微分方程数值解法 常微分方程是描述运动、变化规律的重要数学方法之一，分为初值问题和边值问题两大类。在高等数学中，介绍了一些典型方程的解析解的求解方法。但是，在实际工程和科学研究中遇到的微分方程往往比较复杂，很难甚至不能求出其解的解析表达式，因而需要研究数值求解方法。本章重点讨论下面的初值问题数值求解方法： 在此基础上推出一阶微分方程组与高阶微分方程初值问题的数值解法。 第七章 常微分方程数值解法 此外，还将对下面边值问题的数值解法做简要介绍： 所谓数值解法，就是寻求微分方程初值问题或边值问题的解y(x)在一系列离散点a=x0x1…xN=b上的近似值y0、y1、…、yN的方法。{yn}称为问题的数值解，数值解所满足的离散方程称为差分格式。hi=xi-xi-1称为步长，常取hi=h为定步长。 显然，只有当微分方程初(边)值问题的解存在且唯一时，使用数值解法才有意义。 第七章 常微分方程数值解法 因此总假定所讨论的问题满足解的存在唯一性定理条件。例如： 定理1 对初值问题(7-1)，当f(x,y)在带形域R={a?x?b, -?y?}上连续，且关于y满足Lipschitz条件： 其中(x, y)与(x, )为R中任意两点；L为Lipschitz常数，则初值问题的解存在且唯一。 7.1.1 Euler方法 一、欧拉(Euler)公式 将求解区间[a,b]分成N等分，令h=(b-a)/N，xn=a+nh (n=0,1,…,N)。在xn点列出方程式(7-1)： 将计算y?(xn)的数值微分公式： 代入式(7-4)可得： 7.1.1 Euler方法 截去 可得y(xn)的近似值yn所满足的递推公式： 此式称为欧拉(Euler)公式。已知y0即可由此式依次求出初值问题(7-1)的数值解y0、y1、…、yN。 欧拉公式的几何意义： 用一条折线近似代替方程的解曲线。 因此欧拉公式(7-7)常称为欧拉折线法。 7.1.1 Euler方法 例1 用欧拉公式解初值问题 7.1.1 Euler方法 表7-1 7.1.2 后退Euler方法 二、后退欧拉公式 在xn+1点处方程(7-1)为： 用向后差商逼近导数 则可得 截去 7.1 .3 局部截断误差 可得另一个数值求解递推公式： 称为后退欧拉公式。 三、局部截断误差与方法的阶 式(7-6)中的T1与(7-9)中的T2分别称为欧拉公式与后退欧拉公式的局部截断误差。 定义 对于一种数值计算公式，记方程的准确解y(x)在xn+1的值为y(xn+1)，用y(xn)代入公式计算出的近似解为yn+1，则称下式为该数值方法的局部截断误差。 7.1.4 改进的Euler方法 若局部截断误差T=O(hp+1)，则称此数值方法具有p阶精度。显然，p越大精度越高。 从本质上讲，局部截断误差是微分方程转化为差分格式时的离散化误差。对显式公式而言，也可以说是由准确解y(xn)出发，用数值方法前进一步时产生的误差。但这不是yn+1的实际误差。 四、改进的欧拉公式 为提高方法的精度，将公式(7-7)、(7-10)两端进行算术平均，得到一个新的公式： 7.1 .4 改进的Euler方法 下面我们分析公式(7-12)的局部截断误差： 将上式右端在xn点泰勒展开 得 7.1 .4 改进的Euler方法 因此公式(7-12)具有二阶精度。该公式称为梯形公式。 与后退欧拉公式一样，梯形公式是隐式公式，一般情况下采用迭代法求解，迭代初值可由欧拉公式提供，即按下式进行计算： 可以证明，当步长h满足h2/L时，上述迭代过程收敛于yn+1，其中L为式(7-3)中的李普希兹常数。 实际计算中通常只迭代一次，从而生成如下的预测-校正型公式： 7.1 .4 改进的Euler方法 该公式称为改进的欧拉公式。局部截断误差为O(h3)，所以它具有二阶精度。 为了编程方便，常将式(7-13)改写为： 7.1 .4 改进的Euler方法 例2 用改进的欧拉公式解例1中的初值问题。 解：取步长n = 0.1，改进欧拉法的具体形式为： 7.1 .4 改进的Euler方法 表7-2 7.2 龙格—库塔法 一、龙格—库塔法的基本思想 设y(x)为初值问题(7-1)的准确解，根据微分中值定理有： f(?n , y(?n))为解曲线y(x)在区间[xn, xn+1]上的平均斜率。可见只要对平均斜率提供一种近似计算方法，就能得到一种对应的差分格式。