电子文档交易市场
安卓APP | ios版本
电子文档交易市场
安卓APP | ios版本

龙格-库塔方法基本原理.ppt

54页
  • 卖家[上传人]:自***
  • 文档编号:127521007
  • 上传时间:2020-04-03
  • 文档格式:PPT
  • 文档大小:798KB
  • / 54 举报 版权申诉 马上下载
  • 文本预览
  • 下载提示
  • 常见问题
    • 1、在连续系统的仿真中 主要的计算工作是求解一阶微分方程y f x y y x0 y0解析法只能用来求解一些特殊类型的方程 实际仿真问题中归结出来的微分方程主要靠数值解法 由于实际运算只能完成有限项或有限步运算 因此要将有些需用极限或无穷过程进行的运算有限化 对无穷过程进行截断 这样产生的误差成为截断误差 根据实际情况建立的数学模型往往难以求解 通常需要通过近似替代 将所求解的数学模型简化为易求解的数值计算问题后再进行求解 数学模型的理论解与数值计算问题的精确解之间的误差称为截断误差 这是计算方法本身带来的误差 所以也成为方法误差 2020 4 3 3 得到高精度方法的一个直接想法是利用Taylor展开假设式y f x y a x b 中的f x y 充分光滑 将y xi 1 在xi点作Taylor展开 若取右端不同的有限项作为y xi 1 的近似值 就可得到计算y xi 1 的各种不同截断误差的数值公式 例如 取前两项可得到 3 2龙格 库塔方法 2020 4 3 4 其中 P阶泰勒方法 若取前三项 可得到截断误差为O h3 的公式 类似地 若取前P 1项作为y xi 1 的近似值 便得

      2、到 2020 4 3 5 显然p 1时 yi 1 yi hf xi yi 它即为我们熟悉的Euler方法 当p 2时 要利用泰勒方法就需要计算f x y 的高阶微商 这个计算量是很大的 尤其当f x y 较复杂时 其高阶导数会很复杂 因此 利用泰勒公式构造高阶公式是不实用的 但是泰勒级数展开法的基本思想是许多数值方法的基础 R K方法不是直接使用Taylor级数 而是利用它的思想 2020 4 3 6 龙格 库塔 R K 法的基本思想 Euler公式可改写成 则yi 1的表达式与y xi 1 的Taylor展开式的前两项完全相同 即局部截断误差为O h2 Runge Kutta方法是一种高精度的单步法 简称R K法 2020 4 3 7 同理 改进Euler公式可改写成 上述两组公式在形式上共同点 都是用f x y 在某些点上值的线性组合得出y xi 1 的近似值yi 1 且增加计算的次数f x y 的次数 可提高截断误差的阶 如欧拉法 每步计算一次f x y 的值 为一阶方法 改进欧拉法需计算两次f x y 的值 为二阶方法 局部截断误差为O h3 2020 4 3 8 于是可考虑用

      3、函数f x y 在若干点上的函数值的线性组合来构造近似公式 构造时要求近似公式在 xi yi 处的Taylor展开式与解y x 在xi处的Taylor展开式的前面几项重合 从而使近似公式达到所需要的阶数 既避免求高阶导数 又提高了计算方法精度的阶数 或者说 在 xi xi 1 这一步内多计算几个点的斜率值 然后将其进行加权平均作为平均斜率 则可构造出更高精度的计算格式 这就是龙格 库塔 Runge Kutta 法的基本思想 一般龙格 库塔方法的形式为 2020 4 3 9 其中ai bij ci为待定参数 要求上式yi 1在点 xi yi 处作Tailor展开 通过相同项的系数确定参数 称为P阶龙格 库塔方法 10 Runge Kutta方法的推导思想 对于常微分方程的初值问题 的解y y x 在区间 xi xi 1 上使用微分中值定理 有 即 2020 4 3 11 引入记号 就可得到相应的Runge Kutta方法 2020 4 3 12 如下图 即 则上式化为 即Euler方法 Euler方法也称为一阶Runge Kutta方法 2020 4 3 二阶龙格 库塔法在 xi xi

      4、1 上取两点xi和xi a2 xi a2h 以该两点处的斜率值K1和K2的加权平均 或称为线性组合 来求取平均斜率k 的近似值K 即 式中 K1为xi点处的切线斜率值K1 hf xi yi hy xi K2为xi a2h点处的切线斜率值 比照改进的欧拉法 将xi a2视为xi 1 即可得 2020 4 3 13 确定系数c1 c2 a2 b21 可得到有2阶精度的算法格式 2020 4 3 14 因此 将y xi 1 在x xi处进行Taylor展开 将在x xi处进行Taylor展开 2020 4 3 15 K1 hf xi yi 2020 4 3 16 这里有4个未知数 3个方程 存在无穷多个解 所有满足上式的格式统称为2阶龙格 库塔格式 令 对应项的系数相等 得到 2020 4 3 17 注意到 就是二阶龙格 库塔公式 也就是改进的欧拉法 因此 凡满足条件式有一簇形如上式的计算格式 这些格式统称为二阶龙格 库塔格式 因此改进的欧拉格式是众多的二阶龙格 库塔法中的一种特殊格式 若取 就是另一种形式的二阶龙格 库塔公式 2020 4 3 18 此计算公式称为变形的二阶龙格 库塔法 式

      5、中为区间的中点 也称中点公式 Q 为获得更高的精度 应该如何进一步推广 2020 4 3 19 二级R K方法是显式单步式 每前进一步需要计算两个函数值 由上面的讨论可知 适当选择四个参数c1 c2 a2 b21 可使每步计算两次函数值的二阶R K方法达到二阶精度 能否在计算函数值次数不变的情况下 通过选择不同的参数值 使得二阶R K方法的精度再提高呢 答案是否定的 无论四个参数怎样选择 都不能使公式的局部截断误差提高到三阶 这说明每一步计算两个函数值的二阶R K方法最高阶为二阶 若要获得更高阶得数值方法 就必须增加计算函数值的次数 三阶龙格 库塔法 2020 4 3 20 为进一步提高精度 在区间 xi xi 1 上除两点xi和xi a2 xi a2h 以外 再增加一点xi a3 xi a3h 用这三点处的斜率值K1 K2和K3的加权平均得出平均斜率K 的近似值K 这时计算格式具有形式 2020 4 3 21 同理推导二阶公式 将y xi 1 和yi 1在x xi处进行Taylor展开 使局部截断误差达到O h4 使对应项的系数相等 得到系数方程组 参数的选择不唯一 从而构成一类不同

      6、的三阶R K公式 下面给出一种常用的三阶R K公式 形似simpson公式 2020 4 3 22 2020 4 3 23 四阶 经典 龙格 库塔法 如果需要再提高精度 用类似上述的处理方法 只需在区间 xi xi 1 上用四个点处的斜率加权平均作为平均斜率K 的近似值 构成一系列四阶龙格 库塔公式 具有四阶精度 即局部截断误差是O h5 推导过程与前面类似 由于过程复杂 这里从略 只介绍最常用的一种四阶经典龙格 库塔公式 2020 4 3 24 K1 hf xi yi K2 hf xi a2h yi b21K1 K3 hf xi a3h yi b31K1 b32K2 K4 hf xi a4h yi b41K1 b42K2 b43K3 其中c1 c2 c3 c4 a2 a3 a4 b21 b31 b32 b41 b42 b43均为待定系数 这里K1 K2 K3 K4为四个不同点上的函数值 分别设其为 设 yi 1 yi c1K1 c2K2 c3K3 c4K4 2020 4 3 25 类似于前面的讨论 把K2 K3 K4分别在xi点展成h的幂级数 代入线性组合式中 将得到的公式与y xi

      7、 1 在xi点上的泰勒展开式比较 使其两式右端直到h4的系数相等 经过较复杂的解方程过程便可得到关于ci ai bij的一组特解a2 a3 b21 b32 1 2b31 b41 b42 0a4 b43 1c1 c4 1 6c2 c3 1 3 26 四阶 经典 Runge Kutta方法 2020 4 3 27 例1 使用高阶R K方法计算初值问题 解 1 使用三阶R K方法 2020 4 3 28 其余结果如下 2 如果使用四阶R K方法 ixik1k2k3yi1 00000 10000 10000 11030 12561 11112 00000 20000 12350 13760 15951 24993 00000 30000 15620 17640 20921 42844 00000 40000 20400 23420 28661 66645 00000 50000 27770 32590 41631 9993 2020 4 3 29 其余结果如下 ixik1k2k3k4yi1 00000 10000 10000 11030 11130 12351 11112 00000 2000

      8、0 12350 13760 13920 15631 25003 00000 30000 15620 17640 17910 20421 42864 00000 40000 20400 23420 23890 27811 66675 00000 50000 27770 32590 33480 40062 0000 2020 4 3 例 已知一阶系统的微分方程为 初始条件 取仿真步长h 0 1 分别用欧拉法 梯形法和龙格 库塔法计算该系统仿真第一步的值 解 原方程可变为 即 数值积分公式应用 1 用欧拉法计算根据欧拉公式 将函数表达式及其初始值代入后 可得该系统仿真第一步的值 数值积分公式应用 2 用梯形法计算 根据预报 校正公式 将函数表达式及其初始值代入后 可得仿真第一步的值 用预报公式求起始值 数值积分公式应用 再用校正公式得到系统仿真第一步的值 数值积分公式应用 二阶龙格 库塔公式 3 用二阶龙格 库塔法计算根据公式先计算出两个系数 再计算仿真第一步的值 数值积分公式应用 则系统仿真第一步的值为 数值积分公式应用 四阶龙格 库塔 Runge Kutta 法 4 用四阶龙格 库塔公式

      9、计算根据公式先计算出4个系数 再计算仿真第一步的值 数值积分公式应用 数值积分公式应用 则系统仿真第一步的值为 数值积分公式应用 从上述结果可以看出 对于同一个系统进行仿真计算时 其值的精度是随着数值积分公式的变化而改变的 其中欧拉法计算精度最低 其次为梯形法和二阶龙格 库塔法 四阶龙格 库塔法计算精度最高 数值积分公式应用 例2 用matlab演示 龙格 库塔法的误差估计 一个高精度的仿真方法必须将步长控制作为手段 实现步长控制涉及局部误差估计和步长控制策略两方面的问题 龙格 库塔法的误差估计 RK方法的误差估计通常是设法找一个低一阶的RK公式 将两个公式计算结果之差作为估计误差 例如Runge Kutta Fehlberg法的计算公式是RKF1 2公式用另一个一阶公式来估计误差 龙格 库塔法的误差估计 RKM3 4公式 误差估计式用3阶 计算公式为4阶 RKM3 4公式误差估计公式 龙格 库塔法的步长控制 龙格 库塔法的误差估计和步长控制策略的基本思想是 每积分一步都设法估计出本步的积分误差en 然后判断是否满足允许误差E 据此选择相应的步长控制策略 每一步的局部误差通常取以下形式

      10、en En yn 1 其中 yn 是利用误差估计式计算出的本步的估计误差 当 yn 较大时 en是相对误差 当 yn 较小时 en是绝对误差 这样作的目的是避免当y的值很小时 en变得过大 仿真模型的运行速度与实际系统运行速度一致 称为实时仿真 一般方法难以满足实时仿真的需要 所得模型的执行速度较慢 机理也不符合实时仿真的需要 假设对一般形式的系统进行仿真 以RK 2为例进行分析 其公式为 实时仿真 假定在h 2的时间内计算机刚好计算一次右端函数f 则计算分为两步 1在tk时刻利用当前的un yn计算K1 2在tn h 2时刻计算K2 此时un 1无法得到 但实时仿真除了要满足执行速度的要求外 还要求实时接收外部输入 并实时得到输出 此种情况下 解决的方法有两个 对un 1进行预报 增大仿真误差 或仿真延迟半个计算步距 后者的计算流程如下可见 后种方法的输出也会延迟半个计算步距 为了克服这个缺陷 人们提出了如下形式的实时二阶RK法 实时RK 2公式 其计算流程 假定在h 2的时间内计算机也刚好计算一次右端函数f 则计算也分为两步 1在tn时刻利用当前的un yn计算K1 2在tn h

      《龙格-库塔方法基本原理.ppt》由会员自***分享,可在线阅读,更多相关《龙格-库塔方法基本原理.ppt》请在金锄头文库上搜索。

      点击阅读更多内容
    最新标签
    发车时刻表 长途客运 入党志愿书填写模板精品 庆祝建党101周年多体裁诗歌朗诵素材汇编10篇唯一微庆祝 智能家居系统本科论文 心得感悟 雁楠中学 20230513224122 2022 公安主题党日 部编版四年级第三单元综合性学习课件 机关事务中心2022年全面依法治区工作总结及来年工作安排 入党积极分子自我推荐 世界水日ppt 关于构建更高水平的全民健身公共服务体系的意见 空气单元分析 哈里德课件 2022年乡村振兴驻村工作计划 空气教材分析 五年级下册科学教材分析 退役军人事务局季度工作总结 集装箱房合同 2021年财务报表 2022年继续教育公需课 2022年公需课 2022年日历每月一张 名词性从句在写作中的应用 局域网技术与局域网组建 施工网格 薪资体系 运维实施方案 硫酸安全技术 柔韧训练 既有居住建筑节能改造技术规程 建筑工地疫情防控 大型工程技术风险 磷酸二氢钾 2022年小学三年级语文下册教学总结例文 少儿美术-小花 2022年环保倡议书模板六篇 2022年监理辞职报告精选 2022年畅想未来记叙文精品 企业信息化建设与管理课程实验指导书范本 草房子读后感-第1篇 小数乘整数教学PPT课件人教版五年级数学上册 2022年教师个人工作计划范本-工作计划 国学小名士经典诵读电视大赛观后感诵读经典传承美德 医疗质量管理制度 2 2022年小学体育教师学期工作总结 2022年家长会心得体会集合15篇
     
    收藏店铺
    关于金锄头网 - 版权申诉 - 免责声明 - 诚邀英才 - 联系我们
    手机版 | 川公网安备 51140202000112号 | 经营许可证(蜀ICP备13022795号)
    ©2008-2016 by Sichuan Goldhoe Inc. All Rights Reserved.