常微分方程(ODE)求解方法总结(续)

news2024/10/7 2:22:28

常微分方程(ODE)求解方法总结(续)

  • 1 隐式方法
    • 2 多步法
    • 2.1 二阶方法
    • 2.1.1 非自启动修恩方法
    • 2.2 开型和闭型积分公式
    • 2.3 高阶多步法

1 隐式方法

常微分方程(ODE)求解方法总结 里面介绍了我称为“正常思路”的方法,就是 y i + 1 = y i − a y i ∗ h y_{i+1}=y_i - ay_{i}*h yi+1=yiayih,从 y i y_i yi推出 y i + 1 y_{i+1} yi+1这种思路。这种思路的方法称之为显式方法,也可以称为前向方法。// (原来形式为 y i + 1 = y i + f ( x i , y i ) ∗ h y_{i+1}=y_i+ f(x_i, y_i)*h yi+1=yi+f(xi,yi)h,为了与下面的隐式对比方便,写成 y i + 1 = y i − a y i ∗ h y_{i+1}=y_i - ay_{i}*h yi+1=yiayih)。

还有另外一种思路,就是, y i + 1 = y i − a y i + 1 ∗ h y_{i+1}=y_i - ay_{i+1}*h yi+1=yiayi+1h这种思路,由此式可以得出:
y i + 1 = y i / ( 1 + a ) h y_{i+1}=y_i /(1+a)h yi+1=yi/(1+a)h (1.1)
这种思路方法称为隐式方法,或者后向方法。

这种方法有很多应用,因为其有一个优点就是,这种方法不会因为步长选取不当而出现不稳定的问题,因此,该方法称作无条件稳定。( 这个先不说具体原因,后面会讲到。)
这种方法的特点就是,将 y i + 1 y_{i+1} yi+1先作为已知的点,来参与计算。

2 多步法

前一篇文章和这一篇文章的上一章节里面主要都是讲了单步法,就是从 y i y_i yi推出 y i + 1 y_{i+1} yi+1,过程中只需要一个知道前一个点 y i y_i yi

常微分方程(ODE)求解方法总结 里面“第2.1.5 欧拉方法的改进思路3”里面有讲到,能不能用取三个点、四个点等用二次曲线和三次曲线来插值逼近使积分结果更加接近真实值。
接下来,就会沿着这个思路前进。

2.1 二阶方法

在讨论高阶形式之前,先给出一个简单的二阶方法,便于理解。
在这里插入图片描述
单步法是通过 x i x_i xi点的值估计 x i + 1 x_{i+1} xi+1处的值。二阶方法,就是用 x i − 2 x_{i-2} xi2 x i − 1 x_{i-1} xi1 x i x_{i} xi点的值通过二次插值来估计 x i + 1 x_{i+1} xi+1处的值。

2.1.1 非自启动修恩方法

2.2 开型和闭型积分公式

2.3 高阶多步法

(未完待续。。。)

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.coloradmin.cn/o/717975.html

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈,一经查实,立即删除!

相关文章

Jira UI Locations及注意事项总结

issue view ui locations : https://developer.atlassian.com/server/jira/platform/issue-view-ui-locations/#issue-operations-bar-locations1.问题操作栏Issue Operations Bar Locations模块分为两部分: opsbar-operationsflopsbar-transitions两个location.共同定义了问题…

【力扣】111、二叉树的最小深度

111、二叉树的最小深度 给定一个二叉树,找出其最小深度。 最小深度是从根节点到最近叶子节点的最短路径上的节点数量。 说明:叶子节点是指没有子节点的节点。 // var minDepth function(root){if(!root) return 0;const stack [ [ root ,1] ];//…

零拷贝详解

目录 一、什么是零拷贝 二、传统的IO执行流程 三、零拷贝相关的知识点回顾 1、内核空间&用户空间 2、用户态&内核态 3、上下文切换 4、虚拟内存 5、DMA技术 四、零拷贝实现的几种方式 1、mmapwrite实现的零拷贝 2、sendfile实现的零拷贝 3、sendfileDMA sc…

MySQL原理探索——23 MySQL是怎么保证数据不丢的

