汽车驱动力与行驶阻力平衡图在Matlab中的程序%轻型货车有关参数
n_min=600;n_max=4000;
m1=2000;m0=1800;m=3880;
r=0.367;
eT=0.85;
f=0.013;
CDA=2.77;
i0=5.83;
If=0.218;Iw1=1.798;Iw2=3.598;
ig=[6.09 3.09 1.71 1.00];
L=3.2;
a=1.947;
b=L-a;
hg=0.9;
g=9.8;
n=(n_min:10:n_max);
N=n./1000;
Tq=-19.313+295.27.*N-165.44.*N.^2+40.874.*N.^3-3.8445.*N.^4;
Ft1=Tq.*ig(1).*i0.*eT./r;
Ft2=Tq.*ig(2).*i0.*eT./r;
Ft3=Tq.*ig(3).*i0.*eT./r;
Ft4=Tq.*ig(4).*i0.*eT./r;
s=2*pi*r*0.001;
Ua1=n./ig(1)./i0.*60.*s;
Ua2=n./ig(2)./i0.*60.*s;
Ua3=n./ig(3)./i0.*60.*s;
Ua4=n./ig(4)./i0.*60.*s;
Ff=m*g*f;
ua=(0:1:120);
Fw=CDA.*ua.^2./21.15;
F=Ff+Fw;
plot(Ua1,Ft1,Ua2,Ft2,'r',Ua3,Ft3,'b',Ua4,Ft4,'g',ua,F,'m'),grid