odeint 相关问题

Odeint是一个现代C ++库,用于数值求解常微分方程。

在多维数组上使用 odeint

我想将 odeint 应用于初始条件数组,并返回与这些初始条件大小相同的导数。我可以循环遍历每个初始条件,但我认为......

回答 1 投票 0

如何停止积分时间相关的 ODE,但时间跨度的上限取决于另一个变量

我正在求解一些与时间相关的方程式,这些方程式处理充电和放电时离子浓度的变化。 (c2、c3、c4、c5 离子的浓度) 例如:dc2dt = ((I/96485.0) - (k2c2S) - (

回答 0 投票 0

使用odeint求解不同阶数的多极子ode

我正在尝试解决三个问题,其中一个是一阶,其他是二阶。 这是方程式: 这是代码: #注意: p[0]=y p[1]=dy/dt p[2]=q p[3]=u p[4]=du/dt #在...

回答 0 投票 0

Odeint 与 openmp 的集成和并行计算

我正在尝试使用 odeint 在 C++ 中实现三体问题集成,但是要进化,系统需要很多时间(系统集成仅 50 年大约需要 10 小时),所以我想...

回答 0 投票 0

我想通过求解运动方程来跟踪一个小时的粒子路径。对于时间步长,t= np。排列 (0, 3600,10)。这是正确的吗?

我想通过求解运动方程来跟踪一个小时的粒子路径。对于时间步长, t= np。排列 (0, 3600,10)。 它是否正确? 持续时间是否取决于时间步长。如果我...

回答 1 投票 0

boost 中的 Runge Kutta k 向量

使用 boost::odeint 求解 ODE 系统时,我可以使用 stepper.do_step 函数及时控制每次迭代。这样,我可以存储 u(t)。 我正在研究一种算法

回答 1 投票 0

如何使scipy.odeint的速度更快?

我目前正在解决一个559个非线性微分方程的综合系统.我必须通过改变常数c1,c2 b和g来拟合一些实验数据的解,我正在使用 ...

回答 1 投票 0

如何使用args向scipy.integrate.solve_ivp传递数组?

我的代码是: import numpy as np from scipy.integration import odeint def func(t, x, u): dydt = (-x + u) 5 return dydt y0 = 0 t = np.linspace(0, 10, 1000) u = 1.0 sol = odeint(fund, y0=y0, t=...)

回答 1 投票 0

在python上用scipy.integrated的Radau方法解决中子动力学方程。

我试着用python中的RADAU方法求解有两个反馈(燃料温度反馈和冷却剂温度反馈)的中子动力学方程组,导入numpy作为 ...

回答 1 投票 0

Eigen :: VectorXd和Boost :: Odeint,不起作用

我正在测试Eigen :: VectorXd作为boost :: odeint的state_type。我使用以下代码:#include #include #include

回答 1 投票 0

Scipy整合:solve_ivp是否总是比odeint更差?

示例ODE的解析解决方案,我正在通过解决以下初始值问题来测试这两种方法之间的差异:y'= 2 * y-t您可以通过考虑...来解决此问题,方法是...

回答 1 投票 1

如何使用args将数组传递给scipy.integrate.solve_ivp?

如何将数组传递给scipy.integate.solve_ivp函数?现在u = 1.0,我想要的是u = np.random.uniform(-1,1,1000)。代码是:从scipy.integrate import numpy as np import numpy ...

回答 1 投票 1

为什么我的Python BZ Oregonator无法编译?

我正在尝试通过求解ODE来模拟Belousov-Zhabotinsky反应溶液中的颜色变化,并生成一个图形,该图形将使用odeint演示振荡。我有一个...

回答 1 投票 0

正在调用基类虚拟方法,而不是子级的重写实现

我正在使用BOOST :: ODEINT编写ODE集成器。它要求带有逻辑的类的operator()带有以下签名“ void operator()(const std :: vector&x,std :: vector&...

回答 1 投票 0

如何将时间作为变量包含在Python Gekko中?

我需要在模型中包括时间以解决一组复杂的微分方程。这是一个简单的问题,演示了常数k = 0.1和初始条件y(0)= 10的问题。 ...

回答 1 投票 1


我不知道为什么使用随附的python代码获得W与PT的多个图

我想得到三个图(W vs. X,W vs. y,W vs PT)。但是我可以得到两个合适的图(W vs. X和W vs. y)。不幸的是,我终于得到了W vs PT(绿线)的多个图形。我不知道如何...

回答 1 投票 0

是否可以在Python的ODEINT中更改我的ODE?

我正在尝试模拟一个具有微小差异的振荡系统:当身体向上运动时,我希望它使用特定的运动方程(ODE),而当身体向上时,我希望它使用另一个运动方程。]]] >

回答 1 投票 0

编写与odeint一起使用的方程式的循环

我有一个初值问题需要解决;微分方程是从一个类似如下的字典派生的:eqs = {'a':array([-1。,2.,4.,0.,...]),'b':array([1。 ,-10 ....

回答 1 投票 -1

odeint:无法根据规则“安全”将数组数据从dtype('complex128')转换为dtype('float64')

以下代码给出错误:无法根据规则'safe'从numpy.fft导入numpy,将数组数据从dtype('complex128')强制转换为dtype('float64'),从scipy.integrate导入fft。 。

回答 1 投票 0

© www.soinside.com 2019 - 2024. All rights reserved.