如何在python中为ODE中的变量添加约束?我需要解一个微分方程组,我有一个返回dx/dt和du/dt的函数,但我需要加上u在区间[0,1]内的约束。我该怎么做?到目前为止,我一直在使用scipy.integrate中的odeint来解决这个问 ...2024-06-26 已阅读: n次
Python中直线法求解偏微分方程的时空参数我希望调整this method of lines based solution of a pde,以便k是空间和时间的函数,或者如果不满足阈值条件,k等于零。这不是我正在解决的实际偏微分方程,但它是 ...2024-06-26 已阅读: n次
如何在求解耦合微分方程时引入突变?我一直在尝试使用odeint和solve_ivp求解以下三个耦合微分方程组。(这不是确切的系统。我修改它只是为了提高问题的可读性。) def model(x,t): # x is a list ...2024-06-26 已阅读: n次
空间离散线法求解des系统我正试图使this example using the method of lines so solve a pde适应des系统。这不是我试图解决的系统,只是一个例子 如何合并第二个de?我已尝试将 ...2024-06-26 已阅读: n次
在python中使用scipy.integrate.solve_ivp时,有没有办法绑定变量的值?我想在python中用两个变量x和u来求解IVP,但我需要u的值在0和1之间。现在它给了我一个关于u的负值的解决方案。这是我的密码 def model_case_3(t, z, Kmax, k, b, ...2024-06-26 已阅读: n次
Pyomo中的时延响应我试图用几个ODE模拟一个DAE系统,其中一个(控制器)显示出与模拟时间尺度相比的巨大时滞。我应该如何在Pyomo中实现这一点(不考虑其他软件包,Gekko已经考虑过了,但由于apm.exe sour ...2024-06-26 已阅读: n次
共解向量函数我把一个修正的FitzHugh–Nagumo model的向量函数lambdiE(没有扩散项): from sympy import symbols, Matrix, solve, Eq from s ...2024-06-26 已阅读: n次
在pythoncapicod中引发异常并停止执行我的代码有两个Python函数,run()和dudt(u,t,…)。run调用一个C函数odepack_odeint,它调用Fortran子例程LSODA,后者调用另一个C函数ode_function ...2024-06-26 已阅读: n次
Python:比较NumPy和SymPy-ufuncified函数的速度刚刚写了一个代码来比较一个函数的计算速度,这个函数写为numpy,另一个函数使用ufuncifyfrom{}: import numpy as np from sympy import symbols ...2024-06-26 已阅读: n次
如何用谱方法解偏微分方程?我想用谱方法来解偏微分方程。像这样的方程,formula,初始条件是u(t=0,x)=(a^2)*sech(x),u''ut(t=0)=0 为了解决这个问题,我使用python和spectrum方法。 ...2024-06-26 已阅读: n次
scipy的solve_ivp未集成tspan的整个范围我试图在Python中使用来自scipy的solve_ivp来解决IVP问题。我将solve_ivp的tspan参数指定为(0,10),如下所示。然而,由于某种原因,我得到的解总是停在t=2.5左右 ...2024-06-26 已阅读: n次