王浩軍

（中國空空導(dǎo)彈研究院洛陽471009）

基于窗函數(shù)的數(shù)字FIR濾波器設(shè)計

王浩軍

（中國空空導(dǎo)彈研究院洛陽471009）

簡要介紹了窗函數(shù)設(shè)計FIR濾波器的理論基礎(chǔ)，運用Matlab語言實現(xiàn)了低通FIR濾波器的設(shè)計與仿真。窗函數(shù)法設(shè)計FIR濾波器具有設(shè)計簡單、方便等優(yōu)點。

窗函數(shù)；FIR濾波器；海明窗；布萊克曼窗

Class NumberTN713.7

1 引言

數(shù)字濾波器用于改變信號在時域或頻域中的信號特性，提取出含在原信號內(nèi)部有用信號供下一級使用。

數(shù)字濾波器是完成信號濾波處理功能的、用有限精度算法實現(xiàn)的離散時間線性非時變系統(tǒng)，是數(shù)字信號處理中使用最為廣泛的一種線性系統(tǒng)，本文主要研究數(shù)字濾波器中的窗函數(shù)設(shè)計法。

2 數(shù)字濾波器類型

根據(jù)單位脈沖響應(yīng)的不同，數(shù)字濾波器可以分為兩大類FIR和IIR，各自具有優(yōu)缺點，F(xiàn)IR的優(yōu)點：1）具有線性相位。由于系數(shù)的對稱性，F(xiàn)IR的頻率響應(yīng)曲線的相位是線性的，在通帶內(nèi)的線性相位只是增加了通帶內(nèi)信號在時域上的延遲，這非常有利于圖像或者通信上對信號處理的要求。2）由于是有限個抽頭，F(xiàn)IR的輸出只是和當前以及過去有限個輸入有關(guān)，傳遞函數(shù)只有零點，沒有極點，所有FIR濾波器永遠是穩(wěn)定的，不會因為系數(shù)精度的影響導(dǎo)致沖擊響應(yīng)發(fā)散。FIR的缺點：由于沖擊響應(yīng)個數(shù)有限，所以對于同樣濾波器的要求，F(xiàn)IR實現(xiàn)的階數(shù)要比IIR大得多。IIR的優(yōu)點在于它的反饋結(jié)構(gòu)，只需要很少的階數(shù)便可以實現(xiàn)在阻帶內(nèi)的幅度衰減，但是IIR缺點也在于它的反饋結(jié)構(gòu)，系數(shù)精度的影響有可能導(dǎo)致濾波器極點移向不穩(wěn)定區(qū)域，導(dǎo)致響應(yīng)發(fā)散。另外，IIR濾波器相位響應(yīng)不是線性的，這限制了它在一些嚴格要求線性相位地方的應(yīng)用。

3 FIR濾波器設(shè)計

FIR濾波器設(shè)計方法的根本思想是逼近所需求的頻率響應(yīng)曲線，目前采用比較多的是窗函數(shù)法和最佳逼近方法。窗函數(shù)法的基本原理是從所要求的頻率響應(yīng)Hd(ejω)經(jīng)過反傅里葉變換導(dǎo)出單位脈沖響應(yīng)hd[n]：

加窗處理主要是因為hd[n]是無限長序列，需要對其進行截斷，得到滿足要求的有限長的單位脈沖響應(yīng)h[n]=w[n]hd[n]。h[n]就是一個逼近理想頻率響應(yīng)的因果FIR濾波器的單位單位脈沖響應(yīng)。

加窗的方式（包括窗形狀和窗長N）決定了濾波器的頻率響應(yīng)對理想頻率響應(yīng)的逼近程度。常用的窗函數(shù)基本參數(shù)如表1所示。

表1 典型窗函數(shù)性功能比較

