汽车理论课后习题Matlab程序

合集下载

汽车理论习题Matlab程序

1.3 确定一轻型货车的动力性能（货车可装用4挡或5挡变速器，任选其中的一种进行整车性能计算）：1）绘制汽车驱动力与行驶阻力平衡图。

2）求汽车最高车速，最大爬坡度及克服该坡度时相应的附着率。

3）绘制汽车行驶加速度倒数曲线，用图解积分法求汽车用2档起步加速行驶至70km/h 的车速－时间曲线，或者用计算机求汽车用2档起步加速行驶至70km/h 的加速时间。

轻型货车的有关数据：汽油发动机使用外特性的Tq -n 曲线的拟合公式为23419.313295.27()165.44()40.874() 3.8445()1000100010001000q n n n n T =-+-+- 式中，Tq 为发动机转矩（N•m ）;n 为发动机转速（r/min ）。

发动机的最低转速n min =600r/min,最高转速n max =4000r/min 。

装载质量 2000kg整车整备质量 1800kg总质量 3880kg车轮半径 0.367m传动系机械效率 ηt =0.85滚动阻力系数 f =0.013空气阻力系数×迎风面积 C D A =2.77m 2主减速器传动比 i 0=5.83飞轮转动惯量 I f =0.218kg•m 2二前轮转动惯量 I w1=1.798kg•m 2四后轮转动惯量 I w2=3.598kg•m 2变速器传动比 ig(数据如下表)轴距L=3.2m质心至前轴距离（满载）a=1.974m质心高（满载）hg=0.9m解：Matlab程序：(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.3006km/h(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)求最大爬坡度相应的附着率和求汽车行驶加速度倒数曲线程序：clearn=[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;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);ad1=1./a1;a2=(Ft2-Ff-Fw2)/(deta2*m);ad2=1./a2;a3=(Ft3-Ff-Fw3)/(deta3*m);ad3=1./a3;a4=(Ft4-Ff-Fw4)/(deta4*m);ad4=1./a4;a5=(Ft5-Ff-Fw5)/(deta5*m);ad5=1./a5;plot(ua1,ad1,ua2,ad2,ua3,ad3,ua4,ad4,ua5,ad5);axis([0 99 0 10]);title('汽车的加速度倒数曲线');xlabel('ua(km/h)');ylabel('1/a');gtext('1/a1');gtext('1/a2');gtext('1/a3');gtext('1/a4');gtext('1/a5'); a=max(a1);af=asin(max(Ft1-Ff-Fw1)/G);C=tan(af)/(a/L+hg*tan(af)/L);disp('假设后轮驱动，最大爬坡度相应的附着率=');disp(C);假设后轮驱动，最大爬坡度相应的附着率=0.4219(4)>>clearnT=0.85;r=0.367;f=0.013;CDA=2.77;i0=5.83;If=0.218;Iw1=1.798;Iw2=3.598;L=3.2;a=1.947;hg=0.9;m=3880;g=9.8; G=m*g; ig=[5.56 2.769 1.644 1.00 0.793];nmin=600;nmax=4000;u1=0.377*r*nmin./ig/i0;u2=0.377*r*nmax./ig/i0;deta=0*ig;for i=1:5deta(i)=1+(Iw1+Iw2)/(m*r^2)+(If*(ig(i))^2*i0^2*nT)/(m*r^2);endua=[6:0.01:99];N=length(ua);n=0;Tq=0;Ft=0;inv_a=0*ua;delta=0*ua;Ff=G*f;Fw=CDA*ua.^2/21.15;for i=1:Nk=i;if ua(i)<=u2(2)n=ua(i)*(ig(2)*i0/r)/0.377;Tq=-19.313+295.27*(n/1000)-165.44*(n/1000)^2+40.874*(n/1000)^3-3.8445*(n/100 0)^4;Ft=Tq*ig(2)*i0*nT/r;inv_a(i)=(deta(2)*m)/(Ft-Ff-Fw(i));delta(i)=0.01*inv_a(i)/3.6;elseif ua(i)<=u2(3)n=ua(i)*(ig(3)*i0/r)/0.377;Tq=-19.313+295.27*(n/1000)-165.44*(n/1000)^2+40.874*(n/1000)^3-3.8445*(n/100 0)^4;Ft=Tq*ig(3)*i0*nT/r;inv_a(i)=(deta(3)*m)/(Ft-Ff-Fw(i));delta(i)=0.01*inv_a(i)/3.6;elseif ua(i)<=u2(4)n=ua(i)*(ig(4)*i0/r)/0.377;Tq=-19.313+295.27*(n/1000)-165.44*(n/1000)^2+40.874*(n/1000)^3-3.8445*(n/100 0)^4;Ft=Tq*ig(4)*i0*nT/r;inv_a(i)=(deta(4)*m)/(Ft-Ff-Fw(i));delta(i)=0.01*inv_a(i)/3.6;elsen=ua(i)*(ig(5)*i0/r)/0.377;Tq=-19.313+295.27*(n/1000)-165.44*(n/1000)^2+40.874*(n/1000)^3-3.8445*(n/100 0)^4;Ft=Tq*ig(5)*i0*nT/r;inv_a(i)=(deta(5)*m)/(Ft-Ff-Fw(i));delta(i)=0.01*inv_a(i)/3.6;enda=delta(1:k);t(i)=sum(a);endplot(t,ua);axis([0 80 0 100]);title('汽车2档原地起步换挡加速时间曲线');xlabel('时间t（s）');ylabel('速度ua（km/h）');>> ginputans =25.8223 70.073725.7467 70.0737所以汽车2档原地起步换挡加速行驶至70km/h 的加速时间约为25.8s2.7已知货车装用汽油发动机的负荷特性与万有特性。

汽车理论习题Matlab程序

1.3 确定一轻型货车的动力性能（货车可装用4挡或5挡变速器，任选其中的一种进行整车性能计算）：1）绘制汽车驱动力和行驶阻力平衡图。

2）求汽车最高车速，最大爬坡度及克服该坡度时相应的附着率。

轻型货车的有关数据：汽油发动机使用外特性的Tq-n 曲线的拟合公式为23419.313295.27()165.44()40.874() 3.8445()1000100010001000q n n n n T =-+-+- 式中，Tq 为发动机转矩（N•m ）;n 为发动机转速（r/min ）。

发动机的最低转速n min =600r/min,最高转速n max =4000r/min 。

装载质量 2000kg整车整备质量 1800kg总质量 3880kg车轮半径 0.367m传动系机械效率 ηt =0.85滚动阻力系数 f =0.013空气阻力系数×迎风面积 C D A =2.77m 2主减速器传动比 i 0=5.83飞轮转动惯量 I f =0.218kg•m 2二前轮转动惯量 I w1=1.798kg•m 2四后轮转动惯量 I w2=3.598kg•m 2变速器传动比 ig(数据如下表)轴距L=3.2m质心至前轴距离（满载）a=1.974m质心高（满载）hg=0.9m解：Matlab程序：(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.3006km/h(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)求最大爬坡度相应的附着率和求汽车行驶加速度倒数曲线程序：clearn=[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;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);ad1=1./a1;a2=(Ft2-Ff-Fw2)/(deta2*m);ad2=1./a2;a3=(Ft3-Ff-Fw3)/(deta3*m);ad3=1./a3;a4=(Ft4-Ff-Fw4)/(deta4*m);ad4=1./a4;a5=(Ft5-Ff-Fw5)/(deta5*m);ad5=1./a5;plot(ua1,ad1,ua2,ad2,ua3,ad3,ua4,ad4,ua5,ad5);axis([0 99 0 10]);title('汽车的加速度倒数曲线');xlabel('ua(km/h)');ylabel('1/a');gtext('1/a1');gtext('1/a2');gtext('1/a3');gtext('1/a4');gtext('1/a5'); a=max(a1);af=asin(max(Ft1-Ff-Fw1)/G);C=tan(af)/(a/L+hg*tan(af)/L);disp('假设后轮驱动，最大爬坡度相应的附着率=');disp(C);假设后轮驱动，最大爬坡度相应的附着率=0.4219(4) >>clearnT=0.85;r=0.367;f=0.013;CDA=2.77;i0=5.83;If=0.218;Iw1=1.798;Iw2=3.598;L=3.2;a=1.947;hg=0.9;m=3880;g=9.8; G=m*g; ig=[5.56 2.769 1.644 1.00 0.793];nmin=600;nmax=4000;u1=0.377*r*nmin./ig/i0;u2=0.377*r*nmax./ig/i0;deta=0*ig;for i=1:5deta(i)=1+(Iw1+Iw2)/(m*r^2)+(If*(ig(i))^2*i0^2*nT)/(m*r^2);endua=[6:0.01:99];N=length(ua);n=0;Tq=0;Ft=0;inv_a=0*ua;delta=0*ua;Ff=G*f;Fw=CDA*ua.^2/21.15;for i=1:Nk=i;if ua(i)<=u2(2)n=ua(i)*(ig(2)*i0/r)/0.377;Tq=-19.313+295.27*(n/1000)-165.44*(n/1000)^2+40.874*(n/1000)^3-3.8445*(n/100 0)^4;Ft=Tq*ig(2)*i0*nT/r;inv_a(i)=(deta(2)*m)/(Ft-Ff-Fw(i));delta(i)=0.01*inv_a(i)/3.6;elseif ua(i)<=u2(3)n=ua(i)*(ig(3)*i0/r)/0.377;Tq=-19.313+295.27*(n/1000)-165.44*(n/1000)^2+40.874*(n/1000)^3-3.8445*(n/100 0)^4;Ft=Tq*ig(3)*i0*nT/r;inv_a(i)=(deta(3)*m)/(Ft-Ff-Fw(i));delta(i)=0.01*inv_a(i)/3.6;elseif ua(i)<=u2(4)n=ua(i)*(ig(4)*i0/r)/0.377;Tq=-19.313+295.27*(n/1000)-165.44*(n/1000)^2+40.874*(n/1000)^3-3.8445*(n/100 0)^4;Ft=Tq*ig(4)*i0*nT/r;inv_a(i)=(deta(4)*m)/(Ft-Ff-Fw(i));delta(i)=0.01*inv_a(i)/3.6;elsen=ua(i)*(ig(5)*i0/r)/0.377;Tq=-19.313+295.27*(n/1000)-165.44*(n/1000)^2+40.874*(n/1000)^3-3.8445*(n/100 0)^4;Ft=Tq*ig(5)*i0*nT/r;inv_a(i)=(deta(5)*m)/(Ft-Ff-Fw(i));delta(i)=0.01*inv_a(i)/3.6;enda=delta(1:k);t(i)=sum(a);endplot(t,ua);axis([0 80 0 100]);title('汽车2档原地起步换挡加速时间曲线');xlabel('时间t（s）');ylabel('速度ua（km/h）');>> ginputans =25.8223 70.073725.7467 70.0737所以汽车2档原地起步换挡加速行驶至70km/h 的加速时间约为25.8s2.7已知货车装用汽油发动机的负荷特性和万有特性。

汽车理论课后题matlab程序

n=600:1:4000; r=; i0=; eff=; f=;m=3880; g=; G=m*g; CdA=; a=; hg=; L=;Iw1=; Iw2=; Iw=Iw1+Iw2; If=;Ttq=+*n/*(n/1000).^2+*(n/1000).^*(n/1000).^4;%驱动力行驶阻力平衡图for ig=[,,,,]Ua=*r*n/ig/i0;Ft=Ttq*ig*i0*eff/r;plot(Ua,Ft);hold on;endFf=G*f;ua=0::max(Ua);Fw=CdA*ua.^2/;plot(ua,(Ff+Fw)); title('驱动力-行驶阻力平衡图');xlabel('Ua/(km/h)');ylabel('Ft/N');gtext('Ft1'),gtext('Ft2'),gtext('Ft3'),gtext('Ft4'),gtext('Ft5'),gtext('Ff+Fw') [x,y]=ginput(1); disp('汽车的最高车速');disp(x);disp('km/h');%最大爬坡度及最大爬坡度时的附着率Ua=*r*n/i0;Ft=Ttq**i0*eff/r;Fw=CdA*Ua.^2/;i=tan(asin((Ft-(Ff+Fw))/G));disp('汽车的最大爬坡度');disp(max(i));C=max(i)/(a/L+hg/L*max(i));disp('克服最大爬坡度时的附着率');disp(C);%加速度倒数曲线figure;for ig=[,,,,]Ua=*r*n/ig/i0;q=1+Iw/(m*r^2)+If*ig^2*i0^2*eff/(m*r^2);Ft=Ttq*ig*i0*eff/r;Fw=CdA*Ua.^2/;as=(Ft-(Ff+Fw))/q/m;plot(Ua,1./as);hold on;endaxis([0 98 0 10]);title('行驶加速度倒数曲线');xlabel('Ua/(km/h)');ylabel('1/a');gtext('1/a1'),gtext('1/a2'),gtext('1/a3'),gtext('1/a4'),gtext('1/a5'); %加速时间曲线u(1)=*r*600/i0/;dt=;w=1;t(1)=0;ig=;while(u(w)<n(w)=u(w)*ig*i0/r;Ttq(w)=+*n(w)/*(n(w)/1000)^2+*(n(w)/1000)^*(n(w)/1000)^4;Ft(w)=Ttq(w)*ig*i0*eff/r;Ff=G*f;Fw(w)=CdA*u(w)^2/;q=1+Iw/(m*r^2)+If*ig^2*i0^2*eff/(m*r^2);a(w)=*(Ft(w)-(Ff+Fw(w)))/q/m;u(w+1)=u(w)+a(w)*dt;if(u(w+1)>=*4000*r/i0)ig=;endif(u(w+1)>=*4000*r/i0) ig=; endif(u(w+1)>=*4000*r/i0) ig=; endt(w+1)=(w+1)*dt; w=w+1; endfigure;plot(t,u);axis([0 60 0 100]);title('2挡原地起步加速时间曲线');xlabel('时间t/s'),ylabel('车速u/(km/h)');[x1,y1]=ginput(1);disp('加速到70km/h 的时间');disp(x1);0204060801001202000400060008000100001200014000驱动力-行驶阻力平衡图Ua/(km/h)F t /N102030405060708090012345678910行驶加速度倒数曲线Ua/(km/h)1/a10203040506001020304050607080901002挡原地起步加速时间曲线时间t/s车速u /(k m /h)Iw1=; Iw2=; Iw=Iw1+Iw2; If=; m=1600; g=;G=m*g; Ttqmax=140; i1=; i0=; eff=; r=; hg=; b=; L=; q=1+Iw/(m*r^2)+If*i1^2*i0^2*eff/(m*r^2); Ft=Ttqmax*i1*i0*eff/r;a=Ft/q/m;q=a/g;C1=q/(b/L-hg/L*q);if(C1>disp('发动机转矩不能充分发挥');elsedisp('发动机转矩能够充分发挥');endb=(q/+hg/L*q)*L;disp('调整后的b值');disp(b);disp('调整后的前轴负荷率'),disp(b/L);Iw=; If=; m=2152; g=; r=; eff=; i4=1; i0=;amax=; Ua=50/; CdA=;q=1+Iw/(m*r^2);f=q*g;disp('汽车在该路面上的滚动阻力系数');disp(f);q=1+Iw/(m*r^2)+If*i4^2*i0^2*eff/(m*r^2);D4max=f+q*amax/g;disp('最大动力因数');disp(D4max);Fw=CdA*Ua^2/;D1max=(D4max*G+Fw)*i1/i4/G;imax=tan(asin((D1max-f*sqrt(1-D1max^2+f^2))/(1+f^2))); disp('最大爬坡度'),disp(imax);n=600:1:4000; r=; i0=; eff=; f=;CdA=;m=3880; g=; G=m*g;Ttq=+*n/*(n/1000).^2+*(n/1000).^*(n/1000).^4;for ig=[,,,,]Ua=*r*n/ig/i0;Pe=Ttq.*n/9550;plot(Ua,Pe);hold on;endUa=0::max(Ua);Pf=G*f*Ua/3600;Pw=CdA*Ua.^3/76140;plot(Ua,(Pf+Pw)/eff);title('汽车的功率平衡图'),xlabel('Ua/(km/h)'),ylabel('P/kw');%等速百公里燃油消耗n1=[815,1207,1614,2012,2603,3006,3403,3804];b0=[,,,,,,,];b1=[,,,,,,,] ;b2=[,,,,,,,];b3=[,,,,,,,];b4=[,,,,,,,];B0=spline(n1,b0,n); B1=spline(n1,b1,n); B2=spline(n1,b2,n); B3=spline(n1,b3,n); B4=spline(n1,b4,n);deng=7;figure;for ig=[,]Ua=*r*n/ig/i0;Pf=G*f*Ua/3600;Pw=CdA*Ua.^3/76140;P=(Pf+Pw)/eff;b=B0+B1.*P+B2.*P.^2+B3.*P.^3+B4.*P.^4; Q=P.*b./*Ua*deng); plot(Ua,Q); hold on; endaxis([0 100 12 28]);title('最高挡与次高挡等速百公里油耗曲线'),xlabel('Ua/(km/h)'),ylabel('百公里油耗Qs/(L/100km)');204060801001200102030405060708090100汽车的功率平衡图Ua/(km/h)P /k w020406080100121416182022242628最高挡与次高挡等速百公里油耗曲线Ua/(km/h)百公里油耗Q s /(L /100k m )beta=; L=; a1=; b1=L-a1; hg1=; a2=; b2=L-a2; hg2=;g=; z=0::1;faif1=beta*z./(1/L*(b1+z*hg1)); fair1=(1-beta)*z./(1/L*(a1-z*hg1)); faif2=beta*z./(1/L*(b2+z*hg2)); fair2=(1-beta)*z./(1/L*(a2-z*hg2));plot(z,faif1,'r',z,fair1,'g',z,faif2,'b',z,fair2,'k'); fai=0::1;Ef1=(b1/L)./(beta-fai*hg1/L); Er1=(a1/L)./((1-beta)+fai*hg1/L); Ef2=(b2/L)./(beta-fai*hg2/L); Er2=(a2/L)./((1-beta)+fai*hg2/L); figure;plot(fai,Ef1,'r',fai,Er1,'g',fai,Ef2,'b',fai,Er2,'k');axis([0 1 0 1]);fai=;Er1=(a1/L)/((1-beta)+fai*hg1/L);Er2=(a2/L)/((1-beta)+fai*hg2/L);abmax1=fai*Er1*g;abmax2=fai*Er2*g;t21=; t22=; u0=30;s1=1/*(t21+t22/2)*u0+u0^2/abmax1;s2=1/*(t21+t22/2)*u0+u0^2/abmax2;disp('空载时的制动距离'),disp(s1);disp('满载时的制动距离'),disp(s2);%前部管路损坏时制动距离z1=a1*fai/(L+fai*hg1); z2=a2*fai/(L+fai*hg2); abmax1=z1*g; abmax2=z2*g;s1=1/*(t21+t22/2)*u0+u0^2/abmax1;s2=1/*(t21+t22/2)*u0+u0^2/abmax2;disp('前部管路损坏时空载制动距离');disp(s1); disp('前部管路损坏时满载制动距离');disp(s2); %后部管路损坏时的制动距离z1=b1*fai/(L-fai*hg1);z2=b2*fai/(L-fai*hg2); abmax1=z1*g; abmax2=z2*g;s1=1/*(t21+t22/2)*u0+u0^2/abmax1;s2=1/*(t21+t22/2)*u0+u0^2/abmax2;disp('后部管路损坏时空载制动距离');disp(s1); disp('后部管路损坏时满载制动距离');disp(s2);hg=; b=; L=; beta=; fai=; g=;fai0=(L*beta-b)/hg;if(fai<fai0)disp('前轮先抱死,制动效率为'); E=(b/L)/(beta-fai*hg/L);elsedisp('后轮先抱死,制动效率为'); E=(a/L)/((1-beta)+fai*hg/L);enddisp('最大制动减速度'),disp(E*fai*g);disp('a)、b)、c)的双回路制动系统制动系增益都为G');disp('1回路失效时，a)、b)、c)的制动系增益分别为,,');disp('1回路失效时，a)、b)、c)的制动系增益分别为,,');disp('a)中1回路失效时最大制动强度');zmax1=a*fai/(L+fai*hg)disp('a)中1回路失效时制动效率');disp(zmax1/fai)disp('a)中1回路失效时最大制动减速度');abmax1=zmax1*gdisp('a)中2回路失效时最大制动强度');zmax2=a*fai/(L+fai*hg)disp('a)中1回路失效时制动效率');disp(zmax1/fai)disp('a)中2回路失效时最大制动减速度');abmax2=zmax2*gdisp('b)或c)中1或2回路失效时最大制动强度');z=fai*b/(2*L*beta-fai*hg)disp('b)或c)中1或2回路失效时最大制动减速度');abmax=z*gdisp('b)或c)中1或2回路失效时制动效率');disp(z/fai)disp('双回路系统a制动系增益最大，一个回路失效时的最大制动减速度也比b、c大，所以其性能较优')disp('b和c制动系增益相同,两者在一个回路失效时的制动效率相同。

汽车理论习题Matlab程序

1.3 确定一轻型货车的动力性能（货车可装用4挡或5挡变速器，任选其中的一种进行整车性能计算）： 1）绘制汽车驱动力与行驶阻力平衡图。

2）求汽车最高车速，最大爬坡度及克服该坡度时相应的附着率。

轻型货车的有关数据：汽油发动机使用外特性的Tq-n 曲线的拟合公式为23419.313295.27()165.44()40.874() 3.8445()1000100010001000q n n n n T =-+-+-式中，Tq 为发动机转矩（N•m ）;n 为发动机转速（r/min ）。

发动机的最低转速n min =600r/min,最高转速n max =4000r/min 。

装载质量 2000kg 整车整备质量 1800kg 总质量 3880kg 车轮半径 0.367m传动系机械效率ηt=0.85滚动阻力系数f=0.013空气阻力系数×迎风面积C D A=2.77m2主减速器传动比i0=5.83飞轮转动惯量I f=0.218kg•m2二前轮转动惯量I w1=1.798kg•m2四后轮转动惯量I w2=3.598kg•m2变速器传动比ig(数据如下表)轴距L=3.2m质心至前轴距离（满载）a=1.974m质心高（满载）hg=0.9m解：Matlab程序：(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/1000) .^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.3006km/h(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/1000) .^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)求最大爬坡度相应的附着率和求汽车行驶加速度倒数曲线程序：clearn=[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.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; 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);ad1=1./a1;a2=(Ft2-Ff-Fw2)/(deta2*m);ad2=1./a2;a3=(Ft3-Ff-Fw3)/(deta3*m);ad3=1./a3;a4=(Ft4-Ff-Fw4)/(deta4*m);ad4=1./a4;a5=(Ft5-Ff-Fw5)/(deta5*m);ad5=1./a5;plot(ua1,ad1,ua2,ad2,ua3,ad3,ua4,ad4,ua5,ad5);axis([0 99 0 10]);title('汽车的加速度倒数曲线');xlabel('ua(km/h)');ylabel('1/a');gtext('1/a1');gtext('1/a2');gtext('1/a3');gtext('1/a4');gtext('1/a5'); a=max(a1);af=asin(max(Ft1-Ff-Fw1)/G);C=tan(af)/(a/L+hg*tan(af)/L);disp('假设后轮驱动，最大爬坡度相应的附着率=');disp(C);假设后轮驱动，最大爬坡度相应的附着率=0.4219(4) >>clearnT=0.85;r=0.367;f=0.013;CDA=2.77;i0=5.83;If=0.218;Iw1=1.798;Iw2=3.598;L=3.2;a=1.947;hg=0.9;m=3880;g=9.8;G=m*g; ig=[5.56 2.769 1.644 1.00 0.793];nmin=600;nmax=4000;u1=0.377*r*nmin./ig/i0;u2=0.377*r*nmax./ig/i0;deta=0*ig;for i=1:5deta(i)=1+(Iw1+Iw2)/(m*r^2)+(If*(ig(i))^2*i0^2*nT)/(m*r^2); endua=[6:0.01:99];N=length(ua);n=0;Tq=0;Ft=0;inv_a=0*ua;delta=0*ua; Ff=G*f;Fw=CDA*ua.^2/21.15;for i=1:Nk=i;if ua(i)<=u2(2)n=ua(i)*(ig(2)*i0/r)/0.377;Tq=-19.313+295.27*(n/1000)-165.44*(n/1000)^2+40.874*(n/1000)^3-3.8445*(n/1000)^4;Ft=Tq*ig(2)*i0*nT/r;inv_a(i)=(deta(2)*m)/(Ft-Ff-Fw(i));delta(i)=0.01*inv_a(i)/3.6;elseif ua(i)<=u2(3)n=ua(i)*(ig(3)*i0/r)/0.377;Tq=-19.313+295.27*(n/1000)-165.44*(n/1000)^2+40.874*(n/1000)^3-3.8445*(n/1000)^4;Ft=Tq*ig(3)*i0*nT/r;inv_a(i)=(deta(3)*m)/(Ft-Ff-Fw(i));delta(i)=0.01*inv_a(i)/3.6;elseif ua(i)<=u2(4)n=ua(i)*(ig(4)*i0/r)/0.377;Tq=-19.313+295.27*(n/1000)-165.44*(n/1000)^2+40.874*(n/1000)^3-3.8445*(n/1000)^4;Ft=Tq*ig(4)*i0*nT/r;inv_a(i)=(deta(4)*m)/(Ft-Ff-Fw(i));delta(i)=0.01*inv_a(i)/3.6;elsen=ua(i)*(ig(5)*i0/r)/0.377;Tq=-19.313+295.27*(n/1000)-165.44*(n/1000)^2+40.874*(n/1000)^3-3.8445*(n/1000)^4;Ft=Tq*ig(5)*i0*nT/r;inv_a(i)=(deta(5)*m)/(Ft-Ff-Fw(i));delta(i)=0.01*inv_a(i)/3.6;enda=delta(1:k);t(i)=sum(a);endplot(t,ua);axis([0 80 0 100]);title('汽车2档原地起步换挡加速时间曲线');xlabel('时间t（s）');ylabel('速度ua（km/h）');>> ginputans =25.8223 70.073725.7467 70.0737所以汽车2档原地起步换挡加速行驶至70km/h的加速时间约为25.8s2.7已知货车装用汽油发动机的负荷特性与万有特性。

汽车理论习题Matlab程序

1.3 确定一轻型货车的动力性能（货车可装用4挡或5挡变速器，任选其中的一种进行整车性能计算）：1）绘制汽车驱动力与行驶阻力平衡图。

2）求汽车最高车速，最大爬坡度及克服该坡度时相应的附着率。

轻型货车的有关数据：汽油发动机使用外特性的Tq-n 曲线的拟合公式为23419.313295.27()165.44()40.874() 3.8445()1000100010001000q n n n n T =-+-+- 式中，Tq 为发动机转矩（N•m ）;n 为发动机转速（r/min ）。

发动机的最低转速n min =600r/min,最高转速n max =4000r/min 。

装载质量 2000kg整车整备质量 1800kg总质量 3880kg车轮半径 0.367m传动系机械效率 ηt =0.85滚动阻力系数 f =0.013空气阻力系数×迎风面积 C D A =2.77m 2主减速器传动比 i 0=5.83飞轮转动惯量 I f =0.218kg•m 2二前轮转动惯量 I w1=1.798kg•m 2四后轮转动惯量 I w2=3.598kg•m 2质心至前轴距离（满载） a=1.974m质心高（满载） hg=0.9m解：Matlab 程序：(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/1000).^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.3006km/h(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)求最大爬坡度相应的附着率和求汽车行驶加速度倒数曲线程序：clearn=[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;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);ad1=1./a1;a2=(Ft2-Ff-Fw2)/(deta2*m);ad2=1./a2;a3=(Ft3-Ff-Fw3)/(deta3*m);ad3=1./a3;a4=(Ft4-Ff-Fw4)/(deta4*m);ad4=1./a4;a5=(Ft5-Ff-Fw5)/(deta5*m);ad5=1./a5;plot(ua1,ad1,ua2,ad2,ua3,ad3,ua4,ad4,ua5,ad5);axis([0 99 0 10]);title('汽车的加速度倒数曲线');xlabel('ua(km/h)');ylabel('1/a');gtext('1/a1');gtext('1/a2');gtext('1/a3');gtext('1/a4');gtext('1/a5'); a=max(a1);af=asin(max(Ft1-Ff-Fw1)/G);C=tan(af)/(a/L+hg*tan(af)/L);disp('假设后轮驱动，最大爬坡度相应的附着率=');disp(C);假设后轮驱动，最大爬坡度相应的附着率=0.4219(4) >>clearnT=0.85;r=0.367;f=0.013;CDA=2.77;i0=5.83;If=0.218;Iw1=1.798;Iw2=3.598;L=3.2;a=1.947;hg=0.9;m=3880;g=9.8;G=m*g; ig=[5.56 2.769 1.644 1.00 0.793];nmin=600;nmax=4000;u1=0.377*r*nmin./ig/i0;u2=0.377*r*nmax./ig/i0;deta=0*ig;for i=1:5deta(i)=1+(Iw1+Iw2)/(m*r^2)+(If*(ig(i))^2*i0^2*nT)/(m*r^2);endua=[6:0.01:99];N=length(ua);n=0;Tq=0;Ft=0;inv_a=0*ua;delta=0*ua;Ff=G*f;Fw=CDA*ua.^2/21.15;for i=1:Nk=i;if ua(i)<=u2(2)n=ua(i)*(ig(2)*i0/r)/0.377;Tq=-19.313+295.27*(n/1000)-165.44*(n/1000)^2+40.874*(n/1000)^3-3.8445*(n/100 0)^4;Ft=Tq*ig(2)*i0*nT/r;inv_a(i)=(deta(2)*m)/(Ft-Ff-Fw(i));delta(i)=0.01*inv_a(i)/3.6;elseif ua(i)<=u2(3)n=ua(i)*(ig(3)*i0/r)/0.377;Tq=-19.313+295.27*(n/1000)-165.44*(n/1000)^2+40.874*(n/1000)^3-3.8445*(n/100 0)^4;Ft=Tq*ig(3)*i0*nT/r;inv_a(i)=(deta(3)*m)/(Ft-Ff-Fw(i));delta(i)=0.01*inv_a(i)/3.6;elseif ua(i)<=u2(4)n=ua(i)*(ig(4)*i0/r)/0.377;Tq=-19.313+295.27*(n/1000)-165.44*(n/1000)^2+40.874*(n/1000)^3-3.8445*(n/100 0)^4;Ft=Tq*ig(4)*i0*nT/r;inv_a(i)=(deta(4)*m)/(Ft-Ff-Fw(i));delta(i)=0.01*inv_a(i)/3.6;elsen=ua(i)*(ig(5)*i0/r)/0.377;Tq=-19.313+295.27*(n/1000)-165.44*(n/1000)^2+40.874*(n/1000)^3-3.8445*(n/1000)^4;Ft=Tq*ig(5)*i0*nT/r;inv_a(i)=(deta(5)*m)/(Ft-Ff-Fw(i));delta(i)=0.01*inv_a(i)/3.6;enda=delta(1:k);t(i)=sum(a);endplot(t,ua);axis([0 80 0 100]);title('汽车2档原地起步换挡加速时间曲线');xlabel('时间t （s ）');ylabel('速度ua （km/h ）');>> ginputans =25.8223 70.073725.7467 70.0737所以汽车2档原地起步换挡加速行驶至70km/h 的加速时间约为25.8s2.7已知货车装用汽油发动机的负荷特性与万有特性。

汽车理论习题Matlab程序

确定一轻型货车的动力性能（货车可装用4挡或5挡变速器，任选其中的一种进行整车性能计算）：1）绘制汽车驱动力与行驶阻力平衡图。

2）求汽车最高车速，最大爬坡度及克服该坡度时相应的附着率。

轻型货车的有关数据：汽油发动机使用外特性的Tq-n 曲线的拟合公式为23419.313295.27()165.44()40.874() 3.8445()1000100010001000q n n n n T =-+-+- 式中，Tq 为发动机转矩（N•m）;n 为发动机转速（r/min ）。

发动机的最低转速n min =600r/min,最高转速n max =4000r/min 。

装载质量 2000kg整车整备质量 1800kg总质量 3880kg车轮半径传动系机械效率 ηt =滚动阻力系数 f =空气阻力系数×迎风面积 C D A =主减速器传动比 i 0=飞轮转动惯量 I f =•m 2二前轮转动惯量 I w1=•m 2四后轮转动惯量 I w2=•m 2变速器传动比 ig(数据如下表)轴距 L=质心至前轴距离（满载） a=质心高（满载） hg=解：Matlab程序：(1) 求汽车驱动力与行驶阻力平衡图和汽车最高车速程序： n=[600:10:4000];Tq=+*(n/1000)*(n/1000).^2+*(n/1000).^*(n/1000).^4; m=3880;g=;nmin=600;nmax=4000;G=m*g;ig=[ ];nT=;r=;f=;CDA=;i0=;L=;a=;hg=;If=;Iw1=;Iw2=;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=*r*n/ig(1)/i0;ua2=*r*n/ig(2)/i0;ua3=*r*n/ig(3)/i0;ua4=*r*n/ig(4)/i0;ua5=*r*n/ig(5)/i0;ua=[0:5:120];Ff=G*f;Fw=CDA*ua.^2/;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'),gtex t('Ff+Fw');zoom on;[x,y]=ginput(1);zoom off;disp('汽车最高车速=');disp(x);disp('km/h');汽车最高车速=km/h(2)求汽车最大爬坡度程序：n=[600:10:4000];Tq=+*(n/1000)*(n/1000).^2+*(n/1000).^*(n/1000).^4;m=3880;g=;nmin=600;nmax=4000;G=m*g;ig=[ ];nT=;r=;f=;CDA=;i0=;L=;a=;hg=;If=;Iw1=;Iw2=;Ft1=Tq*ig(1)*i0*nT/r;ua1=*r*n/ig(1)/i0;Ff=G*f;Fw1=CDA*ua1.^2/;Fz1=Ff+Fw1;Fi1=Ft1-Fz1;Zoom on;imax=100*tan(asin(max(Fi1/G)));disp('汽车最大爬坡度=');disp(imax);disp('%');汽车最大爬坡度=%(3)求最大爬坡度相应的附着率和求汽车行驶加速度倒数曲线程序：clearn=[600:10:4000];Tq=+*(n/1000)*(n/1000).^2+*(n/1000).^*(n/1000).^4;m=3880;g=;nmin=600;nmax=4000;G=m*g;ig=[ ];nT=;r=;f=;CDA=;i0=;L=;a=;hg=;If=;Iw1=;Iw2=;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=*r*n/ig(1)/i0;ua2=*r*n/ig(2)/i0;ua3=*r*n/ig(3)/i0;ua4=*r*n/ig(4)/i0;ua5=*r*n/ig(5)/i0;Fw1=CDA*ua1.^2/;Fw2=CDA*ua2.^2/;Fw3=CDA*ua3.^2/;Fw4=CDA*ua4.^2/;Fw5=CDA*ua5.^2/;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);ad1=1./a1;a2=(Ft2-Ff-Fw2)/(deta2*m);ad2=1./a2;a3=(Ft3-Ff-Fw3)/(deta3*m);ad3=1./a3;a4=(Ft4-Ff-Fw4)/(deta4*m);ad4=1./a4;a5=(Ft5-Ff-Fw5)/(deta5*m);ad5=1./a5;plot(ua1,ad1,ua2,ad2,ua3,ad3,ua4,ad4,ua5,ad5);axis([0 99 0 10]);title('汽车的加速度倒数曲线');xlabel('ua(km/h)');ylabel('1/a');gtext('1/a1');gtext('1/a2');gtext('1/a3');gtext('1/a4');gtext('1/a5') ;a=max(a1);af=asin(max(Ft1-Ff-Fw1)/G);C=tan(af)/(a/L+hg*tan(af)/L);disp('假设后轮驱动，最大爬坡度相应的附着率=');disp(C);假设后轮驱动，最大爬坡度相应的附着率=(4) >>clearnT=;r=;f=;CDA=;i0=;If=;Iw1=;Iw2=;L=;a=;hg=;m=3880;g=;G=m*g; ig=[ ];nmin=600;nmax=4000;u1=*r*nmin./ig/i0;u2=*r*nmax./ig/i0;deta=0*ig;for i=1:5deta(i)=1+(Iw1+Iw2)/(m*r^2)+(If*(ig(i))^2*i0^2*nT)/(m*r^2);endua=[6::99];N=length(ua);n=0;Tq=0;Ft=0;inv_a=0*ua;delta=0*ua;Ff=G*f;Fw=CDA*ua.^2/;for i=1:Nk=i;if ua(i)<=u2(2)n=ua(i)*(ig(2)*i0/r)/;Tq=+*(n/1000)*(n/1000)^2+*(n/1000)^*(n/1000)^4; Ft=Tq*ig(2)*i0*nT/r;inv_a(i)=(deta(2)*m)/(Ft-Ff-Fw(i));delta(i)=*inv_a(i)/;elseif ua(i)<=u2(3)n=ua(i)*(ig(3)*i0/r)/;Tq=+*(n/1000)*(n/1000)^2+*(n/1000)^*(n/1000)^4; Ft=Tq*ig(3)*i0*nT/r;inv_a(i)=(deta(3)*m)/(Ft-Ff-Fw(i));delta(i)=*inv_a(i)/;elseif ua(i)<=u2(4)n=ua(i)*(ig(4)*i0/r)/;Tq=+*(n/1000)*(n/1000)^2+*(n/1000)^*(n/1000)^4; Ft=Tq*ig(4)*i0*nT/r;inv_a(i)=(deta(4)*m)/(Ft-Ff-Fw(i));delta(i)=*inv_a(i)/;elsen=ua(i)*(ig(5)*i0/r)/;Tq=+*(n/1000)*(n/1000)^2+*(n/1000)^*(n/1000)^4;Ft=Tq*ig(5)*i0*nT/r;inv_a(i)=(deta(5)*m)/(Ft-Ff-Fw(i));delta(i)=*inv_a(i)/;enda=delta(1:k);t(i)=sum(a);endplot(t,ua);axis([0 80 0 100]);title('汽车2档原地起步换挡加速时间曲线');xlabel('时间t （s ）');ylabel('速度ua （km/h ）');>> ginputans =所以汽车2档原地起步换挡加速行驶至70km/h 的加速时间约为已知货车装用汽油发动机的负荷特性与万有特性。

汽车理论习题Matlab程序

1.3 确定一轻型货车的动力性能〔货车可装用4挡或5挡变速器，任选其中的一种进行整车性能计算〕：1〕绘制汽车驱动力与行驶阻力平衡图。

2〕求汽车最高车速，最大爬坡度及克服该坡度时相应的附着率。

3〕绘制汽车行驶加速度倒数曲线，用图解积分法求汽车用2档起步加速行驶至70km/h 的车速－时间曲线，或者用计算机求汽车用2档起步加速行驶至70km/h 的加速时间。

轻型货车的有关数据：汽油发动机使用外特性的Tq-n 曲线的拟合公式为23419.313295.27()165.44()40.874() 3.8445()1000100010001000q n n n n T =-+-+- 式中，Tq 为发动机转矩〔N•m 〕;n 为发动机转速〔r/min 〕。

发动机的最低转速n min =600r/min,最高转速n max =4000r/min 。

装载质量 2000kg整车整备质量 1800kg总质量 3880kg传动系机械效率 ηt滚动阻力系数 f空气阻力系数×迎风面积 C D A 2主减速器传动比 i 0飞轮转动惯量 I f =0.218kg•m 2二前轮转动惯量 I w1=1.798kg•m 2四后轮转动惯量 I w2=3.598kg•m 2变速器传动比 ig(数据如下表)解：Matlab程序：(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');汽车最高车速=km/h(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)求最大爬坡度相应的附着率和求汽车行驶加速度倒数曲线程序：clearn=[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;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);ad1=1./a1;a2=(Ft2-Ff-Fw2)/(deta2*m);ad2=1./a2;a3=(Ft3-Ff-Fw3)/(deta3*m);ad3=1./a3;a4=(Ft4-Ff-Fw4)/(deta4*m);ad4=1./a4;a5=(Ft5-Ff-Fw5)/(deta5*m);ad5=1./a5;plot(ua1,ad1,ua2,ad2,ua3,ad3,ua4,ad4,ua5,ad5);axis([0 99 0 10]);title('汽车的加速度倒数曲线');xlabel('ua(km/h)');ylabel('1/a');gtext('1/a1');gtext('1/a2');gtext('1/a3');gtext('1/a4');gtext('1/a5'); a=max(a1);af=asin(max(Ft1-Ff-Fw1)/G);C=tan(af)/(a/L+hg*tan(af)/L);disp('假设后轮驱动，最大爬坡度相应的附着率=');disp(C);假设后轮驱动，最大爬坡度相应的附着率=(4) >>clearnT=0.85;r=0.367;f=0.013;CDA=2.77;i0=5.83;If=0.218;Iw1=1.798;Iw2=3.598;L=3.2;a=1.947;hg=0.9;m=3880;g=9.8; G=m*g; ig=[5.56 2.769 1.644 1.00 0.793];nmin=600;nmax=4000;u1=0.377*r*nmin./ig/i0;u2=0.377*r*nmax./ig/i0;deta=0*ig;deta(i)=1+(Iw1+Iw2)/(m*r^2)+(If*(ig(i))^2*i0^2*nT)/(m*r^2);endua=[6:0.01:99];N=length(ua);n=0;Tq=0;Ft=0;inv_a=0*ua;delta=0*ua;Ff=G*f;Fw=CDA*ua.^2/21.15;for i=1:Nk=i;if ua(i)<=u2(2)n=ua(i)*(ig(2)*i0/r)/0.377;Tq=-19.313+295.27*(n/1000)-165.44*(n/1000)^2+40.874*(n/1000)^3-3.8445*(n/100 0)^4;Ft=Tq*ig(2)*i0*nT/r;inv_a(i)=(deta(2)*m)/(Ft-Ff-Fw(i));delta(i)=0.01*inv_a(i)/3.6;elseif ua(i)<=u2(3)n=ua(i)*(ig(3)*i0/r)/0.377;Tq=-19.313+295.27*(n/1000)-165.44*(n/1000)^2+40.874*(n/1000)^3-3.8445*(n/100 0)^4;Ft=Tq*ig(3)*i0*nT/r;inv_a(i)=(deta(3)*m)/(Ft-Ff-Fw(i));delta(i)=0.01*inv_a(i)/3.6;elseif ua(i)<=u2(4)n=ua(i)*(ig(4)*i0/r)/0.377;Tq=-19.313+295.27*(n/1000)-165.44*(n/1000)^2+40.874*(n/1000)^3-3.8445*(n/100 0)^4;Ft=Tq*ig(4)*i0*nT/r;inv_a(i)=(deta(4)*m)/(Ft-Ff-Fw(i));delta(i)=0.01*inv_a(i)/3.6;elsen=ua(i)*(ig(5)*i0/r)/0.377;Tq=-19.313+295.27*(n/1000)-165.44*(n/1000)^2+40.874*(n/1000)^3-3.8445*(n/100 0)^4;Ft=Tq*ig(5)*i0*nT/r;inv_a(i)=(deta(5)*m)/(Ft-Ff-Fw(i));delta(i)=0.01*inv_a(i)/3.6;enda=delta(1:k);t(i)=sum(a);endplot(t,ua);axis([0 80 0 100]);title('汽车2档原地起步换挡加速时间曲线');xlabel('时间t〔s〕');ylabel('速度ua〔km/h〕');>> ginputans =所以汽车2档原地起步换挡加速行驶至70km/h 的加速时间约为s货车装用汽油发动机的负荷特性与万有特性。

汽车理论习题Matlab程序68411

1.3 确定一轻型货车的动力性能（货车可装用4挡或5挡变速器，任选其中的一种进行整车性能计算）：1）绘制汽车驱动力与行驶阻力平衡图。

2）求汽车最高车速，最大爬坡度及克服该坡度时相应的附着率。

轻型货车的有关数据：汽油发动机使用外特性的Tq-n 曲线的拟合公式为23419.313295.27()165.44()40.874() 3.8445()1000100010001000q n n n n T =-+-+-式中，Tq 为发动机转矩（N •m ）;n 为发动机转速（r/min ）。

发动机的最低转速n min =600r/min,最高转速n max =4000r/min 。

装载质量 2000kg 整车整备质量 1800kg 总质量 3880kg 车轮半径 0.367m 传动系机械效率 ηt =0.85 滚动阻力系数 f =0.013 空气阻力系数×迎风面积 C D A =2.77m 2 主减速器传动比 i 0=5.83 飞轮转动惯量 I f =0.218kg •m 2 二前轮转动惯量 I w1=1.798kg •m 2四后轮转动惯量I w2=3.598kg•m2变速器传动比ig(数据如下表)轴距L=3.2m质心至前轴距离（满载）a=1.974m质心高（满载）hg=0.9m解：Matlab程序：(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.3006km/h(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)求最大爬坡度相应的附着率和求汽车行驶加速度倒数曲线程序：clearn=[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;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);ad1=1./a1;a2=(Ft2-Ff-Fw2)/(deta2*m);ad2=1./a2;a3=(Ft3-Ff-Fw3)/(deta3*m);ad3=1./a3;a4=(Ft4-Ff-Fw4)/(deta4*m);ad4=1./a4;a5=(Ft5-Ff-Fw5)/(deta5*m);ad5=1./a5;plot(ua1,ad1,ua2,ad2,ua3,ad3,ua4,ad4,ua5,ad5);axis([0 99 0 10]);title('汽车的加速度倒数曲线');xlabel('ua(km/h)');ylabel('1/a');gtext('1/a1');gtext('1/a2');gtext('1/a3');gtext('1/a4');gtext('1/a5'); a=max(a1);af=asin(max(Ft1-Ff-Fw1)/G);C=tan(af)/(a/L+hg*tan(af)/L);disp('假设后轮驱动，最大爬坡度相应的附着率=');disp(C);假设后轮驱动，最大爬坡度相应的附着率=0.4219(4)>>clearnT=0.85;r=0.367;f=0.013;CDA=2.77;i0=5.83;If=0.218;Iw1=1.798;Iw2=3.598;L=3.2;a=1.947;hg=0.9;m=3880;g=9.8;G=m*g; ig=[5.56 2.769 1.644 1.00 0.793];nmin=600;nmax=4000;u1=0.377*r*nmin./ig/i0;u2=0.377*r*nmax./ig/i0;deta=0*ig;for i=1:5deta(i)=1+(Iw1+Iw2)/(m*r^2)+(If*(ig(i))^2*i0^2*nT)/(m*r^2); endua=[6:0.01:99];N=length(ua);n=0;Tq=0;Ft=0;inv_a=0*ua;delta=0*ua;Fw=CDA*ua.^2/21.15;for i=1:Nk=i;if ua(i)<=u2(2)n=ua(i)*(ig(2)*i0/r)/0.377;Tq=-19.313+295.27*(n/1000)-165.44*(n/1000)^2+40.874*(n/1000)^3-3.8445*(n/100 0)^4;Ft=Tq*ig(2)*i0*nT/r;inv_a(i)=(deta(2)*m)/(Ft-Ff-Fw(i));delta(i)=0.01*inv_a(i)/3.6;elseif ua(i)<=u2(3)n=ua(i)*(ig(3)*i0/r)/0.377;Tq=-19.313+295.27*(n/1000)-165.44*(n/1000)^2+40.874*(n/1000)^3-3.8445*(n/100 0)^4;Ft=Tq*ig(3)*i0*nT/r;inv_a(i)=(deta(3)*m)/(Ft-Ff-Fw(i));delta(i)=0.01*inv_a(i)/3.6;elseif ua(i)<=u2(4)n=ua(i)*(ig(4)*i0/r)/0.377;Tq=-19.313+295.27*(n/1000)-165.44*(n/1000)^2+40.874*(n/1000)^3-3.8445*(n/100 0)^4;Ft=Tq*ig(4)*i0*nT/r;inv_a(i)=(deta(4)*m)/(Ft-Ff-Fw(i));delta(i)=0.01*inv_a(i)/3.6;elsen=ua(i)*(ig(5)*i0/r)/0.377;Tq=-19.313+295.27*(n/1000)-165.44*(n/1000)^2+40.874*(n/1000)^3-3.8445*(n/100 0)^4;Ft=Tq*ig(5)*i0*nT/r;inv_a(i)=(deta(5)*m)/(Ft-Ff-Fw(i));delta(i)=0.01*inv_a(i)/3.6;enda=delta(1:k);t(i)=sum(a);endplot(t,ua);axis([0 80 0 100]);title('汽车2档原地起步换挡加速时间曲线');xlabel('时间t（s）');ylabel('速度ua （km/h ）'); >> ginput ans =25.8223 70.0737 25.7467 70.0737所以汽车2档原地起步换挡加速行驶至70km/h 的加速时间约为25.8s2.7已知货车装用汽油发动机的负荷特性与万有特性。

汽车理论课后作业matlab编程详解(带注释)

1.3matlab程序：（1）%驱动力-行驶阻力平衡图%货车相关参数。

m=3880;g=9.8;nmin=600;nmax=4000;G=m*g;ig=[5.56 2.769 1.644 1.000.793];y=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;n=600:10:4000;%发动机转数转换成汽车行驶速度。

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;%计算各档位驱动力。

Tq=-19.313+295.27*(n/1000)-165.44*(n/1000).^2+40.874*(n/1000).^3-3.84 45*(n/1000).^4;Ft1=Tq*ig(1)*i0*y/r;Ft2=Tq*ig(2)*i0*y/r;Ft3=Tq*ig(3)*i0*y/r;Ft4=Tq*ig(4)*i0*y/r;Ft5=Tq*ig(5)*i0*y/r;%计算行驶阻力。

Fz1=m*g*f+2.77*ua1.^2/21.15;Fz2=m*g*f+2.77*ua2.^2/21.15;Fz3=m*g*f+2.77*ua3.^2/21.15;Fz4=m*g*f+2.77*ua4.^2/21.15;Fz5=m*g*f+2.77*ua5.^2/21.15;%驱动力-行驶阻力平衡图。

plot(ua1,Ft1,ua2,Ft2,ua3,Ft3,ua4,Ft4,ua5,Ft5,ua1,Fz1,'k',ua2,Fz2,'k', ua3,Fz3,'k',ua4,Fz4,'k',ua5,Fz5,'k');title('驱动力-行驶阻力平衡图');xlabel('ua(km/s)');ylabel('Ft(N)');legend('Ft1','Ft2','Ft3','Ft4','Ft5','Ff+Fw')%求最高车速。

汽车理论课后习题MATLAB编程-武汉理工版

汽车理论课后习题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.3006km/h(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)求最大爬坡度相应的附着率和求汽车行驶加速度倒数曲线程序：clearn=[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;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);ad1=1./a1;a2=(Ft2-Ff-Fw2)/(deta2*m);ad2=1./a2;a3=(Ft3-Ff-Fw3)/(deta3*m);ad3=1./a3;a4=(Ft4-Ff-Fw4)/(deta4*m);ad4=1./a4;a5=(Ft5-Ff-Fw5)/(deta5*m);ad5=1./a5;plot(ua1,ad1,ua2,ad2,ua3,ad3,ua4,ad4,ua5,ad5);axis([0 99 0 10]);title('汽车的加速度倒数曲线');xlabel('ua(km/h)');ylabel('1/a');gtext('1/a1');gtext('1/a2');gtext('1/a3');gtext('1/a4');gtext('1/a5');a=max(a1);af=asin(max(Ft1-Ff-Fw1)/G);C=tan(af)/(a/L+hg*tan(af)/L);disp('假设后轮驱动，最大爬坡度相应的附着率=');disp(C);假设后轮驱动，最大爬坡度相应的附着率=0.4219(4) >>clearnT=0.85;r=0.367;f=0.013;CDA=2.77;i0=5.83;If=0.218;Iw1=1.798;Iw2=3.598;L=3.2;a=1.947;hg=0.9;m=3880;g=9.8;G=m*g; ig=[5.56 2.769 1.644 1.00 0.793];nmin=600;nmax=4000;u1=0.377*r*nmin./ig/i0;u2=0.377*r*nmax./ig/i0;deta=0*ig;for i=1:5deta(i)=1+(Iw1+Iw2)/(m*r^2)+(If*(ig(i))^2*i0^2*nT)/(m*r^2);endua=[6:0.01:99];N=length(ua);n=0;Tq=0;Ft=0;inv_a=0*ua;delta=0*ua;Ff=G*f;Fw=CDA*ua.^2/21.15;for i=1:Nk=i;if ua(i)<=u2(2)n=ua(i)*(ig(2)*i0/r)/0.377;Tq=-19.313+295.27*(n/1000)-165.44*(n/1000)^2+40.874*(n/1000)^3-3.8445*(n/100 0)^4;Ft=Tq*ig(2)*i0*nT/r;inv_a(i)=(deta(2)*m)/(Ft-Ff-Fw(i));delta(i)=0.01*inv_a(i)/3.6;elseif ua(i)<=u2(3)n=ua(i)*(ig(3)*i0/r)/0.377;Tq=-19.313+295.27*(n/1000)-165.44*(n/1000)^2+40.874*(n/1000)^3-3.8445*(n/100 0)^4;Ft=Tq*ig(3)*i0*nT/r;inv_a(i)=(deta(3)*m)/(Ft-Ff-Fw(i));delta(i)=0.01*inv_a(i)/3.6;elseif ua(i)<=u2(4)n=ua(i)*(ig(4)*i0/r)/0.377;Tq=-19.313+295.27*(n/1000)-165.44*(n/1000)^2+40.874*(n/1000)^3-3.8445*(n/100 0)^4;Ft=Tq*ig(4)*i0*nT/r;inv_a(i)=(deta(4)*m)/(Ft-Ff-Fw(i));delta(i)=0.01*inv_a(i)/3.6;elsen=ua(i)*(ig(5)*i0/r)/0.377;Tq=-19.313+295.27*(n/1000)-165.44*(n/1000)^2+40.874*(n/1000)^3-3.8445*(n/100 0)^4;Ft=Tq*ig(5)*i0*nT/r;inv_a(i)=(deta(5)*m)/(Ft-Ff-Fw(i));delta(i)=0.01*inv_a(i)/3.6;enda=delta(1:k);t(i)=sum(a);endplot(t,ua);axis([0 80 0 100]);title('汽车2档原地起步换挡加速时间曲线');xlabel('时间t（s）');ylabel('速度ua（km/h）');>> ginputans =25.8223 70.073725.7467 70.0737所以汽车2档原地起步换挡加速行驶至70km/h的加速时间约为25.8s2.7已知货车装用汽油发动机的负荷特性与万有特性。

汽车理论习题Matlab程序

1.3 确定一轻型货车的动力性能（货车可装用4挡或5挡变速器，任选其中的一种进行整车性能计算）：1）绘制汽车驱动力与行驶阻力平衡图。

2）求汽车最高车速，最大爬坡度及克服该坡度时相应的附着率。

轻型货车的有关数据：汽油发动机使用外特性的Tq-n 曲线的拟合公式为23419.313295.27()165.44()40.874() 3.8445()1000100010001000q n n n n T =-+-+-式中，Tq 为发动机转矩（N•m ）;n 为发动机转速（r/min ）。

发动机的最低转速n min =600r/min,最高转速n max =4000r/min 。

装载质量 2000kg 整车整备质量 1800kg 总质量 3880kg 车轮半径 0.367m 传动系机械效率 ηt =0.85 滚动阻力系数 f =0.013 空气阻力系数×迎风面积 C D A =2.77m 2 主减速器传动比 i 0=5.83飞轮转动惯量 I f =0.218kg•m 2 二前轮转动惯量 I w1=1.798kg•m 2 四后轮转动惯量 I w2=3.598kg•m 2质心至前轴距离（满载） a=1.974m 质心高（满载） hg=0.9m解：Matlab 程序：(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/1000).^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.3006km/h- (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)求最大爬坡度相应的附着率和求汽车行驶加速度倒数曲线程序：clearn=[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;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);ad1=1./a1;a2=(Ft2-Ff-Fw2)/(deta2*m);ad2=1./a2;a3=(Ft3-Ff-Fw3)/(deta3*m);ad3=1./a3;a4=(Ft4-Ff-Fw4)/(deta4*m);ad4=1./a4;a5=(Ft5-Ff-Fw5)/(deta5*m);ad5=1./a5;plot(ua1,ad1,ua2,ad2,ua3,ad3,ua4,ad4,ua5,ad5);axis([0 99 0 10]);title('汽车的加速度倒数曲线');xlabel('ua(km/h)');ylabel('1/a');gtext('1/a1');gtext('1/a2');gtext('1/a3');gtext('1/a4');gtext('1/a5');a=max(a1);af=asin(max(Ft1-Ff-Fw1)/G);C=tan(af)/(a/L+hg*tan(af)/L);disp('假设后轮驱动，最大爬坡度相应的附着率=');disp(C);假设后轮驱动，最大爬坡度相应的附着率=0.4219(4) >>clearnT=0.85;r=0.367;f=0.013;CDA=2.77;i0=5.83;If=0.218;Iw1=1.798;Iw2=3.598;L=3.2;a=1.947;hg=0.9;m=3880;g=9.8;G=m*g; ig=[5.56 2.769 1.644 1.00 0.793];nmin=600;nmax=4000;u1=0.377*r*nmin./ig/i0;u2=0.377*r*nmax./ig/i0;deta=0*ig;for i=1:5deta(i)=1+(Iw1+Iw2)/(m*r^2)+(If*(ig(i))^2*i0^2*nT)/(m*r^2);endua=[6:0.01:99];N=length(ua);n=0;Tq=0;Ft=0;inv_a=0*ua;delta=0*ua;Ff=G*f;Fw=CDA*ua.^2/21.15;for i=1:Nk=i;if ua(i)<=u2(2)n=ua(i)*(ig(2)*i0/r)/0.377;Tq=-19.313+295.27*(n/1000)-165.44*(n/1000)^2+40.874*(n/1000)^3-3.8445*(n/100 0)^4;Ft=Tq*ig(2)*i0*nT/r;inv_a(i)=(deta(2)*m)/(Ft-Ff-Fw(i));delta(i)=0.01*inv_a(i)/3.6;elseif ua(i)<=u2(3)n=ua(i)*(ig(3)*i0/r)/0.377;Tq=-19.313+295.27*(n/1000)-165.44*(n/1000)^2+40.874*(n/1000)^3-3.8445*(n/100 0)^4;Ft=Tq*ig(3)*i0*nT/r;inv_a(i)=(deta(3)*m)/(Ft-Ff-Fw(i));delta(i)=0.01*inv_a(i)/3.6;elseif ua(i)<=u2(4)n=ua(i)*(ig(4)*i0/r)/0.377;Tq=-19.313+295.27*(n/1000)-165.44*(n/1000)^2+40.874*(n/1000)^3-3.8445*(n/100 0)^4;Ft=Tq*ig(4)*i0*nT/r;inv_a(i)=(deta(4)*m)/(Ft-Ff-Fw(i));delta(i)=0.01*inv_a(i)/3.6;elsen=ua(i)*(ig(5)*i0/r)/0.377;Tq=-19.313+295.27*(n/1000)-165.44*(n/1000)^2+40.874*(n/1000)^3-3.8445*(n/1000)^4;Ft=Tq*ig(5)*i0*nT/r;inv_a(i)=(deta(5)*m)/(Ft-Ff-Fw(i)); delta(i)=0.01*inv_a(i)/3.6; enda=delta(1:k); t(i)=sum(a); endplot(t,ua);axis([0 80 0 100]);title('汽车2档原地起步换挡加速时间曲线'); xlabel('时间t （s ）');ylabel('速度ua （km/h ）'); >> ginputans =25.8223 70.0737 25.7467 70.0737所以汽车2档原地起步换挡加速行驶至70km/h 的加速时间约为25.8s2.7已知货车装用汽油发动机的负荷特性与万有特性。

汽车理论习题Maab程序

确定一轻型货车的动力性能（货车可装用4挡或5挡变速器，任选其中的一种进行整车性能计算）：1）绘制汽车驱动力与行驶阻力平衡图。

2）求汽车最高车速，最大爬坡度及克服该坡度时相应的附着率。

轻型货车的有关数据：汽油发动机使用外特性的Tq-n曲线的拟合公式为式中，Tq为发动机转矩（N?m）;n为发动机转速（r/min）。

发动机的最低转速nmin =600r/min,最高转速nmax=4000r/min。

装载质量 2000kg整车整备质量 1800kg总质量 3880kg车轮半径 0.367m传动系机械效率ηt=0.85滚动阻力系数f=0.013空气阻力系数×迎风面积C D A=2.77m2主减速器传动比i0=5.83飞轮转动惯量 I f=0.218kg?m2二前轮转动惯量I w1=1.798kg?m2四后轮转动惯量I w2=3.598kg?m2变速器传动比 ig(数据如下表)轴距 L=3.2m质心至前轴距离（满载） a=1.974m质心高（满载） hg=0.9m解：Matlab程序：(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/100 0).^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.3006km/h(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/1000).^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)求最大爬坡度相应的附着率和求汽车行驶加速度倒数曲线程序：clearn=[600:10:4000];Tq=-19.313+295.27*(n/1000)-165.44*(n/1000).^2+40.874*(n/1000).^3-3.8445*(n/100 0).^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;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);ad1=1./a1;a2=(Ft2-Ff-Fw2)/(deta2*m);ad2=1./a2;a3=(Ft3-Ff-Fw3)/(deta3*m);ad3=1./a3;a4=(Ft4-Ff-Fw4)/(deta4*m);ad4=1./a4;a5=(Ft5-Ff-Fw5)/(deta5*m);ad5=1./a5;plot(ua1,ad1,ua2,ad2,ua3,ad3,ua4,ad4,ua5,ad5);axis([0 99 0 10]);title('汽车的加速度倒数曲线');xlabel('ua(km/h)');ylabel('1/a');gtext('1/a1');gtext('1/a2');gtext('1/a3');gtext('1/a4');gtext('1/a5'); a=max(a1);af=asin(max(Ft1-Ff-Fw1)/G);C=tan(af)/(a/L+hg*tan(af)/L);disp('假设后轮驱动，最大爬坡度相应的附着率=');disp(C);假设后轮驱动，最大爬坡度相应的附着率=0.4219(4) >>clearnT=0.85;r=0.367;f=0.013;CDA=2.77;i0=5.83;If=0.218;Iw1=1.798;Iw2=3.598;L=3.2;a=1.947;hg=0.9;m=3880;g=9.8; G=m*g; ig=[5.56 2.769 1.644 1.00 0.793];nmin=600;nmax=4000;u1=0.377*r*nmin./ig/i0;u2=0.377*r*nmax./ig/i0;deta=0*ig;for i=1:5deta(i)=1+(Iw1+Iw2)/(m*r^2)+(If*(ig(i))^2*i0^2*nT)/(m*r^2);endua=[6:0.01:99];N=length(ua);n=0;Tq=0;Ft=0;inv_a=0*ua;delta=0*ua;Ff=G*f;Fw=CDA*ua.^2/21.15;for i=1:Nk=i;if ua(i)<=u2(2)n=ua(i)*(ig(2)*i0/r)/0.377;Tq=-19.313+295.27*(n/1000)-165.44*(n/1000)^2+40.874*(n/1000)^3-3.8445*(n/1000) ^4;Ft=Tq*ig(2)*i0*nT/r;inv_a(i)=(deta(2)*m)/(Ft-Ff-Fw(i));delta(i)=0.01*inv_a(i)/3.6;elseif ua(i)<=u2(3)n=ua(i)*(ig(3)*i0/r)/0.377;Tq=-19.313+295.27*(n/1000)-165.44*(n/1000)^2+40.874*(n/1000)^3-3.8445*(n/1000) ^4;Ft=Tq*ig(3)*i0*nT/r;inv_a(i)=(deta(3)*m)/(Ft-Ff-Fw(i));delta(i)=0.01*inv_a(i)/3.6;elseif ua(i)<=u2(4)n=ua(i)*(ig(4)*i0/r)/0.377;Tq=-19.313+295.27*(n/1000)-165.44*(n/1000)^2+40.874*(n/1000)^3-3.8445*(n/1000) ^4;Ft=Tq*ig(4)*i0*nT/r;inv_a(i)=(deta(4)*m)/(Ft-Ff-Fw(i));delta(i)=0.01*inv_a(i)/3.6;elsen=ua(i)*(ig(5)*i0/r)/0.377;Tq=-19.313+295.27*(n/1000)-165.44*(n/1000)^2+40.874*(n/1000)^3-3.8445*(n/1000)^4;Ft=Tq*ig(5)*i0*nT/r;inv_a(i)=(deta(5)*m)/(Ft-Ff-Fw(i)); delta(i)=0.01*inv_a(i)/3.6;enda=delta(1:k);t(i)=sum(a);endplot(t,ua);axis([0 80 0 100]);title('汽车2档原地起步换挡加速时间曲线'); xlabel('时间t（s）');ylabel('速度ua（km/h）');>> ginputans =25.8223 70.073725.7467 70.0737所以汽车2档原地起步换挡加速行驶至70km/h 的加速时间约为25.8s2.7已知货车装用汽油发动机的负荷特性与万有特性。

汽车理论习题Matlab程序

确信一轻型货车的动力性能（货车可装用4挡或5挡变速器，任选其中的一种进行整车性能计算）：1）绘制汽车驱动力与行驶阻力平稳图。

2）求汽车最高车速，最大爬坡度及克服该坡度时相应的附着率。

3）绘制汽车行驶加速度倒数曲线，用图解积分法求汽车用2档起步加速行驶至70km/h 的车速－时刻曲线，或用运算机求汽车用2档起步加速行驶至70km/h 的加速时刻。

轻型货车的有关数据：汽油发动机利用外特性的Tq-n 曲线的拟合公式为23419.313295.27()165.44()40.874() 3.8445()1000100010001000q n n n n T =-+-+-式中，Tq 为发动机转矩（N•m）;n 为发动机转速（r/min ）。

发动机的最低转速n min =600r/min,最高转速n max =4000r/min 。

装载质量 2000kg 整车整备质量 1800kg 总质量 3880kg 车轮半径传动系机械效率 ηt = 转动阻力系数 f = 空气阻力系数×迎风面积 C D A = 主减速器传动比 i 0= 飞轮转动惯量 I f =•m 2 二前轮转动惯量 I w1=•m 2 四后轮转动惯量 I w2=•m 2变速器传动比 ig(数据如下表)轴距 L=质心至前轴距离（满载） a=质心高（满载） hg=解：Matlab程序：(1) 求汽车驱动力与行驶阻力平稳图和汽车最高车速程序： n=[600:10:4000];Tq=+*(n/1000)*(n/1000).^2+*(n/1000).^*(n/1000).^4; m=3880;g=;nmin=600;nmax=4000;G=m*g;ig=[ ];nT=;r=;f=;CDA=;i0=;L=;a=;hg=;If=;Iw1=;Iw2=;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=*r*n/ig(1)/i0;ua2=*r*n/ig(2)/i0;ua3=*r*n/ig(3)/i0;ua4=*r*n/ig(4)/i0;ua5=*r*n/ig(5)/i0;ua=[0:5:120];Ff=G*f;Fw=CDA*ua.^2/;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');汽车最高车速=km/h(2)求汽车最大爬坡度程序：n=[600:10:4000];Tq=+*(n/1000)*(n/1000).^2+*(n/1000).^*(n/1000).^4;m=3880;g=;nmin=600;nmax=4000;G=m*g;ig=[ ];nT=;r=;f=;CDA=;i0=;L=;a=;hg=;If=;Iw1=;Iw2=;Ft1=Tq*ig(1)*i0*nT/r;ua1=*r*n/ig(1)/i0;Ff=G*f;Fw1=CDA*ua1.^2/;Fz1=Ff+Fw1;Fi1=Ft1-Fz1;Zoom on;imax=100*tan(asin(max(Fi1/G)));disp('汽车最大爬坡度=');disp(imax);disp('%');汽车最大爬坡度=%(3)求最大爬坡度相应的附着率和求汽车行驶加速度倒数曲线程序：clearn=[600:10:4000];Tq=+*(n/1000)*(n/1000).^2+*(n/1000).^*(n/1000).^4;m=3880;g=;nmin=600;nmax=4000;G=m*g;ig=[ ];nT=;r=;f=;CDA=;i0=;L=;a=;hg=;If=;Iw1=;Iw2=;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=*r*n/ig(1)/i0;ua2=*r*n/ig(2)/i0;ua3=*r*n/ig(3)/i0;ua4=*r*n/ig(4)/i0;ua5=*r*n/ig(5)/i0;Fw1=CDA*ua1.^2/;Fw2=CDA*ua2.^2/;Fw3=CDA*ua3.^2/;Fw4=CDA*ua4.^2/;Fw5=CDA*ua5.^2/;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);ad1=1./a1;a2=(Ft2-Ff-Fw2)/(deta2*m);ad2=1./a2;a3=(Ft3-Ff-Fw3)/(deta3*m);ad3=1./a3;a4=(Ft4-Ff-Fw4)/(deta4*m);ad4=1./a4;a5=(Ft5-Ff-Fw5)/(deta5*m);ad5=1./a5;plot(ua1,ad1,ua2,ad2,ua3,ad3,ua4,ad4,ua5,ad5);axis([0 99 0 10]);title('汽车的加速度倒数曲线');xlabel('ua(km/h)');ylabel('1/a');gtext('1/a1');gtext('1/a2');gtext('1/a3');gtext('1/a4');gtext('1/a5'); a=max(a1);af=asin(max(Ft1-Ff-Fw1)/G);C=tan(af)/(a/L+hg*tan(af)/L);disp('假设后轮驱动，最大爬坡度相应的附着率=');disp(C);假设后轮驱动，最大爬坡度相应的附着率=(4) >>clearnT=;r=;f=;CDA=;i0=;If=;Iw1=;Iw2=;L=;a=;hg=;m=3880;g=;G=m*g; ig=[ ];nmin=600;nmax=4000;u1=*r*nmin./ig/i0;u2=*r*nmax./ig/i0;deta=0*ig;for i=1:5deta(i)=1+(Iw1+Iw2)/(m*r^2)+(If*(ig(i))^2*i0^2*nT)/(m*r^2); endua=[6::99];N=length(ua);n=0;Tq=0;Ft=0;inv_a=0*ua;delta=0*ua; Ff=G*f;Fw=CDA*ua.^2/;for i=1:Nif ua(i)<=u2(2)n=ua(i)*(ig(2)*i0/r)/;Tq=+*(n/1000)*(n/1000)^2+*(n/1000)^*(n/1000)^4;Ft=Tq*ig(2)*i0*nT/r;inv_a(i)=(deta(2)*m)/(Ft-Ff-Fw(i));delta(i)=*inv_a(i)/;elseif ua(i)<=u2(3)n=ua(i)*(ig(3)*i0/r)/;Tq=+*(n/1000)*(n/1000)^2+*(n/1000)^*(n/1000)^4;Ft=Tq*ig(3)*i0*nT/r;inv_a(i)=(deta(3)*m)/(Ft-Ff-Fw(i));delta(i)=*inv_a(i)/;elseif ua(i)<=u2(4)n=ua(i)*(ig(4)*i0/r)/;Tq=+*(n/1000)*(n/1000)^2+*(n/1000)^*(n/1000)^4;Ft=Tq*ig(4)*i0*nT/r;inv_a(i)=(deta(4)*m)/(Ft-Ff-Fw(i));delta(i)=*inv_a(i)/;elsen=ua(i)*(ig(5)*i0/r)/;Tq=+*(n/1000)*(n/1000)^2+*(n/1000)^*(n/1000)^4;Ft=Tq*ig(5)*i0*nT/r;inv_a(i)=(deta(5)*m)/(Ft-Ff-Fw(i));delta(i)=*inv_a(i)/;a=delta(1:k);t(i)=sum(a);endplot(t,ua);axis([0 80 0 100]);title('汽车2档原地起步换挡加速时刻曲线');xlabel('时刻t（s）');ylabel('速度ua（km/h）');>> ginputans =因此汽车2档原地起步换挡加速行驶至70km/h的加速时刻约为已知货车装用汽油发动机的负荷特性与万有特性。

汽车理论习题Matlab程序

１、3确定一轻型货车得动力性能（货车可装用４挡或５挡变速器,任选其中得一种进行整车性能计算）:1)绘制汽车驱动力与行驶阻力平衡图。

2)求汽车最高车速，最大爬坡度及克服该坡度时相应得附着率.3）绘制汽车行驶加速度倒数曲线，用图解积分法求汽车用2档起步加速行驶至70km／h得车速—时间曲线,或者用计算机求汽车用2档起步加速行驶至7０ｋm/h得加速时间。

轻型货车得有关数据：汽油发动机使用外特性得Ｔｑ－n曲线得拟合公式为式中，Tq为发动机转矩（N•ｍ);n为发动机转速（r/miｎ）。

发动机得最低转速n miｎ=600r/min，最高转速n max=4000ｒ/mｉn。

装载质量200０kｇ整车整备质量1８００kｇ总质量3880kg车轮半径０、367m传动系机械效率ηt=0、85滚动阻力系数f=0、０13空气阻力系数×迎风面积C DＡ=２、7７m2主减速器传动比ｉ0=5、83飞轮转动惯量Ｉｆ=0、２1８kg•ｍ2二前轮转动惯量Ｉｗ1＝1、７9８kg•ｍ2四后轮转动惯量I w2=３、５98kg•m ２变速器传动比ig（数据如下表)轴距Ｌ=3、2m质心至前轴距离（满载）a=1、9７4m质心高（满载）hg＝0、9ｍ解：Matlａb程序：（１) 求汽车驱动力与行驶阻力平衡图与汽车最高车速程序：n=［600:１0:４００0］;Tq＝－19、313+295、27*（n/100０）-1６5、４4*(ｎ/1000）、＾２+４0、874*（n/1０00）、＾3-３、8445＊(n/１0０0）、^4;m＝３８8０；g=9、8;ｎmin=6０0;nｍax=４000;G=ｍ*g；ig=[5、56 2、769 1、644 １、０0 0、7９3］；nT=0、85；r=0、367;f＝0、０１３；CDA=2、７７；i0=5、８３；L=3、2；a=1、947;ｈg=０、9;If=０、218;Ｉw1＝1、798；Iw2＝３、598；Ft１＝Tq*iｇ（１）＊i0*nT/r;Ft2=Tq*ig(2）＊i0*nT/r；Fｔ３=Ｔq*ｉg（3）＊i0*nT/r；Ft４=Ｔq＊ig(４)*i0*nＴ/r;Ｆｔ5=Tq*ig（5)*i0*nT／r；ua1=0、37７＊ｒ＊n/ig（１)／i0;ua2=0、377*r＊n／ig(2）/ｉ0;ua3=０、37７*r*n/iｇ(３）／i０;ua4＝0、37７＊r＊ｎ／ig（4)/i0;ua5=0、377＊r＊ｎ/iｇ(5）／i０；ua＝［０：5:12０］；Ff=G*ｆ;Fｗ=CＤA*ua、^２/２１、15;Fz=Ff+Fw；ｐlot(ua1,Ft1,ｕa２，Ft２,uａ3,Fｔ3，ua４，Ｆｔ４,uａ5，Fｔ５，ua,Ｆz);tｉｔlｅ（'驱动力-行驶阻力平衡图＇）;xlａｂel('ua(kｍ/ｓ）')；ylabel('Ft（N）＇）;gｔｅxt（'Ft1＇），gｔext（’Ft２'),gtext（’Ft3＇)，gtexｔ(＇Ｆt4＇)，gｔeｘt（'Ｆt５')，gtexｔ(’Ff+Fw’)；zoｏｍon；［x，ｙ］=gｉnpｕｔ(1）；zoｏm oｆf;dｉｓp（'汽车最高车速=’）；ｄisp(x）;disｐ（’kｍ/h’)；汽车最高车速=9９、３006kｍ/h（2）求汽车最大爬坡度程序:n=[６00:10:4000］;Tq＝—19、313+295、27*(n／1000)—165、４4*（n/１0０0)、^２+40、8７４*（n/100０)、^３-3、8４45＊（n/10０0）、＾4；m=３880；g=9、8；nmin=600;nmax=40０0;G=ｍ*g；ig=［5、５62、7691、6441、０0 0、793]；nT=0、８５;r=0、３67；f=0、0１3;CDA=2、77;i０=５、83；L＝3、2;a＝１、94７；hｇ＝0、9;If=0、21８;Iw1＝1、798;Iw2=３、598；Ft１=Tq*ｉg（１）＊ｉ0＊nT／r；ua1=0、3７7＊r＊n/ｉg(1）/i0；Ｆｆ＝Ｇ*f；Fw1=CDA*ｕa1、^2/21、１5;Ｆｚ１=Ff＋Fw1;Ｆi1=Ft1—Fz１;Zｏoｍon；imax=1００＊ｔaｎ（aｓｉn（max(Fi1/G）)）；ｄｉsp（'汽车最大爬坡度＝’）;ｄisp（imax);dｉsp（’％’);汽车最大爬坡度=35、２１９7％(3)求最大爬坡度相应得附着率与求汽车行驶加速度倒数曲线程序:cｌearn=［600：10:40０0］;Tｑ=－１9、３１3＋２９5、27＊(n/1000）-165、44*（ｎ/1000）、＾2+40、874*（ｎ／10０0）、^３－3、8445＊（ｎ/1000）、^４;ｍ=3880;g=9、８；ｎmｉｎ＝6０0；nmax＝4０00;Ｇ=m＊g；ig=［5、５6 2、76９１、644 1、00 0、793]；nT=０、85;ｒ=0、３６7;f=0、013；ＣDＡ=2、77;i0=5、83;L=3、2;a=1、9４７;ｈg=0、9；If=0、218；Iw1=1、798;Iw2＝3、５98；Ft1=Tｑ＊iｇ（1）*ｉ0＊nＴ/r;Ｆt2＝Tq＊ｉg(２）＊i0＊nT／r；Ft3=Ｔq*ig（３)＊i０＊nT/r；Ｆt4=Tq*ｉg（4)*i0＊nT/r;Ｆt5=Tｑ＊ｉg（５）*i0*nT/r;ua１=0、3７７*ｒ＊ｎ/ｉg(1)/i０；ua２＝0、377＊ｒ＊n/ｉｇ（2）/i0;ua3=0、377*r＊ｎ/ig(3）/i0;ua4=0、3７７*r*n/ig（4)／i0;ua5＝0、3７７*r*n/ｉg(5)/i0;Ｆw1=ＣDA＊ｕa1、^2/２１、１5；Fw２=CDA*ｕa2、＾２／21、15；Fｗ3=CDA*ua3、＾２/2１、15;Ｆw4=CDA＊uａ4、^2／21、１５;Fw5＝CDＡ＊ｕa5、^２/21、15；Ff＝G＊f;deta１=１＋（Iw1+Ｉw2）/(ｍ＊r＾2）+(If＊ig（1）^２*i0^2＊nT）／(m＊r^２）; deta2=1+（Iw1+Iｗ２)/（m＊r^２)+(If*ｉg(2)^2＊ｉ0^2＊ｎＴ）/（m*r^2);deta３＝1+（Iw1+Iw2）/(m*ｒ＾2）+(Ｉf*ig（３）^2＊i0^2＊nT)/(m＊ｒ＾2）; dｅtａ４=1+(Iw１＋Iｗ2）/（ｍ*r^２）+(If*ｉg（4)^2＊i0^２*ｎT）/（ｍ＊r^2)；ｄeta5=１+（Ｉw１+Iw2）/(m*r＾２)+（Ｉｆ*iｇ(5）^2＊i0＾2*nT）／（m＊r^２);ａ１=（Ft1－Ff－Fw1)/(deta1＊m）；ａd１＝1、/a1；a２=（Ft2-Ｆf—Fw2）/（deta2＊m);ad２=1、/a2;a3=（Fｔ3-Ff-Ｆw３）/（deta3＊ｍ);aｄ3=1、/a3；a4=（Ft４-Ｆf—Fｗ4）/(deｔa4＊m）；ad４=1、／ａ4;a５=（Ft5-Ff－Ｆw5）/（ｄeta5＊m)；ａd5=1、/a5;pｌｏt（ｕa1，ad1,ua2,ａd２，ua３，ad3,ｕa4,ａd4，ｕa5,aｄ5）；axiｓ(［０９9 ０1０］）；title（＇汽车得加速度倒数曲线＇)；ｘｌabel（＇ｕa(km/h)’）；ylaｂel（'１/ａ');gteｘｔ(＇1/a1');gteｘt（＇1/a2’）；gtexｔ（’1/a３＇）；gtｅxｔ（＇1/ａ４'）；gｔeｘｔ(＇１/a5'）；ａ=max(ａ1）;af=asin（max（Ft1—Ff—Fw1）/Ｇ）；Ｃ＝tan(af）/(a/L+hg＊tａn（ａｆ）/L);disp（’假设后轮驱动，最大爬坡度相应得附着率＝＇)；dｉsｐ(Ｃ);假设后轮驱动，最大爬坡度相应得附着率=0、4219（４) >>cleaｒnＴ=0、８5；r=0、３6７；f=0、0１3;CＤA＝2、77；i０＝5、83;Iｆ＝0、２18;Iw1=1、７9８;Iｗ2＝3、598；L＝3、2；a=1、９47；hg=0、9;m=38８0；ｇ＝９、8；G＝m＊ｇ;ig=[５、56 ２、7６91、644 1、00 0、793］;nmin=６０0；nmax=40０0；ｕ1=0、377*r*nmin、／iｇ／ｉ0;u2=０、３77*ｒ＊nｍax、/ig/i0;ｄeta=0＊iｇ；fｏrｉ=1:5deta（i)=1+(Iw1+Ｉw２）/（m＊r＾2)+（Iｆ*(ig(i))^2*ｉ0＾2＊nT)/(m*ｒ^2)；enｄua＝[6：0、０1：９９]；N=lｅnｇth(uａ）；ｎ=0；Tq=0;Ft=0；inｖ＿a＝０*ua;d ｅｌｔａ=0*ua；Ff=Ｇ＊f；Fw=CDA＊ｕa、^2/21、15;fｏr i＝1:Ｎk=i；ｉｆua（i)＜=u2（２）ｎ＝ｕa（ｉ)*（iｇ(2）＊i0/r)/０、3７7；Tq＝－１9、3１3+２9５、27＊（n/100０)-165、44*（n/1０00）^2+40、874*(n/10０0）^3—３、8４4５＊（n/1000）^4;Fｔ=Tq＊ｉg（2)＊i0＊ｎT／ｒ;inv＿a（ｉ）=（deta（２）＊ｍ)／（Ft-Ｆf—Ｆw（i））；deltａ（i）＝0、01＊ｉnv_a(ｉ)／３、6;elｓeiｆua(i)＜＝ｕ2（3）n=ua（i)*（iｇ（3)*i0/r）/０、３77；Tq＝—１9、3１３+2９5、２７＊（n/１00０）—1６５、44＊（n/1000)＾2+４0、87４＊(n/1０00）^3－3、84４５*(ｎ/10００）^4;Ft=Tq＊ig(３)＊i０＊ｎT／ｒ；inv_a(i）=(deta(3）*ｍ）／(Ft—Ff-Ｆw(i））；deｌtａ(i）=0、01＊inｖ_ａ（ｉ)／3、６；ｅｌseif ｕa(i)〈＝u2（4)n＝ｕa(i）*(ig(4)*ｉ0/r)/0、377；Tq=—19、31３+29５、2７＊（n/1０00）—165、44*（n／1000）^2+４０、874*(n/１000）^3-3、８445＊(ｎ/1０00)^4；Fｔ=Tｑ*ｉg(4）*i0*ｎＴ/ｒ；inv_a（i)=(deｔa（4)*m）/（Ｆt－Fｆ－Fw（ｉ)）;ｄeｌｔａ（i）＝0、01＊inv_a(i）/3、６；elｓen=ua(i）＊(ig(5)＊ｉ０/ｒ）／0、37７;Ｔq=-19、313+2９5、2７*（ｎ/100０）－１６5、44*（ｎ/１00０)^2+40、８7４*(n/100０)＾３－３、84４5＊(ｎ/1000）^4;Ft=Tq＊ｉg(５)＊i０＊nT/r；inv_a(i)=（deｔa(5)*m）/(Ft－Fｆ—Ｆw（i））;ｄelｔa（i）=０、01＊iｎv_a(i）/3、6;enda=deltａ(1：k);t(i)＝sum（a)；endplot（t，ｕa)；axiｓ（[０80 ０100］）;tｉtle('汽车2档原地起步换挡加速时间曲线’）；xlabｅｌ(’时间t（s)’）；ylabｅl（’速度ｕａ（kｍ/h)')；＞〉ｇｉnputans ＝25、82２3 7０、073725、746７70、0737所以汽车2档原地起步换挡加速行驶至7０kｍ/ｈ得加速时间约为25、8s2、7已知货车装用汽油发动机得负荷特性与万有特性。

汽车理论课后习题Matlab程序详解

陕西理工学院杨飞华1 / 341.3 确定一轻型货车的动力性能（货车可装用4挡或5挡变速器，任选其中的一种进行整车性能计算）：1）绘制汽车驱动力与行驶阻力平衡图。

2）求汽车最高车速，最大爬坡度及克服该坡度时相应的附着率。

3）绘制汽车行驶加速度倒数曲线，用图解积分法求汽车用2档起步加速行驶至70km/h的车速－时间曲线，或者用计算机求汽车用2档起步加速行驶至70km/h的加速时间。

轻型货车的有关数据：汽油发动机使用外特性的Tq-n曲线的拟合公式为23419.313295.27()165.44()40.874()3.8445()1000100010001000qnnnnT=?+?+?式中，Tq为发动机转矩（N?m）;n为发动机转速（r/min）。

发动机的最低转速nmin=600r/min,最高转速nmax=4000r/min。

装载质量 2000kg 整车整备质量1800kg 总质量3880kg 车轮半径0.367m 传动系机械效率ηt=0.85 滚动阻力系数 f=0.013 空气阻力系数×迎风面积CDA=2.77m2 主减速器传动比i0=5.83 飞轮转动惯量 If=0.218kg?m2 二前轮转动惯量Iw1=1.798kg?m2 四后轮转动惯量 Iw2=3.598kg?m2 变速器传动比 ig(数据如下表)Ⅰ档Ⅱ档Ⅲ档Ⅳ档Ⅴ档四档变速器 6.09 3.09 1.71 1.00 - 五档变速器 5.56 2.769 1.644 1.00 0.793 轴距 L=3.2m 质心至前轴距离（满载）a=1.974m 质心高（满载）hg=0.9m解：Matlab程序：(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/1000).^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;2 / 34Ft4=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.3006km/h3 / 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)求最大爬坡度相应的附着率和求汽车行驶加速度倒数曲线程序：clearn=[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;Fw1=CDA*ua1.^2/21.15;Fw2=CDA*ua2.^2/21.15;Fw3=CDA*ua3.^2/21.15;4 / 34Fw4=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);ad1=1./a1;a2=(Ft2-Ff-Fw2)/(deta2*m);ad2=1./a2;a3=(Ft3-Ff-Fw3)/(deta3*m);ad3=1./a3;a4=(Ft4-Ff-Fw4)/(deta4*m);ad4=1./a4;a5=(Ft5-Ff-Fw5)/(deta5*m);ad5=1./a5;plot(ua1,ad1,ua2,ad2,ua3,ad3,ua4,ad4,ua5,ad5);axis([0 99 0 10]);title('汽车的加速度倒数曲线');xlabel('ua(km/h)');ylabel('1/a');gtext('1/a1');gtext('1/a2');gtext('1/a3');gtext('1/a4');gtext('1/a5'); a=max(a1);af=asin(max(Ft1-Ff-Fw1)/G);C=tan(af)/(a/L+hg*tan(af)/L);disp('假设后轮驱动，最大爬坡度相应的附着率=');disp(C);假设后轮驱动，最大爬坡度相应的附着率=0.42195 / 34(4) >>clearnT=0.85;r=0.367;f=0.013;CDA=2.77;i0=5.83;If=0.218;Iw1=1.798;Iw2=3.598;L=3.2;a=1.947;hg=0.9;m=3880;g=9.8;G=m*g; ig=[5.56 2.769 1.644 1.00 0.793];nmin=600;nmax=4000;u1=0.377*r*nmin./ig/i0;u2=0.377*r*nmax./ig/i0;deta=0*ig;for i=1:5deta(i)=1+(Iw1+Iw2)/(m*r^2)+(If*(ig(i))^2*i0^2*nT)/(m*r^2);endua=[6:0.01:99];N=length(ua);n=0;Tq=0;Ft=0;inv_a=0*ua;delta=0*ua;Ff=G*f;Fw=CDA*ua.^2/21.15;for i=1:Nk=i;if ua(i)<=u2(2)n=ua(i)*(ig(2)*i0/r)/0.377;Tq=-19.313+295.27*(n/1000)-165.44*(n/1000)^2+40.874*(n/1000)^3-3.8445*(n/1000)^4;Ft=Tq*ig(2)*i0*nT/r;inv_a(i)=(deta(2)*m)/(Ft-Ff-Fw(i));delta(i)=0.01*inv_a(i)/3.6;elseif ua(i)<=u2(3)n=ua(i)*(ig(3)*i0/r)/0.377;Tq=-19.313+295.27*(n/1000)-165.44*(n/1000)^2+40.874*(n/1000)^3-3.8445*(n/100 0)^4;Ft=Tq*ig(3)*i0*nT/r;inv_a(i)=(deta(3)*m)/(Ft-Ff-Fw(i));delta(i)=0.01*inv_a(i)/3.6;elseif ua(i)<=u2(4)n=ua(i)*(ig(4)*i0/r)/0.377;Tq=-19.313+295.27*(n/1000)-165.44*(n/1000)^2+40.874*(n/1000)^3-3.8445*(n/100 0)^4;Ft=Tq*ig(4)*i0*nT/r;inv_a(i)=(deta(4)*m)/(Ft-Ff-Fw(i));delta(i)=0.01*inv_a(i)/3.6;elsen=ua(i)*(ig(5)*i0/r)/0.377;Tq=-19.313+295.27*(n/1000)-165.44*(n/1000)^2+40.874*(n/1000)^3-3.8445*(n/100武汉理工大学6 / 340)^4;Ft=Tq*ig(5)*i0*nT/r;inv_a(i)=(deta(5)*m)/(Ft-Ff-Fw(i));delta(i)=0.01*inv_a(i)/3.6;enda=delta(1:k);t(i)=sum(a);endplot(t,ua);axis([0 80 0 100]);title('汽车2档原地起步换挡加速时间曲线');xlabel('时间t（s）');ylabel('速度ua（km/h）');>> ginputans =25.8223 70.073725.7467 70.0737所以汽车2档原地起步换挡加速行驶至70km/h的加速时间约为25.8s2.7已知货车装用汽油发动机的负荷特性与万有特性。

汽车理论习题Matlab程序文件

1.3 确定一轻型货车的动力性能（货车可装用4挡或5挡变速器，任选其中的一种进行整车性能计算）：1）绘制汽车驱动力与行驶阻力平衡图。

2）求汽车最高车速，最大爬坡度及克服该坡度时相应的附着率。

发动机的最低转速n min =600r/min,最高转速n max =4000r/min 。

(完整word版)汽车理论课后作业matlab编程详解(带注释)(word文档良心出品)

1.3matlab程序：（1）%驱动力-行驶阻力平衡图%货车相关参数。

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;%计算各档位驱动力。

Fz1=m*g*f+2.77*ua1.^2/21.15;Fz2=m*g*f+2.77*ua2.^2/21.15;Fz3=m*g*f+2.77*ua3.^2/21.15;Fz4=m*g*f+2.77*ua4.^2/21.15;Fz5=m*g*f+2.77*ua5.^2/21.15;%驱动力-行驶阻力平衡图。

汽车理论习题Matlab程序

!确定一轻型货车的动力性能（货车可装用4挡或5挡变速器，任选其中的一种进行整车性能计算）：1）绘制汽车驱动力与行驶阻力平衡图。

2）求汽车最高车速，最大爬坡度及克服该坡度时相应的附着率。

轻型货车的有关数据：汽油发动机使用外特性的Tq-n 曲线的拟合公式为23419.313295.27()165.44()40.874() 3.8445()1000100010001000q n n n n T =-+-+-式中，Tq 为发动机转矩（N •m ）;n 为发动机转速（r/min ）。

发动机的最低转速n min =600r/min,最高转速n max =4000r/min 。

装载质量 2000kg整车整备质量 1800kg总质量 3880kg车轮半径传动系机械效率 ηt =滚动阻力系数 f =:空气阻力系数×迎风面积 C D A =主减速器传动比 i 0=飞轮转动惯量 I f =•m 2二前轮转动惯量 I w1=•m 2四后轮转动惯量I w2=•m2变速器传动比 ig(数据如下表)^轴距 L=质心至前轴距离（满载） a=质心高（满载） hg=解：Matlab程序：(1) 求汽车驱动力与行驶阻力平衡图和汽车最高车速程序：n=[600:10:4000];Tq=+*(n/1000)*(n/1000).^2+*(n/1000).^*(n/1000).^4;【m=3880;g=;nmin=600;nmax=4000;G=m*g;ig=[ ];nT=;r=;f=;CDA=;i0=;L=;a=;hg=;If=;Iw1=;Iw2=;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=*r*n/ig(1)/i0;ua2=*r*n/ig(2)/i0;ua3=*r*n/ig(3)/i0;ua4=*r*n/ig(4)/i0;ua5=*r*n/ig(5)/i0;ua=[0:5:120];Ff=G*f;"Fw=CDA*ua.^2/;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'),gtex t('Ff+Fw');zoom on;》[x,y]=ginput(1);zoom off;disp('汽车最高车速=');disp(x);disp('km/h');汽车最高车速=km/h-(2)求汽车最大爬坡度程序：n=[600:10:4000];Tq=+*(n/1000)*(n/1000).^2+*(n/1000).^*(n/1000).^4; m=3880;g=;nmin=600;nmax=4000;G=m*g;ig=[ ];nT=;r=;f=;CDA=;i0=;L=;a=;hg=;If=;Iw1=;Iw2=;！Ft1=Tq*ig(1)*i0*nT/r;ua1=*r*n/ig(1)/i0;Ff=G*f;Fw1=CDA*ua1.^2/;Fz1=Ff+Fw1;Fi1=Ft1-Fz1;Zoom on;imax=100*tan(asin(max(Fi1/G)));：disp('汽车最大爬坡度=');disp(imax);disp('%');汽车最大爬坡度=%(3)求最大爬坡度相应的附着率和求汽车行驶加速度倒数曲线程序：clearn=[600:10:4000];Tq=+*(n/1000)*(n/1000).^2+*(n/1000).^*(n/1000).^4;m=3880;g=;nmin=600;nmax=4000;G=m*g;ig=[ ];nT=;r=;f=;CDA=;i0=;L=;a=;hg=;If=;Iw1=;Iw2=;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=*r*n/ig(1)/i0;ua2=*r*n/ig(2)/i0;ua3=*r*n/ig(3)/i0;ua4=*r*n/ig(4)/i0;ua5=*r*n/ig(5)/i0;Fw1=CDA*ua1.^2/;Fw2=CDA*ua2.^2/;Fw3=CDA*ua3.^2/;Fw4=CDA*ua4.^2/;Fw5=CDA*ua5.^2/;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);ad1=1./a1;a2=(Ft2-Ff-Fw2)/(deta2*m);ad2=1./a2;a3=(Ft3-Ff-Fw3)/(deta3*m);ad3=1./a3;a4=(Ft4-Ff-Fw4)/(deta4*m);ad4=1./a4;a5=(Ft5-Ff-Fw5)/(deta5*m);ad5=1./a5;^plot(ua1,ad1,ua2,ad2,ua3,ad3,ua4,ad4,ua5,ad5);axis([0 99 0 10]);title('汽车的加速度倒数曲线');xlabel('ua(km/h)');ylabel('1/a');gtext('1/a1');gtext('1/a2');gtext('1/a3');gtext('1/a4');gtext('1/a5') ;a=max(a1);af=asin(max(Ft1-Ff-Fw1)/G);"C=tan(af)/(a/L+hg*tan(af)/L);disp('假设后轮驱动，最大爬坡度相应的附着率=');disp(C);假设后轮驱动，最大爬坡度相应的附着率=(4) >>clearnT=;r=;f=;CDA=;i0=;If=;Iw1=;Iw2=;L=;a=;hg=;m=3880;g=;^G=m*g; ig=[ ];nmin=600;nmax=4000;u1=*r*nmin./ig/i0;u2=*r*nmax./ig/i0;deta=0*ig;for i=1:5deta(i)=1+(Iw1+Iw2)/(m*r^2)+(If*(ig(i))^2*i0^2*nT)/(m*r^2); end！ua=[6::99];N=length(ua);n=0;Tq=0;Ft=0;inv_a=0*ua;delta=0*ua;Ff=G*f;Fw=CDA*ua.^2/;for i=1:Nk=i;if ua(i)<=u2(2)n=ua(i)*(ig(2)*i0/r)/;Tq=+*(n/1000)*(n/1000)^2+*(n/1000)^*(n/1000)^4; @Ft=Tq*ig(2)*i0*nT/r;inv_a(i)=(deta(2)*m)/(Ft-Ff-Fw(i));delta(i)=*inv_a(i)/;elseif ua(i)<=u2(3)n=ua(i)*(ig(3)*i0/r)/;Tq=+*(n/1000)*(n/1000)^2+*(n/1000)^*(n/1000)^4; Ft=Tq*ig(3)*i0*nT/r;inv_a(i)=(deta(3)*m)/(Ft-Ff-Fw(i));*delta(i)=*inv_a(i)/;elseif ua(i)<=u2(4)n=ua(i)*(ig(4)*i0/r)/;Tq=+*(n/1000)*(n/1000)^2+*(n/1000)^*(n/1000)^4; Ft=Tq*ig(4)*i0*nT/r;inv_a(i)=(deta(4)*m)/(Ft-Ff-Fw(i));delta(i)=*inv_a(i)/;else.n=ua(i)*(ig(5)*i0/r)/;Tq=+*(n/1000)*(n/1000)^2+*(n/1000)^*(n/1000)^4; Ft=Tq*ig(5)*i0*nT/r;inv_a(i)=(deta(5)*m)/(Ft-Ff-Fw(i));delta(i)=*inv_a(i)/;enda=delta(1:k);t(i)=sum(a);~endplot(t,ua);axis([0 80 0 100]);title('汽车2档原地起步换挡加速时间曲线');xlabel('时间t（s）');ylabel('速度ua（km/h）');>> ginput，ans =所以汽车2档原地起步换挡加速行驶至70km/h的加速时间约为已知货车装用汽油发动机的负荷特性与万有特性。