【正文】 在論文的撰寫過程中老師們給予我很大的幫助，幫助解決了不少的難點，使得論文能夠及時完成，這里一并表示真誠的感謝。老師們認(rèn)真負(fù)責(zé)的工作態(tài)度，嚴(yán)謹(jǐn)?shù)闹螌W(xué)精神和深厚的理論水平都使我收益匪淺。最后，我要特別感謝我的導(dǎo)師趙達睿老師、和研究生助教熊偉麗老師。四年的風(fēng)風(fēng)雨雨，我們一同走過，充滿著關(guān)愛，給我留下了值得珍藏的最美好的記憶。感謝老師四年來對我孜孜不倦的教誨，對我成長的關(guān)心和愛護。從這里走出，對我的人生來說，將是踏上一個新的征程，要把所學(xué)的知識應(yīng)用到實際工作中去。最后，我要感謝我的父母對我的關(guān)系和理解，如果沒有他們在我的學(xué)習(xí)生涯中的無私奉獻和默默支持，我將無法順利完成今天的學(xué)業(yè)。其次，我要感謝大學(xué)四年中所有的任課老師和輔導(dǎo)員在學(xué)習(xí)期間對我的嚴(yán)格要求，感謝他們對我學(xué)習(xí)上和生活上的幫助，使我了解了許多專業(yè)知識和為人的道理，能夠在今后的生活道路上有繼續(xù)奮斗的力量。從他身上，我學(xué)到了許多能受益終生的東西。首先，我要特別感謝我的知道郭謙功老師對我的悉心指導(dǎo)，在我的論文書寫及設(shè)計過程中給了我大量的幫助和指導(dǎo)，為我理清了設(shè)計思路和操作方法，并對我所做的課題提出了有效的改進方案。這期間凝聚了很多人的心血，在此我表示由衷的感謝。本次畢業(yè)設(shè)計是對我大學(xué)四年學(xué)習(xí)下來最好的檢驗。首先非常感謝學(xué)校開設(shè)這個課題，為本人日后從事計算機方面的工作提供了經(jīng)驗，奠定了基礎(chǔ)。（保密論文在解密后遵守此規(guī)定）畢業(yè)設(shè)計（論文）使用授權(quán)聲明本人完全了解濱州學(xué)院關(guān)于收集、保存、使用畢業(yè)設(shè)計（論文）的規(guī)定。 對本文的研究做出重要貢獻的個人和集體均已在文中以明確方式標(biāo)明。論文密級：□公開 □保密（___年__月至__年__月）(保密的學(xué)位論文在解密后應(yīng)遵守此協(xié)議)作者簽名：_______ 導(dǎo)師簽名：______________年_____月_____日 _______年_____月_____日獨 創(chuàng) 聲 明本人鄭重聲明：所呈交的畢業(yè)設(shè)計(論文)，是本人在指導(dǎo)老師的指導(dǎo)下，獨立進行研究工作所取得的成果，成果不存在知識產(chǎn)權(quán)爭議。本人完全意識到本聲明的法律結(jié)果由本人承擔(dān)。盡我所知，除文中已經(jīng)特別注明引用的內(nèi)容和致謝的地方外，本論文不包含任何其他個人或集體已經(jīng)發(fā)表或撰寫過的研究成果。law9=20*log10(aw9/a0)。endaw9=sqrt(aw9)。z=50 w_k=。 end if =zamp。amp。z4 w_k=z/4。 end if 2=zamp。amp。for i=1:21 z=f_z9(i,1)。aw9=0。三分之一倍頻程的加速度均方值39。)。ylabel(39。頻率39。),grid。loglog(f_z9,ai9,39。 end ai9(i,1)=sqrt(ai9(i,1))。 end else ai=0。(abs(f_uf0)delta_f)) ai=Pya9(j)*(f_uf_l)。 %頻率值距1/3倍頻帶上限、下限均小于delta_f elseif((abs(f0f_l)delta_f)amp。amp。(abs(f_uf0)delta_f)) ai=Pya9(j)*abs(f0f_l)。(f0=f_u)) %頻率值距1/3倍頻帶下限小于delta_f，但距上限大于delta_f if((abs(f0f_l)delta_f)amp。 if((f0=f_l)amp。 f_z9(i,1)=f_l*2^(1/6)。%中心頻率for i=1:21 f_l=**^(i1)。ai9=zeros(21,1)。hold on。9DOF Model 人體座椅垂直振動39。)。ylabel(39。頻率[Hz]39。off39。YMinorGrid39。off39。XMinorGrid39。off39。YMinorTick39。off39。XMinorTick39。grid on。Linewidth39。k39。subplot(2,1,1)。endPya9=Pya9/(2*n_n1)。)。for i=2:(2*n_n1) [Py09,f]=periodogram(xac9((i1)*nb+1:((i1)*nb+nf)),hann(nf),nf,fs,39。)。 %功率譜密度分析——人體座椅俯系統(tǒng)[Py09,f]=periodogram(xac9(1:nf),hann(nf),nf,fs,39。 %每5s一段，單個數(shù)據(jù)樣本個數(shù)nb=nf/2。fs=100。9DOF人體座椅垂直加速度時域變化39。)。ylabel(39。時間(S)39。axis auto。k39。figure(3)。 x(:,i)=x(:,n)+h/2*((kt.*Q(:,n)+W*x(:,n))+(kt.*Q(:,i)+W*xL))。 xL=x(:,n)+h*(kt.*Q(:,n)+W*x(:,n))。 %定義初始時間%改進的歐拉算法for i=2:nt %nt=9953 n=i1。 %步長，時間間隔x(:,1)=zeros(18,1)。 %數(shù)據(jù)個數(shù)，即9953nf=500。K C]。%p為8*8的矩陣，其中每個數(shù)都是0q=eye(9)。 %將qf的值賦給第9行，：指9953Q(11,:)=pr。Q=zeros(18,size(pf,2))。kt(10,1)=k1/m1。 C(j,i)=C(i,j)endC=C./m。C(9,9)=c9。C(8,8)=a^2*c7+b^2*c8。C(7,8)=a*c7+b*c8。C(6,6)=c^2*c5+d^2*c6。C(5,5)=c5+c6。C(4,7)=l3*c7+l6*c8。C(4,5)=l4*c5+l5*c6。C(3,8)=a*c7b*c8。C(3,6)=c*c5d*c6。C(3,4)=l1*c3+l2*c4l3*c7l4*c5l5*c6l6*c8。C(2,4)=l2*c4。C(2,2)=c4。C(1,3)=c3。C=zeros(9,9)。for i=1:9 j=i:9。k(8,8)=a^2*k7+b^2*k8。k(7,8)=a*k7+b*k8。k(6,6)=c^2*k5+d^2*k6。k(5,5)=k5+k6。k(4,7)=l3*k7+l6*k8。k(4,5)=l4*k5+l5*k6。k(3,8)=a*k7b*k8。k(3,6)=c*k5d*k6。k(3,4)=l1*k3+l2*k4l3*k7l4*k5l5*k6l6*k8。k(2,4)=l2*k4。k(2,2)=k2+k4。k(1,3)=k3。k=zeros(9,9)。m(8,1:9)=i8。m(6,1:9)=i6。m(4,1:9)=i4。m(2,1:9)=m2。m=zeros(9,9)。c=(+)。a=()。l5=()。l3=(+)。l1=。c8=+3。c6=+3。c4=+5。k9=18121。k7=+4。k5=+6。k3=+5。k1=+6。i6=270。m9=51。m5=1653。m2=2140。 pf=p((l+1):10000)。 %計算前后輪時間差l=round(td/)。 %采用對數(shù)坐標(biāo)作圖%根據(jù)軸距l(xiāng)w和車速V（m/s）生成后輪輸入時域數(shù)據(jù)lw=。39。 Sq(i)=Gq_n0**u/(fb(i)^2)。)。 %標(biāo)注ylabel(39。頻率(Hz)39。xlim([,50])。)。loglog(f(1:500/2),Pya(1:500/2),39。 %求平均figure(2)。 %依次每隔250個點（）求一次功率譜 Pya=Pya+Py0。onesided39。 %利用周期圖法求前500個點（即前5s）的功率譜Pya=Py0。onesided39。)。title(39。路面不平度位移(m)39。)。xlabel(39。auto39。 %給圖形加上柵格線xlim([0,100])。plot(a,p)。q=0。 %將時間變量賦值給變量ap(j)=q。 i=i+1。 a_k=sqrt(2*Gq_f*delta_f)。 for n=1:1:1000 f_k=f_l+()*delta_f。 %速度換算md=2*pi*rand(1001,1)。 %選取B級路面不平度等級n0=。%delta_f=(f_uf_l)/ns。%給最小頻率賦值f_u=50。 %定義初始變量j=1。clf。最后，向文中引用到其學(xué)術(shù)論著及研究成果的學(xué)術(shù)前輩與同行們致謝！再次向敬愛的老師表示衷心的感謝！參考文獻[1]余志生．汽車?yán)碚摚ǖ谒陌妫M]．北京：機械工業(yè)出版社，2008[2]樓順天，姚若玉，（第二版）[M].西安：西安電子科技大學(xué)出版社，2007[3]（第二版）[M].武漢：華中科技大學(xué)出版社，2004[4]呂同富，康兆敏，（第二版）.北京：清華大學(xué)出版社，2008[5]王立秋，魏煥彩，:山東大學(xué)出版社，2004[6]徐中明，張志飛，黃澤好，肖建伯，[J].汽車工程，2007[7][J].科技資訊，2009[8] 李朝峰，王得剛，劉杰，聞邦椿 . 11自由度非線性車輛平順性模型的研究及應(yīng)用. 東北大學(xué)學(xué)報，2009（6）[9] 穆立茂，黃海英，張靖. LabVIEW下貨車平順性仿真分析[N].汽車電子報,2009(8)[10] .吉林大學(xué)博士學(xué)位論文，2009（6）[11] 徐中明，譚海偉，張志飛，陳雪梅，蘇周成. 全地形車平順性客觀評價，2008[12] 陳珂，胡修富，徐內(nèi)恒. HFC3081KR1T型自卸汽車抖動問題的研究. 合肥工業(yè)大學(xué)學(xué)報,2009（11）[13] 陳玲，張文明，楊耀東. ，2004（9）[14] 盧海峰，李明，李軍. ，2005[15] 王秉剛. ，2007[16] 李玉生. 我國汽車技術(shù)現(xiàn)狀分析研究[N]. 汽車,2008(5)[17]賀巖松，[J].中南林學(xué)院學(xué)報，2006(6) [18][M].北京：高等教育出版社，1990[19][M].北京：電子工業(yè)出版社，2004[20][M].重慶：重慶大學(xué)出版社，2002[21]張志飛. ，2008[22]徐中明，崔錦羽，張志飛，牟笑靜. 全地域車平順性仿真分析[J].重慶大學(xué)學(xué)報，2009[23]羅衛(wèi)東..汽車“人體—座椅”系統(tǒng)模型與動力學(xué)參數(shù)選擇[J].山地農(nóng)業(yè)生物學(xué)報，1999[24]，2009[25]，2006[26]李孟良,聶彥鑫,過學(xué)迅,，2009[27]張朝杰，郭志軍,[J].河南科技大學(xué)學(xué)報，2009[28]徐中明,盧少波,文俐,張志飛. 基于MATLAB的摩托車舒適性仿真分析[J].計算機仿真，2006附錄A：路面激勵的matlab程序clear。至此，這次畢業(yè)設(shè)計也將告以段落，但老師的教誨卻讓人終生難忘，通過這次畢業(yè)設(shè)計，不但使我學(xué)到了知識，也讓我學(xué)到了許多的道理，總之是受益匪淺。致謝首先向王軍老師表示衷心的感謝，在這做畢業(yè)設(shè)計的時間里，他為我們的成長和進步做出了貢獻。4．在建立矩陣形式的微分方程的基礎(chǔ)上運用改進的歐拉算法，然后根據(jù)平順性評價方法依次求出各種速度下、各種路面等級下的加權(quán)加速度均方根值以及加權(quán)振級。3．根據(jù)1/2汽車模型建立9自由度動力學(xué)模型，模型中包括了駕駛室、發(fā)動機、人體座椅，能夠方便對平順性進行評價。本文重點研究了汽車的時域路面激勵模型，采用三角級數(shù)進行模擬，以離散譜逼近目標(biāo)隨機路面模型。本文主要內(nèi)容可以概括為以下幾點：1．參考國內(nèi)外大量文獻，對汽車平順性的發(fā)展及研究方法進行了論述，并對MATLAB軟件做了基本介紹，針對汽車平順性仿真的關(guān)鍵問題：路面激勵建模、振動模型建模和仿真方法進行了綜述和分析。綜上所述，提高平順性的最好方法是減少前懸架剛度、后懸架剛度、人體座椅系統(tǒng)剛度，增加人體座椅系統(tǒng)的阻尼。表44到表410分別列舉了人體座椅系統(tǒng)加權(quán)加速度均方根值隨前懸架剛度、后懸架剛度、前懸架阻尼、后懸架阻尼、駕駛室前懸架剛度、人座椅系統(tǒng)剛度、人座椅系統(tǒng)阻尼變化的變化率。通過依次調(diào)節(jié)車輛前后懸架剛度和阻尼及駕駛室前懸剛度參數(shù)，利用MATLAB對其進行多次仿真，研究車輛參數(shù)變化后人體座