用matlab绘制汽车的行驶加速度曲线

用matlab绘制汽车的行驶加速度曲线n=600:10:4000;Tq=-19.313+295.27*(n/1000)-165.44*(n/1000).^2+40.874*(n/1000).^3-3.8445*(n/1000).^4;m=3880;g=9.8;nmin=600;nmax=4000;G=m*g;ig=[5.562.7691.6441.000.793];nT=0.85;r=0.367;f=0.013;CDA=2.77;i0=5.83;L=3.2;a=1.947;hg=0.9;If=0.218;Iw1=1.798;Iw2=3.598;Ft1=Tq*ig(1)*i0*nT/r;Ft2=Tq*ig(2)*i0*nT/r;Ft3=Tq*ig(3)*i0*nT/r;Ft4=Tq*ig(4)*i0*nT/r;Ft5=Tq*ig(5)*i0*nT/r;ua1=0.377*r*n/ig(1)/i0;ua2=0.377*r*n/ig(2)/i0;ua3=0.377*r*n/ig(3)/i0;ua4=0.377*r*n/ig(4)/i0;ua5=0.377*r*n/ig(5)/i0;Fw1=CDA*ua1.^2/21.15;Fw2=CDA*ua2.^2/21.15;Fw3=CDA*ua3.^2/21.15;Fw4=CDA*ua4.^2/21.15;Fw5=CDA*ua5.^2/21.15;Ff=G*f;deta1=1+(Iw1+Iw2)/(m*r^2)+(If*ig(1)^2*i0^2*nT)/(m*r^2);deta2=1+(Iw1+Iw2)/(m*r^2)+(If*ig(2)^2*i0^2*nT)/(m*r^2);deta3=1+(Iw1+Iw2)/(m*r^2)+(If*ig(3)^2*i0^2*nT)/(m*r^2);deta4=1+(Iw1+Iw2)/(m*r^2)+(If*ig(4)^2*i0^2*nT)/(m*r^2);deta5=1+(Iw1+Iw2)/(m*r^2)+(If*ig(5)^2*i0^2*nT)/(m*r^2);a1=(Ft1-Ff-Fw1)/(deta1*m);a2=(Ft2-Ff-Fw2)/(deta2*m);a3=(Ft3-Ff-Fw3)/(deta3*m);a4=(Ft4-Ff-Fw4)/(deta4*m);a5=(Ft5-Ff-Fw5)/(deta5*m);plot(ua1,a1,ua2,a2,ua3,a3,ua4,a4,ua5,a5);title('汽车的行驶加速度曲线');xlabel('ua/(km*h^-1)');ylabel('a/￡¨m*s^-2￡?');gtext('￠?'),gtext('￠ò'),gtext('￠ó'),gtext('￠?'),gtext('￠?');（资料素材和资料部分来自网络，供参考。可复制、编制，期待你的好评与关注）