第8章 常微分方程数值解法.pptVIP

第8章 常微分方程数值解法; 所谓数值解法,就是设法将常微分方程离散化,建立差分方程,给出解在一些离散点上的近似值.;对右边的积分应用左矩形公式，则有;梯形公式;对右边的积分应用左矩形公式，则有; 利用Euler方法求初值问题;分别取步长h=0.2 ,0.1 ,0.05,计算结果如下;h;Euler中点公式则不然, 计算yn+1时需用到前两步的值yn , yn-1 ,称其为两步方法,两步以上的方法统称为多步法.; 从数值积分的角度来看,梯形公式; 由迭代法收敛的角度看，当; y?=y-2x/y , 0?x?1;计算结果如下:;n; 在节点xn+1的误差y(xn+1)-yn+1 ,不仅与yn+1这一步计算有关,而且与前n步计算值yn,yn-1,…,y1都有关. ;另外,在yn=y(xn)的条件下,考虑到y?(x)=?(x,y(x)),则有; yn+1=yn+h?(xn,yn);所以, 改进的Euler方法是二阶方法.;称之为p阶Taylor展开方法.;所以,此差分公式是p阶方法. ; 构造差分公式;由于; 一般地, 参数由(8.4)确定的一族差分公式(8.3)统称为二阶R-K方法.; 四阶标准R-K公式; 解 四阶标准R-K公式为;n;是二级二阶隐式R-K方法,也就是梯形公式.但是p级隐式R-K方法的阶可以大于p,例如,一级隐式中点公式为; 如果将步长减半,取h/2为步长, 从xn经两步计算得到y(xn+1)的近似值记为; 求解初值问题; 设y(x)是初值问题(8.1)的解 ,yn是单步法 (8.5)产生的近似解.如果对任意固定的点xn,均有;其中,局部截断误差|Rn(h)|?C1hp+1, 记en=y(xn)-yn ,则有;由于e0=y(a)-y0=0, 所以有;条件,因此改进的Euler方法是收敛的. ; 将Euler方法应用于方程y?=?y, 得到 ;类似前面分析,可知绝对稳定区域为; 解 因y0=1,计算得y10=1024,而y(1)=9.357623?10-14.; 单步显式方法的稳定性与步长密切相关, 在一种步长下是稳定的差分公式,取大一点步长就可能是不稳定的.;当?-1?0时,公式为隐式公式,反之为显式公式.参数??i,?i?的选择原则是使方法的局部截断误差为; ?n=?(xn,y(xn))=y?(xn) ; ?=1, ?0+?1+?2=1, ?1+2?2=-1/2 , ?1+4?2=1/3; 选取不同的插值多项式pr(x),就可导出不同的差分公式.下面介绍常用的Adams公式.;由此,可建立差分公式 ;下面列出几个带有局部截断误差主项的Adams显式公式 ;其中系数为 ; 3.Adams预估-校正公式 ; y?=y-2x/y , 0?x?1;练习题