ode 相关问题

ODE(常微分方程,与偏微分方程相反)是包含一个或多个函数的微分方程。

四阶非线性微分方程的稳定解

我已经在python中使用bvp求解器解决了以下bvp问题。 d4y / dx4 = 0.00033 * V /(0.000001-y)^(2),y(0)= y'(0)= y(1)= y'(1)= 0在上面的等式'V'是一个使用...

回答 1 投票 0

初始条件非零的Python3 ode求解器

我正在尝试模拟具有单个自由度的质量弹簧系统。对于时间积分,我使用scipy中的ode函数。另外,我正在将数值与解析式进行比较...

回答 1 投票 0

如何在MATLAB的ode15s中为选项结构编写自定义“输出函数”以存储解矢量的历史记录?

对于我的代码选项= odeset('OutputFcn',@ myfun)[T,Y] = ode15s(@(T,X)sys(T,X,vin),[t0 .0005],X(:,1 ),options)函数Xdot = sys(T,X,vin)Xdot = A * X + B * U结束我不知道该怎么写...

回答 1 投票 0

我的odeint实现是否存在逻辑错误?

我得到了错误的解决方案,不确定odeint是否是解决此ODE系统的正确工具。我正在尝试通过求解ODE系统来对简单的一阶化学反应进行建模。 ...

回答 1 投票 0

Scipy中的matlab“ OutputFcn”的等效项?

类似于MATLAB ODE求解器具有Outputfcn, 其中ODE求解器在每个成功的时间步长后都会调用函数options = odeset('OutputFcn',@ odeprint)[T,Y] = ode15s(@(T,X)sys(T,X ,vin),[t0 .0005],X(...

回答 1 投票 0

为隐式Runge-Kutta方法编写函数(第四阶)

我正在尝试使用第4阶隐式Runge-Kutta方法(IRK)编写一个能够解决ODES系统的函数,但是我无法正确定义循环。在这里,我们通过...

回答 1 投票 2

改进矢量化的Julia代码或如何对它们进行矢量化处理?

我正在尝试将一些python代码转换为julia。它们旨在解决一些微分方程。由于它们包含方程式,以提高可读性,因此我给出了它们的图像,如下所示。我是...

回答 1 投票 0


使用Matlab使用欧拉方法求解ODE系统

我创建了一个函数Euler.m,以使用Euler方法求解ODE系统。我希望使用此函数来解决由匿名函数func = @(t)([x(t)+ 4 * y(t)-...

回答 1 投票 1

使用scipy ode求解器时的尺寸值误差

当前正在尝试从scipy.interpolate从scipy.interpolate导入interp1d到np来运行以下代码import numpy。 >

回答 1 投票 0

GNU Octave(Matlab)COVID-19 SIR-X模型

我正在尝试执行与本文所述相同的操作:https://science.sciencemag.org/content/sci/early/2020/04/07/science.abb4557.full.pdf我使用GNU Octave。这是带有ODE的功能文件...

回答 1 投票 0

使用fsolve和Euler 3BDF的SIR模型

嗨,我被要求使用MATLAB中的fsolve命令来求解SIR模型,而Euler 3向后指向。我对如何继续感到很困惑,请帮忙。这是我到目前为止所拥有的。我创建了一个函数...

回答 1 投票 1

是否可以在Python的ODEINT中更改我的ODE?

我正在尝试模拟一个具有微小差异的振荡系统:当身体向上运动时,我希望它使用特定的运动方程(ODE),而当身体向上时,我希望它使用另一个运动方程。]]] >

回答 1 投票 0

从odeModel中提取结果?

我没有尝试获得df对象作为我的输出(因此我可以对每个图进行子集/自定义)。我还读到可以使用函数“ out”提取结果,...

回答 1 投票 1

矢量化/并行化GEKKO方程

使用Phython的GEKKO库,并按照该示例从其原始论文中求解波动方程,需要创建大量方程。 m.Equations([[u [i] .dt()== v [i] for i in ...

回答 1 投票 1

GEKKO优化陷入困境(APM.exe)

我正在尝试使用Python 3(通过Jupyter笔记本)使用GEKKO运行微分方程系统求解器。对于错误的初始参数,它当然会立即停止并说未找到解决方案。或...

回答 1 投票 1

在GNU Octave中实施Euler方法

我正在阅读Chapra和Canale撰写的“工程师的数值方法”。在其中,他们为实现Euler方法(用于求解常微分方程)提供了伪代码。这是...

回答 2 投票 0

解决来自sympy的二阶微分方程组

我正在使用二阶拉格朗日方程组进行多自由度动力学问题。我用sympy来了解运动方程。现在这些方程在计算导数后变得相当长,...

回答 1 投票 0

Scipy odeint-每个时间步tx中参数的不同值

我正在尝试使用SciPy模拟一些动态模型。我有模型定义:def model(y,t,control_signal):dy / dt = some_function_of_time_and_y return dy我定义了时间戳列表...

回答 1 投票 0

您如何在deSolve(R)中实现系统动力学样式的管线延迟?

我正在尝试使用R中的deSolve对流水线延迟进行建模。我有一个库存(worktodo)具有恒定的输入(work_arrival),并且我希望流水线延迟执行(work_rate)在库存中...

回答 1 投票 0

© www.soinside.com 2019 - 2024. All rights reserved.