FIR數(shù)字濾波器設計和實現(xiàn)課件

1、概述線性相位 FIR DF 約束條件和頻率響應窗函數(shù)法頻率取樣法FIR數(shù)字濾波器的實現(xiàn)結構 第五章FIR 數(shù)字濾波器設計和實現(xiàn)1概述：IIR 和 FIR 比較IIR與FIR性能比較IIR數(shù)字濾波器幅頻特性較好，但相頻特性較差 FIR數(shù)字濾波器可以嚴格線性相位，又可任意幅度特性因果穩(wěn)定系統(tǒng)可用 FFT 計算 (計算兩個有限長序列的線性卷積)但階次比 IIR 濾波器要高得多2概述：IIR 和 FIR 比較IIR 與 FIR 設計方法比較IIR DF無限沖激響應，H(Z) 是 z-1 的有理分式，借助于模擬濾波器設計方法，階數(shù)低（同樣性能要求）。其優(yōu)異的幅頻特性是以非線性相位為代價的。缺點：只能設計

2、特定類型的濾波器，不能逼近任意的頻響。 FIR DF有限沖激響應，系統(tǒng)函數(shù) H(Z) 是 z-1 的多項式，采用直接逼近要求的頻率響應。設計靈活性強。缺點： 設計方法復雜； 延遲大； 階數(shù)高。 (運算量比較大，因而在實現(xiàn)上需要比較多的運算單元和存儲單元) FIR DF 的技術要求通帶頻率p，阻帶頻率s 及最大衰減p，最小衰減s很重要的一條是保證 H(z) 具有線性相位。3概述：FIR DF 設計方法FIR 數(shù)字濾波器設計 FIR 濾波器的任務給定要求的頻率特性，按一定的最佳逼近準則，選定 h(n) 及階數(shù) N。三種設計方法 窗函數(shù)加權法 頻率采樣法 FIR DF 的 CAD - 切比雪夫等波紋

3、逼近法4概述：FIR DF 零極點 FIR濾波器的I/O 關系： FIR 濾波器的系統(tǒng)傳遞函數(shù)： 在 Z 平面上有 N-1 個零點；在原點處有一個（N-1）階極點，永遠穩(wěn)定。FIR 系統(tǒng)定義：一個數(shù)字濾波器 DF 的輸出 y(n)，如果僅取決于有限個過去的輸入和現(xiàn)在的輸入x(n), x(n-1),. ., x(n-N+1)，則稱之為 FIR DF。 FIR 濾波器的單位沖激響應：5 FIR DF 的頻率響應為：FIR 濾波器的最重要特點是能實現(xiàn)線性相位。具有線性相移特性的 FIR 濾波器是 FIR 濾波器中應用最廣泛的一種。Hr()：振幅響應，它是一個取值可正可負的實函數(shù)。 () = arg

4、H(ejw) 為數(shù)字濾波器的相位響應。 概述：FIR DF 頻率響應6信號通過線性濾波器時，其幅度和相位可能會發(fā)生改變，濾波器幅頻特性 |H()|和相頻特性 () 可能會隨頻率的變化而改變。如：輸入正弦信號 Acos(n0) 則：輸出為 |H(0)| Acos(n0)，其中相移(0) 輸出頻率和輸入頻率相同，但幅度和相位都發(fā)生了變化輸出信號比輸入信號滯后的樣點數(shù) n (位移) 可由下式求得: 設：n00 濾波器在數(shù)字頻率0 處的相位延遲（位移） 由于相位延遲 n 的不同，最終產(chǎn)生了相位失真。 確保不產(chǎn)生相位失真的辦法：使不同頻率的信號通過濾波器時有相同的延遲 n。 概述：相位失真7對不同的頻率

5、有恒定的相移，不同的相位延遲 n，會產(chǎn)生相位失真. 如：方波 y(t) 可以用無數(shù)奇次諧波的正弦波的疊加來得到： 若每個正弦波相移/2 弧度： 確保所有頻率具有相同相位延遲的簡單方法： 隨著頻率的變化而改變相位，使濾波器具有線性相位特性，即使所有頻率的相位延遲保持恒定，這種方法可通過使系統(tǒng)的相位函數(shù)()為頻率的線性函數(shù)來實現(xiàn)。 概述：相位失真可見相移之后正弦波之和已不再是方波。8線性相移FIR DF 約束條件和頻率響應三個內(nèi)容 約束條件恒延時濾波h(n) 偶對稱：恒相延時和恒群延時同時成立 h(n) 奇對稱：僅恒群延時成立 頻率響應Type I：h(n) 偶對稱、N 為奇數(shù)Type II：h(

6、n) 偶對稱、N 為偶數(shù)Type III：h(n) 奇對稱、N 為奇數(shù)Type IV：h(n) 奇對稱、N 為偶數(shù) FIR DF 零極點分布9相延時： 群延時： 線性相移FIR DF 約束條件：恒延時濾波恒延時濾波 濾波器的延時分為相延時和群延時兩種令恒延時濾波器：p() 或g() 是不隨變化的常量，這時濾波器具有線性相位特性。10（負號是因為系統(tǒng)必有時延） 由于 FIR 濾波器的傳遞函數(shù)為 :w(w)0故: 線性相移FIR DF 約束條件：恒延時恒相延時和恒群延時同時成立要使p、g 都不隨 變化,() 必須是一條過原點直線11于是: 線性相移FIR DF 約束條件：恒延時12可以證明，當 線

7、性相移FIR DF 約束條件：恒延時上式成立，此時恒相延時和恒群延時同時成立時，線性相位濾波器的必要條件是： 不管 N 為偶數(shù)，還是 N 為奇數(shù)，系統(tǒng)沖激響應 h(n) 都關于中心點 (N-1)/2 偶對稱。當 N 為奇數(shù)時對稱中心軸位于整數(shù)樣點上; 當 N 為偶數(shù)時對稱中心軸位于非整數(shù)樣點上。h(n) 為偶對稱，N 為偶數(shù)07nh(n)h(n) 為偶對稱，N 為奇數(shù)06nh(n)13于是有：線性相移FIR DF 約束條件：恒群延時只要求恒群延時成立 若只要求群延時g() 為一常數(shù)，則相移特性為不過原點的直線。0()故14可以證明，當 上式成立，此時故線性相移FIR DF 約束條件：恒群延時1

