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 的参数拟合到二维数据,根据拟合多维数据集的示例,这通常应该是可能的。 这是我失败的尝试所以...
如何在 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 个差值(在通过壁进行传质的逆流热交换器中),但只能调用每个单独的差值。问题是...