窗函數(shù)的奇數(shù)階約束問題探討

歐巧鳳　熊邦書

摘要：本文明確了經典FIR濾波器的窗函數(shù)法中窗函數(shù)的奇數(shù)階約束問題：在定采樣率數(shù)字信號處理系統(tǒng)中，窗函數(shù)法設計的線性相位FIR數(shù)字濾波器必然為第一類線性相位數(shù)字濾波器，階數(shù)N必須為奇數(shù)，并闡明了該結論適用于高維濾波器設計。

關鍵詞：數(shù)字濾波器；窗函數(shù)；奇數(shù)階約束；線性相位

中圖分類號：TN911.72 文獻標志碼：A 文章編號：1674-9324（2018）34-0200-02

引言：通過多年的教學研究和多套教材的對比研究發(fā)現(xiàn)，窗函數(shù)法設計線性相位FIR數(shù)字濾波器理論中，對于濾波器階數(shù)N，即窗函數(shù)的窗口長度N，國內外的數(shù)字信號處理教材均給出了N的計算公式N的計算公式N≥ 。由N≥ 取整得到的N可能為奇數(shù)，也可能為偶數(shù)，例如N≥65.8時，N究竟取66還是67，國內外的相關教材中，并沒有說明N能否取偶數(shù)。部分習題參考書中有這樣的說明：“為了方便，N取奇數(shù)”。部分參考書中對分別取奇數(shù)和偶數(shù)進行了計算和解答。

1.窗函數(shù)濾波器的奇數(shù)階約束。窗函數(shù)法也稱傅氏級數(shù)法，其設計線性相位FIR數(shù)字濾波器的思想是利用時域截短序列h（n）近似理想濾波器的無限長脈沖響應h （n），從而在頻域H（e ）近似H （e ）。所設計的線性相位FIR濾波器滿足h（n）=h（N-n-1），當N為奇數(shù)時，可以實現(xiàn)低通、高通、帶通和帶阻濾波器；N為偶數(shù)時，不可以實現(xiàn)高通和帶阻濾波器。因此，高通和帶阻的線性相位FIR濾波器已經存在奇數(shù)階約束。下面證明窗函數(shù)法設計的低通和帶通線性相位FIR濾波器也存在奇數(shù)階約束，即窗函數(shù)存在奇數(shù)階約束。理想低通濾波器LPF幅度特性和相位特性公式表達為：H （e ）=1，0，0≤ω？搖≤ωω <ω？搖≤π θ（ω）=-jω·α。對應的特性曲線如圖1（a）所示。理想帶通濾波器BPF幅度特性和相位特性公式表達為：H （e ）=1，0，ω ≤ω？搖≤ω 其它，θ（ω）=-jω·α。對應的特性曲線如圖1（b），圖中帶通濾波器等效于截止頻率為ω 的低通濾波器函數(shù)減去截止頻率為ω 的低通濾波器函數(shù)。

低通濾波器頻域響應函數(shù)統(tǒng)一表示為：H （e ）=e ，00≤ω？搖≤ω ω <ω？搖≤π，其中α為低通特性的群延時。對H （e ）進行傅立葉逆變換，可求得理想低通濾波器的時域響應序列為：h （n）= H （e ）e dn= e e dw= · （1）

理想帶通濾波器的時域響應序列為：

由式（1）和式（2）可知，低通和帶通的h （n）均具有偶對稱特性，且對稱中心是序列最高點對應的下標n=α。理想低通濾波器h （n）如圖2所示，h （n）是無限長非因果序列，物理上不可實現(xiàn)。實際中，忽略一定階次以上的旁瓣能量，截取有限長序列h（n）作為h （n）的近似逼近，方法是取窗函數(shù)w（n）與h （n）相乘h（n）=h （n）·w（n），物理上表示對h （n）進行截取，并根據(jù)窗函數(shù)的系數(shù)進行旁瓣衰減。幾種常用的窗函數(shù)的時域曲線如圖3所示。

