考虑周期性的python的ODE求解器

文章/答案/技术大牛

发布

1回答

python、scipy、integration、ode

我想集成一个控制问题，即形式为dx/dt = A.f(t)的常微分方程，其中x(t)是R^3中的函数，f(t)是R^4中的函数，A是矩阵3x4。在我的特殊情况下，f(t) = F'(t)，即函数F的时间导数。此外，F是1周期的。因此，在区间0，1上积分ODE应该再次得到起始位置。然而，像scipy.integrate的solve_ivp这样的方法根本不考虑这种周期性(我已经尝试了

浏览 19提问于2021-01-22得票数 0

回答已采纳

1回答

带条件停止ODE求解器

r、ode

当因变量超过特定值时，我想终止ODE求解器。考虑一个简单的ODE模型dModel<- function(t, y, parms) {list(c(y))}我想在y>100时停止求解器并显示输出。我知道有一个ODE事件，但在阅读了文档后，我不知道如

浏览 1提问于2018-07-10得票数 0

1回答

matlab如何解微分方程？

matlab、differential-equations、numerical-analysis

这个问题不是关于如何使用matlab，而是试图找出当matab求解微分方程时发生了什么。当使用不同的数值方法时，结果通常是不同的。我想知道matlab中使用的是哪种数值方法。

浏览 5提问于2020-07-11得票数 0

1回答

dopri5“需要更大的nmax”

python、scipy

在使用scipy 0.13.0，ode(f).set_integrator('dopri5')时，我得到了错误消息-我在nmax中查找ode.py，但看不到变量我猜，要求集成的数字超过了允许的默认值。如何增加nmax值？

浏览 2提问于2014-09-29得票数 3

1回答

在NetLogo中求解ODEs，欧拉vs R-K vs R求解器

netlogo、runge-kutta

在我的模型中，每个智能体在每个节拍上求解一个ODE系统。我采用了欧拉法(类似于NetLogo中的系统动力学模型)来求解这些一阶常微分方程。然而，为了获得稳定的解决方案，我被迫使用非常小的时间步长(dt)，这意味着使用这种方法进行模拟的速度非常慢。我很好奇，有没有人有建议可以更快地解决这些常微分方程？我正在考虑实现Runge-Kutta (具有更大的时间步长？)正如这里所做的()。我也会考虑<

浏览 0提问于2015-04-10得票数 1

1回答

Scipy:挑选一个ode求解器实例

python、scipy、pickle、ode

我正在运行长时间的模拟，这涉及到一个ode的集成。由于运行时间长，我决定产生后援，允许从检查点继续模拟。对于ode集成，我决定选择实例。这种方法适用于可重入的求解器“dopri5”和“vode 853”，但对于“vode”和“lsoda”则失败。对于我的特殊问题，'vode‘algo似乎是最快的，所以我感兴趣的是，是否可以绕过这个问题。我希望这是可能的，因为在重新启动时也只需要一个

浏览 1提问于2018-01-04得票数 2

1回答

更改函数中的参数以与ODE求解器一起使用

matlab、if-statement、parameters、ode

是否有可能使用ODE求解器，例如ode45，并且仍然能够“更改”被调用函数中的参数值？但是，我希望有一种方法可以轻松地更改参数值，同时仍然使用一个函数和一个脚本运行ODE求解器的多次迭代。但是，似乎我不能简单地向ODE求解器中添加更多的术语，例如：a = P(1); b =

浏览 5提问于2012-10-16得票数 0

回答已采纳

1回答

Java编程，尝试创建一个ODE解算器库(面向对象编程的新手！)

java、ode

