这是我遍的,但是有问题;
Z=[];
for tau=4:10
sol=dde23('V',tau,[1;0.2;1.6e-6;10;0.1],[0,50]);
for t=30:50
Z=[Z,tau+(sol.y(:,1)*i];
end
end
plot(Z,'.','markersize',1)
function dx=V(t,x,z);
xlag=z(:,1);
dx(1,1)=-x(5)*x(1)+x(2)*x(3)^x(4)*xlag(1)/(x(3)^x(4)+xlag(1));
dx(2,1)=x(2);
dx(3,1)=x(3);
dx(4,1)=x(4);
dx(5,1)=x(5);