如何在 matplotlib 中制作两个滑块
..
odeint相关内容
为solve_ivp 传递参数(新的SciPy ODE API)
为了使用 SciPy 解决简单的 ODE，我曾经使用 odeint 函数，形式为: scipy.integrate.odeint(func, y0, t, args=(), Dfun=None, col_deriv=0, full_output=0, ml=None, mu=None, rtol=None, atol=None, tcrit=None, h0=0.0, hmax=0.0, hmi
..
在 Python 中优化微分方程中的常量
好的，那么我将如何编写代码来优化微分方程中的常量 a 和 b，例如 dy/dt = a*y^2 + b，使用curve_fit?我将使用 odeint 来求解 ODE，然后使用 curve_fit 来优化 a 和 b.如果您能就这种情况提供意见，我将不胜感激！ 解决方案 你绝对可以这样做: 将 numpy 导入为 np从 scipy.integrate 导入 odeint从 scip
..
多个 scipy.integrate.ode 实例
我想在多线程(每个CPU核心一个)中使用scipy.integrate.ode(或scipy.integrate.odeint)实例来解决一次多个 IVP.但是文档说:“这个集成器不是可重入的.你不能有两个 ode 实例同时使用“vode"集成器." (虽然文档没有说明，但如果多次实例化，odeint 也会导致内部错误.) 知道可以做什么吗? 解决方案 一种选择是使用 mult
..
在 PYTHON 中求解矩阵耦合微分方程时如何绘制特征值?
假设我们有三个复矩阵和一个带有这些矩阵的耦合微分方程组. 导入numpy、scipy从 numpy 导入(real、imag、matrix、linspace、array)从 scipy.integrate 导入 odeint导入 matplotlib.pyplot 作为 plt定义系统(x，t):a1= x[0];a3= x[1];a5= x[2];a7= x[3];a2= x[4];a4=
..
如何使用 scipy.odeint 求解微分方程组
我想使用 odeint 求解方程组，但出现以下错误: 文件“C:"，第 45 行，在 中C_B = odeint(dC_Bdt,C_B0,t)文件“C:\Anaconda3\envs\ChemEng\lib\site-packages\scipy\integrate\odepack.py"，第 233 行，在 odeint 中int(bool(tfirst)))RuntimeError: f
..
为什么在 SciPy 中使用集成.odeint 时不调用 Dfun(gradient)?
谁能提供一个为 SciPy 中的 integrate.odeint 函数提供雅可比行列式的示例?我尝试从 SciPy 教程中运行此代码 odeint 示例 但似乎Dfun()(雅可比函数)从未被调用. from numpy import * # 添加从 scipy.integrate 导入 odeint从 scipy.special 进口伽玛，通风y1_0 = 1.0/3**(2.0/3.0
..
具有多个参数的 ODEINT(时间相关)
我正在尝试使用 ODEINT 解决单个一阶 ODE.以下是代码.我希望在 3 个时间点获得 y 的 3 个值.我正在努力解决的问题是能够传递 mt 和 nt 的第 n 个值来计算 dydt.我认为 ODEINT 会传递 mt 和 nt 的所有 3 个值，而不是仅传递 0th、1st 或 2nd，具体取决于迭代.因此，我收到此错误: RuntimeError: func (4) 返回的数组大小
..
Python odeint 与微分方程中的数组
我有下一个一阶微分方程(示例): dn/dt=A*n;n(0)=28 当A为常数时，用python odeint完美解决.但是我有一个来自 .txt 文件的不同 A 值的数组[不是函数，只是一个值数组] A = [0.1,0.2,0.3,-0.4,0.7,...,0.0028] 而且我希望在求解 ode A 的每次迭代(或在每个时刻 t 中)是数组中的一个新值.我的意思是:第一次迭
..
odeint 中的参数数组
我正在尝试用 odeint 求解微分方程.这里一些常量参数是固定的，一些在列表中. from scipy.integrate import odeint导入 matplotlib.pyplot 作为 plt将 numpy 导入为 np从 scipy.interpolate 导入 LinearNDInterpolator#x 方向的运动方程 ===== ### d^2x/dt^2 = q[Ex
..
python中ODE的求解系统；在此调用上完成了多余的工作
我正在尝试在 python 中解决不同电位的耦合 ODE 系统.它适用于特定类型的势(指数)，但是一旦势由幂律描述，python 生成的图就完全不连贯，它经常只为所有参数分配零值.我的编码适用于: kr1 = 8*np.pi#rho_m = a**(-3)#V = np.e**(-kr1*x_1)#dVdx = -kr1*np.e**(-kr1*x_1)def quintessence (x,
..
使用 Scipy/Numpy-NDSolve 等效求解 Python 中的矩阵微分方程?
我有两个 numpy 数组:9x9 和 9x1.我想在离散时间点求解微分方程，但无法让 ODEInt 工作.我不确定我是否在做正确的事情. 使用 Mathematica，等式是: Solution = {A[t]}/.NDSolve[{A'[t] == Ab.A[t] &&A[0] == A0}, {A[t]}, {t, 0, .5}, MaxSteps ->\[无限]];时间 = 0.2
..
绘制一阶ODE的解及其导数
我有这段代码可以使用 odeint 解决一个简单的一阶 ODE.我设法绘制了解y(r)，但我也想绘制导数y'= dy/dr.我知道y'由f(y，r)给出，但是我不确定如何用积分的输出来调用函数.预先谢谢你. 来自数学导入sqrt的 从 numpy 导入零、linspace、数组从scipy.integrate导入odeint导入 matplotlib.pylab 作为 pltdef f(y，r
..
如何在matplotlib中制作两个滑块
我想在matplotlib中制作两个滑块来手动更改捕食者-猎物模型中的N和P值: 将 numpy 导入为 np导入matplotlib.pyplot作为plt从scipy.integrate导入odeintdef lotka(x,t,params):N, P = xalpha, beta, gamma, delta = 参数derivs = [alpha*N - beta*N*P, gamma*
..
如何用Odeint求解状态空间模型?
我正在尝试使用Eigen和Odeint对状态空间模型进行数值模拟.我的麻烦是，我需要引用控制数据 U (集成前预定义)，以便正确解决状态的 Ax + Bu 部分空间模型.我试图通过使用计数器来跟踪当前时间步长来实现此目的，但是由于某种原因，每次Odeint调用系统功能时，它都会重置为零. 我该如何解决?我对状态空间系统建模的方法有缺陷吗? 我的系统 struct Eigen_SS_
..
对耦合的ODE使用python内置函数
这部分只是您需要的背景 我正在为二阶仓本模型开发数值求解器.下面列出了我用于查找theta和omega导数的函数. #Ω中的n维变化def d_theta(omega):返回欧米茄#欧米茄的n维变化def d_omega(K，A，P，alpha，mask，n):def layer1(theta，omega):T = theta [:，None]-thetaA [掩码] = K [掩码] *
..
如何监视SciPy.odeint的过程?
SciPy可以通过scipy.integrate.odeint或其他软件包来求解ode方程，但是在函数完全求解后才能给出结果.但是，如果ode函数非常复杂，则该程序将花费大量时间(一到两天)才能得出整个结果.那么，如何控制求解方程的步骤(当方程尚未完全求解时打印出结果)? 解决方案 当我在寻找答案时，找不到满意的答案.因此，我使用 tqdm 项目.希望对您有帮助. 编辑:主持人请我解
..
耦合的OIDENT求解器中方程的顺序重要吗?
我的代码运行，但是，如果我更改了公式中定义的方程式的顺序，我的图形也会更改.有人可以告诉我为什么吗?现在，我不知道该系统中的正确图形是什么. def myFunction(r，t):g = 9.81L_L = 20#吃水L_r = 20#吃水L_d = 4 #ukcu_s = 0.08w_d = 60#宽度的容器rho = 1030b_R = 1.0b_L = 1.0b_D = 1.0A_s
..
多个ODE函数中的曲线拟合参数
我将不实施“>诊断延迟对传播的影响的SEIR模型使用Python 我不知道这是否是您想要的最佳选择，但我希望它能帮助您入门.
..
当类的属性包含另一个类实例时，如何指定numba jitclass?
我正在尝试使用numba来提高scipy.integrate.odeint的python性能. 为此，我必须对定义ODE系统的函数使用@ nb.jit(nopython = True).但是，此函数必须在我的程序中使用另一个python类实例作为参数.我还必须使用具有适当规格的@ nb.jitclass(spec)来初始化该类.这很好用，直到当我发现一个严重的问题时，类的规范包括另一种类型的类
..