汽车理论图形MATLAB程序

1、汽车理论图形MATL程序功率平衡图=0.31;yt=0.9;f=0.017;CD=0.31;A绘制汽车驱动力与行驶阻力平衡图汽车驱动力与行驶阻力平衡图ua(k)=0.377*r*n(k)/(ig(4)*i0);Fz(k)=m*g*f+CD*A*(ua(kF2)/21.15;E(k)=abs(Ft(k)-Fz(k);endEmin,kmin=min(E);Umax=ua(kmin)加速度曲线图m=1230;g=9.8;ig=3.6152.0531.3931.0310.837;i0=3.75;r=0.31;yt=0.9;f=0.017;CD=0.31;A=2.2;np=6000;Pemax=83;

2、deta=1.331.121.071.051.04;fori=1:56;n=500:100:6000;Pe(i)=Pemax*(n(i)/np+(n(i)/np)A2-(n(i)/npF3);Tq(i)=9549*Pe(i)/n(i);endforj=1:5fori=1:56Ft(i,j)=Tq(i)*ig(j)*i0*yt/r;ua(i,j)=0.377*r*n(i)/(ig(j)*i0);Fz(i,j)=m*g*f+CD*A*(ua(i,j)A2)/21.15;a(i,j)=(Ft(i,j)-Fz(i,j)/(deta(j)*m);b(i,j)=1./a(i,j);endendx1=ua(

3、:,1);y1=b(:,1);x2=ua(:,2);y2=b(:,2);x3=ua(:,3);y3=b(:,3);x4=ua(:,4);y4=b(:,4);x5=ua(:,5);y5=b(:,5);figureplot(x1,y1,x2,y2,x3,y3,x4,y4,x5,y5);axis(0185015);title(汽车行驶加速度倒数曲线)；xlabel(ua(km/h);ylabel(1/a);nT=0.9;r=0.31;f=0.0177;CD=0.31;A=2.2;i0=3.75;m=3880;g=9.8;Pemax=83;np=6000;G=m*g;ig=3.6152.0531.39

4、31.0310.837;nmin=600;nmax=6000;u1=0.377*r*nmin./ig/i0;u2=0.377*r*nmax./ig/i0;deta=1.331.121.071.051.04;ua=1:0.01:180;N=length(ua);n=0;Tq=0;Ft=0;inv_a=0*ua;delta=0*ua;Ff=G*f;Fw=CD*A*uaS2/21.15;fori=1:Nk=i;ifua(i)=u2(2)n=ua(i)*(ig(2)*i0/r)/0.377;Tq=9549*(Pemax*(n/np+(n/npF2-(n/npF3)/n;Ft=Tq*ig(2)*i0*n

5、T/r;inv_a(i)=(deta(2)*m)/(Ft-Ff-Fw(i);delta(i)=0.01*inv_a(i)/3.6;elseifua(i)=u2(3)n=ua(i)*(ig(3)*i0/r)/0.377;Tq=9549*(Pemax*(n/np+(n/np)A2-(n/np)A3)/n;Ft=Tq*ig(3)*i0*nT/r;inv_a(i)=(deta(3)*m)/(Ft-Ff-Fw(i);delta(i)=0.01*inv_a(i)/3.6;elseifua(i)=u2(4)n=ua(i)*(ig(4)*i0/r)/0.377;Tq=9549*(Pemax*(n/np+(n/

6、np)A2-(n/np)A3)/n;Ft=Tq*ig(4)*i0*nT/r;inv_a(i)=(deta(4)*m)/(Ft-Ff-Fw(i);delta(i)=0.01*inv_a(i)/3.6;elsen=ua(i)*(ig(5)*i0/r)/0.377;Tq=9549*(Pemax*(n/np+(n/np)A2-(n/np)A3)/n;Ft=Tq*ig(5)*i0*nT/r;inv_a(i)=(deta(5)*m)/(Ft-Ff-Fw(i);delta(i)=0.01*inv_a(i)/3.6;enda=delta(1:k);t(i)=sum(a);endplot(t,ua);axis(

7、0800100);title(汽车2档原地起步换挡加速时间曲线)；xlabel(时间t(s);ylabel(速度ua(km/h);功率平衡图m=1230;g=9.8;ig=3.6152.0531.3931.0310.837;i0=3.75;r=0.31;yt=0.9;f=0.017;CD=0.31;A=2.2;np=6000;Pemax=83;deta=1.331.121.071.051.04;forj=1:5fori=1:56;n=500:100:6000;Pe(i)=Pemax*(n(i)/np+(n(i)/np)A2-(n(i)/npF3);Tq(i)=9549*Pe(i)/n(i);u

8、a(i,j)=0.377*r*n(i)/(ig(j)*i0);Pz(i,j)=(m*g*f*ua(i,j)/3600+CD*A*ua(i,jF3/76140)/yt;endendplot(ua,Pe,ua,Pz);title(功率平衡图)；xlabel(ua(km/h);ylabel(P(kW);text(40,90,档Pe);text(80,90,二档Pe);text(110,90,三档Pe);text(150,90,四档Pe);text(200,90,五档Pe);text(100,25,Pz);%n0=815,1207,1614,2012,2603,3006,3403,3804;forj=

9、1:5fori=1:8Pd(i)=Pemax*(nO(i)/np+(nO(i)/np)A2-(nO(i)/npF3);Td(i)=9549*Pd(i)/n0(i);u(i,j)=O.377*r*nO(i)/(ig(j)*iO);endendb(1)=1326.8-416.46*Pd(1)+72.379*Pd7-5.8629*Pd(1)A3+0.17768*Pd人4;b(2)=1354.7-303.98*Pd(2)+36.657*Pd(2F2-2.0553*Pd(2)A3+0.043072*Pd(2)A4;b(3)=1284.4-189.75*Pd(3)+14.524*Pd(3)A2-0.511

10、84*Pd(3)A3+0.0068164*Pd(3)A4;b(4)=1122.9-121.59*Pd(4)+7.0035*Pd(4)A2-0.18517*Pd(4)A3+0.0018555*Pd(4)A4;b(5)=1141.0-98.893*Pd(5)+4.4763*Pd(5)A2-0.091077*Pd(5)A3+0.00068906*Pd(5)A4;b(6)=1051.2-73.714*Pd(6)+2.8593*Pd(6)A2-0.05138*Pd(6)A3+0.00035032*Pd(6)A4;b(7)=1233.9-84.478*Pd(7)+2.9788*Pd(7)A2-0.0474

11、49*Pd(7F3+0.00028230*Pd(7F4;b(8)=1129.7-45.291*Pd(8)+0.71113*Pd(8)A2+0.00075215*Pd(8)A3-0.000038568*Pd(8)A4;u1=u(:,1);u2=u(:,2);u3=u(:,3);u4=u(:,4);u5=u(:,5);B1=polyfit(u1,b,3);B2=polyfit(u2,b,3);B3=polyfit(u3,b,3);B4=polyfit(u4,b,3);B5=polyfit(u5,b,3);fork=1:56bh(k,1)=polyval(B1,ua(k,1);bh(k,2)=pol

12、yval(B2,ua(k,2);bh(k,3)=polyval(B3,ua(k,3);bh(k,4)=polyval(B4,ua(k,4);bh(k,5)=polyval(B5,ua(k,5);endfori=1:5fork=1:56Q(k,i)=Pz(k,i)*bh(k,i)/(1.02*ua(k,i)*7.10);endendplot(ua(:,4),Q(:,4),ua(:,5),Q(:,5);axis(0185020);title(最高档(五档)和次高档(四档)等速百公里油耗图);xlabel(ua(km/h);ylabel(Qs(L/100km);text(90,10,四档)；text

13、(125,10,五档)；爬坡度图m=1230;g=9.8;G=m*g;ig=3.6152.0531.3931.0310.837;i0=3.75;r=0.31;yt=0.9;f=0.017;CD=0.31;A=2.2;np=6000;Pemax=83;fori=1:56;n=500:100:6000;Pe(i)=Pemax*(n(i)/np+(n(i)/np)A2-(n(i)/npF3);Tq(i)=9549*Pe(i)/n(i);endforj=1:5fori=1:56Ft(i,j)=Tq(i)*ig(j)*i0*yt/r;ua(i,j)=0.377*r*n(i)/(ig(j)*i0);Fz(

14、i,j)=m*g*f+CD*A*(ua(i,jF2)/21.15;F=Ft-Fz;I=tan(asin(Ft/G)*100;endendplot(ua,I);holdon;title(汽车的爬坡度图)；xlabel(ua(km/h);ylabel(i(%);text(50,55,档)；text(50,35,二档)；text(100,22,三档)；text(150,15,四档)；text(200,13,五档);动力因数图m=1230；g=9.8；G=m*g；ig=3.6152.0531.3931.0310.837；i0=3.75；r=0.31；yt=0.9；f=0.017；CD=0.31；A=2.2；np=6000；Pemax=83；%绘制汽车驱动力与行驶阻力平衡图fori=1:56；n=500:100:6000；Pe(i)=Pemax*(n(i)/np+(n(i)/np)A2-(n(i)/npF3);Tq(i)=9549*Pe(i)/n(i)；endforj=1:5fori=1:56Ft(i,j)=Tq(i)*ig(j)*i0*yt/r；ua(i,j)=0.377