ODE(常微分方程,与偏微分方程相反)是包含一个或多个函数的微分方程。
我正在尝试解决耦合的ODES,但我在这里提出了一小部分问题。我试图沿整个网格求解R。网格大小太大而无法运行常规方法,因此我试图在Julia中并行处理,但是我仍然无法在Mac环境中运行代码。我将M1与8GB RAM一起使用,并尝试了多种方法来解决该问题。如果我做错了什么,请让我知道,如果我做错了什么?
我正在 R 中对传染病的传播进行建模,我需要通过多组参数循环模型。到目前为止,我有一个工作模型(下面是一个虚拟版本),但只是为了罪恶......
在Python中使用solve_ivp()求解耦合微分方程时,我的代码中存在无法解释的增益
我正在尝试使用 scipy.integrate 中的solve:ivp() 函数求解python 中的两个耦合微分方程,但我在其中一项项上得到了无法解释的增益,从而破坏了方程。它我...
我编写了以下ODE文件: 函数 [dndt,dmdt,dhdt,dVdt]=hh(n,m,h,V) C=1; GK=36; GNa=120; GL=0.3; EK=-72; ENa=55; EL=-49.4; dndt=(1-n)*alphan(V)-n*betan(V); dmdt=(1-m)*alpham(V...
让我们考虑以下代码: 将 sympy 导入为 sy 将 scipy 导入为 sc 从 sympy.physical.units.quantities 导入数量 从 sympy.solvers.ode.systems 导入 dsolve_system 来自 IPython.display
我正在尝试求解弹性基础上梁的四阶 ODE。我将这个问题应用于嵌入多个土层中的横向加载桩,这些土层可以在...
我正在使用solve_ivp,但根据问题的设置,我得到了奇怪的结果,尽管我认为这更多是求解器的实现问题而不是编码问题,但我想...
在 NDSolve 中使用 NIntegrate 函数的问题
我试图将使用 NDSolveValue 求解的偏微分方程生成的插值函数传递给对 NDSolveValue 的第二次调用,以获得简单 ODE 的解。这是我的代码
我是编码新手,正在尝试求解简单的一维色散方程。 方程和边界条件: adC/dx = bd^2C/dx^2 x = hj,C = C0; x = - inf,C = 0 解析解为 C = C0 *...
我这样做是为了使用射击方法解决BVP功能。 我不明白为什么当我的步长 h = 0.1 时,限制 x_end 超出了我在函数solve下的 while 循环中设置的值...
我正在尝试为由方程组控制的优化系统找到多个参数。 在我的代码中,我需要能够使用 4 或 5 个方程并确定
我正在使用这个优化脚本: 使用微分方程、优化、优化多元算法、OptimizationOptimJL、绘图、SciMLSensitivity 参数=\[\] 函数lotka_vol...
我正在尝试模拟具有两个入口流和一个出口的水箱的液位。这个想法是,之后将用于控制问题。然而,我无法让它与 gekko 一起工作,但它......
我正在尝试使用前向欧拉方法解决 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 的参数拟合到二维数据,根据拟合多维数据集的示例,这通常应该是可能的。 这是我失败的尝试所以...