今天这篇文章,我会继续介绍在业务高峰期临时提升性能的方法。从文章标题“MySQL 是怎么保证数据不丢的”,你就可以看出来,今天我介绍的方法,跟数据的可靠性有关。 在前面文章,我都着重介绍了 WAL 机制(你可…

ElementUI plus框架Table表格cell-style属性的使用

官方文档说明: 例:设置单元格文字居中 Object方式: function方式:

安全 --- http报文包详解及burp简单使用

HTTP HTTP(超文本传输协议)是今天所有web应用程序使用的通信协议。最初HTTP只是一个为了获取基本文本的静态资源而开发的简单协议,后来对其进行扩展和利用,使其发展为能够支持如今常见的复杂分布式应用程序。 (1&…

PADS-LAYOUT菜单及工具使用

目录 1菜单栏 1.1文件菜单 1.2编辑菜单 1.3查看菜单 1.4设置菜单 1.5工具菜单 1.6帮助菜单 2工具栏 2.1标准工具栏 2.2绘图工具栏 2.3设计工具栏 2.4尺寸标注工具栏 2.5ECO工具栏 3系统配置 3.1全局选项 3.2设计选项 3.3栅格和捕获选项 3.4显示选项 3.5布线选…

【UnityDOTS 六】预制实例化成Entity

预制实例化成Entity 前言 从Authoring模式中,如何通过预制件来实例化一个对应的Entity对象到DOTS系统中。 一、Authoring模式与Runtime模式 Authoring创作模式:即我们熟悉的方便操作的创建预制的模式 Runtime模式:运行模式,即在…

Three.js教程:网格模型

推荐:将 NSDT场景编辑器 加到你的3D工具链 工具集: NSDT简石数字孪生 网格模型(三角形概念) 本节课给大家演示网格模型Mesh渲染自定义几何体BufferGeometry的顶点坐标,通过这样一个例子帮助大家建立**三角形(面)**的概念 三角形(面) 网格模型Mesh其实…

Spring Boot 中的 XA 事务

Spring Boot 中的 XA 事务 在现代化的应用程序开发中,事务管理是一个重要的话题。事务管理可以确保数据的一致性和完整性,同时也可以避免数据丢失和冲突等问题。在分布式环境中,XA 事务是一种常用的事务管理方式。在本文中,我们将…

基于麦克风阵列模块I2s6路slot数字音频信号的ADSP/STM32F4处理

hezkz17进数字音频系统研究开发交流答疑 1麦克风阵列 2 ADAU1452 DSP 输入接口 3 PCM数据算法处理

uni-app获取节点的相关信息

获取单个节点: selectorQuery.select(selector) 在当前页面下选择第一个匹配选择器 selector 的节点,返回一个 NodesRef 对象实例,可以用于获取节点信息。 selector 说明: selector 类似于 CSS 的选择器,但仅支持下列…

Jenkins 参数化构建

目录 前言: 以下是Jenkins官方的介绍: 首先新建一个自由风格的job,然后勾选“This project is parameterized”,如下图: 勾选参数构建 choice parameterized: 选项: 使用参数:…

linux驱动struct file_operations`中设置`owner`成员变量。

在Linux内核中,struct file_operations结构体是用来定义一个文件操作的集合。该结构体中包含了一组函数指针,这些函数指针指向了与文件操作相关的函数,例如文件的打开、关闭、读写、定位等操作。同时,该结构体中还包含了一些元数据…

Opencv项目实战:00 专栏内容介绍

目录 Opencv项目实战专栏介绍: 01 文字检测OCR 02 角度探测器 03 扫描二维码&条形码 04 全景图片拼接 05 物体检测 06 文档扫描仪 07 人脸识别和考勤系统 08 Yolov3更高精度的检测物体 09 物体尺寸测量 10 面部特征提取及添加滤镜 11 使用Opencv高亮…

接近完美的监控系统—普罗米修斯

内容 隐藏 1.怎么采集监控数据? 2.采集的数据结构与指标类型 2.1 数据结构 2.2 指标类型 3.实例概念 4.数据可视化 5.应用前景 普罗米修斯(Prometheus)是一个SoundCloud公司开源的监控和告警系统。当年,由于SoundCloud公司生产了太多的服务&#…

ant vue 自定义首页菜单实现

如图,我把 views/dataScreening/DataScreening 页面,设置成首页,登录以后自动跳转到这个页面 实现代码: 文件:/src/config/router.config.js import { UserLayout, TabLayout, RouteView, BlankLayout, PageView } from @/components/layouts/*** 走菜单,走权限控制* @t…

【深度学习】AIGC ,ControlNet 论文,原理,训练,部署,实战,教程

论文:https://arxiv.53yu.com/pdf/2302.05543 代码:https://github.com/lllyasviel/ControlNet 得分几个博客完成这个事情的记录了,此篇是第一篇,摘录了一些论文内容。ControlNet 的原理极为朴实无华(对每个block添加…

c++ connect函数连接失败

bool Connect() {//初始化网络WSADATA wsadata;WSAStartup(MAKEWORD(2, 2), &wsadata);if (0 ! WSAStartup(MAKEWORD(2, 2), &wsadata))return false;if (LOBYTE(wsadata.wVersion) ! 2 ||HIBYTE(wsadata.wVersion) ! 2) {printf("请求协议版本失败!\n");ret…

Ghatgpt正式登录苹果手机应用商城,并支持Siri和快捷指令

根据最新信息,OpenAI 发布的 ChatGPT 官方 iOS 应用程序迎来了重大更新。该应用程序已经在上个月登陆了美国、英国、法国、德国和韩国等 App Store,并且成为该市场上最受欢迎的免费应用程序。 作为生产力类应用的领导者,该应用程序完全免费&a…