在动态模拟中直接包括时滞微分方程
以下展示一个微分时滞逻辑方程
相较于一个普通微分逻辑方程
的动态模拟.
Manipulate[
Module[{T = 50, sol, sold},
sol = First[y /. NDSolve[{y'[t] == a y[t](1 - y[t]), y[0] == 0.1}, y, {t, 0, T}]];sold = First[y /. NDSolve[{y'[t] == a y[t](1 - y[t - 1]), y[t /; t ≤ 0] == 0.1}, y, {t, 0, T}]];
ParametricPlot[{{sol[t], sol[t - 1]}, {sold[t], sold[t - 1]}}, {t, 1, T}, PlotRange -> {{0, 3}, {0, 3}}, PlotStyle -> {Gray, Blue}]], {{a, 1}, 0, 2}]