梁瑞剛　姚熾偉

摘 要：針對(duì)數(shù)字濾波器設(shè)計(jì)問(wèn)題，利用諧波小波在頻域具有良好的盒形特性，以及傅里葉變換的尺度、線性、頻移特性，研究了一種通過(guò)疊加多個(gè)具有不同中心頻率的諧波小波，在頻域構(gòu)建出平頂濾波器的方法；結(jié)果表明：基于諧波小波設(shè)計(jì)的梳狀濾波器設(shè)計(jì)方法易于理解和編程實(shí)現(xiàn)，通過(guò)合理選擇帶寬參數(shù)即可設(shè)計(jì)出工程應(yīng)用所需的高性能多通帶濾波器。

關(guān)鍵詞：信號(hào)處理；梳狀濾波；濾波器設(shè)計(jì)；諧波小波

濾波是信號(hào)分析中較為常用的手段之一，它的目的主要在于信號(hào)選頻，將所需要的頻率選取出來(lái)，而將不需要的頻率成分衰減掉。當(dāng)信號(hào)在復(fù)雜系統(tǒng)中傳輸時(shí)，每通過(guò)其中的一個(gè)環(huán)節(jié)，都會(huì)受到該環(huán)節(jié)傳輸特性的影響，使信號(hào)有所變化（衰減、放大、延遲等），這就形成了更為廣泛的濾波和濾波器的概念。

按照信號(hào)處理的性質(zhì)，濾波器可分為模擬濾波器和數(shù)字濾波器。數(shù)字濾波器的實(shí)現(xiàn)不但比模擬濾波器容易的多，而且還能獲得較理想的濾波器性能，在數(shù)字信號(hào)處理中應(yīng)用非常廣泛。在經(jīng)典數(shù)字濾波器設(shè)計(jì)中包括無(wú)限沖激響應(yīng)濾波器（IIR，Infinite Impulse Response）設(shè)計(jì)和有限沖激響應(yīng)濾波器（FIR，F(xiàn)inite Impulse Response）設(shè)計(jì)兩大類。FIR濾波器可實(shí)現(xiàn)線性相位，為獲得較好的性能，常需要較高的階次，若對(duì)信號(hào)處理的實(shí)時(shí)性沒(méi)有特殊要求，F(xiàn)IR濾波器是較好的選擇。

經(jīng)典濾波器的設(shè)計(jì)原理和方法都很成熟。筆者利用諧波小波在時(shí)域、頻域都具有良好集中性的特點(diǎn)，以及傅里葉變換的尺度、線性、頻移特性，研究了一種通過(guò)多個(gè)具有不同中心頻率的諧波小波時(shí)域疊加，在頻域構(gòu)建出平頂濾波器的設(shè)計(jì)方法。利用MATLAB軟件編寫了應(yīng)用程序，并結(jié)合仿真信號(hào)進(jìn)行了驗(yàn)證。

1 諧波小波

1993年由英國(guó)劍橋大學(xué)D.E.Newland教授首先提出的諧波小波概念，這種小波有優(yōu)秀的緊支性以及完全"盒形"的頻域特性，具有更為廣泛意義的正交性，可以更為靈活的實(shí)現(xiàn)時(shí)頻分解，沒(méi)有二進(jìn)限制。

2 梳狀濾波

廣義諧波小波實(shí)質(zhì)上等價(jià)于一個(gè)可以任意調(diào)節(jié)通帶位置的理想帶通濾波器，在實(shí)際應(yīng)用中，由于存在時(shí)域截?cái)?，所以其頻譜會(huì)有嚴(yán)重的Gibbs現(xiàn)象，即紋波現(xiàn)象。為了壓制紋波，常用的做法是疊加邊瓣較小的窗函數(shù)，比如漢寧窗、海明窗或高斯窗。其中，漢寧窗和海明窗由于無(wú)法調(diào)節(jié)衰減頻率，因此在實(shí)際應(yīng)用中較為受限，而高斯窗可以通過(guò)調(diào)節(jié)帶寬參數(shù)來(lái)調(diào)節(jié)衰減頻率從而可以適用于各種場(chǎng)合。引入高斯因子而構(gòu)造的復(fù)解析帶通濾波器的時(shí)域表達(dá)式為：

其中，a為高斯窗的帶寬參數(shù)。通過(guò)設(shè)置不同的帶寬參數(shù)a就可以達(dá)到調(diào)節(jié)衰減速率的目的，使得其頻率響應(yīng)具有近似“盒形”特性，因而可以適用不同的應(yīng)用場(chǎng)合。

考慮到故障頻率往往是以倍頻的形式出現(xiàn)，因此將不同中心頻率的帶通濾波器相互疊加構(gòu)成—梳狀濾波器就可以更好的完成濾波功能，組合帶通濾波器的時(shí)域表達(dá)式為：

