制动系统有关曲线的matlab程序

clcclearcloseallm=4080;g=9.8;G=m*g;h=0.845;L=3.950;a=2.100;b=L-a;m1=9324;g=9.8;G1=m1*g;h1=1.324;L1=3.950;a1=3.024;b1=L1-a1;%空载I曲线k=0:0.1:1;Fz1=G*(b+k*h)/L;Fz2=G*(a-k*h)/L;Fu1=Fz1.*k;Fu2=k.*Fz2;figure(1)plot(Fu1,Fu2,'r')gridonxlabel('{\itF}_{\rmxu1}/(N)')ylabel('{\itF}_{\rmxu2}/(N)')text(17300,9000,'\uparrow{\it空载}','FontSize',18)holdon%满载I曲线k=0:0.1:1;Fz1=G1*(b1+k*h1)/L1;Fz2=G1*(a1-k*h1)/L1;Fu1=Fz1.*k;Fu2=k.*Fz2;figure(1)plot(Fu1,Fu2,'b')gridonxlabel('{\itF}_{\rmxu1}/(N)')ylabel('{\itF}_{\rmxu2}/(N)')text(23879,30000,'\uparrow{\it满载}','FontSize',18)holdon%空载β线B=0.38;Fu=Fu1/B;Fu2=(1-B)*Fu;figure(2)plot(Fu1,Fu2,'g')gridonxlabel('{\itF}_{\rmxu1}/(N)')ylabel('{\itF}_{\rmxu2}/(N)')text(20456,41000,'\downarrow{\it空载}','FontSize',18)holdon%满载β线B1=0.4;Fu=Fu1/B1;Fu2=(1-B1)*Fu;figure(2);plot(Fu1,Fu2,'r')gridonxlabel('{\itF}_{\rmxu1}/(N)')ylabel('{\itF}_{\rmxu2}/(N)')text(36739,52000,'\uparrow{\it满载}','FontSize',18)holdon%空载f和r线fork=0:0.1:1;Fxb1=1:10:100000;Fxb2=Fxb1.*(L-k.*h)./(k.*h)-G.*b/h;figure(3)holdonplot(Fxb1,Fxb2,'c')D=L+k*h;Fxb2=k*G*a./D-k*h*Fxb1./D;figure(3)holdonplot(Fxb1,Fxb2,'r')gridonholdonlegend('f线','r线',10)xlabel('{\itF}_{\rmxb1}/(N)')ylabel('{\itF}_{\rmxb2}/(N)')endholdon