matlab ddesd,具有状态相关时滞的 DDE

编写时滞代码

首先，编写一个函数来定义方程组中的时滞。此方程组中唯一存在的时滞位于项 -y2(e1-y2(t)) 中。

function d = dely(t,y)

d = exp(1 - y(2));

end

注意：所有函数都作为局部函数包含在示例的末尾。

编写方程代码

现在，创建一个函数来编写方程的代码。此函数应具有签名 dydt = ddefun(t,y,Z)，其中：

t 是时间(自变量)。

y 是解(因变量)。

Z(n,j) 对时滞 yn(d(j)) 求近似值，其中时滞 d(j) 由 dely(t,y) 的分量 j 给出。

求解器会自动将这些输入传递给该函数，但是变量名称决定如何编写方程代码。在这种情况下：

Z(2,1)→y2(e1-y2(t))

function dydt = ddefun(t,y,Z)

dydt = [y(2);

-Z(2,1)*y(2)^2*exp(1 - y(2))];

end

编写历史解代码

接下来，创建一个函数来定义历史解。历史解是时间 t≤t0 的解。

function v = history(t) % history function for t < t0

v = [log(t);

1./t];

end

求解方程

最后，定义积分区间 [t0tf] 并使用 ddesd 求解器对 DDE 求解。

tspan = [0.1 5];

sol = ddesd(@ddefun, @