8、5FIR濾波器單獨滿足恒定群延時的必要條件為： 沖激響應 h(n) 對中心點 (N-1)/2 成奇對稱。此時，無論 N 為奇數(shù)或偶數(shù)，濾波器的相頻特性均為線性，并包含有/2 的固定相移： 因此，信號通過此類濾波器時不僅產(chǎn)生 (N-1)/2 個取樣點的延遲，還將產(chǎn)生 90o 的相移，通常這類濾波器又被稱為 90o 移相器，并具有很好的應用價值。當 N 為奇數(shù)時，故07h(n) 為奇對稱，N 為偶數(shù)nh(n)06h(n) 為奇對稱，N 為奇數(shù)nh(n)線性相移FIR DF 約束條件：恒群延時16 奇對稱：() 對所有的頻率成分都有一個 90相移。因此，有四種類型的 FIR DF： 線性相移 FIR

9、 DF 約束條件線性相位約束條件對于任意給定的值 N，當 FIR 濾波器的 h(n) 相對其中心點 (N-1)/2 是對稱時，不管是偶對稱還是奇對稱，此時濾波器的相移特性是線性的，且群延時都是 = (N-1)/2 。偶對稱 ： () 為過原點的，斜率為 - 的一條直線17線性相移 FIR DF 頻率響應：Type Ih(n) 偶對稱，N 為奇數(shù)（恒相時延、恒群時延此時，由于 h(n) 序列的長度為奇數(shù)，因此濾波器的頻率響應函數(shù)可進行以下拆分（前后對稱部分、中心點）：h(n) 為偶對稱，N 為奇數(shù)06nh(n)對上式的第二和式作變量替換（n=N-1-m) 后得到:由對稱條件則 H(ej) 表示為

10、：18線性相移 FIR DF 頻率響應：Type I令 則上式為 19 由此可以看出其線性相位特性。由于 cos(n) 對于 =0、2都是偶對稱，所以幅度響應 Hr() 對=0、2也是偶對稱。線性相移 FIR DF 頻率響應：Type I其中振幅響應：相頻響應：N=9Hr (w)20h(n) 偶對稱，N 為偶數(shù)（恒相時延、恒群時延由于h(n) 序列的長度為偶數(shù)，因此濾波器的頻率響應函數(shù)可拆分成如下兩部分（前后對稱部分，中心點處無值）：線性相移 FIR DF 頻率響應：Type IIh(n) 為偶對稱，N 為偶數(shù)07nh(n)對上式的第二和式作變量替換（n=N-1-m) 后得到:由對稱條件則 H

11、(ej) 表示為：21 線性相移 FIR DF 頻率響應：Type II令 ，則上式為： 其中（注意 n 從1 開始，即 b(0)=0，或沒有定義） 22線性相移 FIR DF 頻率響應：Type II 與所設計的 b(n) 或 h(n) 無關，恒為 0。這種類型（即 h(n) 偶對稱，N為偶數(shù)）不能用于高通或帶阻濾波器。 （2）由于 cos(n-1/2) 對于 =是奇對稱，所以，Hr(w) 對 =也是奇對稱；以 =0、2為偶對稱。 振幅響應：相頻響應：N=8n 從1開始Hr (w)注意： （1） 在 = 處，有：23h(n) 奇對稱，N 為奇數(shù)（恒群時延h(n) 長度為奇數(shù)，拆分成前后兩部分

12、：線性相移 FIR DF 頻率響應：Type III對上式的第二和式作變量替換，并利用對稱條件 h(n)=-h(N-1-n)，得:06h(n) 為奇對稱，N 為奇數(shù)nh(n)24Hr (w)線性相移 FIR DF 頻率響應：Type III，則上式為： 其中令振幅響應：相頻響應：n 從1開始25與 c(n) 或 h(n) 的值無關，因此，這種類型的濾波器不適用于低通、帶阻或高通濾波器設計，而且 ，這說明 jHr(w) 是純虛數(shù)，對于逼近理想數(shù)字希爾伯特變換和微分器，它是很有用的。理想的希爾伯特變換是一個全通濾波器，它對輸入信號產(chǎn)生 90 度的相移，它頻繁用于通信系統(tǒng)中的調(diào)制。微分器廣泛用于模擬

13、和數(shù)字系統(tǒng)中對信號求導。 (2)由于 sin(n) 對于 =0、2 都是奇對稱，所以，Hr(w) 以 =0、 、2為奇對稱。 注意： (1) 在 = 0 和 處，有：線性相移 FIR DF 頻率響應：Type III26線性相移 FIR DF 頻率響應：Type IVh(n) 奇對稱，N 為偶數(shù)（恒群時延07h(n) 為奇對稱，N 為偶數(shù)nh(n)其中27線性相移 FIR DF 頻率響應：Type IVHr (w)與 d(n) 或 h(n) 的取值無關，因此傳輸函數(shù) H(z) 在 z = 1 處為零點。顯然，這種類型不能用于實現(xiàn)低通濾波器。又有 ，所以這類濾波器適用于設計希爾伯特變換和微分器。

14、注意： (1) 在 =0 處，有：(2)由于 sin(n-1/2) 在 =處偶對稱，在0、2 是奇對稱，所以，Hr(w) 以 = 偶對稱，0、2為奇對稱。 28一般形式： 偶對稱： 奇對稱： （兩個恒時延條件） （一個恒時延條件） ( Hr() 為 的實函數(shù) )線性相移 FIR DF 頻率響應：小結29四類線性相位FIR濾波器第一類FIR 系統(tǒng)是 的線性組合，在 時，易取得最大值，因此這一類濾波器易體現(xiàn)低通特性，且是偶函數(shù)。通過頻率移位，又可體現(xiàn)高通、帶通、帶阻特性。所以，經(jīng)典的低通、高通、帶通和帶阻濾波器的 都是偶對稱的。第三、四類FIR 系統(tǒng)是 的線性組合，在 時， 的值為零，且是奇函數(shù)。