圖2中的α為延時單元的個數(shù)。數(shù)字信號中根據(jù)采樣間隔來確定時間t=n*Δt= ，在定采樣率信號處理系統(tǒng)中，小數(shù)個延時單元沒有意義，因此，數(shù)字序列下標α必須為整數(shù)。

h （n）對稱中心上有一個最大序列值點h（α）。待設計的FIR濾波器需要線性相位，因此h（n）是偶對稱或奇對稱的序列。由圖3可知，常規(guī)窗函數(shù)是有限長序列，具有偶對稱性：w（n）=w（N-n-1）。根據(jù)序列運算規(guī)律可知，濾波器函數(shù)h（n）=h （n）·w（n）依然是偶對稱序列，h（n）對應的數(shù)字濾波器為第一類線性相位。根據(jù)序列相乘須對應點相乘，h（α）=h （α）·w（α），w（α）是序列中的一個點，且必須是序列對稱中心，所以α必須為整數(shù)，h（n）總長度N=2α+1，必然為奇數(shù)。反之，假設窗函數(shù)階數(shù)N為偶數(shù)，要保證線性相位，則h（n）必須偶對稱，N=6的h（n）如圖4（a）所示，其中延時單元的個數(shù)α= =2.5為小數(shù)，則h （n）=S （n-α）=S （n-2-0.5），對于定采樣率數(shù)字信號處理，小數(shù)下標無意義，α=2.5，N=6。對于變采樣率數(shù)字信號處理，在N之間插值N-1個點，在2倍重采樣條件下h（n）長度為N+（N-1）=11，仍然存在奇數(shù)階約束。如果h（n）非偶對稱，h（n）如圖4（b）所示，此時設計的FIR濾波器不是線性相位，不屬于設計線性相位FIR數(shù)字濾波器設計范疇。

2.窗函數(shù)法設計FIR濾波器理論中的奇數(shù)階約束引申到二維圖像濾波，同樣適用。對二維圖像數(shù)據(jù)I（x，y）進行空間濾波，數(shù)字圖像處理教材中很早就會講到均值濾波、高斯濾波等算法，所采用的算子模板尺寸有3×3和5×5等，模板尺寸為奇數(shù)的原因，教材中一般解釋為考慮到關于當前像素點對稱，本質原因是定采樣率濾波器存在奇數(shù)階約束，下面以二維高斯濾波器為例進行說明。高斯函數(shù)的傅立葉變換仍是高斯函數(shù)，因此，頻域的二維高斯低通濾波器w （u，v），其時域響應函數(shù)也為高斯函數(shù)：

H（x，y）為無限長的二維函數(shù)，通過尺寸為N×N的二維矩形窗函數(shù)ω （x，y）進行截斷，則得到N階二維高斯數(shù)字濾波器。其中取平滑系數(shù)σ=2、N=5的二維矩形窗函數(shù)進行截斷，可得到平滑系數(shù)為2的5×5高斯算子模板H（x，y），如圖5所示。采用二維時域卷積，可得到高斯濾波后的圖像F（x，y）=I（x，y）？鄢H（x，y）。

二維高斯函數(shù)具有可分離性：H（x，y）=e =

e ·e =H（x）·H（y），因此對輸入圖像I（x，y）先進行x方向的一維高斯濾波，再進行y方向的一維高斯濾波，可達到二維高斯濾波的相同效果。一維加窗高斯濾波器存在奇數(shù)階約束，因此，二維窗函數(shù)在兩個維度上存在奇數(shù)階約束。

結論：從上述推導可以看出，采用窗函數(shù)法設計定采樣率FIR數(shù)字濾波器，濾波器階數(shù)N必須為奇數(shù)，且設計出來的濾波器必然為第一類線性相位數(shù)字濾波器。

參考文獻：

[1]門愛東.數(shù)字信號處理系統(tǒng)分析與設計[M].楊波，譯.北京：電子工業(yè)出版社，2004：147-183.（原著：Paulo S.R.Diniz Eduardo A.B.da Silva.Digital Signal Processing System Analysis and Design.ISBN 7-121-00063-6）

[2]高西全，丁玉美.數(shù)字信號處理[M].西安：西安電子科技大學出版社，2008：196-213.