常微分方程的数值解法.pptxVIP

1第九章常微分方程数值解法第一节Euler方法第三节单步法的收敛性和稳定性第二节Runge-Kutta方法上一页下一页返回

本章介绍求解微分方程数值解的基本思想和方法.含有自变量、未知函数和它的一阶导数和高阶导数的方程.常微分方程它是描述运动、变化规律的重要数学方法之一，分为两类：初值问题即给出未知函数及导数在初始点的值；边值问题即给出未知函数及（或）它的某些导数在区间两个端点的值。

?考虑一阶常微分方程的初值问题:3只要f(x,y)在[a,b]?R1上连续，且关于y满足Lipschitz条件，即存在与x,y无关的常数L使对任意定义在[a,b]上的y1(x)和y2(x)都成立，则上述问题解存在唯一解。所谓数值解法就是要计算出初值问题的解函数y(x)在一系列离散点a=x0x1…xN=b上的近似值：y0,y1,……yN.节点间距为步长，通常采用等距节点，即取hi=h(常数)。{yn}称为问题的数值解.数值解所满足的离散方程统称为差分格式.上一页下一页返回

第一节欧拉方法4上一页下一页返回一、欧拉公式令yn为y(xn)的近似值，将上式代入(*)式可得此式称为欧拉(Euler)公式.为Euler方法的局部截断误差.

例1用欧拉公式解初值问题5上一页下一页返回解：取步长h=0.1,欧拉公式的具体形式为：依次计算可得………y3，y4，y5，y6，y7，y8，y9，y10

其部分结果见下表6上一页下一页返回可见Euler方法的计算结果精度不太高。

?欧拉公式的几何意义：7上一页下一页返回x0P0x1P1x2P2xnPn几何意义：用折线近似代替方程的解曲线,因而也称Euler方法为折线法.

二、后退的欧拉公式8上一页下一页返回令yn+1为y(xn+1)的近似值，则可得也用一阶差商逼近导数已知yn时,必须通过解方程才能求出yn＋1,这样的公式称为隐式公式,而Euler公式为显式公式.Euler公式和后退Euler公式都是由yn去计算yn+1，因此，称它们为单步法。称为后退Euler公式

三、局部截断误差与方法的阶上一页下一页返回9定义在假设yi=y(xi)，即第i步计算是精确的前提下，考虑的截断误差Ti+1=y(xi+1)?yi+1称为局部截断误差。定义若某算法的局部截断误差为O(hp+1)，则称该算法有p阶精度。显然,p越大,精度越高.（将准确解代入公式的左、右两端，其左端与右端之差）Euler方法的精度其中：

1所以，Euler方法具有1阶精度。2将3在点4处一阶Taylor展开

所以，后退的Euler方法也具有1阶精度。将在点处一阶Taylor展开隐式Euler方法的精度

?欧拉公式的改进12上一页下一页返回—显、隐式两种算法的平均其局部误差为：此公式具有2阶精度.称平均公式或梯形公式梯形公式可由下迭代式计算：其中迭代初值是Euler公式提供.

四、改进的欧拉公式Step1:先用显式欧拉公式作预测，算出),(1iiiiyxfhyy+=+Step2:再将代入隐式梯形公式的右边作校正，得到1+iy)],(),([2111+++++=iiiiiiyxfyxfhyy注：此法亦称为预测-校正法。可以证明该算法具有2阶精度，同时可以看到它是个单步递推格式，比隐式公式的迭代求解过程简单。它的精度高于显式欧拉法。

为了便于编程,常将改进的欧拉公式写为：

例2用改进的欧拉法解例1中的初值问题.解：取步长h=0.1,改进欧拉法的具体形式为具体计算过程如下

16xn改进的欧拉法误差xn改进的欧拉法误差0100.61.4859560.0027160.21.1840960.0000880.81.6164760.0040240.41.3433600.0017191.01.7378690.005818依次计算可得………y3，y4，y5，y6，y7，y8，y9，y10其部分结果见下表上一页下一页返回

例3对下面的初值问题17上一页下一页返回解（1）取步长h=0.1,欧拉方法的具体公式为取步长h=0.1,改进的欧拉方法的具体公式为取步长h=0.1，分别用Euler方法、改进的Euler方法求数值解。

计算结果见下表