matlab画多个垂直的线段,新手求助，图形最后多出一条垂直线

clear

clc

z1=140;%齿轮齿数

z2=140;

mn=4;%法向模数

an=20/180*pi;%法向压力角

B=15.0908/180*pi;

b=75;b1=75;b2=75;%齿宽

mt=mn/cos(B);

at=atan(tan(an)/cos(B));%端面压力角,弧度

Bb=atan(tan(B)*cos(at));%基圆螺旋角

zv1=z1/cos(B)/cos(B)/cos(B);

zv2=z2/cos(B)/cos(B)/cos(B);

xn1=0.115;

xn2=0.03391032344401/cos(B);

xt1=xn1*cos(B);

xt2=0.03391032344401;

L12=580.6;

a12=mt*(z1+z2)/2;

y12=(L12-a12)/mt;

ha1=(1*cos(B)+y12-xt2)*mt;     %齿顶高,m

ha2=(1*cos(B)+y12-xt1)*mt;

hf1=(1*cos(B)+0.25*cos(B)-xt1)*mt;%齿根高,m

hf2=(1*cos(B)+0.25*cos(B)-xt2)*mt;

ra1=mt*z1/2+ha1;%齿顶圆半径

ra2=mt*z2/2+ha2;

rf1=mt*z1/2-hf1;%齿根圆角半径，m

rf2=mt*z2/2-hf2;

rb1=mt*z1/2*cos(at);

rb2=mt*z2/2*cos(at);

am12=acos(cos(at)/(1+2*y12/(z1+z2)));%由中心距变动系数y求得啮合角

a_a1=acos(rb1/ra1);

a_a2=acos(rb2/ra2);

n1=7500;  %r/h

tz=3600/z1/n1; %s

e=1/2/pi*(z1*(tan(a_a1)-tan(am12))+z2*(tan(a_a2)-tan(am12)))+b1*sin(B)/pi/mn;

et=1/2/pi*(z1*(tan(a_a1)-tan(am12))+z2*(tan(a_a2)-tan(am12)));

%H=ha1+hf1;

H=2.25*mn;

rm1=rb1/cos(am12);

r1=mt*z1/2;r2=mt*z2/2;

Bm=atan(tan(B)*(rm1/r1));%Bm是节圆柱上的螺旋角

PB1=r1*cos(at)*(tan(a_a1)-tan(am12));

PB2=r2*cos(at)*(tan(a_a2)-tan(am12));

B1B2=PB1+PB2;

B2B1=B1B2;

pb=mt*pi*cos(at);%基圆齿距

f=B1B2;

pt=pi*mt;

f2=4*pt-f-b*tan(Bb);

syms t

v=7500*2*pi*r1/3600;

L1=v*t/sin(Bb).*(t>=0&t

% fplot(L1,[0,tz]);

% hold on;

L2=(v*t+pt)/sin(Bb).*(t>=0&t<=(b*tan(Bb)-pt)/v)+b/cos(Bb).*(t>=(b*tan(Bb)-pt)/v&t<=(f-pt)/v)+(f-(pt+v*t-b*tan(Bb)))/sin(Bb).*(t>=(f-pt)/v&t

% fplot(L2,[0,tz]);

% hold on;

L3=(f-(2*pt-b*tan(Bb))-v*t)/sin(Bb).*(t>=0&t

% fplot(L3,[0,tz]);

% hold on;

L4=(f-(3*pt-b*tan(Bb))-v*t)/sin(Bb).*(t>=0&t<=(f+tan(Bb)*b-3*pt)/v)+0.*(t>=(f+tan(Bb)*b-3*pt)/v&t

% fplot(L4,[0,tz]);

L=L1+L2+L3+L4;

fplot(L,[0,tz]);