我正在尝试编写一个ODE求解器(类似于MATLAB中的ODE45 )。我希望这个求解器是“独立的”，即我可以传递给它我想要它求解的颂歌的名称，然后它运行那个特定微分方程的代码。到目前为止，我有我的主类(程序的主要部分)和一个名为ODEsolver的类( ODE求解器所在<em

浏览 1提问于2011-10-08得票数 0

回答已采纳

2回答

在寻找一个好的Java ODE解算器

java、matlab、solver、ode

我正在做一个使用ODE求解器(ode45)为MATLAB中的算法创建图形用户界面的项目。所以我必须把MATLAB代码翻译成Java。问题是ode45求解器。Java似乎没有现成的求解器，而且ODE也不是我的专长。是我看起来不太好，还是真的没有为Java实现的ODE解算器？谢谢

浏览 4提问于2012-07-17得票数 6

1回答

用SciPy求解一组常微分方程

python、scipy、ode

我正在尝试解决一组模拟淀粉酶(酶)水解淀粉的常微分方程。/usr/bin/env pythonimport math当我在Matlab上工作的时候，我经常用ode

浏览 0提问于2014-07-17得票数 2

4回答

将刚性代码与Python集成

python、scipy、pygsl

我正在寻找一个好的库，它将集成Python中的刚性代码。问题是，and‘s odeint有时会给我很好的解决方案，但是初始条件的微小变化会导致它倒下并放弃。同样的问题由MATLAB的刚性求解器(ode15s和ode23s)很高兴地解决了，但是我不能使用它(即使是来自Python，因为MATLAB的Python绑定没有一个实现回调，而且我需要将一个

浏览 2提问于2010-01-18得票数 21

回答已采纳

1回答

使用scipy.integrate.ode检索其他返回值

python、scipy、ode

我正在设置一个相当复杂的函数，我希望将其作为一个ODE来解决。我最好使用在大多数python环境中都可以访问的求解器，所以我尝试使用scipy.integrate.ode。问题是我的函数生成了很多我想在每次迭代中提取的值。下面是一个例子:我想在每次完成的r.integrate()迭代中获取someValues，并将它们保存在某个地方：

浏览 16提问于2021-03-02得票数 1

1回答

如何了解python中ode求解器的中间轨迹

python、math、scipy、ode、differential-equations

我正在使用Python包来求解ODE方程。但是，我需要知道在python中的中间状态，或者换句话说，ode求解器的轨迹。from scipy.integrate import odeint在这里，输出只是给我最后的状态，而不是中间的步骤，我如何才能得到中间的步骤？我如何知道初值和最终输出

浏览 7提问于2022-11-15得票数 1

回答已采纳

1回答

simulink中的采样时间计数器

simulink、sampling

如何从simulink中的“时钟”块中提取simulink时间，然后在每个采样周期Ts下增加matlab内嵌函数中的计数器(k)？假设simulink时钟时间不是，那么我们可以考虑 k=k+1;但这是行不通的，因为simulink时间是可变的。有什么想法吗？谢谢。

浏览 1提问于2014-08-15得票数 0

回答已采纳

1回答

如何测试ODE集成是否达到均衡？

matlab、ode

我介绍了对ode45求解器的一些修改。我有时使用多达64个组件，所有组件都在0,1间隔内，组件总和为1。在某些时间间隔，我会暂停集成过程，以便快速检查是否需要进一步集成，我正在寻找一些聪明的方法来有效地解决这个问题。2:三个或更多的组件以周期性的方式剧烈波动。 3:一个或两个分量变化非常快，振幅低，频率短。 4:以上均不正确，必须继续集成。给出一个想法:我发现将ode45求解器生成的</em

浏览 25提问于2020-11-04得票数 1

1回答

ODE45:与稀疏矩阵中的`expm`结果不同

matlab、matrix、ode、differential-equations

在从t1到t1+dt的特定时间步长dtK = ...Method2：K = ...[T, Y] = ode45(@(t,y)dy(y,K),[t1 t1+dt],y0);其中：ydot = -1i*K*y; 对于大型稀疏矩阵，这两种方法给出了不同的结果。哪一

浏览 1提问于2016-08-22得票数 1

1回答

用于DAE的隐式/向后Euler

matlab、differential-equations、implicit-methods

我有一个关于隐式Euler的问题。我知道如何计算隐式欧拉法，但我的问题是如何将其用于DAE (微分代数方程)。在对原始DAE应用索引缩减后，我获得了正确的解决方案，因此我获得了ODE，然后应用了隐式Euler。但是，任务是在DAE上部署隐式Euler。有没有人能给我一个提示，告诉我如何改进我的代码，使其也适用于DAE？非常感谢，请看我的代码附件。[t,y]=beul('system','dsyste

浏览 11提问于2017-10-19得票数 0

1回答

求解一阶ODE，其中包含另一个ODE* (Python中的odeint/ solve_ivp )*

python、numpy、ode、odeint

我试图用Python建立一个快速的数值求解器，以解决这个形式的微分问题：其中r是一个常数。我想在一段时间内整合A，这是我感兴趣的。然而，这是一个复杂的事实，因为dA/dt方程包含另一个变量B，它本身是由ODE /dt描述的。B实际上是一个向量，但我简化了表达式，试图更清楚地突出我的问题。目前，我有一个使用手动欧拉方法的解决方案:即计算dB/dt (然后使用

浏览 0提问于2018-11-04得票数 1

回答已采纳

2回答

矢量化SciPy ode解算器

python、scipy、ode

我的问题是关于当前的赋值颂歌求解者。在中，它们的用法是：def jac(t, y, arg1):r = ode:它使用了很多python循环(

浏览 5提问于2014-01-21得票数 2

回答已采纳

1回答

python中的Ode求解器

python、scipy、ode

我正试图用来求解一个ODE (dx^2/dt^2 = -4(x^2+y^2)^(3/2))，但我无法让它工作。这是我的代码：from scipy.integrate import odeint x1 = x[0]我得到了这个错误：但是，我将额外

浏览 0提问于2019-11-16得票数 0

回答已采纳

点击加载更多