%计算5个档的驱动力
for k=1:5
u=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的每次循环的图形都保留显示
end
title('每一档的驱动力和速度之间的关系')
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; break
end
end
disp '====汽车动力性能仿真计算结果===='
fprintf('汽车理论最高车速(驱动力与行驶阻力曲线交点)Vmax=%3.1fkm/h\n ',umax)