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.09 3.09 1.71 1.00];
for i=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
for i=1:3401;
for j=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')
grid on
axis(0,120,0,14000)