汽车理论课后习题matlab程序.docx
《汽车理论课后习题matlab程序.docx》由会员分享，可在线阅读，更多相关《汽车理论课后习题matlab程序.docx（7页珍藏版）》请在冰豆网上搜索。
汽车理论课后习题matlab程序
汽车理论课后习题matlab程序
1.3.1驱动力阻力平衡图
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;
end
fori=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;
end
end
plot（ua,Ft,ua,F）
xlabel（'ua/（km/h）'）;
ylabel（'F/N'）
title（'汽车驱动力-行驶阻力平衡图'）;
gtext（'Ft1'）
gtext（'Ft2'）
gtext（'Ft3'）
gtext（'Ft4'）
gtext（'Ff+Fw'）
gridon
axis（0,120,0,14000）
1.3.3加速度倒数曲线
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;
end
fori=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）;%加速度倒数%
end
end
plot（ua,b）;
xlabel（'ua/（km/h）'）;
ylabel（'1/a/（s^2/m）'）;
title（'加速度倒数曲线'）
gtext（'1/a1'）;
gtext（'1/a2'）;
gtext（'1/a3'）;
gtext（'1/a4'）;
gridon;
axis（[0,50,0,5]）
2.7功率平衡曲线
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;
end
fori=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）;
Pe（i,j）=Ft（i,j）*ua（i,j）/3600;%PeEXPRESSION%
c（i,j）=7.647*ua（i,j）/1000+3.638*ua（i,j）^3/100000;%（Pf+Pw）/gtEXPRESSION%
end;
end;
plot（ua,Pe,ua,c）;
xlabel（'ua/（km/h）'）;
ylabel（'Pe/kW'）;
title（'功率平衡图'）;
gtext（'Pe1'）;
gtext（'Pe2'）;
gtext（'Pe3'）;
gtext（'Pe4'）;
gtext（'P阻力'）;
gridon;
axis（[0,90,0,20]）;
最高档等速百公里油耗曲线
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）/（ig4*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:
100;
f=polyval（b,ua4）;
figure;
plot（ua,Qs,'ko',ua4,f,'LineWidth',2）;
title（'汽车最高档等速百公里油耗曲线'）;
xlabel（'ua/（km/h）'）;
ylabel（'Qs/[1/L（100km）]'）;
次高档等速百公里油耗曲线
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,Qs,'ko',ua4,f,'LineWidth',2）;
title（'汽车次高档等速百公里油耗曲线'）;
xlabel（'ua/（km/h）'）;
ylabel（'Qs/[1/L（100km）]'）;
升级会员 