1、、 3n=600:1 ： 4000; r= 、 7； i0=5 、 83； e f=0 、85 ； f 0、 01；m 3 80; g、 8； = *g ； Cd=2 、77;a=1 、 4 ; hg 0、9;L=3 、2； Iw1 1、798 ； Iw =3 、 5 8; Iw=I 1+Iw2 ； If 0、 18 ；Ttq=-19 、31 + 95、 7 / 000 、 44*p we（ n ，2)/1000000+ 0、 74*po r（， ) 000 00000-3 、 84 powe （ n ， )/1 ;fig= 、 5 2、 6 1、 644 1 、 、 7 3 ; Ua 0、

2、7 r*n ig/i ；Ft=Ttq*ig 0*ef / ； m g f; w= dA o e（ Ua， )/2 、 1; s f+Fw;pl ( a， Ft,Ua， s）；hon； nd x（ 5， 100 , t1);text( 0,70， F_；）te (40， 4000, F _3；）te t( 5,2 0 ,F t_ ）； t（ 70,1900,F_ 5 ）； xt(105 ， 1 0, F f F_）;xab l( _a km/h ） ;ylbel （ F_tN）;tit e（ 驅(qū)動力行駛阻力平衡圖）； , = npu （ 1）； di p( 汽車得最高車速） ;disp（ x）

3、;di p( m/h ）； Ua=0 、 377*r n ig/i0 ； t= t 5 、 * 0 eff/ ; =m g* ;Fw=CdA*powe （ a， )/21 、 12； Ff+Fw ；i=atan （ Ft-Fs)/ ） ;di p( 汽車得最大爬坡度 ）； i p（ma (i） ;C a (i)/（ a/L+h /L max(i） );d sp（ 克服最大爬坡度時(shí)得附著率fi ur ; ）; isp （）；for ig=5 、 6 、 76 1 、 44Ua=0 、 377 r* i /i0;Ft= q g i0ef /r ；1、 00、 793 ;q=1+Iw/ （ m 2)

4、+If i 0 e / （ m*r ); f m g f;F=C A p wer （ Ua, )/ 1、 1 ;Fs Ff Fw；Ac= （ F s)/ / ； lot(Ua ， 1、 /Ac ） ; hold ；endaxi (8 ）； ex ( ,0、 5， 1); text （ 1， 1， 2 ) ；te t(20,1 、 5， ）； ext （ 30， 2、 5， 4）；text( 0， 4， 5 ）；t tl（行駛加速度倒數(shù)曲線);xlabel （ Ua/ （k / ） ); labe （ 1/a );figure ；u( ) 0、 77 r*6 0/i0/ 、 7 9； dt= 、

5、 01； w=1 ；（ 1)=0;i =2 、 7 9; wh le（ u（ w) 99、 0）n(w ） =u(w ） *ig i0/0 、 37 r;T q(w)=-1 、 3+ 5 、2 *n （w ） 10 0-165 、 44*( (w） / 0 0) + 、87 (（ w） / 0 0） 3 3 、84（ n （w)/1000 ） 4 ；Ft（ w)=Ttq(w ） ig*i0*eff ； G f；Fw（ w） =C *u(w)2 1、 15 ； 1+I / （ m r ） +I g 2 i0 e f（ r 2); a（） =3 、 *（ t（）（ Ff+F（ w ）)） / m;

6、 u（ w+1) （ w） +a （ w) t；i (u（ +1 ）= 、 3 7*40 *r/ 、 7 /i0 ）ig=1 、 644；end f(u （ +1 ） =0 、 377*400 r/1 、 / ) i = 、；eni（ u(w 1） = 、 377 400 r/1 、 0/i0 ）ig=0 、 793 ；ent( +1)= （ w )*d ；ww+1 ；ed lo （ ,u） ; old n u 7；lot(t,u ）；axis（ 0 60 0 100 ）;t tle （ 2 擋原地起步加速時(shí)間曲線）；xla (時(shí)間 t/s) ， ylabel(車u速（ km/h ） ） ;，

7、 y1 ginpu （ 1);disp （加速到70km/h 得時(shí)間 ） ; is（ x1） ;2、7=600 ：1 ：4000;r 0、3；i0=5 、83； ef =0 、8 ；f、 0 3；m=3880 ；=9 、 8； Gm*g ；Cd =2 、 7； a= 、94; hg=0 、 9； L=3 、； w1 、798; I2 3、5 8；Iw=I + w2; I =0 、 2 8; q 19、 313+295、 2 n/ 0 0 165、 44 pow (n,2) 100 00 +40 、84*po e（ n，3） /10 0 00 3、 8445*power （n ， 4） 10；f

8、o ig=5 、 56， 2、 769 ， 1、 44， 1、 00， 0、 Ua=0 、 77*r n/ig/i0;P q 、 * /95 0；plot （Ua ，P）；hold n; dPmax=m x(Pe）；lot(Ua ， Pmx);Ua=0 ：、 1: ax（ Ua） ;P G*f* a/3 00；Pw=C A a、3/76140 ；plot(U ,(Pf+Pw ）/ef );tex ( 0， 50， ) ；tex （ 15,35 ， 2)； ext （ 30,40, 3；）t xt（ 4， 30， 4 ）；text （ 7， 4， 5 ) ；tex (11 ,80， ( _f+P

9、_w)/ )； l el（ u a k /h ) ；ylabel( Ft );title （驅(qū)動力行駛阻力平衡圖） ;t le（汽車得功率平衡圖） , l bel(U / （ m/h ）， ylab （ P/k )；fi ure ；n1 815,1 0， 161 , 01， 6 3,3006 ， 3403 ， 04 ;b0=1326 、8，1354 、7，1284 、4 ，1122 、9， 14 、， 105 、 ,1 3、 9, 9、；b1= 416 、46， -30 、 98,- 8 、 7 ,-121 、 59, 98、8 ， -73 、714 ， 8、 478 ，4 、 291 ； 2

10、= 7、 79， 6、657,14 、 24，7 、00 5，4、 63，2 、8 93,2、 88，0 、7 11；b3 5、8 9，-2 、05 ,、 5 1 4, 0、 517, 0、091 7， 0、 5138 ， 0 、047449 ，、 0 752 5 ;b4= 0、 17 68，0 、0 02 ，0、 008 6 ,0、 0018 55， 0、 068，、 0003503 ， 0、 00 3 ,- 、 00 8568 ;B0=spl ne（ n1,， n)；B1spl n ( 1， b1, )；B2=spline(n1 ， b2,n);3= pl n （ n， b3，）；B4=s

