81常微分方程定解问题数值解得概念82初值问题的Euler方法局部截断误差.docxVIP

第八章常微分方程数值解 8.1引言（基本求解公式） 在工程和科学技术的实际问题中，常需要求解微分方程 只有简单的和典型的微分方程可以求出解析解 而在实际问题中的微分方程往往无法求出解析解 在高等数学中我们见过以下常微分方程： = axb 侦。）=% ⑴ (2)卩 = f（x 双?/） axb (2) \y（a） = yQ/y\a） = a = 以?/) a xb [y伝)=%，对)=义 (1),(2)式称为初值问题,(3)式称为边值问题 另外，在实际应用中还经常需要求解常微分方程组: (?Vi =/[(人/力,光)1(玉))=力0 (4) 1 义=认X, ylfy2) J2(XO) = y20 本课程主要研究问题⑴的数值解法，对(2)?⑷只作简单介绍 我们首先介绍初值问题(1)的解存在的条件 定理1?如果连续函麴Xx,’)满足￡如祯该条件，即 HE数乙，使得Vrc 口0],均有 I / (x, y 1) -/(X, y2)\L\yi-y2\ 则初值问题(1)的解存在且唯一. 对于问题(1),要求它的数值解 就是求未知函麴(X)在区间劣用上的一系列离散成节点) a = x0 x1x2 - 上函数值心D的近似值Va (k = 12 ? ? ?, 〃) 而= 就是问题：1)的数值解 从(1)的表达式 W = f^r W = f^ry) axb ［火。)=No (1) 可以看出，求它的数值解的关键在于 兴尤)数值计算问题 或者它的等价的积分方程y(x) = % +「/?(/,y(t))dt中 Ja 积分「/(3(，))山的数值计算问题 J a 求解微分方程的数值方法数值微分 求解微分方程的数值方法 数值微分 数值积分 而数值积分问题我们已经学习过,下考虑数值微分方法 微积分中，关于导数的定义如下： TOC \o 1-5 \h \z 5 h /I。 h /i° 2h 自然而又简单的方法就是，取极限的近似值，即差商! 向前差商 由Taylor展开 f(气 + 力)=f(x°) + hf\x.) + §尸G)/。Mg%x° + h 因此，有误差 R(x)=广(W-EE—E= *%)= O(h) ■向后差商 fM-了 fM- 了（X。） - f（Xo - h）h 误差. R（x）=广 3。）_13。）一。_”）=￡广怎）=o（/i） 中心差商 2hf3o）a 2h h2 h2 R（x） =—[fy^） + 厂 w）] = 7 0）二。（层） 12 o 佑徼分方糕散0解的基本懲您 基本问题：对于微分方程 = axb [y(G = % 要在区间口，国上的若干离散点。=X。V玉 ? ? ? X=人 处计算解函数y(x)的近似值y°双,…，义 实际应用中通常取求解区间[心仞的等分点作为离散点， b —(2 艮卩xk=a+kh. * = 0,1, ???,〃，其中步长* = 一 推导初值问题的数值方法的途径: 求初值问题数值解的方法是步进法，即从已知的初值％ 出发，通过一定的计算求加然后由》或％和以求出力， 依次计算到义,即,在计算出爪? ＜化）后计算以十，这时 Taylor 展开, Taylor 展开, 利用差商离散导数，利用数值积分方法 单步法：只利用以来计算％+】； 多步法:穽f* 1.泰勒展开的求解方法思路： 可将)，(也+/7)=贝％])按泰勒级数展开为 火叫+1)=)3)+人土务 e[xt,xt+i] 略去力2得,yg)? y(也)+ Mm) =y(也)+M (Xp %) 将近似号改等号，则得到数值解序列刃计算公式: Jy(o) =〉o i处+1 =处+你％处)，k=O，l，??『T 2.化导数为差商的求解方法思路： 若在点从处的导数用差商来近似代替，如向前差商 yM yg）/「也 则微分方程初值问题化为 ）（*瓦_ - /（叫，＞（也）），# =0,1,…,〃_1 ? h W）= % 将近似号改等号，精确解y（他）改为近似解序列以满足， i %+1 =光 + hf（xqk），* = o, 1,…- I 3.数值积分的求解方法思路: 如果将微分方程y = f(x, y)化成dy = /(x, y)dx, 然后在各小区间［x,,x,+1］上对其两边进行积分, 即 「「dy = \Xk+l f(x,y)dx,上=0,1,??一1 J * J Xk k =0』,???，〃 一 k =0』,???，〃 一1 yg) 一 y(xk) = f (x, y(x))dx, ffl yk和yM近似y(xQ和y(3),右边用数值积分公式, 如用矩形数值积分公式可得， f 贝。)二 ％ 侦+1 — 乂二母皿損)，SO,1,…,〃T 以上三种方法推导出同一个数值求解公式: \ y（。）= % 侦+i 二 % + 你如刃），s。, 1,? ?