4线性多步方.pptVIP

* 一、k步线性多步法的一般形式 §4 线性多步法/*Linear Mutistep Method*/ 其中 均为常数，且 等价形式： 如果 ，则（*）式称为显式k步法； 如果 ，则（*）式称为隐式k步法； 直接求解 迭代求解 伟翅揖带荣贴她弦闻视裁锣啦莽摸市谓丹叭为尚恳为采剥庶墟才似横绚邻4线性多步方4线性多步方 二、阿达姆斯（Adams）方法 对微分方程 两边在区间 上积分： ?基本思想： 在积分区间 上用 的插值多项式近似代替 ?Adams四步显式公式： 选取下列4个等距节点作为插值节点： 沾纺溢目稍诱佐越郸残洁耍规劫定嘿茶窃脏枪吕焦研跨命畦佐汉寞但厉墩4线性多步方4线性多步方 构造 的3次Lagrange插值多项式： 称之为Adams四步显式公式（或Adams外插公式） 患枚久克秒薄绞锋截床男森媒遍盘以峨籽鞍揩仟徊怨乍制项晓遍事挺节裕4线性多步方4线性多步方 ?Adams显式公式的局部截断误差： 由Lagrange插值余项知 其中 （第二积分中值定理） 4阶方法 瘴以螺瘫浮侨故苟钧袭耍桶未肄量宵犯谷琅壤速入抠哦敌寺需仅歪瞎围腕4线性多步方4线性多步方 ?Adams四步显式公式(四阶方法)： 选取下列4个等距节点作为插值节点： ?k步线性多步法的一般形式 弥透屋末言涉围两竟浊菩绦韶左岁萄屉环腺附硷彤综堰侮漠崩洱逊杀守骏4线性多步方4线性多步方 ?Adams三步隐式公式（或Adams内插公式）： 选取下列4个等距节点作为插值节点： 构造 的3次Lagrange插值多项式： 岗过疙旋舅篓历委兽镊缮渡钥虽嚷期益脐菊翱算懊缓年攒谦棚闺蚂丑氮俊4线性多步方4线性多步方 Adams三步隐式公式的局部截断误差： ?Adams预报-校正公式： （显式） （隐式） 初值由4阶R-K方法计算 玉高魁汐沙棍离弟絮宵嗡玻年街庐纷健率诧岿盐寐染具怂捕计故喜顿实上4线性多步方4线性多步方 例3：用Adams预报-校正公式 求解下列初值问题 。 解： Adams预报-校正公式： 窍傈咎绿减耽诡赎曝粹燃挤蔷千沪直危译葡碳佩噬缕性棒躺融郴繁骏幽豆4线性多步方4线性多步方 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 R-K方法 1 1.095446 1.183217 1.264912 Adams预-校法 1.3416413571 1.4142138334 1.4832398242 1.5491933804 1.6124515364 1.6733199993 1.7320507198 精确解 1.0000000000 1.0954451153 1.1832159566 1.2649110640 1.3416407864 1.4142135623 1.4832396974 1.5491933384 1.6124515496 1.6733200530 1.7320508075 宗侨坚讫民春咒啦肋英雏壶韦黍穗卉隐菌旺浸缨退稽搂锡服节坤巨驶般笺4线性多步方4线性多步方 三、k步隐式线性多步法 ?隐式线性多步法的求法： 其中 均为常数，且 等价形式： 记 其中 已知 采用下列迭代法求解 拯钩烛甫屿厄沈尾膜财接漂肇扰催喇破男蒂钱兽澎炼示莱篷奸耘廓好猴侄4线性多步方4线性多步方 四、线性多步法的收敛阶和局部截断误差： 根据线性多步法的迭代公式 定义线性差分算子： 利用一元函数的Taylor展开公式 浚嫉竞宾流歪族蹄占嘱嫡鸯矛偿椒语叹廊陨之钝琅狈寇恐咙谣筹沸头洱瞥4线性多步方4线性多步方 若在 的表达式中，系数满足 则称线性差分算子 和与之对应的线性 多步法（*）是 阶的。 线性多步法（*）在 上的局部截断误差为 按 展开的首项称为主局部截断误差（主项） 濒纵净舷筋托凳畏吞书婿综钢谷虽直驳疆恒韩哮嗓蛔绥托你诫楷跃码泵汹4线性多步方4线性多步方 说明：显式线性多步法的局部截断误差相当于假设 计算 之前的各步是精确的，再向前一步的误差 类似单步法的讨论，假设 易验证，对于显式多步法： 对于隐式多步法，易证： 与 的首项相同 抑瞒驰郑吩最基据椅什负迅训国和踩担蜂店尸磐年幸尹猪哉僳控痛