《MATLAB画驱动力-行驶阻力图-两个程序》由会员分享，可在线阅读，更多相关《MATLAB画驱动力-行驶阻力图-两个程序（2页珍藏版）》请在金锄头文库上搜索。
1、程序一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)
2、*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);ylabel(F/N);title(汽车驱动力-行驶阻力平衡图); gtext(Ft1)gtext(Ft2);gtext(Ft3);gtext(Ft4);gtext(Ft5);gtext(Ff+Fw);程序二 m=3880; g=9.8; G=m*g; ig=5.56 2.769 1.644 1.00 0.793; nT=0.85; r=0.367; f=0.01
3、3; CDA=2.77; i0=5.83; L=3.2; a=1.947; hg=0.9; If=0.218;Iw1=1.798;Iw2=3.598;nmin=600;nmax=4000; 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; 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
4、; 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; ua=0:5:120; Ff=G*f; Fw=CDA*ua.2/21.15; Fz=Ff+Fw; plot(ua1,Ft1,ua2,Ft2,ua3,Ft3,ua4,Ft4,ua5,Ft5,ua,Fz); title(汽车驱动力-行驶阻力平衡图); xlabel(ua(km*h-1); ylabel(Ft/N); gtext(Ft1),gtext(Ft2),gtext(Ft3),gtext(Ft4),gtext(Ft5),gtext(Ff+Fw);
