驱动力-行驶阻力平衡图matlab编程

汽车驱动力与行驶阻力平衡图在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