c++ - 使用 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。

如果我有时间，我会写一个更具体的答案，也许还有一些代码片段。