【正文】 ylabel(39。 subplot(3,2,6)。)。ylabel(39。 subplot(3,2,4)。)。ylabel(39。 subplot(3,2,2)。)。 % 加上網(wǎng)格 y_fil=fftfilt(h_bs,y)。)。ylabel(39。 subplot(2,2,3)。)。ylabel(39。 % 調(diào)用自編函數(shù)計(jì)算濾波器的頻率特性 figure(2) subplot(2,2,1)。 wsu=fsu/fs*2*pi。As=36。Y39。plot(f,Y)。 xlabel(39。X39。plot(f,X)。 xlabel(39。 Y=Y(1:length(Y)/2)。 % 按指定的采樣率和每樣本編碼位數(shù)回放 N=length(x)。 整一個(gè)流程下來，使自己更好地掌握課本的知識，理論結(jié)合實(shí)際，學(xué)以致用。 5. 在做課程設(shè)計(jì)的文檔時(shí)，發(fā)現(xiàn)流程圖和公式不會繪制。用漢寧窗設(shè)計(jì) FIR 濾波器濾掉了語音中的噪聲，但是也把原始語音的很小的一部分也濾掉了，所以回放語音的時(shí)候聽起來比以前的更加平滑，說明這段程序設(shè)計(jì)是成功的。 15 圖 35 FIR 濾波器的頻率響應(yīng) 信號濾波處理 用自己設(shè)計(jì)的各濾波器分別對采集的信號進(jìn)行濾波，在 Matlab 中， FIR 濾波器利用函數(shù) fftfilt 對信號進(jìn)行濾波， IIR 濾波器利用函數(shù) filter 對信號進(jìn)行濾波。 M=ceil(*pi/dw)+1。 % 帶阻濾波器設(shè)計(jì)指標(biāo) fcd=(fpd+fsd)/2。其中 b 為系統(tǒng)函數(shù)的分子系數(shù)， a為系統(tǒng)函數(shù)分母系數(shù)。 % 對原始信號和加噪信號進(jìn)行 fft 變換 X=X(1:length(X)/2)。 % 按指定的采樣率和每樣本編碼位數(shù)回放 N=length(x)。采樣頻率越高，即采樣的間隔時(shí)間越短，則在單位時(shí)間內(nèi)計(jì)算機(jī)得到的聲音樣本數(shù)據(jù)就越多，對聲音波形的表示也越精確。 漢寧窗 漢寧窗 （ Hanning Window） 又稱升余弦窗，漢寧窗可以看作是 3 個(gè)矩形時(shí)間窗的頻譜之和，或者說是 3 個(gè) ??tsin 型函數(shù)之和，而括號中的兩項(xiàng)相對于第一個(gè)譜窗向左、右各移動了 π/T ，從而使旁瓣互相抵消，消去高頻干擾和漏能。因此 ??nhd 是無限長非因果序列，不能直接作為 FIRDF 的單位脈沖響應(yīng)。與前面的表示形式，即)()()( ??? ? jgj eHeH ? 相比， )(?gH 與 ? 不同。 MATLAB 可以進(jìn)行 矩陣 運(yùn)算、繪制函數(shù)和數(shù)據(jù)、實(shí)現(xiàn)算法、創(chuàng)建用戶界面、連接其他編程語言的程序等，主要應(yīng)用于工程計(jì)算、控制設(shè)計(jì)、信號處理與通訊、圖像處理 、 信號檢測 、金融建模設(shè)計(jì)與分析等領(lǐng)域 。此課程設(shè)計(jì)比較簡單， 主要是將書本中的知識運(yùn)用到現(xiàn)實(shí)中，并且根據(jù)自己對設(shè)計(jì)題目的理解，運(yùn)用軟件編寫出程序?qū)崿F(xiàn)這一設(shè)計(jì)，也是我們對數(shù)字信號處理的原理進(jìn)行驗(yàn)證的一個(gè)過程。程序經(jīng)過多次的調(diào)試運(yùn)行，基本實(shí)現(xiàn)了設(shè)計(jì)的目標(biāo)，并且經(jīng)過適當(dāng)?shù)男薷暮?，從波形可以看出基本可以濾完全除噪聲而且達(dá)到語音不失真的效果。 關(guān)鍵詞 DSP 應(yīng)用實(shí)例； 濾波去噪； FIR濾波器；漢寧窗； MATLAB 使用 3 目 錄 引言 ．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．． 4 ．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．． ．．． 4. 課程設(shè)計(jì)的要求 ．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．． 4 設(shè)計(jì)平臺 ．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．． 5 設(shè)計(jì)原理 ．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．． 6 FIR濾波器 ．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．． 6 窗口設(shè)計(jì)法 ．．．．．．．．．．．．．．．．．．．．．．． ．．．．．．．．．．．．．．．．．．．． 7 漢寧窗 ．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．． 8 設(shè)計(jì)步驟 ．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．． 9 設(shè)計(jì)流程圖 ．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．． 9 錄制語音信號 ．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．． 10 濾波器設(shè)計(jì) ．．．．．． ．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．． 12 信號濾波處理 ．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．． 14 結(jié)果分析 ．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．． 15 出現(xiàn)的問題及解決方法 ．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．． 16 項(xiàng)目完成 總 結(jié) ．．．．．． ． ．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．． 17 4 1. 引 言 本課程設(shè)計(jì)主要是對一段語音信號，加入噪聲后，用某種函數(shù)法設(shè)計(jì)出的FIR 濾波器對加入噪聲后的語音信號進(jìn)行濾波去噪處理，并且分析對比前后時(shí)域和頻域波形的程序設(shè)計(jì)。對此，也可以加深我們對所學(xué)知識的理解，培養(yǎng)我們的動手能力。 6 圖 11 MATLAB 的設(shè)計(jì)平臺 MATLAB 的基本數(shù)據(jù)單位是矩陣，它的指令表達(dá)式與數(shù)學(xué)、工程中常用的形式十分 相似，故用 MATLAB 來解算問題要比用 C， FORTRAN 等語言完成相同的事情簡捷得多，并且 MathWork也吸收了像 Maple 等軟件的優(yōu)點(diǎn) ,使 MATLAB 成為一個(gè)強(qiáng)大的數(shù)學(xué)軟件。 )(?? 與 )(?? 不同。窗口設(shè)計(jì)法就是截取 ??nhd 為有限長 8 的一段因果序列，并用合適的窗口函數(shù)進(jìn)行加權(quán)作為 FIRDF 的單位脈沖響應(yīng)??nh 。可以看，漢寧窗主瓣加寬并降低，旁瓣則顯著減小，從減小泄漏觀點(diǎn)出發(fā)，漢寧窗優(yōu)于矩形窗．但漢寧窗主瓣加寬，相當(dāng)于分析帶寬加寬，頻率分辨力下