11、p ne（ n1,b ,n);pg=7 ；forg=1 、 00， 0、 a= 、 377* *n/i /i0 ；PfG f* a/3600 ；P = d a、 3/ 614 ;P=(Pf+P ） /eff;b= 0+B 、 *P+B2 、 *P、 B3、 P、 3 B、 *P、 4; Q= 、 b、 /(1、 0 *Ua*p )；plo （ U， Q）；ho d n;edax（ 0 00 12 28)tex （ ,14， 四擋text(80 ， 16,五檔）；);title （最高擋與次高擋等速百公里油耗曲線百公里油耗Qs/(L km） ） ;） , abe (U /(km/h） ）， ye

12、l（3主函數(shù)：i0= 、 17，5 、 3，、 8 ,6、 17， 6、 3； r 1:1:5y（ i ） =ji s shiji n（ i0(）； dy；fo i= :：b(i） youh o （ i0（i） )； d；plo （ b,y,+ )； old n；b =l nspac (b(1 ），b(5),100) ；y1=s li e（ b,y， b )； ot （ b1， y）；te t （13、， 2、 , i0=5 、 17 );tex （ 13、， 7、 6， _0= 、 4 )；text(13 、， 27、， i 0= 、 3）；te（ 4、， 26、 6，6、 17） ;te

13、t （ 4、 2， 26、 2,i_0=6 、 33）；ti le（燃油經(jīng)濟(jì)性加速時(shí)間曲線)；xabe ( 百公里油耗（L/100km ） ）；yla el(加速時(shí)間s）;子函數(shù) 1f nc ioy=jiasushi ian （ )n =li ace(0,5000 ）；nm x=4000;nmin=6 0; =0 、 36； yit 0、85; DA=2 、77； f、 013 ；G= （ 3 8 )*9 、； ig= 6、 9,3、09,、 71 ， 1、00 ；for=1: ： 4uamax(i)= hes（ nmax ， r， ig（） ,i )；uami (i） hesu （ nmin

14、 ， r， ig(i ）,i0）; a（i， :) l nspace(uami (i）， uama （ )， 100 ）；n( ,:） = h an（ u（ i，： ),r,ig(i), 0）；Ttq（ i,： )=zhuanju （ n （i，： )；（，： )=q ngli(T q( ，：）,ig（），i0,yi a， )； (i,： )= G+CDA*(ua （ i,:）、 2) 21 、 5 ；delta（i ，：) +（、7 8 3、59 +0 、2 8 (ig(i) 2)( 2) y a)/(3880*r）；a（ i,：)=1 、 /(de a（ i， :)、 388 、 / （

