使用 boost odeint 的 DDE

Question

是否可以使用 C++ Boost 求解时间延迟微分方程 - odeint 图书馆 ？例如下面的等式：

x'(t) = r*x(t)*(1 - x(t-tau)), 

其中 tau 是时间延迟的常数值。

Answer 1

是的，你可以。但是 odeint 并不是专门为 DDE 设计的。使用 odeint 解决 DDE 有两种可能性：

1. 您将 x 及其离散化历史视为因变量并直接使用步进器。
2. 您仅将 x 视为因变量，并使用系统函数（您的 r.h.s.）传递历史。但在这种情况下，您应该只使用在多个时间步上评估状态的步进器，例如 Euler 或 RK2。

如果我有时间，我会写一个更具体的答案，也许用一些代码 sn-ps。