常微分方程数值解法本章主要内容欧拉法改进欧拉法2

1、第8章常微分方程数值解法本章主要内容：1 .欧拉法、改进欧拉法.2龙格-库塔法。3.单步法的收敛性与稳定性。重点、难点一、微分方程的数值解法在工程技术或自然科学中，我们会遇到的许多微分方程的问题，而我们只能对其中具有较简单形式的微分方程才能够求出它们的精确解。对于大量的微分方程问题我们需要考虑求它们的满足一定精度要求的近似解的方法，称为微分方程的数值解法。本章我们主要乎二 f (x, y)讨论常微分方程初值问题的数值解法。dxIy(x) = yo数值解法的基本思想是：在常微分方程初值问题解的存在区间a,b内，取n+1个节点a=xo X! VV XN=b (其中差hn= Xn -Cn-1称为步长

2、，一般取 h为常数，即等步长)，在这 些节点上把常微分方程的初值问题离散化为差分方程的相应问题，再求出这些点的上的差 分方程值作为相应的微分方程的近似值(满足精度要求)。、欧拉法与改进欧拉法欧拉法与改进欧拉法是用数值积分方法对微分方程进行离散化的一种方法。将常微分方程 y =f(x, y)变为 y(XnJ = y(Xn)Xn1 f(t,y(t)dtxn1.欧拉法(欧拉折线法)欧拉法是求解常微分方程初值问题的一种最简单的数值解法。欧拉法的基本思想：用左矩阵公式计算(*)式右端积分，则得欧拉法的计算公式为：b a yn 1 =ynhf(Xn,yn) (n =0,1，N-1) h =N欧拉法局部截断

3、误差h2Rn 1 y ( n 1) Xn咗n 1咗X. i或简记为0()。我们在计算时应注意欧拉法是一阶方法，计算误差较大。欧拉法的几何意义：过点 Ao (Xo, yo), Al (Xi, yi),，A n (x n, y n )，斜率分 别为f (x0, y0), f (x1, y1),，f (x n, y n)所连接的一条折线，所以欧拉法亦称为 欧拉折线法。例1用欧拉法解初值问题dy2xydx(0) =1(0 _ x _1)在x = 0 (0.2) 1处的近似解。(计算过程保留4位小数)。【思路】用欧拉法求解常微分方程的初值问题时，首先熟练掌握欧拉公式的一般形式, 根据具体题目写出找出欧拉

4、公式的迭代式，并根据初始条件和所给步长进行迭代求解。解/ f(x, y) = - 2xy , h= 0.2,欧拉公式为：yn 1 二 yn hf (Xn,n )二n 0.2( Xn) = (1 - 0.4Xn)yn(n =0,1,2,3,4,5)列表计算如下:nXnyny(xn)y(xn)-yn0011010.210.9608-0.039220.40.920.8521-0.067930.60.77280.6977-0.075140.80.58730.5273-0.06510.39940.3679-0.03152.改进欧拉法改进欧拉法比欧拉法的计算准确，是对欧拉法的改进。改进欧拉法的基本思想：用

5、梯形公式计算(*)式右端积分，则得改进欧拉法的计算公式为：h |1b - ayn ”詁 2 f(Xnn) 5 ln 1)5 =0,1，,N -1) h 禾U用改进欧拉法计算常微分方程初值问题时，我们应注意此公式为隐式表达式，需要对它进行迭代求解。计算时可以采用一次迭代和多次迭代，因此，就有改进欧拉法预估-校正法公式和反复迭代的改进欧拉法预估-校正法公式。改进欧拉法预估-校正法公式：=yn +hf(Xn, yn)*hyn=yn 蔦 f (Xn,yn) +n =0,1，N 1)反复迭代的改进欧拉法预估-校正法公式：=yn +hf(Xn, yn)*yff)=ynf 区皿)+ f (x1, yn；】(

6、n= 0,n -1, m = 0,1,.,)h?改进欧拉法的局部截断误差Rn 1二一 y ( nd)冷乞；1乞召1或简记为0( h)。12从局部截断误差的形式看，改进欧拉法是二阶方法，因此，它比欧拉法更精确。例2用预估-校正法求初值问题y = y xy2(0 乞 x)y(0)二 1在 x=0 (0.2) 1 的解。【思路】掌握预估-校正法的计算公式，根据已知条件迭代求解。解步长h=0.2，将f (x, y)二-y -xy2代入预估-校正公式，整理得yn0)1 =0.8yn -0.2Xny；% 1 =0.9yn 0.1Xny； -0.1心1 *.1$0)1)2列表计算如下:nX ny ny(X

7、n)001110.20.807200.8046320.40.636900.6314530.60.490480.4891840.80.377800.37720510.291030.29100例3用改进欧拉法求解例 1的初值问题，要求ynm) - ynm4 10。【思路】掌握改进欧拉法的计算公式，根据已知条件迭代求解，并检验迭代解是否满 足精度要求，若满足则确定此解为常微分方程在某点的近似解。解 将f(x, y) = -2xy代入改进欧拉法的计算公式得：y = yn +hf (Xn, yn) = y. + 0.2( -2x. y.) = (1 - Og) y.yJV yn +f (Xn, yn)

