【正文】 。C(3,5)=c5c6。C(3,6)=c*c5d*c6。C(3,7)=c7c8。C(3,8)=a*c7b*c8。C(4,4)=l1^2*c3+l2^2*c4+l3^2*c7+l4^2*c5+l5^2*c6+l6^2*c8。C(4,5)=l4*c5+l5*c6。C(4,6)=c*l4*c5+d*l5*c6。C(4,7)=l3*c7+l6*c8。C(4,8)=a*l3*c7+b*l6*c8。C(5,5)=c5+c6。C(5,6)=c*c5+d*c6。C(6,6)=c^2*c5+d^2*c6。C(7,7)=c7+c8+c9。C(7,8)=a*c7+b*c8。C(7,9)=c9。C(8,8)=a^2*c7+b^2*c8。C(9,9)=c9。C(9,9)=c9。for i=1:9 j=i:9。 C(j,i)=C(i,j)endC=C./m。kt=zeros(18,1)。kt(10,1)=k1/m1。kt(11,1)=k2/m2。Q=zeros(18,size(pf,2))。 %產(chǎn)生(16*size)全零矩陣，其中size為二維矩陣qf（1*9953）第二維的尺寸（即列數(shù)）,Q為（16*9553）的矩陣Q(10,:)=pf。 %將qf的值賦給第9行，：指9953Q(11,:)=pr。 %將qr的值賦給第10行，：指9953p=zeros(9)。%p為8*8的矩陣，其中每個數(shù)都是0q=eye(9)。 %q為8*8的單位矩陣，其中對角線上都是1W=[p q。K C]。 %p q k c均為9*9的矩陣，w為18*18的矩陣nt=size(pf,2)。 %數(shù)據(jù)個數(shù)，即9953nf=500。h=1/100。 %步長，時間間隔x(:,1)=zeros(18,1)。%x為18*1的矩陣，且每個數(shù)都是0t(1)=0。 %定義初始時間%改進的歐拉算法for i=2:nt %nt=9953 n=i1。 t(i)=h*i。 xL=x(:,n)+h*(kt.*Q(:,n)+W*x(:,n))。 %當n=1時，x（：，n）表示x矩陣（18*9553）第n列的列矩陣，kt為16*16的矩陣，Q（；，n）表示矩陣（16*9553）的第n列的列矩陣。 x(:,i)=x(:,n)+h/2*((kt.*Q(:,n)+W*x(:,n))+(kt.*Q(:,i)+W*xL))。 endxac9=W(18,:)*x。figure(3)。plot(t,xac9,39。k39。),grid。axis auto。xlabel(39。時間(S)39。)。ylabel(39。垂直加速度 [rad/s^2]39。)。title(39。9DOF人體座椅垂直加速度時域變化39。)。fs=100。 %分析頻率=1/，nf=500。 %每5s一段，單個數(shù)據(jù)樣本個數(shù)nb=nf/2。 %自譜驗算時的參數(shù) %重復點數(shù)n_n=19。 %功率譜密度分析——人體座椅俯系統(tǒng)[Py09,f]=periodogram(xac9(1:nf),hann(nf),nf,fs,39。onesided39。)。Pya9=Py09。for i=2:(2*n_n1) [Py09,f]=periodogram(xac9((i1)*nb+1:((i1)*nb+nf)),hann(nf),nf,fs,39。onesided39。)。 Pya9=Pya9+Py09。endPya9=Pya9/(2*n_n1)。figure(4)。subplot(2,1,1)。loglog(f(1:nf/2),Pya9(1:nf/2),39。k39。,39。Linewidth39。,2)。grid on。set(gca,39。XMinorTick39。,39。off39。,39。YMinorTick39。,39。off39。,39。XMinorGrid39。,39。off39。,39。YMinorGrid39。,39。off39。)xlim([ 50])xlabel(39。頻率[Hz]39。)。ylabel(39。加速度功率譜密度(m^2/s^3)39。)。title(39。9DOF Model 人體座椅垂直振動39。)。hold on。delta_f=。ai9=zeros(21,1)。f_z9=zeros(21,1)。%中心頻率for i=1:21 f_l=**^(i1)。 f_u=f_l*。 f_z9(i,1)=f_l*2^(1/6)。%求中心頻率 %若頻率大于50HZ,循環(huán)結(jié)束 %在一個三分之一倍頻程內(nèi)積分 for j=1:251 f0=f(j)。 if((f0=f_l)amp。amp。(f0=f_u)) %頻率值距1/3倍頻帶下限小于delta_f，但距上限大于delta_f if((abs(f0f_l)delta_f)amp。amp。(abs(f_uf0)delta_f)) ai=Pya9(j)*abs(f0f_l)。 %頻率值距1/3倍頻帶上限小于delta_f，但距下限大于delta_f elseif((abs(f_uf0)delta_f)amp。amp。(abs(f_lf0)delta_f)) ai=Pya9(j)*(abs(f0f_u)+delta_f)。 %頻率值距1/3倍頻帶上限、下限均小于delta_f elseif((abs(f0f_l)delta_f)amp。amp。(abs(f_uf0)delta_f)) ai=Pya9(j)*(f_uf_l)。 %頻率值距1/3倍頻帶上限、下限均大于delta_f else ai=Pya9(j)*delta_f。 end else ai=0。 end ai9(i,1)=ai9(i,1)+ai。 end ai9(i,1)=sqrt(ai9(i,1))。 end subplot(2,1,2)。loglog(f_z9,ai9,39。k39。),grid。xlim([ 50])xlabel(39。頻率39。)。ylabel(39。加速度均方根值39。)。title(39。三分之一倍頻程的加速度均方值39。)。aw9=0。w_k=0。for i=1:21 z=f_z9(i,1)。 if =zamp。amp。z2 w_k=。 end if 2=zamp。amp。z4 w_k=z/4。 end if 4=zamp。amp。z w_k=1。 end if =zamp。amp。z=50 w_k=。 end aw9=aw9+ai9(i,1)^2*w_k^2。endaw9=sqrt(aw9)。a0=。law9=20*log10(aw9/a0)。學位論文原創(chuàng)性聲明 本人鄭重聲明：所呈交的學位論文，是本人在導師的指導下進行的研究工作所取得的成果。盡我所知，除文中已經(jīng)特別注明引用的內(nèi)容和致謝的地方外，本論文不包含任何其他個人或集體已經(jīng)發(fā)表或撰寫過的研究成果。對本文的研究做出重要貢獻的個人和集體，均已在文中以明確方式注明并表示感謝。本人完全意識到本聲明的法律結(jié)果由本人承擔。學位論文作者（本人簽名）： 年 月 日學位論文出版授權(quán)書本人及導師完全同意《中國博士學位論文全文數(shù)據(jù)庫出版章程》、《中國優(yōu)秀碩士學位論文全文數(shù)據(jù)庫出版章程》(以下簡稱“章程”)，愿意將本人的學位論文提交“中國學術(shù)期刊（光盤版）電子雜志社”在《中國博士學位論文全文數(shù)據(jù)庫》、《中國優(yōu)秀碩士學位論文全文數(shù)據(jù)庫》中全文發(fā)表和以電子、網(wǎng)絡(luò)形式公開出版，并同意編入CNKI《中國知識資源總庫》，在《中國博碩士學位論文評價數(shù)據(jù)庫》中使用和在互聯(lián)網(wǎng)上傳播，同意按“章程”規(guī)定享受相關(guān)權(quán)益。論文密級：□公開 □保密（___年__月至__年__月）(保密的學位論文在解密后應(yīng)遵守此協(xié)議)作者簽名：_______ 導師簽名：______________年_____月_____日 _______年_____月_____日獨 創(chuàng) 聲 明本人鄭重聲明：所呈交的畢業(yè)設(shè)計(論文)，是本人在指導老師的指導下，獨立進行研究工作所取得的成果，成果不存在知識產(chǎn)權(quán)爭議。盡我所知，除文中已經(jīng)注明引用的內(nèi)容外，本設(shè)計（論文）不含任何其他個人或集體已經(jīng)發(fā)表或撰寫過的作品成果。對本文的研究做出重要貢獻的個人和集體均已在文中以明確方式標明。本聲明的法律后果由本人承擔。 作者簽名: 二〇一〇年九月二十日畢業(yè)設(shè)計（論文）使用授權(quán)聲明本人完全了解濱州學院關(guān)于收集、保存、使用畢業(yè)設(shè)計（論文）的規(guī)定。本人愿意按照學校要求提交學位論文的印刷本和電子版，同意學校保存學位論文的印刷本和電子版，或采用影印、數(shù)字化或其它復制手段保存設(shè)計（論文）；同意學校在不以營利為目的的前提下，建立目錄檢索與閱覽服務(wù)系統(tǒng)，公布設(shè)計（論文）的部分或全部內(nèi)容，允許他人依法合理使用。（保密論文在解密后遵守此規(guī)定）作者簽名: 二〇一〇年九月二十日致 謝時間飛逝，大學的學習生活很快就要過去，在這四年的學習生活中，收獲了很多，而這些成績的取得是和一直關(guān)心幫助我的人分不開的。首先非常感謝學校開設(shè)這個課題，為本人日后從事計算機方面的工作提供了經(jīng)驗，奠定了基礎(chǔ)。本次畢業(yè)設(shè)計大概持續(xù)了半年，現(xiàn)在終于到結(jié)尾了。本次畢業(yè)設(shè)計是對我大學四年學習下來最好的檢驗。經(jīng)過這次畢業(yè)設(shè)計，我的能力有了很大的提高，比如操作能力、分析問題的能力、合作精神、嚴謹?shù)墓ぷ髯黠L等方方面面都有很大的進步。這期間凝聚了很多人的心血，在此我表示由衷的感謝。沒有他們的幫助，我將無法順利完成這次設(shè)計。首先，我要特別感謝我的知道郭謙功老師對我的悉心指導，在我的論文書寫及設(shè)計過程中給了我大量的幫助和指導，為我理清了設(shè)計思路和操作方法，并對我所做的課題提出了有效的改進方案。郭謙功老師淵博的知識、嚴謹?shù)淖黠L和誨人不倦的態(tài)度給我留下了深刻的印象。從他身上，我學到了許多能受益終生的東西。再次對周巍老師表示衷心的感謝。其次，我要感謝大學四年中所有的任課老師和輔導員在學習期間對我的嚴格要求，感謝他們對我學習上和生活上的幫助，使我了解了許多專業(yè)知識和為人的道理，能夠在今后的生活道路上有繼續(xù)奮斗的力量。另外，我還要感謝大學四年和我一起走過的同學朋友對我的關(guān)心與支持，與他們一起學習、生活，讓我在大學期間生活的很充實，給我留下了很多難忘的回憶。最后，我要感謝我的父母對我的關(guān)系和理解，如果沒有他們在我的學習生涯中的無私奉獻和默默支持，我將無法順利完成今天的學業(yè)。四年的大學生活就快走入尾聲，我們的校園生活就要劃上句號，心中是無盡的難舍與眷戀。從這里走出，對我的人生來說，將是踏上一個新的征程，要把所學的知識應(yīng)用到實際工作中去?；厥姿哪辏〉昧诵┰S成績，生活中有快樂也有艱辛。感謝老師四年來對我孜孜不倦的教誨，對我成長的關(guān)心和愛護。學友情深，情同兄妹。四年的風風雨雨，我們一同走過，充滿著關(guān)愛，給我留下了值得珍藏的最美好的記憶。在我的十幾年求學歷程里，離不開父母的鼓勵和支持，是他們辛勤的勞作，無私的付出，為我創(chuàng)造良好的學習條件，我才能順利完成完成學業(yè)，感激他們一直以來對我的撫養(yǎng)與培育。最后，我要特別感謝我的導師趙達睿老師、和研究生助教熊偉麗老師。是他們在我畢業(yè)的最后關(guān)頭給了我們巨大的幫助與鼓勵，給了我很多解決問題的思路，在此表示衷心的感激。老師們認真負責的工作態(tài)度，嚴謹?shù)闹螌W精神和深厚的理論水平都使我收益匪淺。他無論在理論上還是在實踐中，都給與我很大的幫助，使我得到不少的提高這對于我以后的工作和學習都有一種巨大的幫助，感謝他耐心的輔導。在論文的撰寫過程中老師們給予我很大的幫助，幫助解決了不少的難點，使得論文能夠及時完成，這里一并表示真誠的感謝。43