matlab刚性与非刚性情境下使用的ode求解器使用指南

最新推荐文章于 2025-03-23 14:02:08 发布
最新推荐文章于 2025-03-23 14:02:08 发布
阅读量5k 收藏 40
点赞数 10
分类专栏： 数学建模 文章标签： matlab
版权声明：本文为博主原创文章，遵循 CC 4.0 BY 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/m0_37149062/article/details/106853333
版权

解决常微分数值解的情况分为刚性问题与非刚性问题。问题不同使用的函数不同，使用错了，迭代时间久了，反而效果不是很明显。

非刚性问题

在这里插入图片描述

刚性问题定义及求解器使用

所谓刚性问题，其数值解变化分量特别诡异，有时很慢，有时很快；变化快的分量很快趋于稳定值，而变化慢的分量缓慢趋于稳定值，变化快时，我们想要积分小步长，慢的趋于稳定的，我们想要步长放大。但是事与愿违，因此称这种在步长选择上有纠结的为刚性问题。在这里插入图片描述

汇总版本ode使用指南

求解器方法描述使用场合
ode232到3阶Runge-Kutta算法，低精度非刚性
ode454-5阶Runge-Kutta算法，中精度非刚性
ode113Adams算法，精度可到10的-3次到10^-4次非刚性，计算时间比ode45快
ode23t梯形算法适度刚性
ode15s反向数值微分算法，中精度刚性
ode23s2阶Rosebrock算法，低精度刚性，当精度较低时，计算时间比ode15s短
ode23tb梯形算法，低精度刚性，当精度较低时，计算时间比ode15s短
ode15i可变秩求法完全隐形微分方程

在这里插入图片描述