15、Ft（ i,：）- （，： )） ; (，：）=F ( ,： )F(i,:); d mp (1,：)=ua(2,: ） 3、6;t mp1( ，：） =1 、/a （ 2,:);n1=1; or j1=1:1 ： 1 0 f ua （， j1)max( (2,：） & （ 3， j1） 0 em 2 （，） =ua( ， j1） /3 、 6;temp （ 2,） 1、 /a （， j1);n1=n1+ ;enenn2=1 ；for j1=1 ： 1:1 0if a（ 4， 1) ma（ ua（ 3,： )& ua（4, 1） =70; 3（ ,n2)=ua （ 4， j1)/3 、 6;te

16、mp3 （2 ，n2） =1 、 /a （， j1）; n2 1;ed nd te 1（ 1,1)* e 1(2，） qiuji （ temp1(1 ，： )， te p （，： )+qiuj （m （ 1,:）， temp2 （2，： )） +qi i（tem 3(1，:)，temp （ ,:）)； end子函數(shù) 2funct on a =ches ( n,， g, 0 )ua=0 、 3 7 r* /i 、 /i ;e子函數(shù) 3fu c ion jiasu( umax,u i， ig,a ， ua,iglobalf DA yi a mrIf Iw1 Iw2 pgB0B1 B2) 3B4 n

17、；ua1 in:1: max;delta=1 (Iw1+Iw2)/(m 2)+ （ i 2*i 2 yi a)/( r )；P0=(G f、* a0、/36 +CDA 、 ua0、3 761 0+(delta m、 a0 36 0） a） /yita;P（ G*f 、* 1/3600+CDA 、 u 、 3/7 10+ （ l a m、 ua /360 )*） /y a;dt=1/ （、 6 a)for=1 ： 1： 401b0（） =B（） + (i)* 0（ i)+B2 （ i）（ P0(i)） 2 B3(i)*(P0(i ） )3+B （） *（ P（ i) ; ndb1 int rp1

18、 （ u ,b0， 1）；Q= 、 *b1 、 / （ 367、 1、 pg) ；i1=size （ t） ;i=i1 （ 2); =Qt （ 2:i 1） ;q= （ Q（） +Q （ i） t 、 /2+sum(Qt1 ） dt ；nd子函數(shù) 4un tio y = huanj (n )y= 9、 33+2 5、27、*（ n 、 1000) 、 44 、（ n、 /1000) 、 +40 、 74、（ n、 / 000 ）、 3 3 、8 45、（、 /1000 ）、 4;ed子函數(shù) 5f nc i nn=ig i、 n = huansu （ua、 / （ 0、 377 r);ua,，

19、i，0 )end子函數(shù) 6f ncti n b ouh o （ 0 )g o al f GCDA y a m IIw1 Iw2 pgig=6 、09,3 、09， 1、 7， 1、 00； r= 、 36 ;B1B3 B4 nyi 0、 85； CD 2 、77;f=0 、 013； %i、 3 ；G (3 80） 9、； If=0 、21； Iw1=1 、7 8;I 2= 、 598;m=38 0;n0=81 1207 1614 20126 3 0 6 34038 ; 0=1326 、8 1、 7 1284 、4 1122 、9 41、051 、2 23、9 129 、7；B10= 416

20、、 6 30、 98 -1 、 7 -121 、59 -98 、 89 - 3、 714 -84 、 7 45、291 ；B20= 2、 379 36 、 6 7、 24 7、 0035、 476 2、 8 93 2 、97 8 0、 7113;30 5、8629 -2 、05530、51184 - 、 18517 0、091 0、0138 -0 、044 9 、 000 5 15 ;B40= 、 1 768 0 、 430720、 00681 40、 001 555、 0 689 60、 003032 0、 002 30 0、 00 0 85 8；=6 ： 400；0=spline （ 0,

21、 00， );B =spline （ n0,B1， n）； B2spl e（ n0,B2， n)；B3=spli e（n0， B30， n ）； B4=sp e（n ， B4， n) ；ua 0、 377r 、 n 、 / （ i0* g （ 4)；F4=f G DA（ ua、 2)/ 、 15；P f =F4 、 ua4 、 /(yi a 3、 1000 ） ;fo i=1: :3 01b4 （ i)= 0(） +B （ i） _fw4(i)+B2 （i） *( _ w4(i） )2+B3 （ )*（ P f 4（ )）3+B4 （ i） *(P_ 4（ i) ;endpg 7、 06；ua

