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

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