matlab ddesd用法例子,求解带有常规时滞的时滞微分方程 (DDE)

ddesd

求解带有常规时滞的时滞微分方程 (DDE)

语法

sol = ddesd(ddefun,delays,history,tspan)

sol = ddesd(ddefun,delays,history,tspan,options)

参数

ddefun用于对微分方程 y′(t) = f(t,y(t),y(d(1),...,y(d(k))) 的右侧进行计算的函数句柄。此函数必须为以下形式：

dydt = ddefun(t,y,Z)

其中 t 对应当前 t，y 是一个求 y(t) 近似值的列向量，Z(:,j) 用于为以 delays(t,y) 的分量 j 形式提供的时滞 d(j) 求 y(d(j)) 近似值。输出是对应 f(t,y(t),y(d(1),...,y(d(k))) 的列向量。

delays返回时滞 d(j) 的列向量的函数句柄。时滞取决于 t 和 y(t) 两者。ddesd 通过使用 min(d(j),t) 施加 d(j) ≤ t 要求。

如果所有时滞函数都采用 d(j) = t – τj 形式，则您可以将参数 delays 设置为常向量 delays(j) = τj。有了这种形式的时滞函数，ddesd 的使用方法与 dde23 完全相同。

history按以下三种方式之一指定 history：

一个 t 函数，要求 y

= history(t) 能够将 t ≤ t0 的解 y(t) 以列向量的形式返回

一个固定列向量(如果 y(t) 为常量)

来自之前积分的解 sol(如果此调用继续该积分)

tspan从 t0=tspan(1) 到 tf=tspan(end) 的积分区间，其中 t0

< tf