8、+ f (Xn*, yfV=y. - 0.2 X.n 十 乂卄丫阳 b列表计算如下:nXnyPyn2)yf)= yny(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.0051三、龙格-库塔法1 龙格-库塔法龙格-库塔法具有精

9、度高、收敛、稳定，不需要计算高阶导数等优点，是求解微分方程初值问题的一组著名的显示单步方法，广泛应用于求解常微分方程的初值问题。本章我们介绍了二、三、四阶龙格 -库塔法。龙格-库塔法的基本思想：在计算初值问题二f(x,y)的数值解时，考虑均差yXG-yX)，则由微分中值y(x0) = yoh定理可得y(xn .J - y(xn)二hy(xn v h) (0 ： : 1)，由初值问题可得公式为：y(xk 勺)=y(xk) hf (xk)h, y(xk h)上式中 f (xk r h, 丫(兀r h)称为区间上的平均斜率。如果给平均斜率一种计算方法，就可得到计算y (x n+1 )的近似值yn+1

10、的公式。如果仅取Xn处的斜率值f(Xn,yn)作为平均斜率的近似值，则得到的n 1的公式为欧拉公式；如果取Xn ,Xn 1处的斜率值f(xn,yn)，f(X. 1,n 1)的平均值作为平均斜率的近似值，则得到的yn 1的公式改进欧拉公式。二阶龙格-库塔法的公式和局部截断误差：yn+ = yn +h(入水1 十人2k2)k1 二 f (Xn, yn)k f (Xn 1, yn hk1)在上式中选择不同的参数，会得到不同的二阶龙格-库塔法公式，所以二阶龙格-库塔法公式不唯一。二阶龙格-库塔法公式的局部截断误差为O ( h3)。常见的二阶龙格-库塔法公式有以下两种改进欧拉法迭代公式yn = yn +

11、 (ki * k2)2也=f (Xn, yn)k2 = f (Xn41,yn +hkj= yn +hk2仏=f (xn, yn)hhk2 = f (Xn 十二，Yn + 二 kJ 2 2三阶龙格-库塔法的公式和局部截断误差：常见的三阶龙格-库塔法公式为khyn = yn + (ki +4k2 +k3)6ki = f (Xn, yn)*k2 = f (Xn 十一，yn + ki)2 2k f (Xn h,yn -hki 2hk2)三阶龙格-库塔法公式的局部截断误差为o( h4)。四阶龙格-库塔法公式yn+ =yhk2* = f(Xn, yn)k2=f(XnJ，ynJkl)2 2通常所说的龙格-库

12、塔法是指四阶龙格-库塔法，也称为标准龙格-库塔法。由于它是一步 法，(即已知yn ,就可以求出yn 1，无需知道ynd, yn/ ,的值)且它的计算精度 高，所以应用较多，但在计算时，因为每一步都需要计算四次f (X, y)的值，计算量较大，所以，一般用来计算前几项的近似值，即“表头”。四阶龙格-库塔法公式为的公式和局部截断误差：yn =yn +( +2k2 +2k3 +k4)6ki = f (Xn, yn)hh也=f (Xn +：,yn + 二 &)22hhk3 = f (Xn+ 二，Yn+ 二k2)22k f (Xn h, yn hk3)四阶龙格-库塔法的局部截断误差为 O ( h5)。四

13、、单步法的收敛性和稳定性1.收敛性如果在无舍入误差且步长 h充分小的情况下，求得的近似值 yn足够精确地逼近真解y(Xn)，即：当h 0时，一致地有头；y(Xn) 欧拉法整体截断误差：5二y(xj - yn其中y(Xn)为真解，yn为在无舍入误差情况下，从yo用欧拉法计算公式求得的近似解。 欧拉法的收敛条件：如果f(x,y)关于y满足Lipschitz条件，且局部截断误差Rn有界，即h 2Rn兰一M 2 ( n = 1,2，N ),则欧拉法收敛。且欧拉法的整体截断误差估计式为：+譽2严)-1)其中L为Lipschitz常数，b-a为求解区间的长度， M 2 = max y(x)。aX亜3稳定性

14、和绝对稳定性稳定性：指初始(或某步)产生的误差在后面的迭代计算中不会再扩大。即存在常数C及ho,OV hW ho时，对任意两个初始值 yo,o满足不等式 n - yn兰Co - yo。 欧拉法稳定性的条件：如果f(x,y)关于y满足Lipschitz条件，则欧拉法稳定。 绝对稳定性：若对固定步长ho及任意两个初始值Yo , o满足不等式n - yn 兰 O - yo。 我们在讨论稳定性时应注意，一般在实际计算中只能取固定步长，它不可能任意缩小。所以绝对稳定性则表示的是对固定步长h。，在初始(或某步)所产生的误差，在以后计算中不会逐步增长。由于绝对稳定性的成立和f(x,y)有关，讨论较为复杂。所以一般fy y地，我们对简单的微分方程yy的绝对稳定