汽车理论编程题.doc

］;
  　Ua＝0、３７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）；
　　　　hoｌｄ　ｏn；　
ｅnd
ｔｅxｔ（１5，1３00０,’Ｆ＿t_1');
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_ｗ’）;
xｌabｅl(’ｕ_a km/h'）;
ylａbel（＇F_t　N’）;
titｌe（'驱动力－行驶阻力平衡图＇）；
［ｘ,ｙ］=ｇｉnpuｔ（1）； diｓp('汽车得最高车速’）;disp（x）;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（'克服最大爬坡度时得附着率'）;ｄisp（Ｃ）；
fiｇurｅ;
for ig=[5、５6 ２、76９ 1、６44　1、０0　0、793］;
   Ua=0、377＊r*ｎ／iｇ/i0;
   Ft=Ｔｔq＊ｉg＊i0＊efｆ/r；
  　q=1+Iw/（m＊ｒ^2)+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 ｏｎ；　
end
axiｓ(［０ ９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ｅ（＇行驶加速度倒数曲线');
xlabel（’Ua/（kｍ/ｈ）');
ｙlabeｌ（'1/a＇);
figure；
u(１)＝0、３77＊r*6０0/i0/２、7６9；dt=０、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／２115；
ｑ＝1+Iｗ/（m＊r＾２）+Iｆ＊ｉg＾2＊i0^２＊eｆf／（ｍ＊r＾2);
a（ｗ）=3、６*（Ｆt（ｗ）—（Ff+Fｗ（w）)）/ｑ／m;
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ｇ=１、００；
enｄ
iｆ（u(w＋1）>=０、377＊400０＊r/1、０0/i0）
 　　ig=0、793；
enｄ
t(ｗ+1)=（w＋１)dｔ；
w＝w+1；
eｎd
ｐloｔ（ｔ,u）;
ｈold ｏn
u＝7０；
ｐlot(t,u）；
axis（［0 60 0 100]）;
tｉtle（＇2挡原地起步加速时间曲线’）；
xlaｂｅｌ('时间t/s')，ylabel(’车速u／（km/h）'）; 　
[ｘ１，y1]＝ginpuｔ（1);disp（’加速到70km/h得时间'）;ｄisｐ（x1）;

2、7
ｎ=600：1：4000;　r＝0、3６７；　i0=5、83； efｆ=0、8５；　f＝０、0１3；  m=3880； ｇ=9、8； G＝m*g； CdＡ=2、７7； a=１、94７; hg=0、9； L=3、２；　Ｉw１＝1、798;　Iｗ2＝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、8７4*poｗeｒ（n，3）/10００0０00００－3、8445*power（n，4）／10；
foｒ　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;
ｅｎd
Pmax=mａx(Pe）；
ｐlot(Ua，Pmａx);
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_w)/η＇)；
ｘlａｂel（＇u＿a kｍ/h’)；
ylabel(’F_t Ｎ＇);
title（’驱动力—行驶阻力平衡图’）;
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=［7２、３79，３6、657,14、５24，7、00３5，4、４７63，2、8５93,2、９７88，0、7１11３］；
b3＝[—5、8６２9，-2、0５5３,—０、5１1８4,—0、１８517,—0、091０7７，—0、０5138，－0、047449，—０、０0０752１5］;
b4=［0、17７68，0、0４３0７2，0、00６8１6４,0、0018５55，0、００068９０６，０、0003503２，0、０00２８２3０,-０、00００３8568］;
B0=splｉne（n1,ｂ０，n)；
B1＝splｉnｅ(ｎ1，b1,ｎ)；
B2=spline(n1，b2,n);
Ｂ3=ｓplｉnｅ（n１，b3，ｎ）；
B4=spｌｉne（n1,b４,n);
pg=7； 　
for ｉg=[1、00，0、７９３]    
　  Ｕa=０、377*ｒ*n/iｇ/i0；  
   Pf＝G＊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;
eｎd
axｉｓ（[0 １00 12 28])；
texｔ（４０,14，'四挡’);
text(80，16,＇五档＇）；
title（＇最高挡与次高挡等速百公里油耗曲线＇）,ｘｌabeｌ('Uａ/(km/h）'），yｌａｂel（’百公里油耗Qs/(L／１００km）'）;

3
主函数：
i0=[５、17，5、４3，５、8３,6、17，63３］；
ｆｏr　ｉ＝1:1:5
   y（i）=jiａsｕshijiａn（i0(ｉ））；
ｅｎd
y；
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７、５,’i＿0=5、17’);
texｔ（13、６，２7、6，＇ｉ_0=５、4３＇)；
text(13、９，27、２，’i＿0=５、８3＇）；
teｘｔ（１4、１，26、6，’ｉ＿０＝6、17'）;
teｘt（１4、2，26、2,'i_0=6、33＇）；
tiｔle（’燃油经济性—加速时间曲线＇)；
xｌabeｌ(’百公里油耗（L/100km）'）；
ylaｂel('加速时间s’）;
子函数1
fｕncｔioｎ　y=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:１：4      　    　             　   　 　　　  　
  　uamax(i)=ｃhesｕ（nmax，r，ig（ｉ）,i０)； 　     　 　
　 　 uamiｎ(i）＝ｃhesu（nmin，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０、/（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+１;　      　
     　　 enｄ
 　  enｄ
 　  n2=1；
　   for j1=1：1:1０0
  　  　　　if ｕa（4，ｊ1)＞maｘ（ua（3,：))&＆ua（4,ｊ1）<=70;  
  　     　　  ｔｅｍｐ3（１,n2)=ua（4，j1)/3、6;  　   　 　　
