【導(dǎo)讀】處理的理論知識(shí)，對(duì)含噪聲的信號(hào)進(jìn)行去噪處理，并對(duì)信號(hào)的時(shí)域和頻域進(jìn)行分析。通過理論推導(dǎo)得出相應(yīng)結(jié)論，再利用MATLAB作為編程工具進(jìn)行計(jì)算機(jī)實(shí)現(xiàn)。應(yīng)用數(shù)字信號(hào)處理的相應(yīng)理論，以MATLAB為實(shí)驗(yàn)工具，分別設(shè)計(jì)了IIR濾。在設(shè)計(jì)實(shí)現(xiàn)的過程中，使用窗函數(shù)。波的方法對(duì)噪聲信號(hào)進(jìn)行處理，并分析了不同去噪方法所適用的信號(hào)類型。了詳細(xì)的闡述，使人們今后在進(jìn)行信號(hào)去噪時(shí)可以更有針對(duì)性的選擇去噪方法。

　　

【正文】 不存在，去噪效果良好。 （ 2）對(duì)實(shí)測(cè)信號(hào) X2 的去噪仿真， 信號(hào) X2 波形 如圖所示： 0 2020 4000 6000 8000 10000 12020 0 . 7 0 . 6 5 0 . 6 0 . 5 5 0 . 5 0 . 4 5 0 . 4 0 . 3 5 圖 實(shí)測(cè)信號(hào)二波形圖 32 經(jīng)過觀察 第二組信號(hào) 的圖像，我們可以發(fā)現(xiàn)，信號(hào) X2 的噪聲比信號(hào) X1的噪聲更加稠密，對(duì)此，我打算 首先 使用最為傳統(tǒng)的 FIR 低通濾波器對(duì)其進(jìn)行去噪。 FIR濾波器參數(shù)如下： Ft=8000， Fp=5， Fs=30； 程序 代碼如下： clf。 Ft=8000。 Fp=5。 Fs=30。 wp=2*Fp/Ft。 ws=2*Fs/Ft。 rp=1。 rs=50。 p=110.^(rp/20)。 %通帶阻帶波紋 s=10.^(rs/20)。 fpts=[wp ws]。 mag=[1 0]。 dev=[p s]。 [n21,wn21,beta,ftype]=kaiserord(fpts,mag,dev)。%kaiserord求階數(shù)截止頻率 b21=fir1(n21,wn21,Kaiser(n21+1,beta))。 %由 fir1設(shè)計(jì)濾波器 [h,w]=freqz(b21,1)。 %得到頻率響應(yīng) plot(w/pi,abs(h))。 title(39。FIR低通濾波器 39。)。 load x2。 S=fft(x2)。 %傅里葉變換 33 z21=fftfilt(b21,x2)。 m21=fft(z21)。 %求濾波后的信號(hào) figure(1)。 subplot(2,1,1)。 plot(x2)。 title(39。濾波前信號(hào)的波形 39。)。 figure(1)。 subplot(2,1,2)。 plot(z21)。 title(39。濾波后的信號(hào)波形 39。)。 0 2020 4000 6000 8000 10000 12020 0 . 7 0 . 6 0 . 5 0 . 4 0 . 3濾波前信號(hào)的波形0 2020 4000 6000 8000 10000 12020 0 . 6 0 . 4 0 . 200 . 2濾波后的信號(hào)波形 圖 FIR低通信號(hào)去噪波形圖 觀察信號(hào)去噪后的圖像，可以看到，去噪后的圖像較為平滑，所以說 FIR 低通濾波器對(duì)信號(hào)的去噪效果很好，但是，對(duì)于此信號(hào)的去噪不止可以用 FIR 濾波器的方式，也可以用均值濾波的方式。 34 均值濾波代碼如下： load x2。 figure(1)。 subplot(2,1,1)。 title(39。原信號(hào) 39。)。 plot(x2) z=ones(68,1)。 x68=conv(x2,z)/68。 figure(1)。 subplot(2,1,2)。 plot(x68) title(39。濾波后信號(hào) 39。)。 去噪后的信號(hào)波形如圖 所示。 0 2020 4000 6000 8000 10000 12020 0 . 7 0 . 6 0 . 5 0 . 4 0 . 30 2020 4000 6000 8000 10000 12020 0 . 8 0 . 6 0 . 4 0 . 20濾波后信號(hào) 圖 中值濾波后信號(hào) （ 3）實(shí)測(cè)信號(hào) X3 去噪仿真， 信號(hào) X3 的波形 如 下 圖所示： 35 0 200 400 600 800 1000 12002 . 4 22 . 4 32 . 4 42 . 4 52 . 4 62 . 4 72 . 4 82 . 4 92 . 52 . 5 1 圖 實(shí)測(cè)信號(hào)三波形圖 對(duì) 第三組信號(hào)我 利用貝塞爾濾波器 進(jìn)行設(shè)計(jì) ， 電子學(xué)和信號(hào)處理中，貝賽爾(Bessel)濾波器是具有最大平坦的群延遲（線性相位響應(yīng)）的線性過濾器。貝賽爾濾波器常用在音頻天橋系統(tǒng)中。模擬貝賽爾濾波器描繪為幾乎橫跨整個(gè)通頻帶的恒定的群延遲，因而在通頻帶上保持了被過濾的信號(hào)波形。貝塞 爾 (Bessel)濾波器具有最平坦的幅度和相位響應(yīng)。帶通（通常為用戶關(guān)注區(qū)域）的相位響應(yīng)近乎呈線性。Bessel濾波器可用于減少所有 IIR濾波器固有的非線性相位失真。其參數(shù)如下： N=5，wn=1000， 程序 代碼如下： clear all。 N=5。 wn=1000。 [b,a]=besself(N,wn)。 [h,w]=freqs(b,a,512)。 magH=abs(h)。 36 phaH=unwrap(angle(h))。 figure。 subplot(2,1,1)。 plot(w/(2*pi),2*log10(magH))。 xlabel(39。頻率 /Hz39。)。 ylabel(39。振幅 /dB39。)。 grid on。 subplot(2,1,2)。 plot(w/2/pi,angle(h)*180/pi)。 grid on。 xlabel(39。頻率 /Hz39。)。 ylabel(39。相位 /dB39。)。 figure。 dt=1/10000。 f1=100。 f2=1000。 t=0:dt:。 load x3。 H=[tf(b,a)]。 [y,t1]=lsim(H,x3,t)。 subplot(2,1,1)。 plot(t,x3)。 title(39。輸入信號(hào) 39。)。 xlabel(39。時(shí)間 /s39。)。 37 subplot(2,1,2)。 plot(t1,y)。 title(39。輸出信號(hào) 39。)。 xlabel(39。時(shí)間 /s39。)。 信號(hào)去噪后波形圖如圖 所示： 0 0 . 0 1 0 . 0 2 0 . 0 3 0 . 0 4 0 . 0 5 0 . 0 6 0 . 0 7 0 . 0 8 0 . 0 9 0 . 12 . 4 22 . 4 42 . 4 62 . 4 82 . 5輸入信號(hào)時(shí)間 /s0 0 . 0 1 0 . 0 2 0 . 0 3 0 . 0 4 0 . 0 5 0 . 0 6 0 . 0 7 0 . 0 8 0 . 0 9 0 . 10123輸出信號(hào)時(shí)間 /s 圖 貝塞爾去噪后信號(hào)波形 根據(jù) 圖 ，我們可以了解到， 經(jīng)過貝塞爾濾波器 去噪 后的實(shí)測(cè)信號(hào)十分平滑，幾乎沒有噪聲的存在 ，信號(hào)去噪 的效果十分好 。 38 第 5 章 總 結(jié) 本文主要 內(nèi)容為信號(hào)的去噪仿真，在文中列舉了幾種不同的去噪方法，也對(duì)不同 方法的使用與濾波效果進(jìn)行了對(duì) 比 。 同時(shí)，也對(duì)幾種 典型 信號(hào)進(jìn)行了頻譜分析，對(duì)頻譜圖的含義進(jìn)行了概述，有助于大家對(duì)頻譜圖的重視和理解。 通 過對(duì) 給定的三組實(shí)測(cè) 信號(hào)進(jìn)行去噪，我們可以對(duì)去噪的方法進(jìn)行簡(jiǎn)單的分析，對(duì)于多種不同含噪聲信號(hào)來說， IIR 與 FIR 濾波器是最為有效也是最為普遍的去噪方 法 ，雖然，個(gè)別時(shí)候這兩種去噪方 法 與其他去噪方 法 相比效果并不十分完美，但其適用廣泛，簡(jiǎn)單，是很好的去噪方 法 。 中值濾波在如今的數(shù)字信號(hào)處理和數(shù)字圖像處理中都是應(yīng)用比較普遍的方法。設(shè)有一個(gè)一維序列 (f1,f2,f3,?? ,fn)。取窗口長(zhǎng)度為 M(M 為奇數(shù) )，對(duì)此一維 信號(hào)進(jìn)行中值濾波，就是從輸入序列中相繼抽出 M 的個(gè)數(shù)(fiv,? ,fi1,fi,fi+1,? ,fi+v),其中 fi 為窗口的中心點(diǎn)值， v=(M1)/2。再將這 M個(gè)點(diǎn)的值按其數(shù)值大小排序，取其序號(hào)為正中間那個(gè)數(shù)作為濾波輸出。中值濾波的窗口形狀和尺寸設(shè)計(jì)對(duì)濾波效果影響很大。不同的圖像要求和不同的應(yīng)用要求，往往采用不同的窗口形狀和尺寸。 均值濾波就是在一定的時(shí)間內(nèi)對(duì)輸入信號(hào)進(jìn)行連續(xù)的采樣然后求出幾次采樣值的平均值。設(shè) x0．．． xi 是采樣值， yk 是對(duì)采樣進(jìn)行平均值濾波的輸出值，其關(guān)系如下： 在工程應(yīng)用中，均值濾 波對(duì)于周期性的干擾有很好的抑制作用，另外，均值濾波對(duì)于熱噪聲產(chǎn)生的干擾也有很好的抑制作用。 n 的取值取決 n 系統(tǒng)的參數(shù)，如果系統(tǒng)對(duì)輸入的靈敏度要求很高，則 n的取值應(yīng)比較小，如果系統(tǒng)的輸入變化速度不快，系統(tǒng)對(duì)輸入的靈敏性要求也不是很高，為了抑制干擾，提高輸入的精度，可以將 n 取得比較大。 39 通過這次設(shè)計(jì)使我懂得了理論與實(shí)際相結(jié)合是很重要的，只有理論知識(shí)是遠(yuǎn)遠(yuǎn)不夠的，只有把所學(xué)的理論知識(shí)與實(shí)踐相結(jié)合起來，從理論中得出結(jié)論，才能真正為社會(huì)服務(wù)，從而提高自己 的實(shí)際動(dòng)手能力和獨(dú)立思考的能力。在設(shè)計(jì)的過程中遇到問題，可以說 是 困難重重，這畢竟第一次做的，難免會(huì)遇到過各種各樣的問題，同時(shí)在設(shè)計(jì)的過程中發(fā)現(xiàn)了自己的不足之處，對(duì)以前所學(xué)過的知識(shí)理解得不夠深刻，掌握得不夠牢固。 但最終我還是克服了各種困難將設(shè)計(jì)完美的完成了。 40 參考文獻(xiàn) [1] 高西全 ,丁玉美 . 數(shù)字信號(hào)處理 [M].西安 :西安電子科技大學(xué)出版社 ,2020 [2] 叢玉良 ,王宏志 . 數(shù)字信號(hào)處理原理 [M].北京 :電子工業(yè)出版社 ,2020 [3] 黃植功 . 基于 MATLAB 的連續(xù)信號(hào)頻譜分析 [D]. 桂林 :廣西師范大學(xué)電子工程學(xué)院 ,2020 [4] 井敏英 ,郭佳林 ,李大偉等 . 連續(xù)時(shí)間信號(hào)頻譜分析研究及 MATLAB 實(shí)現(xiàn) [J]. 大眾科技 , 2020(4):4243 [5] 張德豐 . 詳解 MATLAB數(shù)字信號(hào)處理 [M]. 北京 : 電子工業(yè)出版社 ,2020 [6] 王彬 . MATLAB數(shù)字信號(hào)處理 [M]. 北京 : 機(jī)械工業(yè)出版社 ,2020 [7] 呂晶晶 ,趙晶亮 ,劉國(guó)鵬 .基于 MATLAB的振動(dòng)信號(hào)去噪研究 [J].電子測(cè)試 , 2020(7): 110113 [8] 周輝 ,董正紅 . 數(shù)字信號(hào)處理基礎(chǔ)及 MATLAB實(shí)現(xiàn) [M]. 北京 : 北京希 望電子出版社 , 2020 [9] 趙淑敏 .基于 MATLAB實(shí)現(xiàn)對(duì)語(yǔ)音信號(hào)頻譜分析 [J].技術(shù)前沿 ,2020(1):1215 [10]劉小群 ,周云波 .基于 MATLAB的 DFT及 FFT頻譜分析 [J].山西電子技術(shù) , 2020 (4): 4750 [11]李輝 . 數(shù)字信號(hào)處理及 MATLAB實(shí)現(xiàn) [M]. 北京 :機(jī)械工業(yè)出版社 , 2020 [12]劉志松 .基于小波分析的信號(hào)去噪方法 [J]. 浙江海洋學(xué)院學(xué)報(bào) (自然科學(xué)版 ), 2020,30(2):150154 [13]林川 . MATLAB與數(shù)字信號(hào)處理實(shí)驗(yàn) [M]. 武漢 :武漢大學(xué)出版社 ,2020 [14] Min Q,Shizhen H. Research on Wavelet Threshold DENOISING Method Based on MATLAB[J]電子器件 ,2020, 35(1):103106 [15]門愛東 .數(shù)字信號(hào)處理 [M].北京 :科學(xué)出版社 ,2020 [16]英格爾 .數(shù)字信號(hào)處理 —— 應(yīng)用 MATLAB[M].北京 :科學(xué)出版社 ,2020 [17]姚天任 .數(shù)字信號(hào)處理 [M].北京 :清華大學(xué)出版社 ,2020 41 致 謝 經(jīng)過這幾個(gè)多月的努力 ，這篇基于 MATLAB 頻譜分析與信號(hào)去噪的論文終于完成了， 值此畢業(yè)論文完成之際，我謹(jǐn)向所有關(guān)心、愛護(hù)、幫助我的人們表示最誠(chéng)摯的感謝與最美好的祝愿。 本論文是在我最敬愛的導(dǎo)師崔暢老師的指導(dǎo)下完成的 ，她嚴(yán)謹(jǐn)細(xì)致、一絲不茍的作風(fēng)一直是我工作、學(xué)習(xí)中的榜樣；導(dǎo)師循循善誘的教導(dǎo)和不拘一格的思路給予我無盡的啟迪。 這 幾 個(gè)月中 ，導(dǎo)師嚴(yán)謹(jǐn)?shù)闹螌W(xué)態(tài)度，精益求精的工作作風(fēng)，誨人不倦的高尚師德，樸實(shí)無華、平易近人的人格魅力對(duì)我影響深遠(yuǎn) 。而且 本論文從選題到完成，幾易其稿，每一步都是在導(dǎo)師的指導(dǎo)下完成的，傾注了導(dǎo)師大量的心血，在此我向 我的導(dǎo)師 崔暢老師 表示深切的謝意與祝福！ 而且我也要感謝我的同學(xué)，在畢業(yè)設(shè)計(jì)期間，我的同學(xué)和朋友給了我很大幫助，沒有他們或許我的畢業(yè)設(shè)計(jì)不會(huì)這么順利。 同時(shí)我也深刻的認(rèn)識(shí)到，在對(duì)待一個(gè)新事物時(shí)，一定要從整體考慮，完成一步之后再作下一步，這樣才能更加有效。這份 畢業(yè)設(shè)計(jì)，也許是我大學(xué)生涯交上的最后一個(gè)作業(yè)了，想借此機(jī)會(huì)感謝四年以來給我?guī)椭乃欣蠋煛⑼瑢W(xué)們，你們的友誼是我人生的財(cái)富，是我生命中不可或缺的一部分， 感謝你們對(duì)我巨大 的 幫助和支持，讓我渡過了難關(guān)，謝謝你們 ！