Matlab：解算刚性 ODE
希望我的博客，能帮上你解决学习中工作中所遇到的问题
02-23 478
Matlab：解算刚性 ODE
刚性 ODE 简介：比较刚性和非刚性系统的行为-matlab开发
05-29
如果绝对稳定性要求比精度要求严格得多，则 ODE 是刚性的，我们在选择 ODE 求解器时需要小心。 这个实时脚本从一个简单的（非僵硬的）示例开始，然后将其与一些僵硬的示例进行比较。
Matlab 刚性问题求解器-ode23s
qingfengxd1的博客
07-21 2360
最后需要注意的是，ode23s 求解器对于一些特定类型的刚性问题可能表现出更好的性能，但对于其他类型的刚性问题可能不如其他求解器（例如 ode15s 或 ode23t）高效。因此，在实际使用中，根据具体问题的性质选择合适的求解器是很重要的。在输出中，te 是事件的时间，ye 是事件发生时的解，ie 是触发的事件的索引。该求解器将返回在给定时间范围内求得的微分方程的解。当使用 MATLAB 的 ode23s 求解器来求解刚性问题时，需要提供一个函数句柄来表示微分方程，并设置初始条件和求解的时间范围。
ODE求解器与欧拉法
最新发布
m0_61544805的博客
03-23 178
数值求解形如 dxdt=f(x,t)dtdx​=f(x,t) 的常微分方程。：将连续时间离散化为步长 ΔtΔt，通过迭代更新近似解。：简单高效，但误差较大（一阶精度）。仅使用当前点的导数进行一步近似。：欧拉法、龙格-库塔法等。
matlab刚性和非刚性的区别,刚性微分方程与非刚性的大概区分（自看，求指点）...
weixin_42470196的博客
03-27 1447
(这里的内容是作为自己的一个粗略的总结，不确定是否正确，希望有大佬能够更明确的指出其中的错误，作出指导)对于刚性和非刚性微分方程的区分，可以简单的转变为在将原方程转换为常微分方程组后，进行一个简单的系数判断：例如：y''' - 3y'' - y'y = 0, y(0) = 0, y'(0) = 0, y''(0) = -1;在这里可以设 y1 = y, y2 = y', y3 = y'', 有y...
关于求解微分方程——初学Matlab里的 ODE求解器
热门推荐
weixin_47723114的博客
03-15 2万+
这篇文章简要总结了一些关于ode求解器的一些知识及运用。然后在例子里面提到了一个好玩的东西——混沌理论。大家看个热闹就行hhh。不过还是希望有相关学习研究的友友能够与我多多交流学习哇！
ODE 求解器：实现几种流行的求解器，用于在 MATLAB 中求解 ODE-matlab开发
05-29
ODE 求解器的集合，形式为： y_dot = f(t,y) 然后求解器以以下形式提供此类 ODE 的解y = f(t,y) 到目前为止已实现的求解器： - 欧拉- 四阶龙格-库塔- Runge-Kutta 3/8 - Dormand-Prince - Runge-Kutta-Fehlberg (RKF45) 请检查 /docs 文件夹中的文档。
掌握MATLAB解决一阶微分方程的RK4方法
标签中所提到的'1storder ODE45 OnTheRun rungekutta matlab'意味着这些代码是与Matlab环境密切相关的，特别是针对第一阶微分方程的求解。ODE45是Matlab中用于求解常微分方程初值问题的一个函数，它基于变步长的...
MATLAB仿真与数学建模：理论与实践的结合，深入剖析！
首先概述了MATLAB仿真与数学建模的基本概念和工具箱功能，接着深入探讨了MATLAB在不同类型数学模型（离散数学模型、连续数学模型、统计模型）中的具体应用和分析方法。第四章着重讲述了高级数值计算方法、仿真模型的...
代码生成与集成：Simulink与CANape最佳实践指南
![代码生成与集成：Simulink与CANape最佳实践指南]... 参考资源链接：[CANape中Matlab Simulink模型的集成与应用](https://wenku.csdn.net/doc/6465
多体系统动力学案例研究：理论应用的深度解读
本文首先介绍了多体系统动力学的基本概念、理论基础以及数学模型，包括牛顿力学与拉格朗日力学的基本定律和原理、多体系统的稳定性和响应分析。接着，详细探讨了多体系统动力学仿真工具的选择、动力学模型建立及仿真...
谐响应案例分析：5步解决复杂工程难题
谐响应分析作为结构动力学中的一个重要分析方法，主要应用于工程领域中，用以评估结构在连续周期性载荷下的稳态响应。本文第一章介绍了谐响应分析的基础理论，包括其数学模型和相关计算方法。第二章探讨了谐响应分析...
matlab中ode指令,在Matlab中使用ODE选择步长
weixin_34870207的博客
03-16 2525
嘿，谢谢您抽出宝贵时间看我的问题。这是我先前在physics.stackexchange.com中发布的问题的更新版本我目前正在研究2D激子旋量玻色-爱因斯坦凝聚物，并对这个系统的基态感到好奇。达到基态的数学方法称为虚时法。其中在量子力学时间由假想一个替换方法是非常简单的这种取代导致在我的系统的高能量粒子衰减比低能量更快。在计算的每个步骤中重新归一化粒子数，我们最终得到了一个最低能量的粒子系统。基...
matlab使用教程(24)—常微分方程(ODE)求解器
配电网和matlab的博客
08-25 7015
′ 的某些分量缺失，则这些方程称为微分代数方程或 DAE，并且 DAE 方程组会包含一些代数变量。术语“刚度”无法精确定义，但一般而言，当问题的某个位置存在标度差异时，就会出现刚度。例如，如果 ODE 包含的两个解分量在时间标度上差异极大，则该方程可能是刚性方程。此表包含可用的 ODE 和 DAE 示例文件及其使用的求解器和选项的列表。您可以指定需要解算的任意数量的 ODE 耦合方程，原则上，方程的数量仅受计算机可用内存的限制。为解算该方程，需要将实部和虚部分解为不同的解分量，最后重新组合相应的结果。
数值计算project——ODE求解器
保持写作习惯，完成知识沉淀
06-15 301
一、介绍 这篇文章介绍一个ODE solver的实现，可用于一阶线性系统的求解。 首先，找到一组ODE方程，可以从网上找算例，或者从comsol等软件中导出。 其次，用matlab来求解，看计算结果。 最后，自己实现一个ODE solver，并将结果与matlab的结果进行对比。 ...
选择 MATLAB 的 ODE 求解器
2301_79325339的博客
09-11 558
根据问题的特点和求解需求，MATLAB 提供了几种不同的求解器。MATLAB 提供的 ODE 求解器有两类：显式求解器和隐式求解器。显式求解器适用于问题具有较简单结构和较小的维度的情况，而隐式求解器适用于问题具有更复杂结构和较高维度的情况。如果问题的维度较高，或者需要进行大量的求解器迭代，那么选择一个计算效率较高的求解器可以显著提高求解速度。根据具体问题的特点，可以选择合适的求解器进行求解。根据具体问题的特点，选择适合的求解器，并根据需要调整求解器的参数，可以获得准确且高效的求解结果。，以控制求解的精度。
ode求解器的事件(Event)属性
daijiangtao的专栏
07-07 1万+
检测事件 matlab微分方程如何设定变量的范围  如 dy1=y2 dy2=y1+1 其中y1的范围为0 【解】 M文件： function [value,isterminal,direction] = events1(t,y) value = y(1)-4; isterminal= 1; direction = 0; 命令窗口： dy = @(t,
matlab求解器是什么,Matlab理解ode求解器
weixin_30877741的博客
03-19 2047
我有一个链接微分方程组，我用ode23求解器求解 . 当达到某个阈值时，其中一个参数会改变，这会改变我的函数的斜率 .我使用调试功能跟踪了ode的行为，并注意到它在这个点附近开始跳回“时间” . 基本上它会生成更多的数据点 . 但是，这些并不都在最终的解决方案向量中表示 .有人可以解释这种行为，特别是为什么不是所有计算值都能找到解决方案向量的方法？//编辑：为了澄清，当v从0更改为任何其他值时，行...
刚性微分方程与非刚性的大概区分（自看，求指点）
weixin_34138255的博客
08-05 1万+
（这里的内容是作为自己的一个粗略的总结，不确定是否正确，希望有大佬能够更明确的指出其中的错误，作出指导） 对于刚性和非刚性微分方程的区分，可以简单的转变为在将原方程转换为常微分方程组后，进行一个简单的系数判断： 例如： y''' - 3y'' - y'y = 0, y(0) = 0, y'(0) = 0, y''(0) = -1; 在这里可以设 y1 = y, y2 = y', y3 =...
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值