2024年4月15日发(作者：程筠)
汽车理论课后习题MATLAB编程
1.3 确定一轻型货车的动力性能（货车可装用4挡或5挡变速器，任选
其中的一种进行整车性能计算）：
1）绘制汽车驱动力与行驶阻力平衡图。
2）求汽车最高车速，最大爬坡度及克服该坡度时相应的附着率。
3）绘制汽车行驶加速度倒数曲线，用图解积分法求汽车用2档起步加速行
驶至70km/h的车速－时间曲线，或者用计算机求汽车用2档起步加速行驶至
70km/h的加速时间。
解：(1) 求汽车驱动力与行驶阻力平衡图和汽车最高车速：
n=[600:10:4000];
Tq=-19.313+295.27*(n/1000)-165.44*(n/1000).^2+40.874*(n/1000).^3-3.8445*(n/10
00).^4;
m=3880;g=9.8;nmin=600;nmax=4000;
G=m*g;
ig=[5.56 2.769 1.644 1.00 0.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;
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/s)');
ylabel('Ft(N)');
gtext('Ft1'),gtext('Ft2'),gtext('Ft3'),gtext('Ft4'),gtext('Ft5'),gtext('Ff+Fw');
zoom on;
[x,y]=ginput(1);
zoom off;
disp('汽车最高车速=');disp(x);disp('km/h');
汽车最高车速=
99.3006
km/h
1 / 34
(2)求汽车最大爬坡度程序：
n=[600:10:4000];
Tq=-19.313+295.27*(n/1000)-165.44*(n/1000).^2+40.874*(n/1000).^3-3.8445*(n/10
00).^4;
m=3880;g=9.8;nmin=600;nmax=4000;
G=m*g;
ig=[5.56 2.769 1.644 1.00 0.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;
ua1=0.377*r*n/ig(1)/i0;
Ff=G*f;
Fw1=CDA*ua1.^2/21.15;
Fz1=Ff+Fw1;
Fi1=Ft1-Fz1;
Zoom on;
imax=100*tan(asin(max(Fi1/G)));
disp('汽车最大爬坡度=');
disp(imax);
disp('%');
汽车最大爬坡度=
35.2197%
(3)求最大爬坡度相应的附着率和求汽车行驶加速度倒数曲线程序：
clear
2 / 34
2024年4月15日发(作者：程筠)
汽车理论课后习题MATLAB编程
1.3 确定一轻型货车的动力性能（货车可装用4挡或5挡变速器，任选
其中的一种进行整车性能计算）：
1）绘制汽车驱动力与行驶阻力平衡图。
2）求汽车最高车速，最大爬坡度及克服该坡度时相应的附着率。
3）绘制汽车行驶加速度倒数曲线，用图解积分法求汽车用2档起步加速行
驶至70km/h的车速－时间曲线，或者用计算机求汽车用2档起步加速行驶至
70km/h的加速时间。
解：(1) 求汽车驱动力与行驶阻力平衡图和汽车最高车速：
n=[600:10:4000];
Tq=-19.313+295.27*(n/1000)-165.44*(n/1000).^2+40.874*(n/1000).^3-3.8445*(n/10
00).^4;
m=3880;g=9.8;nmin=600;nmax=4000;
G=m*g;
ig=[5.56 2.769 1.644 1.00 0.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;
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/s)');
ylabel('Ft(N)');
gtext('Ft1'),gtext('Ft2'),gtext('Ft3'),gtext('Ft4'),gtext('Ft5'),gtext('Ff+Fw');
zoom on;
[x,y]=ginput(1);
zoom off;
disp('汽车最高车速=');disp(x);disp('km/h');
汽车最高车速=
99.3006
km/h
1 / 34
(2)求汽车最大爬坡度程序：
n=[600:10:4000];
Tq=-19.313+295.27*(n/1000)-165.44*(n/1000).^2+40.874*(n/1000).^3-3.8445*(n/10
00).^4;
m=3880;g=9.8;nmin=600;nmax=4000;
G=m*g;
ig=[5.56 2.769 1.644 1.00 0.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;
ua1=0.377*r*n/ig(1)/i0;
Ff=G*f;
Fw1=CDA*ua1.^2/21.15;
Fz1=Ff+Fw1;
Fi1=Ft1-Fz1;
Zoom on;
imax=100*tan(asin(max(Fi1/G)));
disp('汽车最大爬坡度=');
disp(imax);
disp('%');
汽车最大爬坡度=
35.2197%
(3)求最大爬坡度相应的附着率和求汽车行驶加速度倒数曲线程序：
clear
2 / 34