15、這一類濾波器都是作為特殊形式的濾波器，如Hilbert變換器、差分器等。請使用時注意。N最好取為奇數(shù)，以便以中心點為對稱。30 一般的 FIR DF 的零、極點： 在z=0處，有一個（N-1）階的極點，故濾波器穩(wěn)定；其零點要求 f(z)=0，根據(jù)代數(shù)理論，它為 N-1階多項式，應有 N-1 個根，所以有 N-1 個零點。如果 h(n) 為實數(shù)值，其根肯定是共軛對稱的。線性相移 FIR DF 零極點分布31令：m=N-1-n于是： 線性相移 FIR DF 零極點分布線性相移 FIR DF 的零極點： 如果 zi 是 H(z) 的零點，即 H(zi) = 0 則 H(z-1) =0，即 zi-1

16、亦為 H(z) 的零點。 32上面提到 zi 肯定是共軛的，故 zi* 亦必為其零點于是零點有：1-1Za1Za21/bb線性相移 FIR DF 零極點分布總結:(1) 一般情況， ，有四個零點： (2) r=1，單位圓上的零點： (共軛對) (3) 位于實軸上的實數(shù)：b, 1/b (實軸上的倒數(shù)對)。(4) zi =1：單零點 33例 ： 設FIR濾波器的系統(tǒng)函數(shù)為： 求出該濾波器的單位取樣響應h(n)，判斷是否具有線性相位，求出其幅度特性和相位特性。解 ：對FIR數(shù)字濾波器，其系統(tǒng)函數(shù)為：34所以，該FIR濾波器具有第一類線性相位特性，設其頻率響應函數(shù)為幅度特性函數(shù)為：相位特性函數(shù)為：35

17、思路：理想數(shù)字濾波器設計的 FIR 數(shù)字濾波器要求：線性相位盡可能降低逼近誤差FIR DF 窗口法（傅里葉級數(shù)法）hd(n) 無限長，且非因果 h(n) 有限長，且因果 36設所要求的 DF 的頻率響應是 Hd(ej)，需要注意：它可能是低通、高通、帶通和帶阻 FIR DF，沒有特指某種類型的數(shù)字濾波器。不管是何種 FIR DF，它的頻率響應是頻域中的周期函數(shù)，周期為 2，所以它可以展開為傅氏級數(shù)形式：窗口法：基本原理 式中 hd(n) 是傅里葉系數(shù)，也是單位取樣響應序列。 由傅里葉級數(shù)理論可得： 37因此，所要求的 DF 的系統(tǒng)函數(shù)便可求得：顯然，Hd(z) 是非因果的，且 hd(n) 的持

18、續(xù)時間為 - +，物理上不可實現(xiàn)。我們可以采用逼近 Hd(ej) 的方法 首先把 hd(n) 先截短為有限項，把 hd(n) 截為2M+1項，得：窗口法：基本原理38然后把截短后的 hd(n) 右移，使之變成因果性的序列。 令 H(z) 等于 H1(z) 乘以 z-M 得：令 h(n)= hd (n-M), n=0, 1, 2, ., 2M，則 頻率響應 z=ej窗口法：基本原理顯然H(z) 是物理可實現(xiàn)的其沖激響應 h(n) 的持續(xù)時間也是有限的選擇 hd(n) = hd(N-1-n)，保證H(z) 具有線性相位。 39對 hd(n) 的截短必然產(chǎn)生誤差，即以 |H(ej)| 近似 |Hd(

19、ej)| 。定義逼近誤差為均方誤差：而 Hd(ej) 可以展開為：式中:窗口法：性能分析 |H(ej)| 對 |Hd(ej)| 的逼近40 因為 |H(ej)| 是對 hd(n) 截短而產(chǎn)生的，假定： 即當 |n|M 時，An = 0, Bn =0。所以把上述兩式代入逼近誤差中，利用三角函數(shù)的正交性可得：由于上式中每一項都是正的，所以，只有當 最小。 窗口法：性能分析41說明： 當用 |H(ej)|Hd(ej)| 時，要使 2 =min， |H(ejw)| 的截短后的單位取樣響應 h(n) 的系數(shù)必須等于所要求的幅頻響應 |Hd(ejw)|展成傅里葉級數(shù)的系數(shù) hd(n)。 有限項傅氏級數(shù)是在

20、最小均方意義上對原信號的最佳逼近 其逼近誤差為： 截短的長度 M 越大，逼近誤差2 愈?。ㄒ驗?hd(n) 值愈小）。 窗口法：性能分析42將 hd(n) 截短： 相當于將 hd(n ) 與一窗函數(shù) wR(n) 相乘，即窗口法：Gibbs 效應其中在一定意義上來看，窗函數(shù)決定了我們能夠 “看到” 多少個原來的沖激響應，“窗” 這個用詞的含義也就在此。43窗函數(shù)的頻譜 窗口法：Gibbs 效應此矩形窗譜為一鐘形偶函數(shù)，在 +2/N 之間為其主瓣，主瓣寬度 = 4/N，在主瓣兩側有無數(shù)幅度逐漸減小的旁瓣, 見圖所示。2/N-2/N主瓣第1個旁瓣第2個旁瓣44 截短，根據(jù)時域相乘映射為頻域卷積，得：

21、窗口法：Gibbs 效應為便于分析，我們假定 |Hd(ejw)| 是理想低通濾波器 LPF。式中積分等于由 c 到 c 區(qū)間內(nèi) WNej(w-) 下的面積，隨著變化，窗函數(shù)的主瓣和不同正負、不同大小的旁瓣移入和移出積分區(qū)間，使得此面積發(fā)生變化， 也即 |H(ejw)| 的大小產(chǎn)生波動。-wc0wc450WR()-c0cHd()0.50.50.08950.08950.04680.0468卷積窗口法：Gibbs 效應-wc0wc-Hd()46現(xiàn)在分析幾個特殊頻率點的濾波器性能 = 0 時： 由于一般情況下都滿足 c 2 / N，因此，H(0) 的值近似等于窗譜函數(shù) WR(ejw) 與軸圍出的整個面

22、積。 窗口法：Gibbs 效應0WR()-wc0wc-Hd() =c 時： 此時窗譜主瓣一半在積分區(qū)間內(nèi)一半在區(qū)間外，因此，窗譜曲線圍出的面積，近似為=0 時所圍面積的一半，即 。-wc0wcHd()w =wcWR(w-)47=c - 2/N 時，正肩峰 此時窗譜主瓣全部處于積分區(qū)間內(nèi)，而其中一個最大負瓣剛好移出積分區(qū)間，這時得到最大值，形成正肩峰。之后，隨著值的不斷增大，H(ejw) 的值迅速減小，此時進入濾波器過渡帶。=c + 2/N 時，負肩峰 此時窗譜主瓣剛好全部移出積分區(qū)間，而其中一個最大負瓣仍全部處于區(qū)間內(nèi)，因此得到最小值，形成負肩峰。之后，隨著值的繼續(xù)增大， H(ejw) 的值振

