四阶显式Adams法求方程组C语言,第五讲第4章线性多步法(续

《第五讲第4章线性多步法(续》由会员分享，可在线阅读，更多相关《第五讲第4章线性多步法(续(15页珍藏版)》请在人人文库网上搜索。

1、第4章 线性多步法4.1 线性多步法的一般公式前面给出了求解初值问题(1.2.1)的单步法，其特点是计算 时只用到 的值，此时 的值均已算出.如果在计算 时除用 的值外，还用到 的值，这就是多步法.若记,h为步长，则线性多步法可表示为(4.1.1)其中为常数，若(即不同时为零)，称(4.1.1)为线性k步法.计算时用到前面已算出的k个值.当时，(4.1.1)为显式多步方法，当则称(4.1.1)为隐式多步法.隐式方法与梯形方法一样，计算时要用迭代法求.多步法(4.1.1)的局部截断误差定义也与单步法类似.举例来说，对于初值问题，步数k=2时，线性多步法表示为当时，格式为显示的：，而时，格式为隐式。

2、的：。定义4.1设y(x)是初值问题(1.2.1)的精确解，线性多步法(4.1.1)在处的局部截断误差定义为(4.1.2)若，则称线性多步法(4.1.1)是p阶的.如果我们希望得到的多步法是p阶的，则可利用Taylor公式展开，将在处展开到阶，它可表示为(4.1.3)注意，(4.1.2)式按Taylor展开可得经整理比较系数可得(4.1.4)若线性多步法(4.1.1)为p阶，则可令于是得局部截断误差(4.1.5)右端第一项称为局部截断误差主项.称为误差常数.要使多步法(4.1.1)逼近初值问题(1.2.1),方法的阶p1,当p=1时，则，由(4.1.4)得(4.1.6)称为相容性条件.公式(4。

3、.1.1)当k=1时即为单步法，若,由(4.1.6)则得式(4.1.1)就是，即为Euler法.此时，方法为p=1阶.若，由得,为确定及，必须令,由(4.1.4)得及此时(4.1.1)就是即为梯形法.由故p=2,方法是二阶的，与3.1节中给出的结果相同.实际上，当k给定后，则可利用(4.1.4)求出公式(4.1.1)中的系数及，并求得的表达式(4.1.5).4.2 Adams显式与隐式方法形如(4.2.1)的k步法称为 Adams 方法，当 时为 Adams 显式方法，当时，称为Adams隐式方法.对初值问题(1.2.1)的方程两端从到积分得显然只要对右端的积分用插值求积公式，求积节点取为即可。

4、推出形如(4.2.1)的多步法，但这里我们仍采用Taylor展开的方法直