matlab汽车理论第一章.doc

m1=2000;m2=1800;mz=3880;g=9.81;r=0.367;CdA=2.77;f=0.013;nT=0.85;ig=5.56 2.769 1.644 1.00 0.793;i0=5.83;If=0.218;Iw1=1.798;Iw2=3.598;L=3.2;a=1.947;hg=0.9;Iw=2*Iw1+4*Iw2;for i=1:171n(i)=(i+29)*20;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)/1000)4;endfor j=1:5for i=1:171Ft(i,j)=Ttq(i)*ig(j)*i0*nT/r;ua(i,j)=0.377*r*n(i)/(ig(j)*i0);Fz(i,j)=CdA*ua(i,j)2/21.15+mz*g*f;endendplot(ua,Ft,ua,Fz)title(汽车驱动力与行驶阻力平衡图);xlabel(ua(km/h);ylabel(Ft(N);for k=1:175n1(k)=3300+k*0.1;Ttq(k)=-19.313+295.27*(n1(k)/1000)-165.44*(n1(k)/1000)2+40.874*(n1(k)/1000)33.8445*(n1(k)/1000)4;Ft(k)=Ttq(k)*ig(5)*i0*nT/r;ua(k)=0.377*r*n1(k)/(ig(5)*i0);Fz(k)=CdA*ua(k)2/21.15+mz*g*f;E(k)=abs(Ft(k)-Fz(k);endfor k=1:175if(E(k)=min(E)disp(汽车最高车速=); disp(ua(k);disp(km/h);endendfor p=1:150n2(p)=2000+p*0.5;Ttq(p)=-19.313+295.27*(n2(p)/1000)-165.44*(n2(p)/1000)2+40.874*(n2(p)/1000)3-3.8445*(n2(p)/1000)4;Ft(p)=Ttq(p)*ig(1)*i0*nT/r;ua(p)=0.377*r*n2(p)/(ig(1)*i0);Fz(p)=CdA*ua(p)2/21.15+mz*g*f; af(p)=asin(Ft(p)-Fz(p)/(mz*g);endfor p=1:150if(af(p)=max(af)i=tan(af(p);disp(汽车最大爬坡度=);disp(i);endendCf2=i/(a/L+hg*i/L);disp(克服该坡度时相应的附着率=);disp(Cf2);for i=1:171n(i)=(i+29)*20;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)/1000)4;endfor j=1:5for i=1:171deta=1+Iw/(mz*r2)+If*ig(j)2*i02*nT/(mz*r2);ua(i,j)=0.377*r*n(i)/(ig(j)*i0);a(i,j)=(Ttq(i)*ig(j)*i0*nT/r-CdA*ua(i,j)2/21.15-mz*g*f)/(deta*mz);if(a(i,j)0.05)b1(i,j)=a(i,j);u1(i,j)=ua(i,j);elseb1(i,j)=a(i-1,j);u1(i,j)=ua(i-1,j);endb(i,j)=1/b1(i,j);endendx1=u1(:,1);y1=b(:,1);x2=u1(:,2);y2=b(:,2);x3=u1(:,3);y3=b(:,3);x4=u1(:,4);y4=b(:,4);x5=u1(:,5);y5=b(:,5);plot(x1,y1,x2,y2,x3,y3,x4,y4,x5,y5);title(加速度倒数曲线);axis(0 120 0 30);xlabel(ua(km/h);ylabel(1/aj);for i=1:171A=ua(i,3)-ua(171,2);if (A0)j=i;endB=ua(i,4)-ua(171,3);if(B0)k=i;endif(ua(i,4)=70)m=i;endendt=ua(1,2)*b(1,2);for p1=2:171t1(p1)=(ua(p1,2)-ua(p1-1,2)*(b(p1,2)+b(p1-1,2)*0.5;t=t+t1(p1);endfor p2=j:171t2(p2)=(ua(p2,3)-ua(p2-1,3)*(b(p2,3)+b(p2-1,3)*0.5;t=t+t2(p2);endfor p3=k:mt3(p3)=(ua(p3,4)-ua(p3-1,4)*(b(p3,4)+b(p3-1,4)*0.5;t=t+t3(p3);endt=t+(ua