输出是对应 f(t,y(t),y(d(1),…,y(d(k))) 的列向量。
delays返回时滞 d(j) 的列向量的函数句柄。时滞取决于 t 和 y(t) 两者。...如果想要 ddesd 计算此事件函数的所有零,则 direction(k) = 0;如果仅计算事件函数呈上升趋势时的零,则 +1,如果仅计算事件函数呈下降趋势时的零,则 -1。...,[0,5]);
使用函数 ddex1delays 指定的时滞以及 ddex1de 计算的微分方程对 [0,5] 区间上的 DDE 求解。...t ≤ 0 条件下的历史记录由函数 ddex1hist 计算求得。...计算时在 [0,5] 区间内放入了 100 个等间距点,以此来求解:
tint = linspace(0,5);
yint = deval(sol,tint);
同时使用以下函数绘图:
plot(tint