矩形窗是不隨時間變量的改變而改變的窗函數(shù)，優(yōu)點是主瓣相對居中，缺點是旁瓣比較高，并且還有負旁瓣，所以在變換中容易帶進高頻率的干擾和泄漏，可能還會有負譜現(xiàn)象出現(xiàn)。它的主瓣寬度最小，頻率分辨率最高，但幅頻識別最低。

漢寧窗屬于余弦窗，主瓣寬度增加并且高度降低，旁瓣則是減小的程度明顯，如果要考慮減小泄漏的因素，那么漢寧窗比矩形窗更加適合，但是漢寧窗主瓣寬度增加，相當于分析帶寬增加，所以漢寧窗的頻率的分辨力會下降。

海明窗也是余弦窗，也被稱為改進的升余弦窗。海明窗與漢寧窗的加權(quán)系數(shù)不同，由于它的加權(quán)系數(shù)，旁瓣達到最小。分析表明，海明窗的第一瓣衰減為-42dB，與漢寧窗相比，旁瓣達到了更小，并且衰減速度也相對于比較慢。

窗函數(shù)的選擇原則：

1）主辧寬度應(yīng)盡可能窄，提高頻域分辨率和減少泄露，以獲得較大的阻帶衰減。

2）具有較低的旁辧幅度，尤其是第一旁辧幅度。旁辧幅度隨頻率盡快衰減。

如這兩條不能同時滿足，往往是增加主瓣寬度以換取對旁瓣的抑制。

利用Matlab仿真，w1=bartlett（64）；

w2=hamming（64）；

wvtool（w1，w2）；

窗口長度64的三角窗和海明窗，在同一張圖中看時域和頻域，如圖1所示。其中，藍色是三角窗，綠色是海明窗。

對于相同窗口長度，布萊克曼窗相比海明窗和漢寧窗主瓣更寬，邊頻泄露小。

4 結(jié)語

不同的窗函數(shù)之間，它們的特點不同，優(yōu)缺點也大不一樣。并且，不同的窗函數(shù)對濾波器的濾波信號影響也不一樣。所以在設(shè)計濾波器的過程中，選擇窗函數(shù)是至關(guān)重要的一步。如果只是要求讀出來主瓣頻率準確值，而不用考慮幅值精密度，那么就可以選擇矩形窗，最好是主瓣寬度稍微窄些，并且矩形窗也便于分辨；如果是分析頻帶較窄的信號，并且干擾噪音較強，那么所選用的窗函數(shù)的旁瓣幅度最好較小，例如漢寧窗、三角窗等窗函數(shù)。通過選擇不同的窗函數(shù)并且對它們的影響進行控制，那么就可以更好地降低設(shè)計要求，同時可以提高效率。

［1］賀超英，王少渝.MATLAB應(yīng)用與試驗教程［M］.北京：電子工業(yè)出版社，2013.

［2］劉興釗，李力利.數(shù)字信號處理［M］.北京：電子工業(yè)出版社，2010.

［3］李在林，陳坤.基于MATLAB的數(shù)字濾波器設(shè)計方法的分析［J］.信息通信，2015（5）：17-18.

［4］徐璐.基于MATLAB的FIR和IIR數(shù)字濾波器的設(shè)計研究［J］.數(shù)字技術(shù)與應(yīng)用，2016（10）：165.

［5］董胡.基于窗函數(shù)與MATLAB的數(shù)字FIR濾波器設(shè)計［J］.微型電腦應(yīng)用，2016（3）：30-32.

Digital FIR Filter Design Based on Window Function

WANG Haojun（China Air-to-Air Missile Academy，Luoyang471009）

The article simply introduces the design of digital FIR filter with window function，and realizes the design and emu?lation of low-pass filter by means of Matlab language.It is simple，convenient and useful to design FIR filter with window function.

window function，F(xiàn)IR Filter，Hamming，Blackman

TN713.7

10.3969/j.issn.1672-9730.2017.08.041

2017年2月5日，

2017年3月24日

王浩軍，男，碩士，工程師，研究方向：雷達接收機設(shè)計，微波設(shè)計。