【正文】 |w|≤wc Hd(ejw)= (22) 0 , wc|w|≤ π 其中 wc 也稱為截止頻率， α稱為樣本延遲 （注意，根據(jù) DTFT 性質(zhì)， ejαw意味者在正 n 方向的位移或延遲 ） 。 這個濾波器 的脈沖響應(yīng)具有無限長， 對 Hd(ejw)做 傅里葉變換獲得理想濾波器的單位脈沖響應(yīng) hd(n)。注意 hd(n)是關(guān)于 α對稱的，這一點對線性相位 FIR 濾波器來說是有用的。 為了從 hd(n)得到一個 FIR 濾波器必須在 hd(n)兩邊將它截斷。為了得到一個長度為 M的因果且線性相位 FIR 濾波器 h(n)，就必須有 hd(n) , 0≤n≤M1 h(n)= 和 α= 21M? (23) 0 , 其余 n 這種運(yùn)算稱為“加窗”。一般來說， h(n)可當(dāng)作是由 hd(n)和某一窗函數(shù) w(n)相乘而得到的。 其中 《 使用 BOHMAN窗設(shè)計的頻率采樣型 FIR濾波器 》 第 ５ 頁 共 23 頁 在 0≤n≤M1 內(nèi)關(guān)于 α的某個對稱函數(shù) w(n)= (24) 0 根據(jù)如何定義上面的 w(n)，可以得到不同的窗函數(shù)設(shè)計。 [5] 用窗口設(shè)計法設(shè)計 FIR 濾波器 的主要設(shè)計步驟為： (1)對 Hd(ejw)做 傅里葉變換獲得理想濾波器的單位脈沖響應(yīng) hd(n)。 hd(n)=)αn(π )]αn(wsin[ c (25) (2)由性能指標(biāo)確定窗函數(shù) w(n)和窗口長度 N。 (3)求得實際 濾波器的單位脈沖響應(yīng) h(n)。 h(n)=hd(n)w(n) (26) 窗口設(shè)計法的基本思想：對于給定的濾波器技術(shù)要求，選擇濾波器長度 M 和具有最窄主瓣寬度和盡可能最小的旁瓣衰減的某個窗函數(shù) w(n)。 窗函數(shù)頻譜中的主瓣寬度 較 窄，獲得的過渡帶較陡；旁瓣衰減 越 大，阻帶的衰減 越快。 表 21 常見窗函數(shù)性能小結(jié) 過渡帶寬 Δw 窗函數(shù)名 最小阻帶衰減 近似值 準(zhǔn)確值 矩形 Mπ4 21dB Bartlett Mπ8 25dB Hanning Mπ8 44dB Hamming Mπ8 53dB Blackman M12? Mπ11 74dB BOHMAN 窗 本課程設(shè)采用 BOHMAN 窗，通過查閱資料和實際調(diào)試， BOHMAN 窗的過渡帶寬定為 Mπ10 ，最小阻帶衰減為 。 當(dāng)運(yùn)用計算機(jī)實現(xiàn)工程 測試信號處理 時，不可能對無限長的信號進(jìn)行測量和運(yùn)算，而 《 使用 BOHMAN窗設(shè)計的頻率采樣型 FIR濾波器 》 第 ６ 頁 共 23 頁 是取其有限的時間片段進(jìn)行分析。做法是從信號中截取一個時間片段，然后用觀察的信號時間片段進(jìn)行周期延拓處 理，得到虛擬的無限長的信號，然后就可以對信號進(jìn)行傅里葉變換、相關(guān)分析等數(shù)學(xué)處理。無線長的信號被截斷以后，其頻譜發(fā)生了畸變，原來集中在 f(0)處的能量被分散到兩個較寬的頻帶中去了（這種現(xiàn)象稱之為頻譜能量泄漏） 。 [6] Bohman 窗對應(yīng)的 Bohmanwin 函數(shù)定義源程序如下： function w = bohmanwin(n) error(nargchk(1,1,nargin))。 [n,w,trivialwin] = check_order(n)。 if trivialwin, return, end。 q = abs(linspace(1,1,n))。 w = [ 0。 ((1q(2:end1)).*cos(pi*q(2:end1))+(1/pi)*sin(pi*q(2:end1)))39。 0]。 當(dāng) n=64 時， Bohman 窗 如圖 21 所示。 圖 21 Bohwin 窗 3 設(shè)計步驟 《 使用 BOHMAN窗設(shè)計的頻率采樣型 FIR濾波器 》 第 ７ 頁 共 23 頁 設(shè)計流程圖 語音信號濾波去噪 —— 使用 BOHMAN 窗設(shè)計的頻率采樣型 FIR 濾波器 的設(shè)計流程圖如圖 31 所示。 開始 錄制一段格式為 .wav 的語音信號 利用函數(shù) wavread 對語音信號進(jìn)行采樣 對語音信號進(jìn)行加噪 繪制 加噪前后 語音 信號的時域和 頻域波形 使用 BOHMAN 窗設(shè)計 FIR 濾波器 繪制濾波器的 頻率響應(yīng) 用濾波器對信號進(jìn)行濾波 比較 原始語音信號和 濾波前后語音信號的波形及頻譜 回放語音信號 計算頻率采樣型 FIR 濾波器系數(shù) 繪制頻率采樣型 FIR 濾波器的 結(jié)構(gòu)圖并填好系數(shù) 結(jié)束 圖 31 設(shè)計流程圖 《 使用 BOHMAN窗設(shè)計的頻率采樣型 FIR濾波器 》 第 ８ 頁 共 23 頁 錄制語音信號 通過點擊 開始 /程序 /附件 /娛樂 /錄音機(jī)，打開 Windows 下的錄音機(jī) 工具 ，如圖 所示，然后 錄制 語音信號 “大家好，我是 ***” ，時間在 23 s 左右 。 點擊文件 /另存為，將文件名命名為 ， 并保存。 然后在 Matlab 軟件平臺下，利用函數(shù) wavread 對語音信號進(jìn)行采樣，記住采樣頻率和采樣點數(shù)。采集完成后在信號中加入一個 頻率為 9000Hz 的單頻噪聲 。 圖 32 錄音機(jī) 為了使后面濾波器的設(shè)計更精確，可對語音信號進(jìn)行一些調(diào)整。 首先通過錄音機(jī)打開語音信號，點擊文件 /屬性，如圖 33 所示 。再點擊 按鈕進(jìn)入聲音選定窗口，并將屬性改為 ，如圖 34 所示，點擊確定并保存調(diào)整后的語音信號。 圖 33 聲音屬性 《 使用 BOHMAN窗設(shè)計的頻率采樣型 FIR濾波器 》 第 ９ 頁 共 23 頁 圖 34 聲音選定 對應(yīng)程序如下： [x,fs,bits]=wavread(39。F:\39。)。 % 輸入?yún)?shù)為文件的全路徑和文件名，輸出的第一個參數(shù)是每個樣本的值， fs 是生成該波形文件時的采樣率， bits 是波形文件每樣本的編碼位數(shù)。 sound(x,fs,bits)。 % 按指定的采樣率和每樣本編碼位數(shù)回放