此函数必须为以下形式:
dydt = ddefun(t,y,Z)
其中 t 对应当前 t,y 是一个求 y(t) 近似值的列向量,Z(:,j) 用于为以 delays(t,y) 的分量 j 形式提供的时滞...输出是对应 f(t,y(t),y(d(1),…,y(d(k))) 的列向量。
delays返回时滞 d(j) 的列向量的函数句柄。时滞取决于 t 和 y(t) 两者。...有了这种形式的时滞函数,ddesd 的使用方法与 dde23 完全相同。...ddesd 以结构体 sol 的形式返回解。使用辅助函数 deval 和输出 sol 来计算区间 tspan = [t0,tf] 中的特定点 tint 的解。...sol.ypsol.x 网格点处的 y(x) 近似值
sol.solver求解器名称 ‘ddesd’
sol = ddesd(ddefun,delays,history,tspan,options) 的解算方法与上述方法相同