ODE(常微分方程,与偏微分方程相反)是包含一个或多个函数的微分方程。
如何在matlab中使用interp1函数对symfun类型的函数?
我在matlab脚本中有一个微分方程组的解决方案。我将我的函数定义如下:syms V_x(t)V_y(t)X(t)Y(t); ode_V_x = diff(V_x,t)== - (B / m)* V_x; ode_V_y = diff(V_y,t)= ...
试图使用solve_ivp的矢量化选项,奇怪的是它会抛出y0必须是1维的错误。 MWE:来自scipy.integrate import solve_ivp import numpy as np import math def f(t,...
我有以下简单的ODE:dx / dt = -1初始条件x(0)= 5,我感兴趣的是x(t)== 1。所以我有以下事件函数:function [value,isterminal,direction] = ...
我正在尝试使用scipy中的odeint解决二阶ODE问题。我遇到的问题是函数隐含地耦合到二阶项,如简化片段中所示(请忽略...
我正在尝试使用python解决Brusselator模型,一个非线性ODE。我以前用MATLAB做这个,但现在正在用python作为后端构建一个应用程序。这就是为什么我要切换ti ...
我试图在距离矩阵上计算一个普通的ODE(常微分方程),但我不知道如何并行化我的代码。来自scipy.integrate import math from math import exp import ...
我希望在时间间隔I = [0,10]的时间内用python求解方程式,初始条件(x_0,y_0)=(1,0),参数值μ∈{-2,-1,0,1, 2}使用函数scipy.integrate ....
代码t0 = 0; tfinal = 1; tspan = [t0 tfinal]; optionsII = odeset('Events',@ myEventsFcnII); y0 = [-0.5972 6.2068]; yp0 = 10000 * [0.0006 1.0178]; [t,y,te,ye,ie] = ode15i(@(t,y,yp)StateII(t,y,yp),...
我在使用我的ode求解器时遇到了一些问题,我正在尝试解决SEIR问题而且我一直在得到相同的错误,这些错误会根据我的代码非常相似地调整代码。我的代码是:...
R simecol和稳态:rootfun无法访问SimObj方程槽
我使用simecol包框架模拟生态模型(作为ODE系统),通过对象类SimObj(见此处)非常容易地使用和共享生态模型。我想要 ...
我试图使用deSolve拟合一个合理的基本ODE模型,包括在模型中一个随时间变化的参数(感染力; FOI)。在没有此参数的情况下运行模型时...
我正在尝试使用Python的solve_ivp来解决ODE。但是,我想基于当前解决方案与先前解决方案之间的比较动态地更改ODE的右侧。 ......
有没有办法强制ode()[deSolve-R包]在ode函数的每个集成步骤中提供输出
我想在数值方案的每一步提取状态变量的值。 deSolve R包中的ode()函数使用一个已实现的ODE求解器来数值求解......
我正在使用包deSolve构建基于R中Lotka-Volterra导数的捕食者 - 食饵模型。我定义了参数,初始状态和时间步长以及模型函数。然后我用...来解决每一个问题
我正在尝试计算大型一阶ODE系统的雅可比行列式。我需要Jacobian才能在我的solve_ivp代码中使用。这是我如何尝试计算雅可比。 C1 = sp.symbols('C [0],C [1],...
我正在尝试为有限差分法设计算法,但我有点困惑。所讨论的ODE是y'' - 5y'+ 10y = 10x,y(0)= 0且y(1)= 100。所以我需要一种方法以某种方式获得...
我正在使用刚性求解器(ode15s)对ODE系统进行时间积分。它工作正常,但我想加快速度。方程组以状态空间形式给出:函数[dx] = ...
我想在python中以数字方式求解ODE,如y'= f(x)(边界条件y(0)= 0)。关键是我不知道这个函数f(x)的分析性表达是什么,相反我有一个集...
我正在对一台高架起重机进行建模并得到以下公式:当谈到Scilab时,我就是noob,到目前为止,我只使用不超过两个自由度的线性系统(使用ODE)进行模拟,...