—常微分方程的数值解法.pptxVIP

  1. 1、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。。
  2. 2、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载
  3. 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
  4. 4、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
  5. 5、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们
  6. 6、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
  7. 7、VIP文档为合作方或网友上传,每下载1次, 网站将根据用户上传文档的质量评分、类型等,对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
第六章 常微分方程的数值解法;本章内容;一. 问题提出 有一个或多个导数及其函数的方程式称为微分方程,在工程中常遇到求解微分方程的问题。;§6.1 引言;;实际问题中还存在初边值混合问题,如梁在横向激励下的弯曲振动。 高阶常微分方程可以化成一阶的常微分方程组; 很多微分方程的解不能用初等函数来表示,有时即使能够用解析式表示其解,但计算量太大而不实用(表达式过于复杂)。 需要用数值方法来求解,一般只要求得到若干个点上的近似值或者解的简单的近似表达式(精度要求满足即可)。;§6.1 引言;;初值 问题 的 常见 解法;§6.2 欧拉方法及其改进 Euler’s Method ; 6.2.1 欧拉公式:/* Euler’s Method */;;§6.2 欧拉方法及其改进;§7.2 欧拉方法;§6.2 欧拉方法及其改进;§7.2 欧拉方法;§7.2 欧拉方法;§6.2 欧拉方法及其改进;§6.2 欧拉方法及其改进;§6.2 欧拉方法及其改进;显、隐式两种算法的平均。需要迭代求解,能否不迭代?; ;§7.2 欧拉方法及其改进;§7.2 欧拉方法;§7.2 欧拉方法;§6.3 龙格 —库塔方法;单步法:即利用前一个节点的函数值yi,计算后一个节点的函数值yi+1。 目的:建立高精度的单步递推格式。 单步递推法的基本思想是从 ( xi , yi ) 点出发,以某一斜率沿直线达到 ( xi+1 , yi+1 ) 点。欧拉法及其各种变形所能达到的最高精度为2 阶。;斜率一定取K1 K2 的平均值吗?;§7.3 龙格 —库塔方法;Step 3: 将 y( xi+1 ) 在 xi 点的泰勒展开并与yi+1作比较;其中?i ( i = 1, …, m ),?i ( i = 2, …, m ) 和 ?ij ( i = 2, …, m; j = 1, …, i?1 ) 均为待定系数,确定这些系数的步骤与前面相似。 ;3阶龙格-库塔法;§6.3 龙格 —库塔方法;注: ? 龙格-库塔法的主要运算在于计算 Ki 的值,即计算 f 的值。Butcher 于1965年给出了计算量与可达到的最高精度阶数的关系:;§7.3 龙格 —库塔方法;§7.3 龙格 —库塔方法;§7.3 龙格 —库塔方法;§7.3 龙格 —库塔方法; 步长过大,达不到精度要求;步长过小,虽然局部截断 误差小,加大了计算工作量,舍入误差的累积增大。 解决途径之一——引入变步长技术,常用的有Richardson外 推法。 从结点 xi出发,先以h为步长,通过一步计算出y(xi+1)的近似值;§7.2 欧拉方法;例:就初值问题 考察欧拉显式格式的收敛性。;§7.2 欧拉方法;§7.2 欧拉方法;一般分析时为简单起见,只考虑试验方程;§6.5 收敛性与稳定性;;§7.3 龙格 —库塔方法;§6.6 一阶常微分方???组与高阶方程 ;设 为节点上的近似解, 则有改进的Euler格式为 ;;;;6.6.1 一阶常微分方程组;例 求解下列二阶微分方程的初值问题 ;;;§6.7 边值问题的数值解法

文档评论(0)

131****9843 + 关注
实名认证
文档贡献者

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

1亿VIP精品文档

相关文档