用matlab绘制汽车功率平衡图.docx

用matlab绘制汽车功率平衡图
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;
m=3880;g=9.8;
nmin=600;nmax=4000;
G=m*g;
ig=[5.562.7691.6441.000.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;

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;
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;
Pe1=Ft1.*ua1/3600;
Pe2=Ft2.*ua2/3600;
Pe3=Ft3.*ua3/3600;
Pe4=Ft4.*ua4/3600;
Pe5=Ft5.*ua5/3600;

ua=0:5:120;
Pf=G*f.*ua/3600;
Pw=CDA.*ua.^3/76140;
Pz=(Pf+Pw)/nT;
Pe=max(Pe1);

plot(ua1,Pe1,ua2,Pe2,ua3,Pe3,ua4,Pe4,ua5,Pe5,ua,Pz,ua,Pe);
title('汽车功率平衡图');
xlabel('ua/(km*h^-1)');
ylabel('Pe/kw');
gtext('¢ñ'),gtext('¢ò'),gtext('¢ó'),gtext('¢ô'),gtext('¢õ'),gtext('(Pf+Pw)/nT'),gtext('Pe');