程序如下：
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 ; % 四后轮转功惯量
Ig5=[5.56 2.769 1.4 1.00 0.793];%变速器传动比
%绘制功率平衡图，以五档为例
for i=1:3401;
n(i)=i+599;
Ttq(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:5;
Ft(i,j)=Ttq(i)*i0*Ig5(j)*gt/r0;%计算驱动力
ua(i,j)=0.377*r0*n(i)./(Ig5(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('Ft5')
gtext('Ff+Fw')
grid
axis([90 100 1500 2000])
ginput(2)
