[理学]数值分析第9章5-6节.pptVIP

  1. 1、本文档共63页,可阅读全部内容。
  2. 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
  3. 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载
  4. 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
  5. 5、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
  6. 6、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们
  7. 7、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
  8. 8、VIP文档为合作方或网友上传,每下载1次, 网站将根据用户上传文档的质量评分、类型等,对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
[理学]数值分析第9章5-6节

9.5 线性多步法 9.5.1 线性多步法的一般公式 9.5.2 阿当姆斯显式与隐式公式 9.5.3 米尔尼方法与辛普森方法 9.5.4 汉明方法 9.5.5 预测-校正方法 9.5.6 构造多步法公式的注记和例 9.6 方程组和高阶方程 9.6.2 化高阶方程为一阶方程组 9.6.3 刚性方程组 利用米尔尼公式(5.11)和汉明公式(5.15)相匹配, 并利用截断误差(5.12),(5.16)改进计算结果,可类 似地建立四阶修正米尔尼-汉明预测-校正格式PMECME): P: M: E: C: M: E: (5.15) (5.11) (5.16) (5.12) 将例6的初值问题用修正的米尔尼-汉明预测-校 正公式计算 及 ,初值 仍用已算出的精 确解, 例7 即 给出计算结果及误差. 解 其中 根据修正的米尔尼-汉明预测-校正公式可得 误差 从结果看,此方法的误差比四阶阿当姆斯隐式法和四阶汉 明方法小,这与理论分析一致. 构造多步法公式有基于数值积分和泰勒展开两种途径, 而用泰勒展开则可构造任意多步法公式,其做法是根 据多步法公式的形式,直接在 处做泰勒展开. 确定多步法(5.1)的系数 及 时不 必套用系数公式(5.4),因为多步法公式不一定如(5.1) 的形式, 数值积分方法只能用于可转化为等价的积分方程的微分方程, 它是有局限性的. 而且套用公式容易记错. (5.1) 例8 用显式二步法 其中 试确定参数 使方法阶数尽可能高,并求局部 截断误差. 解 解初值问题 根据局部截断误差定义,用泰勒展开确定参数满足的方程. 由于 为求参数 使方法阶数尽量高,可令 即得方程组 解得 ,此时公式为三阶的, 即为所求局部截断误差. 而所得二步法为 而且 例9 证明存在 的一个值,使线性多步法 是四阶的. 证明 只要证明局部截断误差 ,则方法 仍用泰勒展开,由于 为四阶. 当 时, ,故方法是四阶的. 9.6.1 一阶方程组 前面研究了单个方程 的数值解法,只要把 和 理解为向量,那么,所提供的各种计算公式即可应用 到一阶方程组的情形. 考察一阶方程组 的初值问题,初始条件给为 则上述方程组的初值问题可表示为 (6.1) 求解这一初值问题的四阶龙格-库塔公式为 式中 若采用向量的记号,记 或表示为 其中 这里 是第 个因变量 在节点 的近似值. 考察两个方程的特殊情形: 这时四阶龙格-库塔公式具有形式 其中 (6.2) (6.3) 这是一步法,利用节点 上的值 ,由(6.3) 式顺序计算 ,然后代入(6.2) 式即可求得节点 上的 . (6.3) 高阶微分方程(或方程组)的初值问题,原则上总可 以归结为一阶方程组来求解. 例如,考察下列 阶微分方程 (6.4) 初始条件为 (6.5) 只要引进新的变量 即可将 阶方程(6.4)化为如下的一阶方程组: (6.6) 初始条件(6.5)则相应地化为 (6.7) 初值问题(6.4),(6.5)和(6.6),(6.7)是彼此等价的. (6.4) (6.5) 特别地,对于下列二阶方程的初值问题: 引进新的变量 ,即可化为下列一阶方程组的初值问题: 在逐步推进的求解过程中,计算 之前事实上已经 求出了一系列的近似值 如果充分利用前面多步的信息来预测 ,则可以期望获得较高的精度. 这就是构造所谓线性多步法的基本思想. 本节主要介绍基于泰勒展开的构造方法. 构造多步法的主要途径是基于数值积分方法和基于泰 勒展开方法,前者可直接由方程 两端积分后 利用插值求积公式得到. 一般的线性多步法公式可表示为 (5.1) 其中 为 的近似, 计算时需先给出前面 个近似值

文档评论(0)

jiupshaieuk12 + 关注
实名认证
文档贡献者

该用户很懒,什么也没介绍

版权声明书
用户编号:6212135231000003

1亿VIP精品文档

相关文档