延迟微分方程
matlab提供了dde23求解非中性微分方程。...dde23的调用格式如下:
sol = dde23(ddefun,lags,history,tspan)
lags是延迟量,比如方程中包含y1(t-0.2)和y2(t-0.3)则可以使用lags=[...y(1)) );
y(2)*( (2+sin(t))*10^(-5) + 9*Z(1,2)/(1+Z(1,2)) – Z(2,1) )];
%y(1)表示x_1(t),因为dde求解的结果中sol...会有个x,为了区别用y(1)表示x_1(t);Z(1,1)表示时滞项x_1(t-0.1);Z(1,2)表示时滞项x_1(t-0.3)
sol = dde23(ddex1dez,[0.1, 0.3],[...,tau,history,tspan);
%[0.1, 0.3]是时滞,[2 2]是初值,[0, 50]是时间范围
figure;
% plot(sol.x,sol.y)
plot(sol.x