汽车加速度倒数曲线编程

If=0.218;Iw1=1.798;Iw2=3.598;
r=0.367;nt=0.85;f=0.013;CdA=2.77;i0=5.83;ig1=5.56;ig2=2.769;ig3=1.644;ig4=1.0;ig5=0.793;
m=3880; g=9.8;
sgm1=1+(Iw1+Iw2)/(m*r.^2)+If*ig1.^2*i0.^2*nt/(m*r.^2);
sgm2=1+(Iw1+Iw2)/(m*r.^2)+If*ig2.^2*i0.^2*nt/(m*r.^2);
sgm3=1+(Iw1+Iw2)/(m*r.^2)+If*ig3.^2*i0.^2*nt/(m*r.^2);
sgm4=1+(Iw1+Iw2)/(m*r.^2)+If*ig4.^2*i0.^2*nt/(m*r.^2);
sgm5=1+(Iw1+Iw2)/(m*r.^2)+If*ig5.^2*i0.^2*nt/(m*r.^2);
n=linspace(600,4000);
tq=-19.313+295.27*n/1000-(n/1000).^2*165.44+(n/1000).^3*40.874-(n/1000).^4*3.8445;
ft1=tq*ig1*i0*nt/r;ft2=tq*ig2*i0*nt/r;ft3=tq*ig3*i0*nt/r;ft4=tq*ig4*i0*nt/r;
ft5=tq*ig5*i0*nt/r;
fw=CdA*ua.^2/21.15;
ff=m*g*f;
a1=(ft1-fw-ff)/(sgm1*m);
a2=(ft2-fw-ff)/(sgm2*m);
a3=(ft3-fw-ff)/(sgm3*m);
a4=(ft4-fw-ff)/(sgm4*m);
a5=(ft5-fw-ff)/(sgm5*m);
ua1=0.377*r*n/(ig1*i0);
ua2=0.377*r*n/(ig2*i0);
ua3=0.377*r*n/(ig3*i0);
ua4=0.377*r*n/(ig4*i0);
ua5=0.377*r*n/(ig5*i0);
y1=1./a1;y2=1./a2;y3=1./a3;y4=1./a4;y5=1./a5;
ua=linspace(0,160);
plot(ua1,y1,ua2,y2,ua3,y3,ua4,y4,ua5,y5)
xlabel('ua/km/h')
ylabel('1/a')
legend('y1','y2','y3','y4','y5')
title('汽车的加速度倒数曲线')