　         　  temp3（2，n2）=1、/a（４，j1）;  　   　
        　　　 ｎ２＝n2＋1;    　　  
        eｎd
 　 　ｅ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
子函数2
functｉon ［ ｕa　］ =　chesｕ( n,ｒ，ｉg,ｉ0 )
ua=0、3７7＊r*ｎ/i０、/iｇ;
eｎｄ
子函数3
fuｎcｔion ［ ｑ　] ＝ jiasu( umax,uｍiｎ，ig,a，ua０,i０　)
global　f Ｇ　ＣDA yiｔa m　r　If Iw1 Iw2 pg　B0　B1 B2 Ｂ3　B4 n；  
ua1＝ｕｍin:1:ｕmax; 　      　  
delta=1＋(Iw1+Iw2)/(m＊ｒ＾2)+（Ｉｆ＊iｇ^2*i０^2＊yiｔa)/(ｍ＊r^２)； 　
P0G＊f、*ｕa0、/36００+CDA、＊ua0、^3／761４0+(delta＊m、＊ｕa0／36０0）＊a）/yita;
P＝（G*f、*ｕａ1/3600+CDA、＊uａ１、^3/7６1４0+（ｄｅlｔa＊m、＊ua１/360０)*ａ）/yｉｔa;
dt=1/（３、6＊a) 　      　 　　  
for ｉ=1：1：３401     　　  
　　  b0（ｉ）=B０（ｉ）+Ｂ１(i)*Ｐ0（i)+B2（i）＊（P0(i)）^2＋B3(i)*(P0(i）)^3+B４（ｉ）*（P０（i))＾４;
ｅnd
b1＝intｅrp1（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
子函数4
ｆunｃtioｎ [ y ］ = ｚhuanjｕ(　n )
y=－１9、3１3+2９5、27、*（n、／1000)－１６５、44、＊（n、/1000)、^２+40８74、＊（n、/１000）、^3－3、8４45、＊（ｎ、/1000）、^4;
eｎd
子函数5
fｕncｔiｏn　[ n ］ =　ｚhuansu（ ua,ｒ，iｇ，ｉ0 )
n=ig＊i０、＊ua、/（0、377＊r);
end
子函数6
fｕnctiｏn [ b　] ＝ ｙouhａo（ ｉ0 )
gｌoｂal f G　CDA yｉｔa m ｒ Iｆ　Iw1 Iw2 pg Ｂ０ B1 Ｂ２ B3 B4 n  
ig=[6、09,3、09，1、7１，1、00］；r=０、36７;
yiｔａ＝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 201２ 26０3 ３0０6 3403 ３8０４］;
Ｂ０0=[1326、8　1３５４、7 1284、4 1122、9 １１41、０　１051、2 １23３、9 １129、7]；
B10=［－416、４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、711１3];

Ｂ30＝[—5、8629 -2、0553 －0、51184　-０、18517 －0、091０７７ —0、0５138　-0、04７4４9　－０、000７5２15］;
B40=［０、1７768　0、０43072 0、00681６4 0、001８555 ０、0００689０6 0、０003５032　0、０002８２30 －0、00０0３85６8］；
ｎ=6００：１：400０；
Ｂ0=spline（ｎ0,Ｂ00，ｎ);　B１=spline（n0,B1０，n）； B2＝splｉｎe（n0,B2０，n)；
B3=spliｎe（n0，B30，n）； B4=spｌｉｎe（n０，B4０，n)；

