精选优质文档-倾情为你奉上程序一m1=2000m2=1800;m=3880;r0=0.367gt=0.85;f=0.013;CDA=2.77;i0=5.83;If=0.218;Iw1=1.798;Iw2=3.598;Ig5=5.56 2.769 1.644 1.00 0.793;for i=1:3401;n(i)=i+599;Ttq(i)=-19.313+295.27*(n(i)/1000)-165.44*(n(i)/1000).2+40.874*(n(i)/1000).3-3.8445*(n(i)/1000).4;end for i=1:3401;for j=1:5; Ft(i,j)=Ttq(i)*i0*Ig5(j)*gt/r0;ua(i,j)=0.377*r0*n(i)./(Ig5(j)*i0);F(i,j)=f*m*9.8+CDA*ua(i,j).2/21.15;endend plot(ua,Ft,ua,F)xlabel(ua/(km/h);y