FIR數(shù)字濾波器的原理與設計.ppt

第七章FIR數(shù)字濾波器的原理與設計 宜春學院理工學院 內容提要 7 1線性相移FIR數(shù)字濾波器的特性7 2窗口法7 3頻率取樣法7 4FIR數(shù)字濾波器的優(yōu)化設計7 5IIR數(shù)字濾波器與FIR數(shù)字濾波器的比較習題及作業(yè) 學習目標 掌握線性相位FIR數(shù)字濾波器的特點掌握窗函數(shù)設計法理解頻率抽樣設計法了解設計FIR濾波器的最優(yōu)化方法理解IIR與FIR數(shù)字濾波器的比較 7 1FIR數(shù)字濾波器的差分方程 沖激響應 系統(tǒng)函數(shù)及其零極點FIR數(shù)字濾波器是非遞歸的線性時不變因果系統(tǒng) 其差分方程為 系統(tǒng)的沖激響應為 可見這個系統(tǒng)的沖激響應是有限長度的 即有限沖激響應 FIR 濾波器 上式兩邊進行Z變換后 可得FIR濾波器的系統(tǒng)函數(shù) 可見 FIR濾波器的系統(tǒng)函數(shù)的極點都位于z 0處 為N 1階極點 與系數(shù)h n 無關 因此FIR濾波器總是穩(wěn)定的 而N 1個零點由沖激響應h n 決定 可以位于有限z平面的任何位置 兩種濾波器的比較一 IIRDF的特點1 DF的設計依托AF的設計 有圖表可查 方便簡單 2 相位的非線性H Z 的頻響 其中 是幅度函數(shù) 是相位函數(shù) 通常 與不是呈線性的 這是IIRfilter 無限長響應濾波器 的一大缺點 因此限制了它的應用 如圖象處理 數(shù)據(jù)傳輸都要求信道具有線性相位特性 3 用全通網(wǎng)絡進行相位校正 可以得線性特性 二 FIRDF的特點1 單位抽樣響應h n 是有限長的 因此FIRDF一定是穩(wěn)定的 2 經延時 h n 總可變成因果序列 所以FIRDF總可以由因果系統(tǒng)實現(xiàn) 3 h n 為有限長 可以用FFT實現(xiàn)FIRDF 4 FIR的系統(tǒng)函數(shù)是Z 1的多項式 故IIR的方法不適用 5 FIR的相位特性可以是線性的 因此 它有更廣泛的應用 非線性的FIR一般不作研究 FIR與IIR數(shù)字濾波器比較 優(yōu)點 1 很容易獲得嚴格的線性相位 避免被處理的信號產生相位失真 這一特點在寬頻帶信號處理 陣列信號處理 數(shù)據(jù)傳輸?shù)认到y(tǒng)中非常重要 2 可得到多帶幅頻特性 3 極點全部在原點 永遠穩(wěn)定 無穩(wěn)定性問題 4 任何一個非因果的有限長序列 總可以通過一定的延時 轉變?yōu)橐蚬蛄?所以因果性總是滿足 5 無反饋運算 運算誤差小 缺點 1 因為無極點 要獲得好的過渡帶特性 需以較高的階數(shù)為代價 2 無法利用模擬濾波器的設計結果 一般無解析設計公式 要借助計算機輔助設計程序完成 7 2 1線性相移FIR數(shù)字濾波器條件 所謂線性相移濾波器 也就是指其相移特性或頻率響應的幅角是頻率的線性函數(shù) FIR數(shù)字濾波器頻率響應為 1 恒時延濾波相延時群延時所謂恒延時濾波就是要求相延時與群延時都是不隨頻率變化的常量 其中有 所謂時延是指信號通過傳輸通道所需要的傳輸時間 它是濾波器平均延遲的一個度量 它是濾波器某一頻率延遲的一個度量 7 2線性相移FIR數(shù)字濾波器 2 要求恒相延時與恒群延時同時成立 如圖7 1 的圖像是一條經過原點的直線 式中H 是正或負的實函數(shù) 等式中間和等式右邊的實部與虛部應當各自相等 同樣實部與虛部的比值應當相等 由上式交叉相乘后利用三角函數(shù)恒等公式得 滿足上式的條件是 上述條件下 就有 即 為一常數(shù) 恒相延時與恒群延時同時成立 如上所述 沖激相應h n 關于中心點偶對稱 由圖7 2可見無論N是偶數(shù)還是奇數(shù) 對稱中心都位于 N 1 2 只是當N為偶數(shù)時 N 1 2不是整數(shù) 圖7 2h n 為偶對稱的情形 3只要求恒群延時成立 相移特性為一條不經過原點的直線 如下圖可見沖激響應關于中心點奇對稱 無論N為奇數(shù)還是偶數(shù) 對稱中心都位于 N 1 2 當N為奇數(shù)時有 圖7 4h n 為奇對稱的情形 總之 線性相移FIR濾波器的必要條件是其沖激響應為偶對稱或奇對稱 7 1 2線性相移FIR濾波器的網(wǎng)絡結構 1 偶對稱的情形 偶對稱時 a N為偶數(shù)時 利用對稱性可作如下化簡 b N為奇數(shù)時 利用對稱性可作如下化簡 可見 以其偶對稱性作這樣的簡化可以使FIR濾波器比一般的直接型結構的乘法器減少一半 2 奇對稱的情形 奇對稱時 a 當N為偶數(shù)時 b 當N為奇數(shù)時 可見 以其奇對稱性作這樣的簡化可以使FIR濾波器比一般的直接型結構的乘法器減少近一半 7 2 3線性相移FIR濾波器的頻率響應 1 偶對稱 N為奇數(shù) 則其頻率響應為 則 圖7 5偶對稱 N為奇數(shù) 該類濾波器適合于設計任何關于為偶對稱特性頻率的濾波器 特點 對皆為偶對稱 所以幅度函數(shù)對也是偶對稱 2 偶對稱 N為偶數(shù) 其頻率響應為 圖7 6偶對稱 N為偶數(shù) 3 奇對稱 N為奇數(shù) 推導方法與前面類似 可得 其幅頻特性為 圖7 7奇對稱 N為奇數(shù) 4 奇對稱 N為偶數(shù) 推導方法與前面類似 可得 其幅頻特性為 圖7 8奇對稱 N為偶數(shù) 這種情況不適合做在處為偶對稱的濾波器 如低通濾波器 特點 當時 相當于在z 1處有一個零點 并且由于對呈奇對稱 對呈偶對稱 因而也對呈奇對稱 對呈偶對稱 以上四種情況可以用統(tǒng)一的形式 即 其中 的實函數(shù) 相移由 決定 而 的線性函數(shù) 當h n 為偶對稱時 當h n 為奇對稱時 圖7 94種類型的線性相位濾波器的相位響應 時域幅度響應和頻域幅度響應的示意圖 總結 第一種情況 偶 奇 四種濾波器都可設計 第二種情況 偶 偶 可設計低 帶通濾波器 不能設計高通和帶阻 第三種情況 奇 奇 只能設計帶通濾波器 其它濾波器都不能設計 第四種情況 奇 偶 可設計高通 帶通濾波器 不能設計低通和帶阻四種FIR數(shù)字濾波器的相位特性只取決于h n 的對稱性 而與h n 的值無關 幅度特性取決于h n 設計FIR數(shù)字濾波器時 在保證h n 對稱的條件下 只要完成幅度特性的逼近即可 例1N 5 h 0 h 1 h 3 h 4 1 2 h 2 2 求幅度函數(shù)H 解 a 0 h 2 2 a 1 2h 3 1 a 2 2h 4 1 H 2 cos cos2 2 cos cos2 分析 為奇數(shù) 并且h n 滿足偶對稱關系 7 2 4線性相移FIR數(shù)字濾波器的零 極點分布 線性相移FIR濾波器有 即 當為實數(shù)時 為實系數(shù)的多項式 此時應是共軛成對的 則也是零點 所以線性相位濾波器的零點必須是互為倒數(shù)的共軛對 這種共軛對共有四種 1 既不在單位圓上 也不在實軸上 有四個互為倒數(shù)的兩組共軛對 如圖zi z i 1 zi 1 z i 2 在單位圓上 但不在實軸上 因倒數(shù)就是自己的共軛 所以有一對共軛零點 zi z i 圖7 10 a 零點分布 3 不在單位圓上 但在實軸上 是實數(shù) 共軛就是自己 所以有一對互為倒數(shù)的零點 zi 1 zi 線性相位濾波器是FIR濾波器中最重要的一種 應用最廣 實際使用時應根據(jù)需用選擇其合適類型 并在設計時遵循其約束條件 圖7 10 b 零點分布 7 3FIR數(shù)字濾波器的設計 窗口法 7 3 0引言設計思路 1 先給定所要求設計的理想濾波器的頻率響應Hd ejw 2 設計一個可實現(xiàn)的FIR濾波器頻率響應H ejw 3 由于設計是在時域中進行 使所設計濾波器的h n 去逼近理想單位取樣響應hd n 如果希望得到的濾波器的理想頻率響應為 那么FIR濾波器的設計就在于尋找一個傳遞函數(shù)去逼近 逼近方法有三種 窗口設計法 時域逼近 頻率采樣法 頻域逼近 最優(yōu)化設計 等波紋逼近 時間窗口設計法是從單位脈沖響應序列著手 使h n 逼近理想的單位脈沖響應序列hd n 我們知道hd n 可以從理想頻響通過付氏反變換獲得 但一般來說 理想頻響是矩形頻率特性 所以 這樣得到的理想單位脈沖響應hd n 往往都是無限長序列 而且是非因果的 但FIR的h n 是有限長的 問題是怎樣用一個有限長的序列去近似無限長的hd n 最簡單的辦法是直接截取一段hd n 代替h n 這種截取等效于在hd n 上施加了一個長度為N的矩形窗 h n 是通過一個 窗口 所看到的一段 因此 h n 也可表達為h n 和一個 窗函數(shù) 的乘積 即h n w n hd n 這一方法通常稱為窗口設計法 設計步驟 1 由定義 3 卷積 插值 7 3 1窗口法的基本思想 1 設計思想在時域 設計逼近理想 下面以理想低通濾波器為例說明其設計過程 為一 以為對稱中心的 偶對稱的 無限長的 非因果序列 圖7 11理想低通濾波器的單位脈沖響應及矩形窗截取 要得到有限長的 最簡單的方法是用一長為的矩形窗截斷 按照線性相位濾波器的要求 必須偶對稱 如上圖 對稱中心必須等于濾波器的延時常數(shù) FIR濾波器的沖擊響應h n 的頻響H ejw 一定與理想的頻響Hd ejw 存在差異 圖7 12理想低通濾波器的頻率響應 7 3 2理論分析 可見 WR ejw 是w的偶函數(shù) 7 13矩形窗的頻譜 7 14矩形窗的卷積過程 正肩峰 負肩峰 圖7 15由 wc到wc區(qū)間曲線WR ej w q 下面積隨w取值變化演示 圖7 16加矩形窗后的頻響與理想頻響的比較 對加矩形窗處理后其頻率響應將產生以下幾點影響 1 當w 0時 主瓣位于積分區(qū)間內 隨著w的移動不同大小的正 負旁瓣移出或移入積分區(qū)間 使得H ejw 的大小產生波動 主瓣附近窗的頻率響應為 隨著N的加大 振蕩變密 主瓣變窄 主瓣與旁瓣的幅度亦有所加大 但主瓣與旁瓣的相對比例不變 吉布斯現(xiàn)象 個負值卻還在此區(qū)間內 使得H ejw 取值最小值 0 0895H ej0 稱為下臂峰 5 w p時 H ejw 隨著區(qū)間內旁瓣的移動而在阻帶內波動 另外 圖7 16表示了0到p范圍內H ejw 變化的情況 0到 p的圖形變化與此對稱 如圖7 17 且以2p為周期 途中假定H ej0 1 在 為過渡帶 3 當w wc時 即主瓣中心移到了wc處 此時 7 17 由圖可見 加矩形框后得到的濾波器的頻響與理想頻響之間存在差異 表現(xiàn)出肩峰 過渡帶及在通帶和阻帶內的波動 只有肩峰和波動盡可能小 而且過渡帶盡可能窄 才能更接近理想特性 綜上 窗口法設計FIR濾波器 h n 長度N增大可使過度帶變窄 而所選窗函數(shù)不僅影響過渡帶的寬度 還能影響肩峰和波動的大小 光束 7 3 3幾種常用窗函數(shù) 2 三角形窗 BartlettWindow 其頻率響應為 主瓣寬度為 1 矩形窗 其頻率響應和幅度響應分別為 是三項矩形窗的幅度響應的移位加權和 它使旁瓣相互抵消 能量更集中在主瓣 但主瓣寬度比矩形窗的主瓣加寬了一倍 為 3 漢寧 Hanning 窗 又稱升余弦窗 下圖為N 31時 矩形窗 三角窗 漢寧窗 漢明窗及布萊克曼這5種窗口函數(shù)的包絡曲線 下圖為N 51時矩形窗 漢寧窗 漢明窗及布萊克曼4種窗口函數(shù)的幅度響應 下圖為N 5時用矩形窗 漢寧窗 漢明窗及布萊克曼設計的低通濾波器的幅度響應 5 凱塞 Kaiser 窗 是一個可選參數(shù) 用來選擇主瓣寬度和旁瓣衰減之間的交換關系 一般說來 越大 過渡帶越寬 阻帶越小衰減也越大 I0 是第一類修正零階貝塞爾函數(shù) 一般取15 25項就可滿足精度要求 若阻帶最小衰減表示為As 20lg s 的確定可采用以下經驗公式 凱澤窗 各種窗函數(shù)的主要性能 7 3 4 窗函數(shù)法的設計1 設計步驟 1 給定頻響函數(shù) 2 求出單位抽樣響應 3 根據(jù)過渡帶寬度和阻帶最小衰減 借助窗函數(shù)基本參數(shù)表 P202表3 確定窗的形式及N的大小 4 最后求及2 設計舉例 例 分別利用矩形窗與漢寧窗設計具有線性相位的FIR低通濾波器 具體要求 其他 并畫出相應的頻響特性 解 1 由于是一理想LF 所以可以得出 2 確定N由于相位函數(shù) 所以呈偶對稱 其對稱中心為 因此 3 加矩形窗 則有 可以求出h n 的數(shù)值 注意偶對稱 對稱中心 由于h n 為偶對稱 N 25為奇數(shù) 所以 例如H 0 0 94789 可以計算的值 畫如下圖 4 加漢寧窗由于可以求出序列的各點值 通過可求出加窗后的h n 相應幅度函數(shù)可用下式求得 如H 0 0 98460 圖如下 2 求hd n 4 確定N值 3 選擇窗函數(shù) 由確定海明窗 53dB 5 確定FIR濾波器的h n 6 求 驗證 若不滿足 則改變N或窗形狀重新設計 5 線性相位FIR高通濾波器的設計 其單位抽樣響應 理想高通的頻響 6 線性相位FIR帶通濾波器的設計 其單位抽樣響應 理想帶通的頻響 7 線性相位FIR帶阻濾波器的設計 其單位抽樣響應 理想帶阻的頻響 4 設計舉例 利用凱澤窗設計一FIR低通filter 要求 解 經驗公式 取38 將N 38 5 653代入表達式 得 0370 01 0000 02040 02 1361 83362 0300 04150 04 2352 55683 3450 07040 07 8294 654819 960 40820 41 3343 0865 2510 10740 11 4333 51117 4410 15220 15 5323 865610 110 20670 21 6314 167813 100 26790 29 7304 428616 440 33620 34 9284 851223 830 48730 49 0 4 8 12 16 18 19 25 29 33 37 21 的圖形如右所示 7 4 頻率抽樣法 一 設計思想 窗函數(shù)設計法是從時域出發(fā) 把理想的用一定 形狀的窗函數(shù)截取成有限長的 以來近似 從而使頻響近似理想頻響 頻率取樣法是從頻域出發(fā) 對理想的頻響 進行等間隔取樣 以有限個頻響采樣去近似理想頻響 即 等間隔取樣 并且 二 利用N個頻域采樣值重構FIR的系統(tǒng)函數(shù)與頻響 1 重構FIR的的單位抽樣響應h n 根據(jù)頻域抽樣理論 由N個頻域采樣點 可以唯一確定h n 即對H k 進行IDFT 2 重構系統(tǒng)函數(shù)H Z 3 FIR的頻響 將代入表達式可得 其中 為大家所知的內插函數(shù) 分析可知 當時 采樣點 有 這說明 重構的頻響 在采樣上嚴格等于H k 而在采樣點之間 頻響則由加權的內插函數(shù)延伸疊加而成 三 線性相位的約束條件 以h n 為偶對稱 N為奇數(shù)的情況進行分析 1 FIR的頻響具有線性相位的一般表達式 當h n 為偶對稱 N為奇數(shù)時 則 而且幅度函數(shù)應為偶對稱 即 2 采樣值H k 具有線性相位的約束 其中 表示采樣值的模 純標量 表示 其相角 因此 在采樣點上具有線性相位的條件應為 而且 必須滿足偶對稱 即 四 設計步驟 1 根據(jù)指標要求 畫出頻率采樣序列的圖形 2 依據(jù)的對稱特點 可以使問題得以簡化 3 根據(jù)線性相位的約束條件 求出 4 將代入FIR的頻響表達式 5 由的表達式畫出實際頻響 四 設計舉例 例 試用頻率采樣法 設計一個具有線性相位 的低通FIR數(shù)字filter 其理想頻率特性為 已知 采樣點N 33 由于h n 為偶對稱 且N 33為奇數(shù) 所以對于 是偶對稱 所以上圖可畫一半 到 截止頻率 即 解 相位約束條件 而為 將代入FIR的頻響 得 考慮到時 所以將負頻部分加進去 有 的圖形如下所示 0 五 IIR和FIR數(shù)字濾波器的比較 IIR濾波器 FIR濾波器 h n 無限長 h n 有限長 極點位于z平面任意位置 濾波器階次低 非線性相位 遞歸結構 不能用FFT計算 可用模擬濾波器設計 用于設計規(guī)格化的選頻濾波器 極點固定在原點 濾波器階次高得多 可嚴格的線性相位 一般采用非遞歸結構 可用FFT計算 設計借助于計算機 可設計各種幅頻特性和相頻特性的濾波器