ODE(常微分方程,与偏微分方程相反)是包含一个或多个函数的微分方程。
我正在尝试使用前向欧拉方法解决 pyomo 框架内的一首颂歌。我明确地写下了离散化方案,但结果并不符合预期。从 pyomo 模型
我正在尝试获取 SIR 模型的以下 ODE 系统的稳态点(稳态点是系统方程为零时的解): 将 sympy 导入为 sp S,...
通过使用 optim() 最小化残差平方和来优化参数,将模型拟合到观测数据
我正在尝试将这个非常简单的 4 种线性 Lotka-Volterra 竞争模型拟合到观察到的数据,但由于某种原因,当我尝试 optim() 函数时,与 deSolve 相关的东西似乎......
我正在尝试用 Javascript 对弹簧上的质量进行 Runge-Kutta 实现,并使用 D3 对其进行可视化。目的是将其与前向欧拉进行比较并评论差异。我的有限元工作...
尝试按照此 PDF 重现 4.1.3 固体材料的排放。 https://www.rivm.nl/bibliotheek/rapporten/2017-0197.pdf 制作了这个版本的模型,这是我得到的最好的......
有没有办法在 MATLAB 中输出 ode 求解器的中间计算?
我有一个刚性 ODE 系统,我用 ode15 求解。时间导数的类型为: yt = A*y - u.*y; 其中 yt、y 和 u 是列向量,A 是稀疏矩阵。 u ha 的计算...
使用 Python symfit ODEModel 拟合多维数据
我正在尝试将 ODE 的参数拟合到二维数据,根据拟合多维数据集的示例,这通常应该是可能的。 这是我失败的尝试所以...
如何在 SymPy 中按预期求解这个微分方程 (1-x^2)y'=x^2-xy-1?
我想用 Sympy 求解 (1-x^2)y'=x^2-xy-1 ,但无法按预期求解。 我期望得到以下解决方案。 -11:y(x)...
Fortran 95 中求解三阶 ODE 的 Runge-Kutta-Fehlberg 方法
我需要帮助构建 Fortran 95 代码,使用 Runge-Kutta-Fehlberg 方法求解三阶 ODE。我试图用 y(0)=3、y'(0)=-2 和 y''(0)=6 求解 y''' = -2y''+y'+2y。确切的解决方案...
获取数据代码: 将 pandas 导入为 pd 进口火炬 数据集 = pd.read_csv('/kaggle/input/fish-bear/population_data.csv') 年= torch.tensor(数据集['年'],dtype = torch.float64) 鱼_p...
我正在尝试解决涉及使用常微分方程组的固定床反应器的问题,但由于某种原因它陷入了无限循环。这是代码: 导入数量...
我正在尝试使用有限差分法作为使用solve_ivp 方法的图像滤波技术来求解二维热方程。显示我的代码,给出错误消息“y0 ...
我对 R 相当陌生,并且拥有这个 ODE 系统 rac{dS}{dT} = - eta(1-\mu)S(t) rac{I(t))}{1-D(t))} rac{dI}{dT} = - eta(1-\mu)S(t) rac{I(t))}{1-D(t))}-\delta I(t))-\phi I (t)) rac{dR}{...
........................................................ ...................................................... ................................................
Scipy_integrate 和 sympy 用于求解 BVP
我有一个小问题,即我有 8 个差值(在通过壁进行传质的逆流热交换器中),但只能调用每个单独的差值。问题是...
有一个基本的 ODE - 微分方程。 y' = 2x ,解为 f(x) = x² 代码: def deriv (y,x) : 返回 2 * x # 时间点 x = np.linspace ( 0 , 10) # 初始条件 y0 = 1 # 溶胶...
我想对 ODE 模型进行参数推断(m.options.IMODE = 5)。我在 Gekko 中有 2 个模型。首先,我从每个模型生成模型数据,并直接在
我有一个与这篇文章非常相似的问题。不幸的是,这个问题没有答案,我还有其他问题。 我正在尝试求解 RLC 电路中的电流...
我正在尝试编辑 odeint 文档中的 Gray-Scott 一维方程示例(页面上的最后一个示例)。 我有下面的代码,它适用于诺伊曼边界条件,但我想要一个
为什么我的基本 Gekko ODE 求解器比 Scipy 慢很多?
在这个最小的例子中,我想求解基本积分器 ODE dT/dt = k[F(t) - T(t)],其中 F(t) 是一个力向量,被选择为方波。我已经实施了两个程序来...