第七章常微分方程数值解.docxVIP

(7.2.2) PAGE (7.2.2) PAGE # 第七章常微分方程数值解 7.1 引言 本章讨论常微分方程初值问题 \y(a) = y0 (7.1.1) 的数值解法，这也是科学与工程计算经常遇到的问题，由于只有很特殊的方 程能用解析方法求解，而用计算机求解常微分方程的初值问题都要采用数值方法. 通常我们假定(7.1.1)中f(x,y)对y满足 i tz条件，即存在常数L0,使 (7.1.2) 则初值问题(7.1.1)的解存在唯一. 假定(7.1.1)的精确解为尸5),求它的数值解就 是要在区间上的一〔④3组离 散点上。=而(为…h求尸⑸的近似外‘无」，以….通常取 /二鼻+叫…)，h称为步长，求(7.1.1)的数值解是按节点的顺餐■ = 12…) 序逐步推进求得巧?》….首先，要对方程做离散逼近，求出数值解的公式，再研 究公式的局部截断误差，计算稳定性 以及数值解的收敛性与整体误差等问题. 7.2 简单的单步法及基本概念 7.2.1 法、后退 er法与梯形法 求初值问题(7.1.1)的一种最简单方法是将节点的导数.小‘(/)用差商代 替必色上幽士 于是(7.1.1)的方程可近似写成 h ) * y(xn) + hAuX4)), 口 =。, L … (7.2.1) 从孙出发M㈤二贝工。尸￥。，由(7.2.1)求得再将代.石尸九十皆(而，%)=乃 乃指尸(勺)入(7.2.1)右端，得到的近似A氏)乃=总+研和为)，一般写成 称为解初值问题的l er法. PAGE PAGE # (7.2.5) PAGE (7.2.5) PAGE # 法的几何意义如图7-1所示.初值问题(7.1.1)的解曲线=y(x)过点 P式殉，尸”，从而出发，以为斜率作 六汽，为)一段直线，与直线交点》 =勺于耳(看,修), 显然有泗=儿十型。，为)，再从写出发，以为斜率作，⑶，力)直线推进到钎 必上 点斗，其余类推，这样得到解 曲线的一条近似曲线，它就是折线此片舄…. (723)y(xn +h) = y(zn) + hy1(xJ + ^-y(^),短日(/,^+1) (723) 略去余项，以为、道「)，就得到近似计算公式(7.2.2). 另外，还可对(7.1.1)的方程两端由到积分得之工田 X^i) - - ￡* X^i) - - ￡* jO))心 (7.2.4) 若右端积分用左矩形公式，用人“,居*1到W%】),则得(7.2.2). 如果在(7.2.4)的积分中用右矩形公式，则得 称为后退(隐式 法.若在(7.2.4)的积分中用梯形公式，则得 jiij - ,二)I 八4丁二一 J - U - (7.2.6) 称为梯形方法. 上述三个公式(7.2.2) ,(7.2.5)及(7.2.6)都是由计算八八十】，这种只用前一 步即可算出的公式称％+i为单步法，其中(7.2.2)可由逐次求K出尸12^…的值， 称为显式方法，而(7.2.5)及(7.2.6)右端含有当料Jf对y非线性时它不 能直接求出八％此时应把它看作一个方程，求解/川，这类方法称为稳式方法. 此时可将(7.2.5)或(7.2.6)写成不动点形式的方程 yn+i = W(^+i-A+i) + ^ 这里对式(7.2.5)有/L乌?瑞，对(7.2.6)则/==八十号力//)， g与八十1无关，可构造迭代法 匚「：士 … (7.2.7) 由于对y满了5足条件(7.1.2)，故有 卜整rj ＜明/a用园J 一丁冗+qz) M歌斓-八」 当瞅＜1或卜 ＜ 圭,迭代法(7.2.7)收敛到九川，因此只要步长h足够小，就 可保证迭代(7.2.7)收敛.对后退 l er法(7.2.5),当时迭代收h ＜ :敛，对梯形 JLr 、， ? 法(7.2.6)，当时迭代序为＜ 7列收敛. 例7.1用 r法、隐式 er法、梯形法解 yi = - y + K + L y(0) = 1 取h=0.1 ,计算到x=0.5,并与精确解 比较. 解 本题可直接用给出公式计算.由于 氏七力=一y +工+ 1』=0」,礴二5yL 1, 法的计算公式为 rs+i =^s + K-yx + / +1) =(1 -增八 + 比［+血=0-9办 + 0 1/ +0.1 n=0时，外?0一9乂 +0.H+0./1 .n=1,2,3,4的计算结果见表 7-1. 对隐式ler法，计算公式为 = Lm + A(-A+l 十五”1 + 1) 解出 八「7^7? +砧+1十/=AM十o X十0 11)1+ ? 1.1 当n=0时，丸■ .5+0 M +0H)= 1009091.其余n=1,2,3,4的计算结果 见表7-1. 表7-1例7.1的三种方法及精确解的计算结果 段 Euler 法 yn 隐式Euler法yn 梯形