qq(i)=i*pi/180;qq1(i)=qq(i)-150*pi/180;
s3=h-h*qq1(i)/(q1-qs);v3=-w*h/(q1-qs);
x(i)=(s0+s3)*sin(qq(i))+e*cos(qq(i));y(i)=(s0+s3)*cos(qq(i))-e*sin(qq(i));
a(i)=(s0+s3)*cos(qq(i))-e*sin(qq(i));b(i)=(s0+s3)*sin(qq(i))-e*cos(qq(i));
xx(i)=x(i)+rr*b(i)/sqrt(a(i)*a(i)+b(i)*b(i));yy(i)=y(i)+rr*a(i)/sqrt(a(i)*a(i)+b(i)*b(i)); end
for i=301:1:360
qq(i)=i*pi/180;
x(i)=(s0+0)*sin(qq(i))+e*cos(qq(i));y(i)=(s0+0)*cos(qq(i))-e*sin(qq(i));
a(i)=(s0+0)*cos(qq(i))-e*sin(qq(i));b(i)=(s0+0)*sin(qq(i))-e*cos(qq(i));
xx(i)=x(i)+rr*b(i)/sqrt(a(i)*a(i)+b(i)*b(i));yy(i)=y(i)+rr*a(i)/sqrt(a(i)*a(i)+b(i)*b(i)); end
plot(x,y,'r',xx,yy,'g')%用plot函数绘制曲线
text(0,20,'理论轮廓线')%理论轮廓线的坐标位于为(0,20)
text(65,40,'实际轮廓线')%实际轮廓线的坐标位于(65,40)
hold on
附图: