汽車理論編程題

1、n=600:1:4000; r=; i0=; eff=; f=; m=3880; g=; G=m*g; CdA=; a=; hg=; L=; Iw1=; Iw2=;Iw=Iw1+Iw2; If=;for ig= ;Ua=*r*n/ig/i0;Ft=Ttq*ig*i0*eff/r;Ff=m*g*f;Fw=CdA*power(Ua,2)/;Fs=Ff+Fw;plot(Ua,Ft,Ua,Fs);hold on;endtext(15,13000,F_t_1);text(20,7000,F_t_2);text(40,4000,F_t_3);text(55,2500,F_t_4);text(70,1900

2、,F_t_5);text(105,2100,F_f+F_w);xlabel(u_a km/h);ylabel(F_t N);title（驅(qū)動力-行駛阻力平衡圖）;x,y=ginput(1); disp( 汽車的最高車速 );disp(x);disp(km/h);Ua=*r*n/ig/i0;Ft=Ttq*i0*eff/r;Ff=m*g*f;Fw=CdA*power(Ua,2)/;Fs=Ff+Fw;i=atan(Ft-Fs)/G);disp(汽車的最大爬坡度);disp(max(i);C=max(i)/(a/L+hg/L*max(i);disp(克服最大爬坡度時的附著率);disp(C);figu

3、re;for ig= ;Ua=*r*n/ig/i0;Ft=Ttq*ig*i0*eff/r;q=1+Iw/(m*rA2)+If*igA2*i0A2*eff/(m*rA2);Ff=m*g*f;Fw=CdA*power(Ua,2)/;Fs=Ff+Fw;Ac=(Ft-Fs)/q/m;plot(Ua,1./Ac);hold on;end axis(0 98 0 10);text(5,1); text(10,1,2);text(20,3);text(30,4);text(40,4,5);title( 行駛加速度倒數(shù)曲線 );xlabel(Ua/(km/h);ylabel(1/a);figure;u(1)=

4、*r*600/i0/;dt=;w=1;t(1)=0;ig=;while(u(w)=*4000*r/i0)ig=;end if(u(w+1)=*4000*r/i0)ig=;endif(u(w+1)=*4000*r/i0)ig=;endt(w+1)=(w+1)*dt;w=w+1;endplot(t,u);hold onu=70;plot(t,u);axis(0 60 0 100);title(2 擋原地起步加速時間曲線 );xlabel(時間 t/s),ylabel(車速 u/(km/h);x1,y1=ginput(1);disp( 加速到 70km/h 的時間 );disp(x1);n=600:

5、1:4000; r=; i0=; eff=; f=;m=3880; g=; G=m*g; CdA=; a=; hg=; L=; Iw1=; Iw2=;Iw=Iw1+Iw2; If=;for ig=,Ua=*r*n/ig/i0;Pe=Ttq.*n/9550;plot(Ua,Pe);hold on;endPmax=max(Pe);plot(Ua,Pmax);Ua=0:max(Ua);Pf=G*f*Ua/3600;Pw=CdA*Ua.A3/76140;plot(Ua,(Pf+Pw)/eff);text(10,50,1);text(15,35,2);text(30,40,3);text(40,30,4

6、);text(70,40,5);text(110,80,(P_f+P_w)/ n )；xlabel(u_a km/h);ylabel(F_t N);title(驅(qū)動力-行駛阻力平衡圖);title(汽車的功率平衡圖),xlabel(Ua/(km/h),ylabel(P/kw);figure;n1=815,1207,1614,2012,2603,3006,3403,3804;b0=,;b1=,;b2=,;b3=,;b4=,;B0=spline(n1,b0,n);B1=spline(n1,b1,n);B2=spline(n1,b2,n);B3=spline(n1,b3,n);B4=spline(n

7、1,b4,n);pg=7;for ig=,Ua=*r*n/ig/i0;Pf=G*f*Ua/3600;Pw=CdA*Ua.A3/76140;P=(Pf+Pw)/eff;b=B0+B1.*P+B2.*P.A2+B3.*P.A3+B4.*P.A4;Q=P.*b./*Ua*pg);plot(Ua,Q);hold on;endaxis(0 100 12 28);text(40,14,四擋);text(80,16,五檔);title(最高擋與次高擋等速百公里油耗曲線),xlabel(Ua/(km/h),ylabel(百公里油耗 Qs/(L/100km);3主函數(shù):i0=,;for i=1:1:5y(i)=

8、jiasushijian(i0(i);endy;for i=1:1:5b(i)=youhao(i0(i);endb;plot(b,y,+r);hold on;b1=linspace(b(1),b(5),100);y1=spline(b,y,b1);plot(b1,y1);text,i_0=);text,i_0=);text,i_0=);text,i_0=);text,i_0=);title( 燃油經(jīng)濟性加速時間曲線 );xlabel(百公里油耗(L/100km);ylabel(加速時間s);子函數(shù) 1function y=jiasushijian(i0)n1=linspace(0,5000);

9、nmax=4000;nmin=600;r=;yita=;CDA=;f=;G=(3880)*;ig=,; for i=1:1:4uamax(i)=chesu(nmax,r,ig(i),i0);uamin(i)=chesu(nmin,r,ig(i),i0);ua(i,:)=linspace(uamin(i),uamax(i),100);n(i,:)=zhuansu(ua(i,:),r,ig(i),i0);Ttq(i,:)=zhuanju(n(i,:);Ft(i,:)=qudongli(Ttq(i,:),ig(i),i0,yita,r);F(i,:)=f*G+CDA*(ua(i,:)A2)/;del

10、ta(i,:)=1+*(ig(i)A2)*(i0A2)*yita)/(3880*rA2);a(i,:)=1./(delta(i,:).*3880./(Ft(i,:)-F(i,:);F2(i,:)=Ft(i,:)-F(i,:);end temp1(1,:)=ua(2,:)/;temp1(2,:)=1./a(2,:);n1=1;for j1=1:1:100if ua(3,j1)max(ua(2,:)&ua(3,j1)max(ua(3,:)&ua(4,j1)=70;temp3(1,n2)=ua(4,j1)/;temp3(2,n2)=1./a(4,j1);n2=n2+1;endendy=temp1(1

11、,1)*temp1(2,1)+qiuji(temp1(1,:),temp1(2,:)+qiuji(temp2(1,:),temp2(2,:)+qiuji (temp3(1,:),temp3(2,:);子函數(shù) 2 function ua = chesu( n,r,ig,i0 )ua=*r*n/i0./ig;end子函數(shù) 3function q = jiasu( umax,umin,ig,a,ua0,i0 )global f G CDA yita m r If Iw1 Iw2 pg B0 B1 B2 B3 B4 n; ua1=umin:1:umax;delta=1+(Iw1+Iw2)/(m*rA2

12、)+(If*igA2*i0A2*yita)/(m*rA2);PO=(G*f*uaO./36OO+CDA.*uaO.A3/7614O+(delta*m*uaO/36OO)*a)/yita;P=(G*f.*ua1/3600+CDA.*ua1.A3/76140+(delta*m.*ua1/3600)*a)/yita; dt=1/*a)for i=1:1:34O1 bO(i)=BO(i)+B1(i)*PO(i)+B2(i)*(PO(i)A2+B3(i)*(PO(i)A3+B4(i)*(PO(i)A4;endb1=interp1(uaO,bO,ua1);Qt=P.*b1./.*pg);i1=size(Q

13、t);i=i1(2);Qt1=Qt(2:i-1);q=(Qt(1)+Qt(i)*dt./2+sum(Qt1)*dt;end子函數(shù) 4 function y = zhuanju( n )y=+.*( n./1000).*(門./1000)八2+*( n./1000).A*( n./IOOO).%;end子函數(shù) 5function n = zhuansu( ua,r,ig,i0 )n=ig*i0.*ua./*r);end子函數(shù) 6function b = youhao( i0 )global f G CDA yita m r If Iw1 Iw2 pg B0 B1 B2 B3 B4 n ig=,;

14、r=;yita=;CDA=;f=;%i0=;G=(3880)*;If=;Iw1=;Iw2=;m=3880;n0=815 1207 1614 2012 2603 3006 3403 3804;B00=;B10=;B20=;B30=;B40=;n=600:1:4000;B0=spline(n0,B00,n); B1=spline(n0,B10,n); B2=spline(n0,B20,n);B3=spline(n0,B30,n); B4=spline(n0,B40,n); ua4=*r.*n./(i0*ig(4);F4=f*G+CDA*(ua4.A2)/;P_fw4=F4.*ua4./(yita*

15、1000);for i=1:1:3401b4(i)=B0(i)+B1(i)*P_fw4(i)+B2(i)*(P_fw4(i)F2+B3(i)*(P_fw4(i)F3+B4(i)*(P_fw4(i)F4; endpg=;ua4_m=25,40,50;s_m=50,250,250;b4_m=spline(ua4,b4,ua4_m);F4_m=f*G+CDA*(ua4_m.A2)/;P_fw4_m=F4_m.*ua4_m./(yita*1000);Q4_m=P_fw4_m.*b4_m.*s_m./(102.*ua4_m.*pg);Q4_a1=jiasu(40,25,ig(4),ua4,i0);Q4_

16、a2=jiasu(50,40,ig(4),ua4,i0);Qid=;tid=;s=1075;Q_i=Qid*tid;Q4all=(sum(Q4_m)+Q4_a1+Q4_a2+Q_i)*100/s;b=Q4all;end子函數(shù) 7 function p = qiuji( x0,y0 )n0=size(x0); n=n0(2);x=linspace(x0(1),x0(n),200) ;y=spline(x0,y0,x);plot(x,y);p=trapz(x,y) ;endbeta=;L=;a1=;b1=L-a1;hg1=;a2=;b2=L-a2;hg2=;g=;z=0:1; faif1=beta

17、*z./(1/L*(b1+z*hg1);fair1=(1-beta)*z./(1/L*(a1-z*hg1);faif2=beta*z./(1/L*(b2+z*hg2);fair2=(1-beta)*z./(1/L*(a2-z*hg2);plot(z,faif1,r,z,fair1,g,z,faif2,b,z,fair2,k);text, _r(空車text, _r(滿載text,1, =z);text, _f()滿);載text, _f()空);車xlabel(制動強度 z/g);ylabel(利用附著系數(shù) );title（利用附著系數(shù)與制動強度的關系曲線）;fai=0:1;Ef1=(b1/L

18、)./(beta-fai*hg1/L); Er1=(a1/L)./(1-beta)+fai*hg1/L);Ef2=(b2/L)./(beta-fai*hg2/L);Er2=(a2/L)./(1-beta)+fai*hg2/L);figure;plot(fai,Ef1,r,fai,Er1,g,fai,Ef2,b,fai,Er2,k); text,E_f);text,滿載);text,E_r);text,E_r);text,空載);xlabel(附著系數(shù))；ylabel(制動效率);title(前、后制動效率曲線);axis(0 1 0 1);fai=;Er1=(a1/L)/(1-beta)+fa

19、i*hg1/L);Er2=(a2/L)/(1-beta)+fai*hg2/L); abmax1=fai*Er1*g;abmax2=fai*Er2*g;t21=;t22=;u0=30;s1=1/*(t21+t22/2)*uO+uOA2/abmax1;s2=1/*(t21+t22/2)*u0+u0A2/abmax2;disp(空載時的制動距離),disp(s1);disp(滿載時的制動距離),disp(s2);z1=a1*fai/(L+fai*hg1);z2=a2*fai/(L+fai*hg2);abmax1=z1*g;abmax2=z2*g;s1=1/*(t21+t22/2)*uO+uOA2/a

20、bmax1;s2=1/*(t21+t22/2)*u0+u0A2/abmax2;disp(前部管路損壞時空載制動距離);disp(s1);disp(前部管路損壞時滿載制動距離 );disp(s2);z1=b1*fai/(L-fai*hg1);z2=b2*fai/(L-fai*hg2);abmax1=z1*g;abmax2=z2*g;s1=1/*(t21+t22/2)*u0+u0A2/abmax1;s2=1/*(t21+t22/2)*u0+u0A2/abmax2;disp(后部管路損壞時空載制動距離);disp(s1);disp(后部管路損壞時滿載制動距離 );disp(s2);m=;Iz=3885;L=;a=;b=;k1=-62618;k2=-110185; i=20;g=;R0=15;u1=;K=m*(a/k2-b/k1)/LA2;Uch=(1/K)A(1/2);disp(穩(wěn)定性因