二阶龙格—库塔公式

第一节 常微分方程 第二节 欧拉方法 第三节 龙格—库塔法 在上一节中,我们得到了一些求微分方程近似解的数值方法,这些方法的局部截断误差较大,精度较低,我们希望得到有更高阶精度的方法。 一阶龙格—库塔方法 如果以y(x)在xi处的斜率作为y(x)在 [xi,xi+1]上的平均斜率k*,即 二阶龙格—库塔方法 在[xi,xi+1]上取两点xi,xi+p(0< p≤1)的斜率值k1,k2的线性组合λ1k1+λ2k2作为k*的近似值(λ1、λ2为待定常数),此公式一般形式可写成 这就是二阶龙格—库塔法公式。 三阶龙格—库塔公式 为了提高精度,考虑在[xi,xi+1]上取三点xi,xi+p,xi+q的斜率值分别为k1,k2,k3,将k1,k2,k3的线性组合作为平均斜率k*的近似值,其中 k* 这就是欧拉法. 则得 其中k1 = f (xi,yi),k2为[xi,xi+1]内任意一点xi+p = xi+ ph (0< p≤1) 的斜率f (xi+p,y(xi+p))。 由于y (xi+p)并没有给出,所以先应该求y (xi+p),仿照改进欧拉公式的构造思想,得到 (8-7) 这样构造出的公式为 k1 k2 k* 公式中含有三个参数λ1,λ2和p,如果我们适当选取参数的值,可以使公式的局部截断误差为O(h3)。 对k1和k2作泰勒展开 代入(8-7)得 (*) 又 y (x)在xi处的二阶泰勒展开式为 当x = xi+1时, ,有 (**) (**) 比较(*)与(**)的系数即可发现, 要使公式(7-7)的局部截断误差满足 ,即要求公式具有二阶 精度只要下列条件成立即可。 (8-8) 满足条件(8-8)的一簇公式统称为 二阶龙格—库塔公式。 特别的,当 塔公式就成为改进欧拉公式。 时,龙格-库   改进欧拉公式就是以y(x)在xi和xi+1 处的斜率k1和k2的算术平均 值作为y(x)在[xi,xi+1]上的 平均斜率k*来进行计算的。 若取 时,龙格-库塔 公式就称为变形的欧拉公式,其形式为 (8-9) 此处的 就是欧拉方法预报出的中点 处的近似解;而 等于中点的斜率值 则近似 ; 所以公式可以看作用中点斜率近似 代替平均斜率k*,因此,公式(8-9)也 称作中点公式。 xi+q=xi+qh (0

二阶龙格库塔公式推导_二阶龙格—库塔公式.PPT相关推荐

  1. 二阶龙格库塔公式推导_[数学]龙格-库塔法

    原理思想 要想求出非常近似的值,有种神器叫做泰勒公式 .泰勒给出了任意一个函数都可以用多项式逼近的方法求出函数值.这与常微分方程的数值方法的思想类似,就是已知初始值,借助导数这个工具,将其近似成求另一 ...

  2. 二阶龙格库塔公式推导_[常微分方程的数值解法系列五] 龙格-库塔(RK4)法

    在惯性导航以及VIO等实际问题中利用IMU求解位姿需要对IMU测量值进行积分得到需要的位置和姿态,其中主要就是求解微分方程.但之前求解微分方程的解析方法主要是应用于一些简单和特殊的微分方程求解中,对于 ...

  3. 二阶龙格库塔公式推导_连续系统数值仿真方法——龙格库塔法

    在决策理论与方法最后一课中,老师讲了连续系统建模与仿真,事实上在本科就已经接触过许多,在这里我把数值仿真方法中的龙格库塔法给详细介绍一下,这是个比较有趣,同时也比较实用的方法.在求解常微分方程初值问题 ...

  4. 二阶龙格库塔公式推导_数值常微分方程-欧拉法与龙格-库塔法

    大三时候在跳蚤市场闲逛,从一位数学院的学长那里买了一些闲书,最近翻出来刚好有李荣华.刘播老师的<微分方程数值解法>和王仁宏老师的<数值逼近>,结合周善贵老师的<计算物理& ...

  5. 二阶龙格库塔公式推导_二阶常系数齐次线性方程通解推导(涉及常数变易法和欧拉公式)...

    欧拉恒等式 二阶微分方程明显比一阶难了很多,下面三图详细地对二阶常系数齐次线性方程的通解进行了推导. 有几下几点需要注意: 1.理解思路. 求二阶常系数齐次线性方程的解,一开始是靠猜的,因为以e为底数 ...

  6. 二阶龙格库塔公式推导_带你走进最美数学公式

    同学们,我们先来跟老师欣赏一下数学中最优美的式子吧? 是什么魔力让以上几个似乎毫不相干的数学中最特殊的数字能如此优美的写在同一个式子呢? 是欧拉,是数学. 0和1--老师就不用介绍啦, e是自然常数( ...

  7. 带通 带阻滤波器 幅频响应_二阶有源带通滤波器设计

    二阶有源带通滤波器设计 1.背景 对于微弱的信号的处理方式一般是:放大和滤波,这个过程中就涉及到放大电路的选取.滤波器的选择以及偏置电路的设计.本例以实例的方式讲解并附带参数计算.仿真.实物测试三个环 ...

  8. matlab编辑二阶线性系统,基于MATLAB的二阶线性系统分析与仿真

    第26卷第5期 河池学院学报 Vol .26No .52006年10月 JOURNAL OF HECH IUN I V ERSI TY Oct .2006基于MAT LAB 的二阶线性系统分析与仿真 ...

  9. 数学公式推导_梯度_gradient

    数学公式推导_梯度_gradient 核心 θt+1=θt−αt∇f(θt)(1)\theta _{t+1}=\theta _t-\alpha _t\nabla f\left( \theta _t \ ...

最新文章

  1. Python元组常用的操作
  2. 价值6.11亿美元的入侵工具无人问津
  3. Linux账号和权限管理详解(超详细示例操作)!
  4. python 百度ai批量识别_Python基于百度AI的文字识别的示例
  5. [react] React什么是有状态组件?
  6. 最完美的xslt数值函数与字符串函数
  7. c语言必背数据结构_严蔚敏数据结构(C语言版)知识点总结笔记课后答案
  8. 联想 m73 黑苹果 软路由 esxi AIO
  9. 雷赛服务器信号er020,伺服与雷赛控制卡配套的小技巧
  10. STM32开发 | 移远4G-Cat.1模组EC200N-CN开发
  11. Windows10下安装MySQL5.0详细教程
  12. NOIP2010导弹拦截
  13. ibmr系列服务器怎么装架子,R440/R540/R640/R740 R820 R930 DELL服务器导轨 滑轨 支架 理线架...
  14. 这100 个网络基础知识,看完成半个网络高手
  15. ES6/ES7/ES8新特性汇总
  16. python png 背景透明_Python - 移除PNG透明图的alpha通道
  17. 听声变位测试软件,刺激战场:听声辩位其实有很大的学问,想了解的朋友请进来...
  18. JAVA的对象list参数的校验
  19. multimap 的使用例子
  20. 视觉伺服控制工具Visual Servoing Platform---VISP(6)----基于4个平面点的姿态估计

热门文章

  1. 写给仿真软件研发的“一篇文章入门”系列(终)
  2. JAVA开发环境配置指南
  3. 嵌入式的选择------
  4. 监控物联网项目:web、flv视频流 、websocket、ESP32
  5. Prometheus 从入门到入土 -----入土
  6. 百度网盘下载失败【1252017】误报违规
  7. 记用Linux半年小结
  8. .bat输出目录下所有文件的绝对路径
  9. 大华城市安防监控系统平台管理存在任意文件下载漏洞
  10. 开源代理软件GAppProxy