DDE с использованием boost odeint

Question

Возможно ли решить дифференциальные уравнения с запаздыванием с использованием C++ Boost - odeint ? Для примера ниже уравнения:

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

где tau - постоянное значение для временной задержки.

Answer 1

Да, вы можете. Но odeint явно не предназначен для DDE. Существует две возможности для решения DDE с odeint:

1. Вы рассматриваете x и его дискретизированную историю как зависимые переменные и используете непосредственно степперы.
2. Вы рассматриваете только переменную x как зависимую переменную и передаете историю с помощью системной функции (ваш r.h.s.). Но в этом случае вы должны использовать только степперы, которые оценивают состояние при умножении ваших временных меток, таких как Euler или RK2.

Если у меня есть время, я напишу более конкретный ответ, возможно, с некоторыми фрагментами кода.