23、蕩并不斷減小，形成濾波器阻帶波動。-wc0wcHd()WR(w-)-wc0wcHd()WR(w-)48 理想濾波器的不連續(xù)點演化為過渡帶 通帶與阻帶內(nèi)出現(xiàn)起伏過渡帶：正負肩峰之間的頻帶。其寬度等于窗口頻譜的主瓣寬度。 對于矩形窗 WR(ejw), 此寬度為 4/N。 肩峰及波動：這是由窗函數(shù)的旁瓣引起的。旁瓣越多，波動越快、越多。相對值越大，波動越厲害，肩峰越強。肩峰和波動與所選窗函數(shù)的形狀有關，要改善阻帶的衰減特性只能通過改變窗函數(shù)的形狀。 加窗處理對理想矩形頻率響應的影響窗口法：Gibbs 效應49 Gibbs 現(xiàn)象在對 hd(n) 截短時，由于窗函數(shù)的頻譜具有旁瓣，這些旁瓣在與 Hd(e

24、jw) 卷積時產(chǎn)生了通帶內(nèi)與阻帶內(nèi)的波動，稱為吉布斯現(xiàn)象。長度 N 的改變只能改變 坐標的比例及窗函數(shù) WR(ejw) 的絕對大小，但不能改變肩峰和波動的相對大?。ㄒ驗椴荒芨淖兇昂瘮?shù)主瓣和旁瓣的相對比例，波動是由旁瓣引起的），即增加 N，只能使通、阻帶內(nèi)振蕩加快，過渡帶減小，但相對振蕩幅度卻不減小。結論：過渡帶寬度與窗的寬度 N 有關，隨之增減而變化。 阻帶最小衰減（與旁瓣的相對幅度有關）只由窗函數(shù) 決定，與 N 無關。 窗口法：Gibbs 效應50Gibbs現(xiàn)象；窗口法：Gibbs 效應51設計FIR DF時，窗函數(shù)不僅可以影響過渡帶寬度，還能影響肩峰和波動的大小(阻帶的衰減)，因此，選擇窗

25、函數(shù)應使其頻譜：主瓣寬度盡量小，以使過渡帶盡量陡。旁瓣相對于主瓣越小越好，這樣可使肩峰和波動減小，即能量盡可能集中于主瓣內(nèi)。對于窗函數(shù)，這兩個要求是相互矛盾的，要根據(jù)需要進行折衷的選擇，窗口法：常用窗函數(shù)w020lg|W(w)/W(0)|B3dBAD (dB/Oct)為了定量地比較各種窗函數(shù)的性能，給出三個頻域指標：3db 帶寬 B，單位為 （最大可能的頻率分辨力）最大旁瓣峰值 A(dB)，A 越小，由旁瓣引起的譜失真越小旁瓣譜峰漸進衰減速度 D（dB/oct）一個好的窗口，應該有最小的 B、A 及最大的 D。52以下介紹的窗函數(shù)均為偶對稱函數(shù)，都具有線性相位特性。 設窗的寬度為N，窗函數(shù)的對

26、稱中心點在(N-1)/2處。因此，均為因果函數(shù)。矩形窗最簡單的窗函數(shù)，從阻帶衰減的角度看，其性能最差。它的頻率響應函數(shù)為：窗口法：基本窗函數(shù)_矩形窗振幅響應53為了對過渡帶和阻帶衰減進行精確分析，對窗振幅響應進行連續(xù)積分（或累積振幅響應），即矩形窗函數(shù) w(n) 以及它的振幅響應、累積振幅響應如下圖所示。窗口法：基本窗函數(shù)_矩形窗性能指標3dB 帶寬 B=0.89最大旁瓣峰值 A= -13dB旁瓣譜峰漸進衰減速度 D=-6dB/oct在 Matlab 中，實現(xiàn)矩形窗的函數(shù)為 w = boxcar(n)。54振幅響應在 = 1 處具有第一個零點：因而主瓣的寬度為 2，所以過渡帶寬也近似為 2。大

27、約在 = 3/N 處，出現(xiàn)第一個旁瓣（即主旁瓣），其幅度為：將它與主瓣振幅 N 比較，則最大旁瓣峰值A(dB) 為 A= -13db。 累積振幅響應第一個旁瓣為 21dB，這個 21dB 的阻帶衰減與窗長度 N 無關。根據(jù)最小阻帶衰減，可以精確地計算出過渡帶寬為：它大約是近似帶寬的一半。 窗口法：基本窗函數(shù)_矩形窗55三角窗（或 巴特利特 Bartlett 窗）由于矩形窗從 0 到 1 （或 1 到 0 ）有一個突變的過渡帶，這造成了吉布斯現(xiàn)象。Bartlett 提出了一種逐漸過渡的三角窗形式，它是兩個矩形窗的卷積。B=1.28, A=-27dB, D= -12dB/oct ，近似過渡帶寬 8

28、/N，精確過渡帶寬 6.1/N，最小阻帶衰減 25dB。與矩形窗來比較，阻帶衰減性能有所改善，但代價是過渡帶的加寬。窗口法：基本窗函數(shù)_三角窗56在 Matlab 中，函數(shù) bartlett(n) 和 triang(n) 用來計算相似的三角窗，但它們有兩個重要的區(qū)別：bartlett 函數(shù)返回的序列兩端總是 0，因此，對于奇數(shù) n，語句 bartlett(n+2) 的中間部分等于 triang(n)；對于偶數(shù) n，bartlett 仍然是兩個矩形序列的卷積，但 n 為偶數(shù)時的三角窗沒有標準定義。5.3.3 窗口法：基本窗函數(shù)_三角窗57余弦窗B=1.2, A=-23dB, D=-12dB/oc

