【正文】 。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的矩陣，其中每個(gè)數(shù)都是0q=eye(9)。 %q為8*8的單位矩陣，其中對(duì)角線上都是1W=[p q。K C]。 %p q k c均為9*9的矩陣，w為18*18的矩陣nt=size(pf,2)。 %數(shù)據(jù)個(gè)數(shù)，即9953nf=500。h=1/100。 %步長(zhǎng)，時(shí)間間隔x(:,1)=zeros(18,1)。%x為18*1的矩陣，且每個(gè)數(shù)都是0t(1)=0。 %定義初始時(shí)間%改進(jìn)的歐拉算法for i=2:nt %nt=9953 n=i1。 t(i)=h*i。 xL=x(:,n)+h*(kt.*Q(:,n)+W*x(:,n))。 %當(dāng)n=1時(shí)，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。時(shí)間(S)39。)。ylabel(39。垂直加速度 [rad/s^2]39。)。title(39。9DOF人體座椅垂直加速度時(shí)域變化39。)。fs=100。 %分析頻率=1/，nf=500。 %每5s一段，單個(gè)數(shù)據(jù)樣本個(gè)數(shù)nb=nf/2。 %自譜驗(yàn)算時(shí)的參數(shù) %重復(fù)點(diǎn)數(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 人體座椅垂直振動(dòng)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é)束 %在一個(gè)三分之一倍頻程內(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)。學(xué)位論文原創(chuàng)性聲明 本人鄭重聲明：所呈交的學(xué)位論文，是本人在導(dǎo)師的指導(dǎo)下進(jìn)行的研究工作所取得的成果。盡我所知，除文中已經(jīng)特別注明引用的內(nèi)容和致謝的地方外，本論文不包含任何其他個(gè)人或集體已經(jīng)發(fā)表或撰寫過的研究成果。對(duì)本文的研究做出重要貢獻(xiàn)的個(gè)人和集體，均已在文中以明確方式注明并表示感謝。本人完全意識(shí)到本聲明的法律結(jié)果由本人承擔(dān)。學(xué)位論文作者（本人簽名）： 年 月 日學(xué)位論文出版授權(quán)書本人及導(dǎo)師完全同意《中國(guó)博士學(xué)位論文全文數(shù)據(jù)庫(kù)出版章程》、《中國(guó)優(yōu)秀碩士學(xué)位論文全文數(shù)據(jù)庫(kù)出版章程》(以下簡(jiǎn)稱“章程”)，愿意將本人的學(xué)位論文提交“中國(guó)學(xué)術(shù)期刊（光盤版）電子雜志社”在《中國(guó)博士學(xué)位論文全文數(shù)據(jù)庫(kù)》、《中國(guó)優(yōu)秀碩士學(xué)位論文全文數(shù)據(jù)庫(kù)》中全文發(fā)表和以電子、網(wǎng)絡(luò)形式公開出版，并同意編入CNKI《中國(guó)知識(shí)資源總庫(kù)》，在《中國(guó)博碩士學(xué)位論文評(píng)價(jià)數(shù)據(jù)庫(kù)》中使用和在互聯(lián)網(wǎng)上傳播，同意按“章程”規(guī)定享受相關(guān)權(quán)益。論文密級(jí)：□公開 □保密（___年__月至__年__月）(保密的學(xué)位論文在解密后應(yīng)遵守此協(xié)議)作者簽名：_______ 導(dǎo)師簽名：______________年_____月_____日 _______年_____月_____日獨(dú) 創(chuàng) 聲 明本人鄭重聲明：所呈交的畢業(yè)設(shè)計(jì)(論文)，是本人在指導(dǎo)老師的指導(dǎo)下，獨(dú)立進(jìn)行研究工作所取得的成果，成果不存在知識(shí)產(chǎn)權(quán)爭(zhēng)議。盡我所知，除文中已經(jīng)注明引用的內(nèi)容外，本設(shè)計(jì)（論文）不含任何其他個(gè)人或集體已經(jīng)發(fā)表或撰寫過的作品成果。對(duì)本文的研究做出重要貢獻(xiàn)的個(gè)人和集體均已在文中以明確方式標(biāo)明。本聲明的法律后果由本人承擔(dān)。 作者簽名: 二〇一〇年九月二十日畢業(yè)設(shè)計(jì)（論文）使用授權(quán)聲明本人完全了解濱州學(xué)院關(guān)于收集、保存、使用畢業(yè)設(shè)計(jì)（論文）的規(guī)定。本人愿意按照學(xué)校要求提交學(xué)位論文的印刷本和電子版，同意學(xué)校保存學(xué)位論文的印刷本和電子版，或采用影印、數(shù)字化或其它復(fù)制手段保存設(shè)計(jì)（論文）；同意學(xué)校在不以營(yíng)利為目的的前提下，建立目錄檢索與閱覽服務(wù)系統(tǒng)，公布設(shè)計(jì)（論文）的部分或全部?jī)?nèi)容，允許他人依法合理使用。（保密論文在解密后遵守此規(guī)定）作者簽名: 二〇一〇年九月二十日致 謝時(shí)間飛逝，大學(xué)的學(xué)習(xí)生活很快就要過去，在這四年的學(xué)習(xí)生活中，收獲了很多，而這些成績(jī)的取得是和一直關(guān)心幫助我的人分不開的。首先非常感謝學(xué)校開設(shè)這個(gè)課題，為本人日后從事計(jì)算機(jī)方面的工作提供了經(jīng)驗(yàn)，奠定了基礎(chǔ)。本次畢業(yè)設(shè)計(jì)大概持續(xù)了半年，現(xiàn)在終于到結(jié)尾了。本次畢業(yè)設(shè)計(jì)是對(duì)我大學(xué)四年學(xué)習(xí)下來(lái)最好的檢驗(yàn)。經(jīng)過這次畢業(yè)設(shè)計(jì)，我的能力有了很大的提高，比如操作能力、分析問題的能力、合作精神、嚴(yán)謹(jǐn)?shù)墓ぷ髯黠L(fēng)等方方面面都有很大的進(jìn)步。這期間凝聚了很多人的心血，在此我表示由衷的感謝。沒有他們的幫助，我將無(wú)法順利完成這次設(shè)計(jì)。首先，我要特別感謝我的知道郭謙功老師對(duì)我的悉心指導(dǎo)，在我的論文書寫及設(shè)計(jì)過程中給了我大量的幫助和指導(dǎo)，為我理清了設(shè)計(jì)思路和操作方法，并對(duì)我所做的課題提出了有效的改進(jìn)方案。郭謙功老師淵博的知識(shí)、嚴(yán)謹(jǐn)?shù)淖黠L(fēng)和誨人不倦的態(tài)度給我留下了深刻的印象。從他身上，我學(xué)到了許多能受益終生的東西。再次對(duì)周巍老師表示衷心的感謝。其次，我要感謝大學(xué)四年中所有的任課老師和輔導(dǎo)員在學(xué)習(xí)期間對(duì)我的嚴(yán)格要求，感謝他們對(duì)我學(xué)習(xí)上和生活上的幫助，使我了解了許多專業(yè)知識(shí)和為人的道理，能夠在今后的生活道路上有繼續(xù)奮斗的力量。另外，我還要感謝大學(xué)四年和我一起走過的同學(xué)朋友對(duì)我的關(guān)心與支持，與他們一起學(xué)習(xí)、生活，讓我在大學(xué)期間生活的很充實(shí)，給我留下了很多難忘的回憶。最后，我要感謝我的父母對(duì)我的關(guān)系和理解，如果沒有他們?cè)谖业膶W(xué)習(xí)生涯中的無(wú)私奉獻(xiàn)和默默支持，我將無(wú)法順利完成今天的學(xué)業(yè)。四年的大學(xué)生活就快走入尾聲，我們的校園生活就要?jiǎng)澤暇涮?hào)，心中是無(wú)盡的難舍與眷戀。從這里走出，對(duì)我的人生來(lái)說(shuō)，將是踏上一個(gè)新的征程，要把所學(xué)的知識(shí)應(yīng)用到實(shí)際工作中去?；厥姿哪?，取得了些許成績(jī)，生活中有快樂也有艱辛。感謝老師四年來(lái)對(duì)我孜孜不倦的教誨，對(duì)我成長(zhǎng)的關(guān)心和愛護(hù)。學(xué)友情深，情同兄妹。四年的風(fēng)風(fēng)雨雨，我們一同走過，充滿著關(guān)愛，給我留下了值得珍藏的最美好的記憶。在我的十幾年求學(xué)歷程里，離不開父母的鼓勵(lì)和支持，是他們辛勤的勞作，無(wú)私的付出，為我創(chuàng)造良好的學(xué)習(xí)條件，我才能順利完成完成學(xué)業(yè)，感激他們一直以來(lái)對(duì)我的撫養(yǎng)與培育。最后，我要特別感謝我的導(dǎo)師趙達(dá)睿老師、和研究生助教熊偉麗老師。是他們?cè)谖耶厴I(yè)的最后關(guān)頭給了我們巨大的幫助與鼓勵(lì)，給了我很多解決問題的思路，在此表示衷心的感激。老師們認(rèn)真負(fù)責(zé)的工作態(tài)度，嚴(yán)謹(jǐn)?shù)闹螌W(xué)精神和深厚的理論水平都使我收益匪淺。他無(wú)論在理論上還是在實(shí)踐中，都給與我很大的幫助，使我得到不少的提高這對(duì)于我以后的工作和學(xué)習(xí)都有一種巨大的幫助，感謝他耐心的輔導(dǎo)。在論文的撰寫過程中老師們給予我很大的幫助，幫助解決了不少的難點(diǎn)，使得論文能夠及時(shí)完成，這里一并表示真誠(chéng)的感謝。43