用Python语言在MATLAB中模拟ode45函数我想知道如何将MATLAB函数ode45导出到python。根据文件规定如下: MATLAB: [t,y]=ode45(@vdp1,[0 20],[2 0]); Python: import ...2024-10-01 已阅读: n次
在表达式中查找和替换函数我需要在SymPy表达式树中找到所有用户生成的函数(以及像用户函数上的导数之类的东西),然后用标准的用户定义符号替换 from sympy import * t=symbols('t') 现在创建用 ...2024-10-01 已阅读: n次
用odeint确定分岔我试图得到分岔图的数据,即触发器[A]的值(由Aspace的linspace指定),在这组微分方程中给出了[N]的稳态(nullcline)值: 由于方程组太长,无法尝试通过重新排列来获得空斜,因此 ...2024-10-01 已阅读: n次
微分方程的数学积分问题我不懂编码,所以请耐心等待。我有一个练习,目的是解一个微分方程,其中有一个函数,它的积分形式为dx cos(x)**2,介于0和t之间 我试着定义一个函数,然后用这个定义在程序结束时得到一个最终值 i ...2024-10-01 已阅读: n次
我不知道无效令牌是什么意思我不知道无效令牌是什么意思。这是错误日志 runfile('C:/Users/evilq/Documents/PBPKmodel2.py', wdir='C:/Users/evilq/Document ...2024-10-01 已阅读: n次
用odein求ODEs的执行时间我想打印出方程组每t+dt的解的执行时间,其中dt是.01秒 import numpy as np from scipy.integrate import odeint import matplotl ...2024-10-01 已阅读: n次
使用神经网络预测常微分方程的初始条件给出的输出假设我们有一个描述X如何影响Y的颂歌系统: dXdt = -k * X dYdt = Kin * (1 - (Vmax * X)/(Km + X)) - kout * Y 我试图用神 ...2024-10-01 已阅读: n次
在函数中访问scipy.integrate.odeint中的当前时间步长有没有办法访问scipy.integrate.odeint中的当前时间步长 我试图解决一个颂歌系统,其中颂歌的形式取决于人口是否会枯竭。基本上,只要{}不低于阈值,我就从总体{}中获取。如果我需要采取 ...2024-10-01 已阅读: n次
微分方程的解在boost::odeint和集成电路我试图将我的快速原型从Python移植到C++。我试着用一个简单的微分方程来测试这个符号,但是对于起始值[2,0],结果却大不相同。 Python正在下降,而C++解决方案正在强劲地增长。你知道吗 它 ...2024-10-01 已阅读: n次
使用西皮·奥德林时间步长不同我正在用scipy.integrate.odeint有不同的时间步长。在 具体来说,我设置了一个一维扩散模拟,粒子在一端扩散到另一端。对于扩散,我使用有限差分二阶导数近似。在 令人惊讶的是,我得到了不 ...2024-10-01 已阅读: n次
如何修改此微分方程解算器以求解大量变量?我想修改下面的代码,以便它可以解决成千上万的变量与成千上万的夫妇微分方程。问题是我希望能够导入变量列表(y),理想情况下是作为numpy数组导入,但可以接受作为常规列表导入。这个列表会很庞大,所以我不 ...2024-10-01 已阅读: n次
将matlab ODE解算器转换为python我试图将MATLAB代码转换为python,但得到的答案完全不同。我尝试使用scipy.ode、solve_ivp和odeint。运行代码时,我得到的值范围为1到0.2,但在MATLAB中,它们的范围 ...2024-10-01 已阅读: n次