求助:对带时滞的方程组仿真
对带时滞的方程组仿真,该程序不报错也不显示图片,是工具有问题还是程序有问题,烦高手指点!
sol=dde23(@twodelay,[0.05 2],[1.4;2.5],[0,300]);
figure (1);
plot(sol.x,sol.y(1,:),sol.x,sol.y(2,:))
figure(2);
plot(sol.y(1,:),sol.y(2,:))
function y=twodelay(t,x,Z)
y=ones(2,1);
xlag1=Z(:,1);
xlag2=Z(:,2);
K=0.02;
k=0.1;
bt=1.77;
gm1=0.1;
gm2=2.4;
A=20;
y(1)=-[K+k*x(2)+bt*x(1)]*x(1)+2*exp(-gm1*0.5)*bt*xlag1(1)*xlag1(1);
y(2)=-gm2*x(2)+A*k*xlag2(2)*xlag2(2);