22、4_m= 25， 0，；s_m 50,25 ， 5；b4_m=spline （ a4， b4， ua4_ )；F _m=f G+CDA*(ua4 、 2 ） /21 、 15；P_ w m= 4 m、 *ua4_m 、 / （yita*3 、 6 1000) ；Q4_ = _fw4_m 、 b m、 *、 /(102 、 * a4、 *p ）Q4 a i u（4 ， 5， ig（ 4 ）， 0、 25， ua4,i0 ）；Q4 a2=j asu（ 0， 0,ig （ 4）， 0、 2， ua4 ， i）；Qi =0 、2 ； tid=19 、 3;s=1075 ；Q_i=Qid*tid ；Q4

23、all=(su ( 4 m） +Q4 1+Q4_a2+ _i) 100/s ；=Q4all ；end子函數(shù) 7funct on p =qi i（ x0,y0 )n =s ze(x0）； n n （ 2）；xl nsp ce（x (1） , 0（ n)， 200)y=spline （ x0，）；plot （ , )；p=tr z(x,y)；en4、3beta= 、 38； L=3 、 95； =2 、； b1 L 1 ； g 、 845 ； a2 2、 95;b2=L 2； hg2=1 、 1； g 9 、8 ； z=0 ： 0 、0： 1;f if1= e a、（ 1/L* （ b1+z* g

24、1) ）；fair1 (1 bet ） *z、（ 1/L*( z h 1)； if2 beta z、 /(1 L*（ b2+z g ） )；fair =(1 beta ） *、 / （ 1/L (a z g );plot （， faif1 ， r , ，f ir1,，gz, aif2 ， b，z， fair2,k) ；te（ 0、， 1、 4， r（空車 )） ;te t （0、 8，1 、2 ， _r( 滿載）);t xt（ 0、 9， =z ）;te t(0、 7， 0、 , _f （滿載）;）te (、 6,0、 3, _（空車）)；la e(制動強(qiáng)度z/ );la l( 利用附著系數(shù)）

25、;titl （ 利用附著系數(shù)與制動強(qiáng)度得關(guān)系曲線fai=0 ： 0 、01： 1;）；Ef1=(b1 ）、 /(be a fai g1/ );Er1= （a /L）、 /(1 beta ） +fai hg1/ );Ef =(b2 ）、 /(beta a g2/L ） ;Er2 (a2/ ）、 /( be a） + ai*hg2/L) ； ig re； ot （ ai,Ef1,ar,Er1， g ， a， f2， b ，f， E 2， k;）t x（ 0、 3，、， E_f）；tex （、 ,0、 9， 滿載） ;te t(、 9,、 8, E_r ) ；te (0、 7,、 7，）;t xt(

26、0 、 7，、 6，空載 );xla (附著系數(shù)）;ylabel( 制動效率 ) ；t le（前、后制動效率曲線 )；ax s（0 10）；ai、；Er1= （a /L） / （ 1 et） fai* g1/L) ；Er2=(a2/L ） /( （ 1-b t ） fai hg /）；a ma fai r1* ; bma 2=fa Er *g; t21=0 、 02；t22=0 、 2； 3；s =1 3 、 6(t21+t 2/ ） 0+ 02/2 、 /abmax1 ； 2 1/3 、 6 ( 1+t22 2） u 02/2 、 92/ab ax2； isp （空載時(shí)得制動距離）， d s

27、p（ s1） ;di（滿載時(shí)得制動距離 ）， i p(s） ; z1=a1*fai （ +f i*hg1) ；z 2 fai（ L+fai hg2 ）；ab 1=z1*g ；a m 2 2 g ；s1=1/3 、 6 (t 1+t 2 2） 0+ 02/25 、 2/a ax；s =1/ 、 6（ 2 t2 / ） *u0+ 2/ 5、 92 abma ; i（前部管路損壞時(shí)空載制動距離 ); isp （ s）；di （前部管路損壞時(shí)滿載制動距離 ）； d s (s2） ;z1 b1 ai/(L-fai* g1）； 2=b2*fai （ L- i*hg2 ）； bmax1=z *g ； abmax2=z2 g ；s1=1/3 、 *（ t21+ 22 2) u+u0 2/25 、 92/abmax1 ；s2=1 3、 6（ t2 22/2 ） u0u0 2/2 、 9 ab ax2；d sp（ 后部管路損壞時(shí)空載制動距