ua４＝0、377＊r、＊n、/（i0*ｉg（4))； 　  
F4=f＊G＋ＣDA＊（ua４、^2)/２１、15；      
P＿fｗ４=F4、＊ua4、/(yiｔa＊3、６＊1000）;    　　
foｒ i=1:１:3４01　    　
b4（i)=Ｂ0(ｉ）+B１（i）＊Ｐ_fw4(i)+B2（i）*(Ｐ_ｆw4(i）)＾2+B3（ｉ)*（P＿fｗ4（ｉ)）^3+B4（i）*(P_ｆｗ4（i))^４;
end　
pg＝7、06； 　  　     &nbsp; 　    
ua4_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； 　  　 　      　 　 　      　
Q4all=(suｍ(Ｑ4＿m）+Q4＿ａ1+Q4_a2+Ｑ_i)＊100/s；
ｂ=Q4all；
end
子函数7
functｉon　［ p ] =　qiｕｊi（ x0,y0 )
n０=sｉze(x0）； n＝n０（2）；　
x＝lｉnspａce（x０(1）,ｘ0（n)，200) ；
y=spline（x0，ｙ０，ｘ）； 　 　　
plot（ｘ,ｙ)；
p=trａｐz(x,y) 　；
enｄ

4、3
beta=０、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*ｈg1)）；
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,’g＇，z,ｆ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,’φ_ｆ（空车）&#39;)；
ｘlaｂeｌ('制动强度z/ｇ’);
ｙlaｂｅl(’利用附着系数φ’）;
titlｅ（'利用附着系数与制动强度得关系曲线＇）；
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,’r’,ｆaｉ,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(0、7，０、6，’空载’);
xlaｂｅｌ('附着系数φ’）;
ylabel(＇制动效率')；
tｉｔle（’前、后制动效率曲线')；
axｉs（[0 1　0　１]）；
ｆ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+ｕ0^2/2５、９２/abmax1；
ｓ2＝1/3、6＊(ｔ２1+t22／2）＊u０＋ｕ0^2/2５、92/abｍax2；
ｄisp（’空载时得制动距离＇），dｉsp（s1）;diｓｐ（＇满载时得制动距离＇），ｄ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+ｕ0^2/25、９2/aｂｍax１；
s２=1/３、6＊（ｔ2１＋t2２/２）*u0+ｕ０^2/２5、92／abmaｘ２;
ｄiｓｐ（’前部管路损坏时空载制动距离＇);ｄisp（s１）；diｓｐ（’前部管路损坏时满载制动距离'）；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）＊u0＋u0＾2/2５、9２／abｍax2；
dｉsp（'后部管路损坏时空载制动距离＇）；disp（s1）；disp(＇后部管路损坏时满载制动距离’）;diｓp(s２)；

5、11
m=1８1８、２;Ｉz=3８8５;L＝３、048；a＝1、463；ｂ=1、5８５；k1=－62618;k2=-11０１85; ｉ=20;g=9、8；R０＝１5;u1=３0、56； K＝m＊(a／k2-b/k1）/L＾２;　
Uch＝(1／K)^(1/2)；
dｉsp(’稳定性因数(s^２/ｍ＾2)K=’）； ｄiｓｐ（Ｋ)；　　
diｓp（’特征车速（ｍ/s)Ｕｃh='）； disp（Uch）；　u=０：０、05：３0；
S=ｕ、/（L＊（１＋K＊u、^2））;
plot(u，S）；
ｔitle(’汽车稳态横摆角速度增益曲线’）;ｘlabel(＇车速u(m/s)')；
ｙｌabel（’稳态横摆角速度增益');
dｉｓｐ（’u=2２、35m/ｓ时,转向灵敏度为’)；disp（S(448));
SM=k2/（k1+k2)—a/L;ａy＝０、4*ｇ;Ａ=K＊ay＊L;B＝L/R0;Ｒ=L/（Ｂ－Ａ);
C=R／Ｒ０；
ｄisp('静态储备系数S、Ｍ、=’）;disp(SM);
disp('侧向加速度为0、４ｇ时前、后轮侧偏角绝对值之差(rad）?a１-a2='）;ｄisp(A)；
diｓp(’侧向加速度为0、4g时转弯半径比值R/R0='）;ｄiｓp(Ｃ)；
W0=L／u1＊(k1*k２/（m＊Iz）*（1+K＊ｕ1^2)）^(1/2);
D=（—ｍ＊（k１＊a^２+k2*b^２）—Iz*（k１+k２)）/（2＊Ｌ＊（ｍ＊Ｉｚ＊k１＊k2＊（1＋Ｋ*u1^2））^（1/2）);
ｔ=ａｔａn（(1－D^2)^（1/2)／(-m＊u1＊a＊W0/(L*k2）－Ｄ))/（W0*(1－D^2）＾(1/２））;
E＝ａtaｎ（（1-D^2）＾（1/2）/D)/(W0＊（1－D＾2）^（１/2）)+ｔ；
ｄｉsｐ（’车速u=３0、56m／ｓ时得瞬态响应参数分别为：')；dｉsp（'横摆角速度波动得固有（圆）频率(ｒaｄ)为？＇）；ｄｉsp（W0）；
disp（＇阻尼比为’);ｄｉsｐ(D）；
diｓp（’反应时间(s）为’）；ｄisp(ｔ）；
diｓp（'峰值反应时间(s）为＇);diｓp(E）;