實(shí)際應(yīng)用中通過(guò)合理設(shè)置不同的參數(shù)就可以實(shí)現(xiàn)不同帶寬及通帶間隔的梳狀濾波器設(shè)計(jì)，達(dá)到多通帶濾波的功能。

3 信號(hào)仿真

利用基于諧波小波設(shè)計(jì)的梳狀濾波器來(lái)進(jìn)行降噪，令高斯窗的帶寬參數(shù)a=2，構(gòu)造中心通帶的中心頻率為400Hz，中心通帶兩側(cè)通帶數(shù)為2，相鄰?fù)◣чg隔為30 Hz，各通帶帶寬為10 Hz，對(duì)信號(hào)進(jìn)行卷積處理后，頻譜如圖2所示。

從圖2中可以看出干擾頻帶已經(jīng)基本濾除，信噪比由之前的4.12dB變?yōu)榱?5.8dB，在處理實(shí)際采集故障信號(hào)時(shí)，由于采集誤差的存在，邊頻帶不可能等間隔分布，因此，梳狀濾波器的通帶帶寬不能太小以確?？梢詫⒐收闲盘?hào)提取出來(lái)。

4 結(jié)束語(yǔ)

基于諧波小波原理的梳狀濾波器設(shè)計(jì)，利用了諧波小波近乎“盒形”的頻域特性，考慮到Gibbs效應(yīng)，在疊加高斯因子后其紋波衰減速率得到調(diào)節(jié)，通過(guò)疊加不同中心頻率的諧波小波就可以設(shè)計(jì)出所要求的梳狀濾波器，并通過(guò)帶寬參數(shù)a可以靈活控制中心帶寬位置和通帶間隔，因而可以利用此方法來(lái)靈活設(shè)計(jì)各種梳狀濾波器來(lái)處理實(shí)際工程應(yīng)用問(wèn)題，此設(shè)計(jì)原理流程易于理解，理論支撐合理，同時(shí)由于諧波小波有明確的函數(shù)表達(dá)式，故易于通過(guò)MATLAB編程實(shí)現(xiàn)。

此方法實(shí)際應(yīng)用中需要注意以下幾個(gè)問(wèn)題：

（1）為了保證濾波器的過(guò)渡帶窄，帶寬參數(shù)a應(yīng)選擇大一些。但此時(shí)為了保證通帶平直，所使用的小波個(gè)數(shù)會(huì)多一些。因此，a和？駐f的取值不唯一，可以有多種不同的組合。

（2）為了簡(jiǎn)化設(shè)計(jì)過(guò)程，。這樣濾波器的通帶完全由起始頻率fL和組合小波的個(gè)數(shù)M確定。

（3）經(jīng)典濾波器中的低通、高通、帶通濾波器均可按單通帶濾波器設(shè)計(jì)，帶阻濾波器可按兩通道帶通濾波器設(shè)計(jì)。

（4）在多通帶濾波器設(shè)計(jì)中優(yōu)勢(shì)明顯，且每個(gè)通帶的起始頻率和帶寬可以靈活選擇，尤其對(duì)于窄帶梳狀濾波器的設(shè)計(jì)可以取得很好的性能。

參考文獻(xiàn)

[1]胡廣書.數(shù)字信號(hào)處理——理論、算法與實(shí)現(xiàn)[M].北京：清華大學(xué)出版社，2003.

[2]胡廣書.現(xiàn)代信號(hào)處理教程[M].北京：清華大學(xué)出版社，2004.

[3]田福慶，羅榮，李萬(wàn)，等.改進(jìn)的諧波小波包峭度圖及其應(yīng)用[J].上海交通大學(xué)學(xué)報(bào)，2014，48（1）：39-42.

[4]徐文明，張梅軍，唐建，等.用Morlet小波進(jìn)行包絡(luò)檢波分析[J].內(nèi)燃機(jī)工程，2002，23（2）：81-84.

[5]唐貴基，鄧飛躍.基于改進(jìn)諧波小波包分解的滾動(dòng)軸承復(fù)合故障特征分離方法[J].儀器儀表學(xué)報(bào)，2015，36（1）：143-149.

[6]張景亭，孫勇軍，周友明，等.Morlet組合小波在顫振試飛中的應(yīng)用[J].飛行力學(xué)，2008，26（4）：62-64.

[7]李遠(yuǎn)軍，許鋼，林園勝，等.基于小波變換的電力系統(tǒng)諧波檢測(cè)的研究[J].重慶工商大學(xué)學(xué)報(bào)（自然科學(xué)版），2014，31（6）：78-81.

作者簡(jiǎn)介：梁瑞剛（1991-），男，碩士研究生，從事機(jī)械振動(dòng)信號(hào)診斷與特征提取研究。