clear
Lh=60;
h=50;
yuliu=3;
w=25/10;
zcwz=57*(1-0.5);
y=(-25+w):w:-w;
x=-sqrt(25^2-y.^2);
x=[-yuliu,x];                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                   
l=x-(-Lh);
k=Lh-yuliu-zcwz;
yo=-k*h/(Lh-yuliu);
xo=-sqrt(k^2-yo^2)+x(1);
d=sqrt((xo-x).^2+yo^2);
xx1=(xo-x).*(l-d)./d+xo;
yy1=yo.*(l-d)./d+yo;
kaikou1=-Lh+(l-d);
kaikou2=(zcwz-Lh)*ones(size(kaikou1));
mutiao=abs(-Lh-x);
kaikoucd=kaikou2-kaikou1;
shuju=[mutiao;kaikou1;kaikou2;kaikoucd;xx1;xo*ones(size(xx1))];
plot([0,-Lh],[-25,-25]);
hold on
for i=1:10
    plot([x(i),-Lh],[-25+w*i,-25+w*i]);hold on;
    plot([x(i),x(i)],[-25+w*i,-25+w*(i-1)]);hold on
end
plot([0,-Lh],[25,25]);
hold on
for i=1:10
    plot([x(i),-Lh],[25-w*i,25-w*i]);hold on;
    plot([x(i),x(i)],[25-w*i,25-w*(i-1)]);hold on;
end
plot([-Lh,-Lh],[-25,25]),hold on;
%===========================================
x=-x;
plot([0,Lh],[-25,-25]);
hold on
for i=1:10
    plot([x(i),Lh],[-25+w*i,-25+w*i]);hold on;
    plot([x(i),x(i)],[-25+w*i,-25+w*(i-1)]);hold on
end
plot([0,Lh],[25,25]);
hold on
for i=1:10
    plot([x(i),Lh],[25-w*i,25-w*i]);hold on;
    plot([x(i),x(i)],[25-w*i,25-w*(i-1)]);hold on;
end
plot([Lh,Lh],[-25,25]),hold on;
axis equal

for i=1:10
    plot([kaikou1(i),kaikou2(i)],[-25-w/2+w*i,-25-w/2+w*i],'r');hold on;
end

for i=1:10
    plot([kaikou1(i),kaikou2(i)],[25+w/2-w*i,25+w/2-w*i],'r');hold on;
end
kaikou1=-kaikou1;
kaikou2=-kaikou2;
for i=1:10
    plot([kaikou1(i),kaikou2(i)],[-25-w/2+w*i,-25-w/2+w*i],'r');hold on;
end

for i=1:10
    plot([kaikou1(i),kaikou2(i)],[25+w/2-w*i,25+w/2-w*i],'r');hold on;
end