汽车理论（第5版）教学课件 非循环.docxVIP

syms T n F Ua; n=0:1:7200; T=1146.*(n=2000)+(443.527-0.165465.*(n-4518.25)+3.60772*10^-5.*(n-4518.25).^2-3.61349*10^-9.*(n-4518.25).^3).*(n2000); F=35.18.*T; Ua=0.00986.*n; Fzu=1266.16+9.3296.*Ua+0.263.*Ua.^2; plot(Ua,F/1000); hold on; plot(Ua,Fzu/1000,--); xlabel(Ua/(km/h)); ylabel(F/kN); title(驱动力-行驶阻力平衡图); legend(驱动力,行驶阻力); grid on; P=F.*Ua/3600; Pzu=Fzu.*Ua/3600; plot(Ua,P); hold on; plot(Ua,Pzu,--); xlabel(Ua/(km/h)); ylabel(P/kW); title(功率平衡图); grid on; D=2.11*10^-4.*T-1.58*10^-6.*Ua.^2; f=0.0076+0.000056.*Ua; H=(D-f.*sqrt(1+f.^2-D.^2))./(1+f.^2); i=tan(asin(H))*100; plot(Ua,i); xlabel(车速/(km/h)); ylabel(i/(%)); title(爬坡度曲线); grid on; i(round(15/0.00986)) ans = 23.9860 a=9.8/1.29.*(D-f); plot(Ua,a); xlabel(车速/(km/h)); ylabel(加速度/(m/s^2)); title(加速度曲线); grid on; plot(Ua,1./a); xlabel(车速/(km/h)); ylabel(加速度倒数/(s^2/m)); title(加速度倒数曲线); grid on; t(1)=Ua(1)/(3.6*a(1)); tt(1)=t(1); for k=2:7201 t(k)=(Ua(k)-Ua(k-1))/(3.6*a(k)); tt(k)=tt(k-1)+t(k); end plot(tt,Ua); xlabel(加速时间/s); ylabel(车速/(km/h)); title(加速时间曲线); grid on; tt(round(50/0.00986)) ans = 12.7598 s(1)=Ua(1)^2/(25.92*a(1)); ss(1)=s(1); for k=2:7201 s(k)=(Ua(k)^2-Ua(k-1)^2)/(25.92*a(k)); ss(k)=ss(k-1)+s(k); end plot(tt,ss) xlabel(加速时间/s); ylabel(加速距离/m); title(加速距离曲线); grid on;