【正文】 平均斜率k*提供一種算法，由式便相應地得到一種微分方程地數(shù)值計算公式。而改進歐拉公式即是利用了與兩個點的斜率與取算術(shù)平均值作為平均斜率k*的近似值：其中是通過已知信息yi利用歐拉公式得到的。因此它啟發(fā)我們，如果設法在［xi，xi＋１］上多預報幾個點的斜率值，然后將他們加權(quán)平均作為k*的近似值，則有可能構(gòu)造出更高精度的計算公式，此則是龍格—庫塔法的基本思想。double t,y[],d[]。}調(diào)用參數(shù)說明：t—雙精度實型變量，積分的起點t0y—雙精度實型一維數(shù)組，長度為n。k—整型變量，積分的步數(shù)h—雙精度實型變量，積分的步長z—雙精度實型二維數(shù)組，體積為n*k以上仿真的基本參數(shù)：M1=1400Kg；M2=70Kg；V1=；V2=；A0=；P1=4Pa；K2=；K1=；6．客車空氣懸架系統(tǒng)的仿真分析空氣懸架的特性在很大程度上取決于空氣彈簧的彈性特性?？諝鈴椈傻妮d荷一位移曲線形狀呈反“S”形，作該曲線上某點的切線便得到了該點的剛度。為了便于研究，簡化了一下，就取曲線上的兩個剛度值來研究，一個是剛度較大的情況下，一個是剛度較小的情況下，來分析車身、車輪、的位移速度加速度圖象汽車在行駛過程中振動時候，其阻尼比一般有三種情況：0ξ1, ξ=1, ξ1,其中：1. 0ξ1時，此時符合汽車小阻尼（欠阻尼）的振動，汽車在大部分情況下都是按照這種情況下振動的。3. ξ1時，此時為過阻尼，汽車振動不明顯。 數(shù)據(jù)的處理利用設計出來的C語言程序計算出來數(shù)據(jù)處理后分別將車輪位移、車輪速度、車身位移、車身速度、車身加速度、車輪加速度轉(zhuǎn)化為test. txt、test2. txt、test3. txt、test4. txt、test5. txt、以供Matlab繪圖時候調(diào)用。在調(diào)用時候，將每個數(shù)據(jù)轉(zhuǎn)化為矩陣形式，與時間參數(shù)一一對應（即t=0::）進行繪圖，從而的出車輪位移圖、車輪速度圖、車身位移圖、車身加速度圖、車身速度圖。通過改變汽車阻尼比進行對比，具體結(jié)果如下：圖1 阻尼比較小的情況，車身位移圖圖2 阻尼比較小的情況下，車身速度圖 圖3 阻尼比較小的情況下，車輪位移圖 圖4 阻尼比較小的情況下，車身加速度圖圖5 阻尼比較小的情況下，車輪加速度圖圖6 阻尼比較大時，車身加速度圖。 仿真結(jié)果分析通過對以上各圖形的相互對比分析，汽車空氣懸架的阻尼比在處于 0ξ1情況下，當車輪受到地面階躍輸入的時候，車身的位移和速度不是隨著地面的階躍輸入突然變化，不是突然變到最大后，車身速度馬上就變?yōu)榱?。在這種情況下，地面的變化情況對車身的沖擊比較小，從而達到汽車平順性和駕駛舒適性的目的。所以汽車在大部分情況下，阻尼比一般處于0ξ1情況下。由于時間倉促和客觀條件限制，再加上作者水平有限，所以在本文研究工作基礎上還有諸多可以完善之處。（2）在空氣懸架控制系統(tǒng)結(jié)構(gòu)及控制算法的設計和研究上，國內(nèi)大多停留在計算機仿真上。 致 謝行文至此，我的這篇論文已接近尾聲；歲月如梭，我四年的大學時光也即將敲響結(jié)束的鐘聲。感謝學院關(guān)心我們畢業(yè)論文的老師和領(lǐng)導們，學院的計算機實驗室給我的課題帶來了很多方便；感謝在班里同學和朋友，感謝你們在我遇到困難的時候幫助我，給我支持和鼓勵，感謝你們。宋老師給我的鼓勵與指引，使我能夠克服重重困難，將論文完成，在此謹向宋宇老師致以誠摯的謝意和崇高的敬意。 Still can reduce vehicle vibration frequency and wheel dynamic load, make its win good ride fort and sex, handling stability and highspeed driving safety, reducing vehicle of automobi damage. This paper automobile air spring suspension and air suspension system characteristic, the influence of vehicle through a simple calculation, and then the calculation results into the puter and the graphic by puter simulation, we will get results. Keywords: the air suspension, air spring, modeling, puter simulation 附錄 include void grkt1(t,y,n,h,k,z) int n,k。 { extern void grkt1f()。 double a[4],tt,*b,*d。 d=malloc(n*sizeof(double))。 a[1]=a[0]。 a[3]=h。 i=n1。 for (l=1。 l++) { grkt1f(t,y,n,d)。 i=n1。 for (j=0。 j++) { for (i=0。 i++) { y[i]=z[i*k+l1]+a[j]*d[i]。 } tt=t+a[j]。 } for (i=0。 i++) y[i]=b[i]+h*d[i]/。 i=n1。 t=t+h。 free(d)。 }汽車車身速度的編程程序：include include include main(){int i,j。double t,h,y[4],z[4][101]。y[1]=。y[3]=。h=。printf(\n)。exit(0)。i=100。 printf(t=%\n,t)。 for (j=0。j++) printf(y(%d)=%e,j,z[j][i])。 } } void grkt1f(t,y,n,d) int n。 {t=t。 d[0]=y[1]。 d[2]=y[3]。 return。FILE *fp。y[0]=。y[2]=。t=。grkt1(t,y,4,h,101,z)。if ((fp=fopen(,w))==NULL){printf (cannot open file\n)。}for (i=0。i++) {t=i*h。fprintf(fp,%\n,(z[3][i]z[3][I1])/)。j=3。 printf(\n)。 double t,y[],d[]。n=n。 d[1]=200*(y[1]y[3])150*(y[0]y[2])。 d[3]=100*(y[3]y[1])+250*(y[2]y[0])4000*(y[2]sin(t))。 }汽車車輪速度的編程include include include main(){int i,j。double t,h,y[4],z[4][101]。y[1]=。y[3]=。h=。printf(\n)。exit(0)。i=100。 printf(t=%\n,t)。 for (j=0。j++