FIR数字滤波器.ppt
第六章 FIR數位濾波器設計,本章重點：本章探討的濾波器架構稱為Finite impulse response(FIR)有別於IIR的設計由於其不具有回授項系統一定穩定因此不必探討極點是否落在Z平面的單位圓內而是探討其線性相位的特性。FIR濾波器的零點的位置與線性相位的關係。各種窗形函數介紹以及其設計方法。最佳化濾波器設計法。採用Parks-McClellan 演算法以MATLAB之remez指令驗證。,FIR數位濾波器結構,FIR為Finite Impulse response的縮寫由於其設計的架構沒有回授的元件所以為有限脈衝響應。其架構一般採用直接法或串接法。,6.1.1 直接形式(Direct form),此種形式如圖6.1-1。,圖6.1-1 FIR直接形式結構圖,6.1.2 串聯形式(Cascade form),此種形式可由二組直接形式串聯而成如圖6.1-2。,圖6.1-2 FIR串聯形式結構圖,6.1.3 線性相位形式(Linear Phase form),所謂線性相位具備線性相位特性的FIR濾波器之差分方程式其係數呈現對稱(Symmetry)的現象例如以長度為5的FIR濾波器為例其乘法次數由5次減為5/2之商為2再加1等於3次。其架構如圖6.1-3。,圖6.1-3 FIR線性相位形式結構圖,線性相位(Linear-phase)FIR特性,所謂線性相位(Linear-phase)FIR指其頻率響應的相位特性呈現直線的關係而不是弧線。,6.2.1 線性相位FIR的對稱與反對稱,對一個長度M的FIR濾波器假如其相位大小具備線性關係則稱為線性相位FIR濾波器。其形式依據脈衝響應的對稱性以及M的奇數或偶數而不同。,圖6.2-1 反對稱奇數脈衝序列與相位圖,圖6.2-2 反對稱偶數脈衝序列與相位圖,圖6.2-3 對稱奇數脈衝序列與相位圖,圖6.2-4 對稱偶數脈衝序列與相位圖,6.2.2 線性相位FIR的零點位置,對一個長度M的線性相位FIR濾波器總共有M-1個零點與M-1個假極點(Trivial poles)。其中的零點具有兩兩對稱的特性由零點的對稱性質可以判定FIR是否為線性相位。,圖6.2-5 反對稱奇數脈衝序列與零點位置圖,圖6.2-6 反對稱偶數脈衝序列與零點位置圖,圖6.2-7 對稱奇數脈衝序列與零點位置圖,圖6.2-8 對稱偶數脈衝序列與零點位置圖,窗形函數(Window Function)設計法,為了去除FIR濾波器在通帶與拒帶的頻率響應震盪有必要把不要的部份予以截除(Truncate)而保留需要的部份使用的方法簡單的說就是把希望的響應與一個適當的窗形函數進行環形摺積使得處理後的響應能夠滿足規格要求。,6.3.1 窗形函數的環形摺積,對於一個希望得到的(Desired)濾波器頻率響應其具有在通帶時為線性相位在拒帶則完全阻絕信號通過。示意圖如圖6.3-1。,圖6.3-1 窗形函數設計法示意圖,6.3.2 窗形函數的種類,常用的窗形函數有Rectangular window、Bartlett window、Hanning window、Hamming window、Blackman window、Kaiser window。,圖6.3-2 窗形函數圖之一,圖6.3-3 窗形函數圖之二,圖6.3-4 窗形函數圖之三,圖6.3-5 窗形函數頻率響應圖之一,圖6.3-6 窗形函數頻率響應圖之二,圖6.3-7 窗形函數頻率響應圖之三,6.3.3 以指令fir1設計FIR濾波器,MATLAB指令fir1只要給定濾波器的長度M、截止頻率w、（或增加窗形函數種類win、濾波器形式f_type）使用它即可輕鬆的滿足設計上的需求。,圖6.3-8 Hamming 低通FIR濾波器,圖6.3-9 Blackman高通FIR濾波器,圖6.3-10 Hanning帶通FIR濾波器,圖6.3-11 Hanning帶拒FIR濾波器,6.3.4 Kaiser window之FIR濾波器設計,Kaiser window是最好用的窗形函數設計過程中必須先算兩個參數其中一個為值目的是要滿足拒帶衰減參數As的需求；另一個則為FIR濾波器的長度M這是我們設計的目的。,圖6.3-12 Kaiser低通FIR濾波器,圖6.3-13 Kaiser高通FIR濾波器,圖6.3-14 Kaiser帶通FIR濾波器,圖6.3-15 Kaiser帶拒FIR濾波器,6.3.5 以指令fir2設計任意形狀(Arbitrary shape)FIR濾波器,fir2類似fir1的指令但是它不用來設計高通、低通、帶通、帶拒濾波器而是設計多頻段(Multilevel)的濾波器。,圖6.3-16 以fir2設計多頻段FIR濾波器,圖6.3-17 以fir2設計多頻段FIR濾波器(階數不夠),最佳化濾波器設計,所謂最佳化就是讓濾波器的頻率響應在衰減帶的起伏(Ripple)等量平均的變化稱為Equiripple。使用的演算法稱為Parks-McClellan algorithm。,6.4.1 最佳化濾波器設計,最佳化濾波器的設計步驟,圖6.4-1 以remez指令設計低通FIR濾波器,圖6.4-2 以remez指令設計高通FIR濾波器,圖6.4-3 以remez指令設計帶通FIR濾波器,圖6.4-4 以remez指令設計帶拒FIR濾波器,6.4.2 微分器設計,remez指令設計頻譜微分器同上一小節介紹的指令格式在尾端加上differentiator的宣告即可。,圖6.4-5 以remez指令設計FIR微分器,6.4.3 Hilbert形式的濾波器設計,同小節介紹的指令格式在尾端加上hilbert的宣告即可執行Hilbert transformer的設計。,圖6.4-6 以remez指令設計帶通FIR濾波器(Hilbert transformer),
