%计算5个档的驱动力for k=1:5u=Umin(k):Umax(k); %每一档的理论速度范围n=ig(k)*i0.*u./r/0.377; %每一档理论速度范围内发动机转速范围Tq=polyval(dt,n); %多项式拟合的发动机转矩和转速关系Ft=Tq.*ig(k)*i0*eta/r;plot(u,Ft) %每一档的驱动力和速度之间的关系hold on, grid on %保证K的每次循环的图形都保留显示endtitle('每一档的驱动力和速度之间的关系')hold on%行驶阻力计算f=0.012; %滚动阻力系数Cd=0.3;A=2.3;u=0:10:180;Ff=G*f; %计算滚动阻力Fw=Cd*A.*u.^2./21.15; %计算空气阻力F=Ff+Fw; %滚动阻力、空气阻力之和plot(u,F,'mo-')title('驱动力阻力图')grid on%求解最高车速k=5;for u=150:0.1:170;n=ig(k)*i0.*u./r/0.377;Tq=polyval(dt,n);Ft=Tq.*ig(k)*i0*eta/r;Fw=Cd*A.*u.^2./21.15;F=Ff+Fw;if abs(Ft-F)<1; %当驱动力与行驶阻力差值小于1N时，近似认为相等，即到达理论最高车速umax=u; breakend enddisp '====汽车动力性能仿真计算结果===='fprintf('汽车理论最高车速(驱动力与行驶阻力曲线交点)Vmax=%3.1fkm/h\n ',umax)