称为绝对稳定域.pptVIP

9.4单步法的收敛性与稳定性 9.4.1收敛性与相容性 数值解法的基本思想是,通过某种离散化手段将微分 方程(1.1)转化为差分方程,如单步法(2.10),即 ym+=yn+hop(xn, yn,h), (4.1) 它在xn处的解为yn,而初值问题(1.1),(1.2)在xn 处的精确解为y(xn),记en=y(xn)-yn称为整体截断 误差 收敛性就是讨论当x=x因定且h=x O时 n→0的问题  定义3若一种数值方法(如单步法(4.1))对于固定 的xn=x+mh,当→O时有,一y(x),其(x) 是(1.1),(1.2)的准确解,则称该方法是收敛的. 显然数值方法收敛是指cn=y(xn)-yn→O对单步 法(4.1)有下述收敛性定理 定理1假设单步法(4.1)具打阶精度,且增量函数 q(x,y,h)关于y满足利普希茨条件 9(x,y,h)-∞(x,y,h)≤L列ly (4.2) 又设初值y是准确的,即y=y(x),则其整体截断误差 y(xn)-y=O(h) 4.3 2  证明设以歹表示取yn=y(xn)用公式(4.1)求 得的结果,即 (n)+hop(n, y(x),h) 则y(xn)一yn为局部截断误差,由于所给方法具有p阶 精度,按定义2,存在定数C,使 V(x ch 又由式(4.4)与(4.1),得 yn-ynH|≤y(x)-yn +ho(xn, y(x,), h)-(,yn, h) 利用假设条件(4.2),有  yn-yn≤(1+hLy(x,)-y 从而有 y(xn+)-uu+i -yn++y(xn+)-y ≤(1+hL)y(xn)-yn|+Ch” 即对整体截断误差en=y(x)-y,成立下列递推关系式 en|≤(1+hL)en|+Ch (4.5) 反复递推,可得 len≤(1+hL)c|+ Chp (1+hLp)-1 (4.6)  再注意到当xn-x0=mh≤7时 (1+ML)≤(e)≤e, 最终得下列估计式 Chl en slope o+ (e-1) (4.7) 由此可以断定,如果初值是准确的,即eo=O,则(4.3) 式成立 依据这一定理,判断单步法(4.1)的收敛性,归结为 验证増量函数φ能否满足利普希茨条件(4.2) 对于欧拉方法,由于其增量函数φ就是∫(x,y),故当 f∫(x,y)关于y满足利普希茨条件时它是收敛的  再考察改进的欧拉方法,其增量函数由(3.2)式 PP(xu, yu, h)==[f(n, y,)+f(x,+h, y, + hf(xn, yu)) 给出,这时有 lp(x,y, h)-(x,y,hs[f(x, y)f(x, y) If(x+,y+hf(x, y))-f(x+h,y+hf(x,y))I 假设∫(xy)关于y满足利普希茨条件,记利普希茨常数为 L,则由上式推得 PP(x,y,h)-(x,y,h)L(+Ly-yl 设限定h≤h0(hn为定数),上式表明卯关于y的利普希茨  常数 L=L(1+L), 因此改进的欧拉方法也是收敛的 类似地,不难验证其他龙格-庳塔方法的收敛性. 定理1表明P≥1时单步法收敛,并且当y(x是初值 问题(1.1),(1.2)的解,(4.1)具有阶精度时,则 有展开式 m+= y(x+h)-y(x)-hop(x,y(x),h) =y(h+y(x), hlp(, y(x), 0)+p(x, y(x), O)h+... hly(x)-(x, y(x),0)]+O(h-)  所以p≥1的充要条件是y(x)-q(x,y(x),O)=0,而 y(x)=f(x,y(x)),于是可给出如下定义: 定义4若单步法(4.1)的增量函麴满足 PP(x,y, 0)=f(x, y), 则称单步法(4.1)与初值问题(1.1),(1.2)相容 以上讨论表明P阶方法(4.1)当P≥1时与(1.1),(1.2) 相容,反之相容方法至少是1阶的 于是由定理1可知方法(4.1)收敛的充分必要条件是此 方法是相容的  9.4.2绝对稳定性与绝对稳定域 定义5若一种数值方法在节点值yn上大小为S的扰动 于以后各节点值yn(mn)上产生的偏差均不超过δ,则称 该方法是稳定的 以欧拉法为例考察计算稳定性 例4考察初值问题 y=-100y, y(O)=1 其准确解y(x)=1∞是一个按指数曲线衰减得很快的函数, 如图9-3所示 用欧拉法解方程y’=-100y得  (1-100h)y 若取h=O.025则欧拉公式的 具体形式为 y 1.5yn, 0.025005000750.100 计算