matlab求解时滞微分方程,dde23调用格式:
sol = dde23(ddefun,lags,history,tspan);
--ddefun函数句柄,求解微分方程y'=f(t,y(t),y(t-τ1),...,y(t-τk))
必须写成下面形式:
dydt =ddefun(t,y,Z);
其中t对应当前时间t,y为列向量,近似于y(t);Z(:,j)近似于y(t-τj)
--lags为延迟时间,为正常数。
例:方程中包含y1(t-0.2)和y2(t-1),则可以表示为lags=[0.2,1]
--history t≤t0 时的状态变量的值