时滞微分方程的求解。许多动力系统随时间的演化不仅依赖于系统当前的状态,而且依赖于系统过去
式中:T≥0为时滞常数。在Matlab中提供了命令dde23来直接求解时滞微分方程。其调用格式为801=dde23(ddefun,lags,history,tspan,options),
其中,ddfun为描述时滞微分方程的函数;lags为时滞常数向量;history为描述t≤to时的状态变量值的函数;tspan为求解的时间区间;options为求解器的参数设置。该函数的返回值sol是结构体数据,其中sol.x成员变量为时间向量l,sol.y成员变量为各个时刻的状态向量构成的矩阵,其每一个行对应着一个状态变量的取值。求解如下时滞微分方程组:
已知,在i≤0时,x(t)=5,x2(t)=0,x(1)=1,试求该方程组在[0,40]上的数值解。