用matlab绘制汽车功率平衡图

用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.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; 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');