数值分析第八章 常微分方程数值解法new.pptVIP

由于e0=y(a)-y0=0, 所以有 则有 设?(x,y)连续且关于y满足Lipschitz条件,对于Euler方法,由于?(x,y,h)=?(x,y),故Euler方法是收敛的. 对于改进的Euler方法,由?(x,y)的Lipschitz条件有 |en|?|e0|eL(b-a)+C1hp/L(eL(b-a)-1) |en|?C1hp/L(eL(b-a)-1)=Chp |?(x,y,h)-?(x,y*,h)|?1/2|?(x,y)-?(x,y*)| +1/2|?(x+h,y+h?(x,y))-?(x+h,y*+h?(x,y*))| ?1/2L(2+hL)|y-y*| 则当h?h0时,?关于y满足常数为1/2L(1+h0L)的Lipschitz 条件,因此改进的Euler方法是收敛的. 可类似验证各阶R-K方法是收敛的. §4.2 单步方法的稳定性 定义8.2 对于初值问题(8.1),取定步长h,用某个差分方法进行计算时,假设只在一个节点值yn上产生计算误差?,即计算值?yn=yn+?, 如果这个误差引起以后各节点值ym(mn)的变化均不超过? ,则称此差分方法是绝对稳定的.