29、t。近似過渡帶寬 8/N，精確過渡帶寬6.5/N，最小阻帶衰減 34dB。窗口法：基本窗函數(shù)_余弦窗或其中頻率響應58窗口法：基本窗函數(shù)_余弦窗59升余弦窗函數(shù)漢寧窗、漢明窗、布萊克曼窗都是升余弦窗的特例。它們都是頻率為 0 2/(N-1) 和 4/(N-1) 的余弦序列的組合。升余弦窗的頻率特性比矩形窗有很大改善。其中 A、B、C 為常數(shù)。當 A = 0.5，B=0.5，C=0 時，為漢寧 (Hanning)窗。 Matlab 中，w = hanning(n)當 A = 0.54，B=0.46，C=0 時，為漢明 (Hamming) 窗。 Matlab 中，w = hamming(n)當 A

30、 = 0.42，B=0.5，C=0.08 時，為布萊克曼窗。 Matlab 中，w = blackman(n)窗口法：升余弦窗函數(shù)60窗口法：升余弦窗函數(shù)61W()0.5u()0.25u(-2/(N-1)0.25u(+2/(N-1)Hanning 窗（升余弦窗）窗口法：升余弦窗函數(shù)_漢寧窗B=1.44, A=-32db, D=-18db/oct，近似過渡帶寬 8/N，精確過渡帶寬 6.2/N，最小阻帶衰減 44dB。與矩形窗來比，最小阻帶衰減性能明顯提高，但過渡帶也明顯增大。 62窗口法：升余弦窗函數(shù)_漢寧窗63Hamming 窗（改進的升余弦窗）窗口法：升余弦窗函數(shù)_漢明窗B=1.3, A=

31、-43dB, D=-6dB/oct，近似過渡帶寬 8/N，精確過渡帶寬6.6/N，最小阻帶衰減 53dB。通過這一系數(shù)調(diào)整，使能量的 99.963% 都集中在了窗譜的主瓣內(nèi)。64窗口法：升余弦窗函數(shù)_漢明窗65Blackman 窗（二階升余弦窗）窗口法：升余弦窗函數(shù)_布萊克曼窗B=1.68, A=-58db, D=-18db/oct，近似過渡帶寬 12/N，精確過渡帶寬11/N，最小阻帶衰減 74dB。通過增加余弦的二次諧波分量，能夠進一步抑制旁瓣，但主瓣寬度卻比矩形窗譜的主瓣寬度大三倍。66窗口法：升余弦窗函數(shù)_布萊克曼窗 比較以上窗函數(shù)，可以看到，矩形窗函數(shù)具有最窄的主瓣B，但也有最大的旁

32、瓣峰值 A 和最慢的衰減速度 D。 漢寧窗主瓣稍寬，但有著較小的旁瓣和較大的衰減速度，因而被認為是較好的窗口。67窗口法：凱瑟窗函數(shù)凱瑟（Kaiser）窗上面討論的幾種窗函數(shù)以犧牲主瓣寬度，換取旁瓣抑制；Kaiser 窗全面反映了這種主瓣和旁瓣衰減之間的互換關系；定義了一組可調(diào)的由零階貝塞爾 Bessel 函數(shù)構成的窗函數(shù)；通過調(diào)整參數(shù)可以在主瓣寬度和旁瓣衰減之間自由選擇它們的比重。從而實現(xiàn)以同一種窗類型來滿足不同窗性能需求的目的。 Kaiser 窗函數(shù)由 J.F. Kaiser 提出，由下式給出： 其中 I0 是修正過的零階貝塞爾 Bessel 函數(shù)是用來調(diào)整窗形狀的參數(shù)，依賴于參數(shù) N。6

33、8對于相同的 N，Kaiser 窗可以提供不同的過渡帶寬，這是其他窗函數(shù)做不到的。通過調(diào)整參數(shù)，就可以方便地完成對過渡帶寬度和阻帶衰減的調(diào)整。參數(shù)越高，其頻譜的旁瓣越小，但主瓣寬度也相應增加。xI0(x)I0()10零階貝塞爾函數(shù) Kaiser 窗函數(shù)依參數(shù) 而變化 Matlab 中，函數(shù) w = kaiser(n, beta) 實現(xiàn) Kaiser 窗 窗口法：凱瑟窗函數(shù)弗里德里希威廉貝塞爾（Friedrich Wilhelm Bessel）1784年1846 年德國天文學家及數(shù)學家69下面圖固定，當窗的長度變化時，相應的旁瓣的高度保持不變。 窗口法：凱瑟窗函數(shù)70 = 5.658，則過渡帶寬

34、等于 7.8pi/N，最小阻帶衰減為 60dB，如下圖所示： 窗口法：凱瑟窗函數(shù)71凱瑟窗的計算由于 Bessel 函數(shù)的復雜性，這種窗的設計公式很難推導，為此，Kaiser 提出了經(jīng)驗公式。給定 p、s、Rp 和 As，參數(shù)定義如下： 對于過渡帶寬 = s - p (rad/s)，濾波器階數(shù)為需要強調(diào)的是：階數(shù)為 N 的濾波器大致能滿足要求，但最后的結果還必須要演算以便證明這一點。在 Matlab 中，函數(shù) w = kaiser(n, beta) 實現(xiàn) Kaiser 窗。5.3.3 窗口法：凱瑟窗函數(shù)72窗函數(shù)旁瓣峰值衰減（dB）窗函數(shù)主瓣寬度加窗后濾波器過渡帶寬（）加窗后濾波器阻帶最小衰減

35、（dB）矩形窗-134/N1.8/N-21漢寧窗(升余弦窗)-318/N6.2/N-44漢明窗(改進升余弦窗)-418/N6.6/N-53布萊克曼窗(二階升余弦窗)-5712/N11/N-74凱塞（=7.865）-5710/N10/N-80窗口法：常用窗函數(shù)的性能指標73FIR DF 窗口法設計步驟 性能要求 Hd(e j)把 Hd(e j) 展成傅里葉級數(shù)，得到 hd(n)；把 hd(n) 自然截短到所需的長度 N=2M+1；將截短后 hd(n) 右移 M 個取樣間隔，得 h(n)；將 h(n) 乘以合適的窗函數(shù)，即得所需的濾波器的沖激響應，這時窗函數(shù)以 n = M對稱（當然窗函數(shù)也可直接加

36、在 hd(n) 上，這時窗函數(shù)以原點為對稱）；利用 h(n)，既可用硬件構成濾波器的系統(tǒng)函數(shù) H(z)，也可直接用計算機軟件實現(xiàn)濾波。窗口法：設計步驟74數(shù)字低通濾波器的設計例： 一個理想低通數(shù)字濾波器的頻率響應如圖所示，為： 假定 c=0.25，分別取 N =11、21、31的線性相位 FIR，觀察加窗后對濾波器幅頻特性的影響。窗口法：設計步驟_ 數(shù)字低通75解: 由于 Hd(e j) 是一個實周期函數(shù)，把它展成為傅氏級數(shù)： 窗口法：設計步驟_ 數(shù)字低通式中將 hd(n) 截短為 N=2M+1，并將截短后的 hd(n) 移位，得：76窗口法：設計步驟_ 數(shù)字低通然后乘以窗函數(shù) wR(n)，得

37、到最后得 h(n)。對于 c = 0.25，由上式得：（當 c=，就得到一個全通濾波器）當 N=11 時，M= 5，求得h(0) = h(10) = -0.045，h(1) = h(9) = 0，h(2) = h(8) = 0.075，h(3) = h(7) = 0.1592，h(4) = h(7) = 0.2251，h(5) = 0.25。當 N=11 時，乘以漢明窗：77窗口法：設計步驟_ 數(shù)字低通78當 N 取不同值時，H(ej) 都不同程度上近似于Hd(ej)。N 過小時，通帶過窄,且阻帶內(nèi)波紋較大,當 N 增加時，通頻帶接近于0.25pi,阻帶內(nèi)波紋減小,但在通帶內(nèi)出現(xiàn)了波紋，隨著

38、N 的繼續(xù)增加，這些波紋并不能消失。由圖中還可以看到，使用漢明窗后，通帶內(nèi)的振蕩基本消失，阻帶內(nèi)的紋波也大大減小，從這一點上來說，濾波器的性能得到了改善，但是，這是以過渡帶的加寬為代價的。 矩形窗漢明窗窗口法：設計步驟_ 數(shù)字低通79適用范圍對于能用解析式表達，且傅里葉級數(shù)的系數(shù)容易求解的濾波器： 此時，窗函數(shù)法是設計 FIR DF 較為方便的一個方法如果 hd(n) 不易求，則使用該方法較為困難窗函數(shù)用窗口法設計 FIR DF，一個重要問題是選用何種窗函數(shù) w(n) 進行截短，以及截短的長度 N。窗函數(shù)的選擇：阻帶衰減指標滿足設計給定的阻帶衰減和其它濾波器性能要求；能量盡量集中于主瓣內(nèi)；個人

39、的經(jīng)驗及喜好有關。窗口法：設計步驟_ 數(shù)字低通80窗函數(shù)長度 N 的選擇：過渡帶寬指標采用試驗方法，即逐漸增大 M，直至 H(ej) 在通帶和阻帶內(nèi)部達到指標要求。若對 |Hd(ej)|的過渡帶提出了具體要求，因為 FIR DF 的過渡帶等于窗函數(shù)的主瓣寬度，那么通過查表計算 N ：例如 5.26，5.86。注意： （1）根據(jù)所要設計線性相位 FIR DF類型來決定最終 N 取奇數(shù)還是偶數(shù)。 （2）一般選擇 N 為奇數(shù)。 5.3.4 窗口法：設計步驟_ 數(shù)字低通81窗口法：設計步驟_ 數(shù)字低通頻率增益(幅度)所要求的通帶邊緣頻率p所要求的阻帶邊緣頻率 s設計中所用的截止頻率 c過渡帶寬度0.5

40、截止頻率 c 的確定：截止頻率 c 對應于明確的 0.5 增益點，而不再標志某個增益點。對于非理想濾波器，其截止頻率 c 不采用通帶邊緣頻率 p 或阻帶邊緣頻率 s， 而使用過渡帶的中點（即通帶邊緣和阻帶邊緣之間的中點。因此，窗函數(shù)法不能精確確定其通帶和阻帶的邊緣頻率）：82例： 根據(jù)下列指標設計一個線性相位FIR低通濾波器 通帶邊緣頻率 fp=2kHz 阻帶邊緣頻率 fstop=3kHz 阻帶衰減 40dB 取樣頻率 fs=10kHz 頻率|H(ejw)|所要求的通帶邊緣頻率p所要求的阻帶邊緣頻率 s設計中所用的截止頻率 c過渡帶寬度-6dB0dB-40dB2kHz3kHz解：(1) 求對應

41、的理想數(shù)字頻率： 過渡帶寬 = 3kHz 2kHz = 1kHz。 轉換為數(shù)字頻率過渡帶： 截止頻率： 數(shù)字截止頻率83(2) 設理想線性相位濾波器為： 由此可得脈沖響應： (4) 由過渡帶寬確定窗口長度： (3) 由阻帶衰減確定窗函數(shù)：因為阻帶衰減 40dB，通過查表知道，可以 選擇 Hanning 窗：(N-1)/2=15則此濾波器的脈沖響應為：8485例: 用矩形窗、漢寧窗和布萊克曼窗設計FIR低通濾波器，設N=11,c=0.2rad。 解: 用理想低通作為逼近濾波器，有86 用漢寧窗設計：用布萊克曼窗設計： 87低通幅度特性88 例：請選擇合適的窗函數(shù)及N來設計一個線性相位低通濾 波器

42、 要求其最小阻帶衰減為-45dB，過渡帶寬為 。 求出 。（設 ） 解：(1) 因為題目要求設計的低通濾波器的最小阻帶衰減 為-45dB，所以應該選擇漢明窗。89 又根據(jù)題目所給低通濾波器的表達式求得： 由此可得： 90例:一段樂曲中夾雜著高頻噪聲，嚴重影響收聽質(zhì)量。下圖給出了噪聲污染后的信號及其頻譜。系統(tǒng)的取樣頻率為 16kHz。設計一個濾波器來提高聲音質(zhì)量。 91從圖中頻譜可以看出，噪聲從 2kHz 的頻率點開始占據(jù)支配地位，因為這個頻率很低，消除噪聲的同時也會損失一部分音樂信息?？梢杂猛◣н吘夘l率為 2kHz 的低通濾波器對此段音樂進行濾波。由于沒有特殊的阻帶衰減要求，任何具有合理特性的

43、窗函數(shù)即可。在此選擇 Hamming 窗。本例中，對過渡帶寬也沒有特殊要求。為了得到合理的陡峭滾降，選擇 N = 101，則過渡帶寬為 92此低通濾波器的參數(shù)總結如下：取樣頻率：16kHz通帶邊緣頻率：2kHz過渡帶寬：545Hz截止頻率： 數(shù)字截止頻率：窗函數(shù)：Hamming 窗，即 93窗口長度：N=101。則此濾波器的脈沖響應為： 得到的脈沖響應、幅頻特性如下圖所示。 94020406080100-0.100.10.20.3h(n)012345678-100-80-60-40-200kHzMagnitudeMagnitude Response95伴有噪聲的音樂經(jīng)過上述濾波器濾除噪聲后的信

44、號及其頻譜如下圖所示。濾波后的信號幾乎沒有噪聲，但歌曲聽起來有點壓抑，這是因為歌曲中高頻分量也隨噪聲一起被濾除的緣故?？梢韵胂螅绻肼暣嬖谟谒械念l率分量上，則不可能在不嚴重降低信號質(zhì)量的情況下濾除噪聲。 96以上設計的是數(shù)字低通濾波器，若希望設計數(shù)字高通、帶通和帶阻濾波器，只需要改變付氏級數(shù)系數(shù)中積分的上、下限即可。數(shù)字高通濾波器窗口法：設計步驟_數(shù)字高通令其時域右移 M 位后的幅頻特性為：（頻域為 e-jM，表示時域右移 M 位） 97窗口法：設計步驟_數(shù)字高通則求得98窗口法：設計步驟_數(shù)字高通 從這個結果可以看出：一個高通濾波器相當于用一個全通濾波器（即c=）減去一個低通濾波器。 傳

45、輸函數(shù)： 脈沖響應：Hap(z)Hlp(z)X(z)Y(z)99數(shù)字帶通濾波器窗口法：設計步驟_數(shù)字帶通令其幅頻特性為：（頻域為 e-jM，表示時域右移 M 位） 則100窗口法：設計步驟_數(shù)字帶通求得 從這個結果可以看出：一個帶通濾波器相當于兩個低通濾波器相減，其中一個截止頻率為 h，另一個為 l。傳輸函數(shù)：脈沖函數(shù)：Hlph(z)Hlpl(z)X(z)Y(z)101窗口法：設計步驟_數(shù)字帶通 或者一個帶通濾波器相當于一個低通濾波器和一個高通濾波器相乘，即先經(jīng)過一個 LP DF，再經(jīng)過一個 HP DF。傳輸函數(shù)：脈沖函數(shù)：（頻域相乘，時域卷積） Hlp (z)Hhp(z)X(z)Y(z)10

46、2 例：用布拉克曼窗設計一個線性相位的理想帶通濾波器 求出 的表達式。 解:可求得此濾波器的時域函數(shù)為 103 采用布拉克曼窗設計時( N =51 )： 104 105 例： 利用窗函數(shù)法完成數(shù)字帶通濾波器的設計，并畫出其線性相位結構的示意圖。低端阻帶邊緣：低端通帶邊緣：高端通帶邊緣：高端阻帶邊緣： 窗函數(shù) 旁瓣峰值衰減（dB）過渡帶（w）阻帶最小衰減（dB）矩形窗 -13 4/N -21漢寧窗 -31 8/N -44漢明窗 -41 8/N -53布萊克曼窗 57 12/N -74解：106 107 因h(n)為偶對稱且N為奇數(shù)，其線性相位結構如圖表示。108 例：為取樣頻率為 22kHz 的

47、系統(tǒng)設計一個 FIR 帶通濾波器，中心頻率為 4kHz，通帶邊緣在 3.5kHz 和 4.5kHz，過渡帶寬為 500Hz，阻帶衰減 50dB。 解：截止頻率： kHzffhhc75.425.05.42/5.0=+=+=數(shù)字截止頻率： 脈沖響應： 109所以窗口長度為 ：窗函數(shù)：因為阻帶衰減 50dB，可以選擇 Hamming 窗，即 則此濾波器的脈沖響應為： 因為過渡帶寬500Hz，轉換為數(shù)字頻率為：110111數(shù)字帶阻濾波器窗口法：設計步驟_數(shù)字帶阻令其幅頻特性為：（頻域為 e-jM，表示時域右移 M 位） 則112窗口法：設計步驟_數(shù)字帶阻求得從這個結果可以看出：一個帶阻濾波器相當于一個

48、低通濾波器加上一個高通濾波器，低通濾波器的截止頻率為 l，高通在 h。傳輸函數(shù)：脈沖函數(shù)：Hlp(z)Hhp(z)+X(z)Y(z)113例：利用窗函數(shù)法設計滿足下列指標的最低階數(shù)的線性相位 FIR數(shù)字帶阻濾波器。該數(shù)字帶阻濾波器的性能指標如下： 低端阻帶邊緣： ；低端通帶邊緣： ； 高端阻帶邊緣： ；高端通帶邊緣： ；阻帶衰減不小于52dB。 窗口法：設計步驟解：由于，可選漢明窗，其過渡帶寬為 所以 取 N=67 (或 N=53)。114窗口法：設計步驟設理想線性相位帶阻濾波器的頻率響應為：其中 = (N-1)/2 =33 (或 = 26)115窗口法：設計步驟則 漢明窗 116FIR DF

49、 Matlab 設計函數(shù)b=fir1(n,wn,options)，單帶 FIR 濾波器b = fir2(n,f,m,options)，多帶 FIR 濾波器兩者可設計低通、高通、帶通、帶阻和通用多帶 FIR 濾波器Fir1 具有以下多種形式：b = fir1(n,wn)b = fir1(n,wn,ftype)b = fir1(n,wn,window)b = fir1(n,wn,ftype,window)b = fir1(.,noscale)參數(shù)向量 b 是 n 階 FIR 濾波器的系數(shù)截止頻率 wn 是從 0 到 1 之間的數(shù)，1 對應著奈氏頻率。對于高通濾波器，ftype 為 high；帶阻濾

50、波器 ftype 為 stop窗口法：Matlab 實現(xiàn)117對于帶通或帶阻濾波器，wn 為包含通頻帶邊帶頻率的一個二元素向量 wn1,wn2 。參量 window 表示所采用的窗函數(shù)類型。window 的長度必須為 n+1 （n 為濾波器的階數(shù)），若 window 卻省，則 fir1 使用漢明窗。注意因為奇數(shù)階的 II 型濾波器（h(n) 為偶對稱，長度N為偶數(shù)）在高頻段的頻率響應為零，所以 fir1 函數(shù)在高通和帶阻情況下不設計 II型濾波器，因此，如果 n 為奇數(shù)時，fir1 將階次加 1 并返回 I 型濾波器。函數(shù) b = fir2(n,f,m) 也可設計加窗 FIR 濾波器，但它針對

51、任意形狀的分段（piece-wise）線性頻率響應。向量 f 由從 0 到 1 的頻率點組成，其中 1 表示奈氏頻率，第一個點必須是 0，最后一個點必須是 1，頻率點必須是遞增的。m 是對應于頻率點 f 處的期望的頻率幅值響應。f 和 m 的長度必須相等。 窗口法：Matlab 實現(xiàn)118Matlab 頻率響應函數(shù)在 Matlab 中提供了一個 freqz 函數(shù)，利用這個函數(shù)開發(fā)一個新的函數(shù) freqz_m，它給出了絕對的和相對的 dB 值幅度響應、相位響應以及群延時響應。窗口法：Matlab 實現(xiàn)function db,mag,pha,grd,w = freqz_m(b,a);% Modif

52、ied version of freqz subroutine% -% db,mag,pha,grd,w = freqz_m(b,a);% db = Relative magnitude in dB computed over 0 to pi radians% mag = absolute magnitude computed over 0 to pi radians % pha = Phase response in radians over 0 to pi radians% grd = Group delay over 0 to pi radians% w = 501 frequency

53、samples between 0 to pi radians% b = numerator polynomial of H(z) (for FIR: b=h)% a = denominator polynomial of H(z) (for FIR: a=1)%H,w = freqz(b,a,1000,whole); H = (H(1:1:501); w = (w(1:1:501); mag = abs(H); db = 20*log10(mag+eps)/max(mag); pha = angle(H); grd = grpdelay(b,a,w);119例5.5 根據(jù)下列技術指標，設計一

54、個數(shù)字 FIR 低通濾波器 wp = 0.2, Rp = 0.25dB ws = 0.3，As = 50dB 采用漢明窗，確定脈沖響應，并給出所設計的濾波器的頻率響應圖。解: 在設計中，沒有使用通帶波動值 Rp=0.25dB，但必須檢查設計的實際波動，驗證它是否確實在給定容限內(nèi)。窗口法：Matlab 實現(xiàn)% Lowpass filter design - Hamming windowwp = 0.2*pi; ws = 0.3*pi;tr_width = ws - wpN = ceil(6.6*pi/tr_width) + 1 n=0:1:N;wc = (ws+wp)/2h = fir1(N,w

55、c/pi);db,mag,pha,grd,w = freqz_m(h,1);delta_w = 2*pi/1000;Rp = -(min(db(1:1:wp/delta_w+1) % Passband RippleAs = -round(max(db(ws/delta_w+1:1:501) % Min Stopband attenuation120運行結果如下： tr_width = 0.3142 （過渡帶寬） N = 67 （濾波器的階數(shù)，長度為 68，Type II 偶對稱偶數(shù) ） wc = 0.7854 （理想 LPF截止頻率） Rp = 0.0364 （實際通帶波動） As = 53

56、（最小阻帶衰減）窗口法：Matlab 實現(xiàn) 從結果看，67 階 Hamming 窗的 FIR 數(shù)字濾波器的實際阻帶衰減為 53dB，通帶波動為 0.0364dB，顯然滿足上面所提的技術要求，其時域和頻域響應曲線如下圖所示。121例5.6 利用例 5.5 給出的設計技術指標，選擇 Kaiser 窗，設計所需的低通濾波器。解：窗口法：Matlab 實現(xiàn)% Lowpass filter design - Kaiser windowwp = 0.2*pi; ws = 0.3*pi; As = 50;tr_width = ws - wp;N = ceil(As-7.95)/(14.36*tr_width

57、/(2*pi)+1) + 1n=0:1:N;beta = 0.1102*(As-8.7)wc = (ws+wp)/2;h = fir1(N,wc/pi,Kaiser(N+1,beta);db,mag,pha,grd,w = freqz_m(h,1);delta_w = 2*pi/1000;As = -round(max(db(ws/delta_w+1:1:501) % Min Stopband Attenuation122運行結果如下：N = 61 （濾波器階數(shù)）beta = 4.5513 As = 51 （實際的阻帶衰減）窗口法：Matlab 實現(xiàn)123例5.7 用 Blackman 窗設計

58、下面的數(shù)字帶通濾波器 阻帶：ws1=0.2，As1 = 60dB；ws2=0.8，As2 = 60dB。 通帶：wp1=0.35，Rp1 = 1dB；wp2=0.65，Rp2 = 1dB；解：窗口法：Matlab 實現(xiàn)% Bandpass filter design - Blackman windowws1 = 0.2*pi; wp1 = 0.35*pi;wp2 = 0.65*pi; ws2 = 0.8*pi;As = 60;tr_width = min(wp1-ws1),(ws2-wp2)N = ceil(11*pi/tr_width) + 1 %; N=75n=0:1:N;wc1 = (w

59、s1+wp1)/2; wc2 = (wp2+ws2)/2;wn = wc1/pi,wc2/pi;h = fir1(N,wn,blackman(N+1);db,mag,pha,grd,w = freqz_m(h,1);delta_w = 2*pi/1000;Rp = -min(db(wp1/delta_w+1:1:wp2/delta_w) % Actua; Passband RippleAs = -round(max(db(ws2/delta_w+1:1:501) % Min Stopband Attenuation124運行結果如下：N = 75 （濾波器的階數(shù)）Rp = 0.0028 （實際

60、的通帶波動）As = 75 （實際的阻帶衰減）窗口法：Matlab 實現(xiàn)125例5.8 理想帶阻濾波器的頻率響應為：窗口法：Matlab 實現(xiàn)用 Kaiser 窗設計一個長度為 N 的帶阻濾波器，它的阻衰減為 60dB。解：% Bandstop filter design - Kaiser windowN = 46; As = 60; n=0:1:N;beta = 0.1102*(As-8.7)+0.3wc1 = pi/3; wc2 = 2*pi/3;wn = wc1/pi,wc2/pi;h = fir1(N,wn,stop,kaiser(N+1,beta);db,mag,pha,grd,w