汽車理論課后習題matlab程序

1.3.1驅(qū)動力阻力平衡圖m1=2000;m2=1800;m=3880;r0=0.367;gt=0.85;f=0.013;CDA=2.77;i0=5.83;If=0.218;Iw1=1.798;Iw2=3.598;Ig4=[6.093.091.711.00];fori=1:3401;n(i)=i+599;Tq(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;endfor i=1:3401;forj=1:4;Ft(i,j)=Tq(i)*i0*Ig4(j)*gt/r0;ua(i,j)=0.377*r0*n(i)./Ig4(j)/i0;F(i,j)=f*m*9.8+CDA*ua(i,j).^2/21.15;endendplot(ua,Ft,ua,F)xlabel('ua/(km/h)');ylabel('F/N')title('汽車驅(qū)動力-行駛阻力平衡圖');gtext('Ft1')gtext('Ft2')gtext('Ft3')gtext('Ft4')gtext('Ff+Fw')gridonaxis(0,120,0,14000)1.3.3加速度倒數(shù)曲線m1=2000;m2=1800;m=3880;r0=0.367;gt=0.85;f=0.013;CDA=2.77;i0=5.83;If=0.218;Iw1=1.798;Iw2=3.598;Ig4=[6.093.091.711.00];q=1.128;fori=1:3401;n(i)=i+599;Tq(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;endfori=1:3401;forj=1:4;Ft(i,j)=Tq(i)*i0*Ig4(j)*gt/r0;ua(i,j)=0.377*r0*n(i)./Ig4(j)/i0;F(i,j)=f*m*9.8+CDA*ua(i,j).^2/21.15;b(i,j)=q*m/(Ft(i,j)-0.131*ua(i,j)^2-9.8*m*f);%加速度倒數(shù)%endendplot(ua,b);xlabel('ua/(km/h)');ylabel('1/a/(s^2/m)');title('加速度倒數(shù)曲線')gtext('1/a1');gtext('1/a2');gtext('1/a3');gtext('1/a4');gridon;axis([0,50,0,5])次高檔等速百公里油耗曲線i0=5.83;r=0.367;t=0.85;m=3880;G=38024;f=0.013;CdA=2.77;pg=7.12;ig3=1.71;ig4=1;n=[8151207161420122603300634033804];B0=[1326.81354.71284.41122.91141.01051.21233.91129.7];B1=[-416.46-303.98-189.75-121.59-98.893-73.714-84.478-45.291];B2=[72.37936.65714.5247.00354.47632.85932.97880.71113];B3=[-5.8629-2.0553-0.51184-0.18517-0.091077-0.05138-0.047449-0.00075215];B4=[0.177680.0430720.00681640.00185550.000689060.000350320.00028230-0.000038568];ua=(0.377*r*n)/(ig3*i0);Tq=-19.313+295.27*n/1000-165.44*((n/1000).^2)+40.874*((n/1000).^3)-3.8445*((n/1000).^4);P=(1/t)*(G*f*ua/3600+CdA*ua.^3/76140);b=B0+B1.*P+B2.*P.^2+B3.*P.^3+B4.*P.^4;Qs=P.*b./(1.02*ua*pg);b=polyfit(ua,Qs,3);ua4=8:1:55;f=polyval(b,ua4);figure;plot(ua