一个关联某个未知函数本身及其各种阶数的衍生物的等式。
在 Sympy 中,我尝试求解这样的微分方程: 从 sympy 导入 * 从 sympy.physicals.vector 导入动态符号 x = 动态符号('x') diffeq = Eq(x(t).diff(t), x(t)) dsolve(差异...
我有一个包含 62500 个联立方程的微分方程组。系统动力学是线性的,其形式为: [dU/dt] = [雅可比]x[U]。 雅可比矩阵是稀疏的(非零值 ~ ...
显然,从 ODE 求解器获得非负解并非易事。在 Matlab 中,某些求解器可以使用 NonNegative 选项来获得非负解。有没有类似的选择...
我想用Python对以下DGL系统进行数值求解: 该过程应该与往常相同。我有 2 个二阶耦合微分方程,我使用替换...
我有以下问题。我想在 python 中对以下 DGL 系统进行数值求解: $$ 开始{对齐*} &g'' + rac{2}{r} g' - rac{2}{r^2} g - rac{3}{r} e g^2 - e^2 g^3 - r...
使用 DeepXDE 包且初始条件中包含可训练变量时,PyTorch 会发生“双重向后错误”
请询问有关如何调试以下 PyTorch 错误的任何想法。这是一个改编自有关反应扩散系统逆向建模的官方示例的示例: “&q...
我正在尝试使用 odeint 在 Python 中求解多个常微分方程。一个简单的版本如下所示: 将 numpy 导入为 np 从 scipy.integrate 导入 odeint deflotka_volterra(xy, ...
所以我一直在尝试对我在论文中找到的耦合一阶微分方程系统进行建模。它们具有生物性质,因此负面解决方案没有任何意义。 ...
我正在使用 MATLAB ODE 套件来求解由 13 个微分方程组组成的系统,这些微分方程组决定了神经元的行为。现在我想添加第二个神经元,它由同一组不同的神经元操作......
我目前正在尝试从以下论文中重新创建分叉图:https://journals.plos.org/plosone/article?id=10.1371/journal.pone.0011238#ack 这就是我的分叉图...
我正在 Julia 中求解随机微分方程组: #定义漂移函数 F!(du,u,t,p)=... #定义噪声函数 G!(du,u,t,p)=... #构建并解决问题 u_0=零...
我正在使用 scipy.integrate.odeint 模拟双摆翻转时间。我的代码的结构方式是,odeint 在固定的时间间隔内求解系统;然后我必须检查结果看看是否...
我正在使用 scipy.odeint() 模拟双摆翻转时间。我的代码的结构方式是,scipy 在给定的时间内解决了整个系统的问题;然后我必须检查数组以查看是否发生翻转...
用 Numerov 方法求解一维薛定谔方程 (python)
我目前正在尝试求解一维薛定谔方程。 (与时间无关)使用 Numerov 方法。该方法的推导对我来说很清楚,但我在实现时遇到了一些问题。我尝试...
我试图在Python中以数值方式求解(x’(t))^2+ax^2(t)-b=0。我知道我可以通过定义如下函数来使用 scipy 的solve_ivp 来求解更简单的方程: def f(t, y): 返回 y ** 2 这会
我需要使用Matlab找到曲面上的一个点(给定相对于起点的角度),其中弧长为给定值。 假设我有一个高阶曲面,其中 z=f(x,y) 是...
Matlab ODE - Simulink 积分器,m 文件中带有求和块,无需 simulink?
是否可以做这样的事情 ---- -------- --->|+ |-------->| 1/s |-----o---------> | -| | | | ---- -------- | /|\ ...
我想对微分方程 dc/dt 进行积分。下面是代码和变量的值。 全部清除; c1=.185;c0=2*10^-6;k3=.1*10^-6; v1=6;v2=.11;v3=.09*10^-6; Ca_ER=10*10^-6;Ca_cyto=1....
我正在尝试用 numpy 求解矩阵方程 $Ax = b$ (实际上我正在使用 scipy.sparse 但我相信问题仍然相同)。在我的设置中,$b$ 是某个函数的导数。 在我的...
我已经浏览了 scipy.integrate.ode 但我无法找到如何实际使用这些集成方法 dorpi5 和 dop853。 我想尝试集成 ode 集成 python 与