列车追踪间隔距离仿真

Vmax=[18:1:160];
Lt=120;
Ls=60; a=1; b=0.9; Tr=1; Td=30;

if (Vmax>=(sqrt(2*a*(Lt+Ls))))
Tz=Tr+(Vmax./3.6)/b+Td+sqrt(2*(Lt+Ls)/a);
else  
Tz=(Tr+(Vmax./3.6)/b+Td+(2*a*(Lt+Ls)+(Vmax./3.6).^2)./(2*a*(Vmax./3.6)));
end
plot(Vmax,Tz);%绘制曲线
title('列车运行间隔图');
xlabel('Vmax(Km/h)');ylabel('Tz(s)');
hold on


Lmin=50;
if (Vmax>=(sqrt(2*a*(Lt+Ls))))
Tz=(Vmax./3.6)/(2*b)+Td+sqrt(2*(Lt+Ls)/a)+(Lmin/(Vmax./3.6)).*(2+floor((Tr*(Vmax./3.6)+(Vmax./3.6).^2/(2*b))/Lmin));
else  
Tz=(Vmax./3.6)/(2*b)+Td+(2*a*(Lt+Ls)+(Vmax./3.6).^2)./(2*a.*(Vmax./3.6))+(Lmin./(Vmax./3.6)).*(2+floor((Tr.*(Vmax./3.6)+(Vmax./3.6).^2/(2*b))/Lmin));
end
plot(Vmax,Tz,'r:');
hold on



Lmin=80;
if (Vmax>=(sqrt(2*a*(Lt+Ls))))
Tz=(Vmax./3.6)/(2*b)+Td+sqrt(2*(Lt+Ls)/a)+(Lmin/(Vmax./3.6)).*(2+floor((Tr*(Vmax./3.6)+(Vmax./3.6).^2/(2*b))/Lmin));
else