logo

您所在位置网站首页 > 海量文档  > 教学课件 > 大学课件

计算方法第5章常微分方程初值问题的数值解法.ppt 27页

本文档一共被下载: ,您可全文免费在线阅读后下载本文档。

  • 支付并下载
  • 收藏该文档
  • 百度一下本文档
  • 修改文档简介
全屏预览

下载提示

1.本站不保证该用户上传的文档完整性,不预览、不比对内容而直接下载产生的反悔问题本站不予受理。
2.该文档所得收入(下载+内容+预览三)归上传者、原创者。
3.登录后可充值,立即自动返金币,充值渠道很便利
特别说明: 下载前务必先预览,自己验证一下是不是你要下载的文档。
  • 上传作者 autohhh(上传创作收益人)
  • 发布时间:2018-07-12
  • 需要金币160(10金币=人民币1元)
  • 浏览人气
  • 下载次数
  • 收藏次数
  • 文件大小:591 KB
下载过该文档的会员
你可能关注的文档:
上一页 下一页 返回首页 第5章 常微分方程初值问题的数值解法 5.1 数值解法的基本思路 5.2 欧拉方法 5.3 龙格-库塔法 5.4 线性多步法 5.1 数值解法的基本思路 所谓初值问题的常微分方程 是指: 为了便于计算,通常在区间[a,b]上作等步长处理,即 思路 从 出发求出 的方法有 1.化导数为差商 2.数值积分法 对于积分 可以用定积分的数值积分的 方法求得,用不同的数值求积公式,即可获得不同的微分方程数值解的不同公式. 若 若 对数值积分若用梯形公式: 3.泰勒(Taylor)展开法 不考虑余项,则有 两边同时求导得 当n=1时 5.2 欧拉方法 5.2.1 基本公式 考虑一阶常微分方程的初值问题 若采用数值积分法,便有 或 和 欧拉公式 后退欧拉公式(隐形) 梯形公式(隐形) 5.2.2 欧拉公式的几何解释 = 5.2.3欧拉方法的误差估计 根据泰勒公式,有 1.欧拉公式 或者写成 2.梯形公式 由泰勒公式有 与 相减 欧拉公式的误差是 梯形公式其每步的截断误差为 比欧拉公式要高,但是梯形公式是隐式的不易计算 5.2.4 改进的欧拉公式 ?先用欧拉公式 对 求出一个初步值 预报值 然后代入式 中,求出校正的 ? 可把它改写成下列形式: 即有 例1 取步长h=0.1,求初值问题 在x=0到x=1各节点上的数值解。 5.3 龙格-库塔法 5.3.1 原理 由微分中值定理 称K为在区间 上的平均斜率 只要平均斜率给出一个算法,就可以得到相应的计算 将改进欧拉法推广为: 相当于 5.3.2 龙格—库塔公式 思路 在区间 上再增加一个或多个预报点,即用三个或更多点上的斜率进行加权平均,作为平均斜率 的近似值,可望得到精度更高的计算公式。 用三个点进行预报: ? ? ? ? ? 由泰勒展开 又有 比较两式 代入前面的公式 每步的截断误差为 令p=1/2(中点),q=1(终点) 令p=1/3,q=2/3 类似的可得 经典龙格—库塔公式(也称为古典形式):其每一步的截断误差为 还有库塔公式 例 用四阶龙格—库塔公式解以下微分方程的初值问题,并与改进的欧拉公式、准确值比较(取步长0.2): 解 对本题的四阶龙格—库塔公式具体形式为 上一页 下一页 返回首页

发表评论

请自觉遵守互联网相关的政策法规,严禁发布色情、暴力、反动的言论。
用户名: 验证码: 点击我更换图片

“原创力文档”前称为“文档投稿赚钱网”,本站为“文档C2C交易模式”,即用户上传的文档直接卖给(下载)用户,本站只是中间服务平台,本站所有文档下载所得的收益归上传人(含作者)所有【成交的100%(原创)】。原创力文档是网络服务平台方,若您的权利被侵害,侵权客服QQ:3005833200 电话:19940600175 欢迎举报,上传者QQ群:784321556