MATLAB画驱动力-行驶阻力图,两个程序

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.644 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)/100 0).^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');

>> m=3880;

>> g=9.8;

>> G=m*g;

>> ig=[5.56 2.769 1.644 1.00 0.793];

>> nT=0.85;

>> r=0.367;

>> f=0.013;

>> CDA=2.77;

>> i0=5.83;

>> L=3.2;

>> a=1.947;

>> hg=0.9;

>> If=0.218;

Iw1=1.798;

Iw2=3.598;

nmin=600;nmax=4000;

>> n=600:10:4000;

>> Tq=-19.313+295.27*(n/1000)-165.44*(n/1000).^2+40.874*(n/1000).^3-3.8445*(n/1000).^4; >> Ft1=Tq*ig(1)*i0*nT/r;

>> Ft2=Tq*ig(2)*i0*nT/r;

>> Ft3=Tq*ig(3)*i0*nT/r;

>> Ft4=Tq*ig(4)*i0*nT/r;

>> Ft5=Tq*ig(5)*i0*nT/r;

>> ua1=0.377*r*n/ig(1)/i0;

>> ua2=0.377*r*n/ig(2)/i0;

>> ua3=0.377*r*n/ig(3)/i0;

>> ua4=0.377*r*n/ig(4)/i0;

>> ua5=0.377*r*n/ig(5)/i0;

>> ua=0:5:120;

>> Ff=G*f;

>> Fw=CDA*ua.^2/21.15;

>> Fz=Ff+Fw;

>> plot(ua1,Ft1,ua2,Ft2,ua3,Ft3,ua4,Ft4,ua5,Ft5,ua,Fz);

>> title('汽车驱动力-行驶阻力平衡图');

>> xlabel('ua(km*h^-1)');

>> ylabel('Ft/N');

>> gtext('Ft1'),gtext('Ft2'),gtext('Ft3'),gtext('Ft4'),gtext('Ft5'),gtext('Ff+Fw');