基于正交匹配追蹤的強脈沖電磁干擾濾波新方法

朱會杰， 王新晴，芮挺,趙洋， 李艷峰

(解放軍理工大學野戰(zhàn)工程學院，南京210007)

摘要：針對強脈沖電磁干擾能量大、規(guī)律性不強，對信號污染嚴重，提出了一種濾波新方法。依據(jù)信號和脈沖電磁干擾結構的不同，構造了單位脈沖原子匹配脈沖電磁干擾，選擇正弦原子、余弦原子和小波原子匹配有用信號。不同于常規(guī)稀疏分解濾波的思路，本文首先將脈沖電磁干擾看作信號成分，使用上述原子利用正交匹配追蹤進行分解。然后，僅選擇正弦原子、余弦原子和小波原子匹配的成分作為無噪信號，并將單位脈沖原子匹配的成分作為干擾濾除。仿真實驗和實例應用證明，該方法能夠有效抑制強脈沖電磁干擾和白噪聲，不僅信噪比高，而且保留了細節(jié)。

關鍵詞：正交匹配追蹤；濾波；脈沖電磁干擾；多原子

中圖分類號:TN911.4文獻標志碼：A

A new denoising method for strong pulse electromagnetic interference signals based on orthogonal matching pursuit

ZHUHui-jie,WANGXin-qing,RUITing,ZHAOYang,LIYan-feng(College of Field Engineering, PLA University of Science and Technology, Nanjing 210007, China)

Abstract:Because strong pulse electromagnetic interferences (PEMIs) have a high energy and a weak regularity, and they pollute signals seriously, a new denoising method to solve this problem was proposed. According to different structures of signals and PEMIs, a unit pulse atom was constructed to match PEMI, and sine atom, cosine atom and wavelet atom were chosen to match signal. Unlike common filtering ways of sparse decomposition, PEMI was firstly taken as a signal component, and all the atoms mentioned above were used to decompose the original signal. Secondly, only the components matched by sine atom, cosine atom and wavelet atom were used to reconstruct the unpolluted signal, and the component matched by the unit pulse atom was eliminated as PEMI. Simulations and applications testified that this technology can be used to filter strong PEMIs, the filtered signals not only have a high signal to noise ratio, but also retain details.

Key words:orthogonal matching pursuit; denoise; pulse electromagnetic interference (PEMI); multi-atom

在信號采集過程中，觀測信號難免會受到白噪聲的污染，而且還常常受到強電磁干擾，使信號淹沒在噪聲中，給信號的分析和處理帶來極大的不便[1]。盡管通過一些昂貴和復雜的硬件設備能夠對部分電磁干擾進行抑制[2]，但數(shù)字濾波的經(jīng)濟性和易用性具有更大的吸引力[3]。對于窄帶或者具有統(tǒng)計規(guī)律的電磁干擾，可以通過小波、高階統(tǒng)計、經(jīng)驗模態(tài)分解(Empirical Mode Decomposition，EMD)等方式進行濾波[4-6]。但對于隨機的強脈沖電磁干擾，瞬間能量很大，在頻域分布的范圍廣，現(xiàn)有方法還有些不足。另外，實測信號往往是多分量的，既有平穩(wěn)成分也有非平穩(wěn)成分，常規(guī)單基、單濾波器方法無法精確地提取有用信號[7]。

稀疏分解限制條件少、自適應性強、靈活高效[8]，為抑制強脈沖電磁干擾提供了一種新思路。常用的稀疏分解方法有匹配追蹤[9](Matching Pursuit，MP)、正交匹配追蹤[10](Orthogonal Matching Pursuit，OMP)和基追蹤[11](Basis Pursuit，BP)。OMP和MP都是貪婪迭代算法，通過每次挑選出與信號內積最大的原子，實現(xiàn)自適應選取原子。但OMP將選用的基進行了正交化處理，比MP有更高的精度，往往有更高的計算效率，得到了許多學者的關注[10, 12-13]，因此本文選用OMP作為稀疏分解的方法。

稀疏分解構建原子庫的方法比較靈活，針對不同的信號成分，可以選用多種原子的組合。由于有用信號中往往既有平穩(wěn)成分也有非平穩(wěn)成分，為了具有代表性，本文選用正弦原子和余弦原子匹配平穩(wěn)成分，小波原子匹配非平穩(wěn)信號成分，這種多原子的方式具有更強的自適應性[7]，提取信號更準確。常規(guī)稀疏分解除噪都是將分解后的殘余當做噪聲[9, 11, 14]，但這種方式不適用于強脈沖電磁干擾。為此，針對強脈沖電磁干擾結構特點，專門構造了單位脈沖原子。對信號進行稀疏分解時，將干擾當作信號成分進行提取，即選用正弦原子、余弦原子、小波原子和單位脈沖原子共同組成的原子庫對信號進行分解，但在重構信號時僅使用正弦原子、余弦原子和小波原子匹配的成分進行重構，將單位脈沖原子匹配的成分和信號分解的殘余剔除。這樣不僅能夠消除白噪聲，而且能有效抑制脈沖電磁干擾。

1基于OMP的強脈沖電磁干擾濾波算法原理

1.1強脈沖電磁干擾濾波的數(shù)學模型

假設觀測信號受到噪聲污染，其數(shù)學模型為：

y=s+σz

(1)

式中：y為觀測信號，是已知的；s為真實信號，是未知的；z為標準的高斯白噪聲，σ為噪聲的幅值，也是未知的。這是常規(guī)的噪聲污染模型。

當觀測信號還同時被強脈沖電磁干擾所污染時，其數(shù)學模型可以表示為：

y=s+p+σz

(2)

式中：p為脈沖電磁干擾，是未知的。因此對脈沖電磁干擾信號進行濾波相當于根據(jù)觀測信號y濾除干擾p和白噪聲σz，提取真實信號s。

1.2基于OMP的強脈沖電磁干擾濾波的算法過程

1.2.1過完備原子庫的構建

構造的原子與有用信號越相似，OMP分解的效果往往就越好[12-13]。強脈沖電磁干擾發(fā)生時間短，瞬時能量大，在結構上與有用信號有明顯的差別，根據(jù)其特點構造單位脈沖原子，假設信號的點數(shù)為N，第i個脈沖原子公式如下：

(3)

式中：1≤n≤N，且1≤i≤N。由公式(3)可知，脈沖原子i僅在第i點幅值為1，其余全為0，這與強脈沖電磁干擾瞬時大幅值的特性相似。

有用信號一般具有一定的規(guī)律，與強脈沖電磁干擾相比，有用信號變化相對平緩。為了具有通適性，選用正弦原子、余弦原子和小波原子匹配有用信號。正弦、余弦原子與周期性的平穩(wěn)信號相似，用于匹配平穩(wěn)成分。小波原子的適應性較強，用于匹配有一定突變的非平穩(wěn)成分，為了更明顯地區(qū)別于單位脈沖原子，僅選用低頻小波基，本文選用coiflets 3小波，分解層數(shù)為4，即僅選用了1～4層的低頻小波基。將正弦原子、余弦原子、小波原子所構成的原子庫用D1表示，使用D1重構的信號為真實信號；脈沖原子構成的原子庫用D2表示，D2所匹配的成分作為脈沖電磁干擾；它們共同構成的原子庫用D表示。

需要說明的是，這種方法不僅對脈沖電磁干擾有效，而且對于可以使用脈沖原子稀疏表示或發(fā)生時間很短的其他類型的電磁干擾也有效，比如發(fā)生時間較短的類方波干擾或三角波干擾等。另外，可以根據(jù)實際的有用信號的特點更換原子庫D1中的原子，可以得到更好的分解效果。

1.2.2OMP計算過程

OMP是一種將信號逐步分解的過程。首先，在原子庫中選取與信號最匹配的原子，然后將匹配的成分從原信號中減去得到殘余；在之后的迭代中，每次選擇與殘余成分最匹配的原子后，都首先將所有選中的原子正交化，之后重新在原始信號上進行投影，提取匹配部分并留下殘余；按照這種方式反復迭代，直到殘余信號的能量小于設定的閾值或達到最大迭代次數(shù)[10]。

OMP能保證殘余與所有已選原子正交，比MP收斂的更快，對原信號的表示也更稀疏[10]；相對于BP，OMP作為一種純貪婪迭代算法所需內存和計算時間都要小許多，對于長度為L，稀疏度為K的信號，K?L，使用OMP計算的復雜度為LK3，因此應用更廣泛[10, 12-13]。

(1)設y為待處理的信號序列，信號殘余R0=y，已選原子的集合ψ-1為空集，迭代次數(shù)n=0。

(4)

(5)

通過以上方式，可以有效地將白噪聲和脈沖電磁干擾濾除，為進一步對信號進行分析和處理奠定了基礎。

2仿真實驗與對比分析

為了驗證上述濾波方法的性能，構造以下仿真信號。原始信號s=s1+s2+s3，其中s1(t)=cos(20πt)，為正弦信號；s2(t)=[1+sin(40πt)]cos(200πt)，為調幅信號；s3(t)=5e-10tcos(60πt)，模擬非平穩(wěn)成分。此外，給原始信號添加白噪聲z和模擬脈沖電磁干擾p1、p2，其中p1為10個隨機發(fā)生的脈沖信號，幅值是[10,11]或者[-10,-11]之間的隨機數(shù)，模擬強脈沖電磁干擾；p2為4個寬度為5的強電磁干擾，幅值是[20,21]之間的隨機數(shù)，模擬類方波電磁干擾；噪聲z為信噪比(Signal to Noise Ratio，SNR)為4的白噪聲。信號采樣頻率為1 024 Hz，采樣點數(shù)為500。

圖1　原始信號與加噪信號的時域圖和頻譜圖 Fig.1 The waveforms in time domain and frequency domain of original signal and polluted signal

圖2　不同濾波方法去噪效果對比圖 Fig.2 The denoising results of different methods

原始信號與加噪信號的波形及頻譜圖如圖1所示。由圖中可以看出，加噪信號在時域和頻域都發(fā)生了畸變，需要將干擾成分濾除才能正確分析信號。為了充分驗證本文方法的有效性，采用1.2節(jié)所述信號重構方案，分別采用MP、OMP和BP對信號分解。此外還采用常規(guī)去噪方法，包括小波[4]、EMD[15-16]以及數(shù)學形態(tài)法[17-18]進行降噪，其中小波選用coiflets 3小波，分解層數(shù)為4，數(shù)學形態(tài)法中的結構元素為幅值為0的平結構元素，長度為6。這些方法的去噪結果如圖2所示，相對應的頻譜如圖3所示，并采用SNR和均方誤差(Mean Square Error，MSE)對這些方法進行評價，如表1所示。

由圖1～圖3和表1可知，得益于本文所提出的原子重構方案以及稀疏分解的自適應性，利用MP、OMP和BP三種方法恢復原始信號的效果較優(yōu)，并且OMP在時域、頻域、SNR和MSE評價指標中都是最優(yōu)的，表現(xiàn)出了優(yōu)良的特性?？陀^來說，BP也是一種優(yōu)異的稀疏分解方法[11]，但使用BP去噪時需要準確估計噪聲的大小，否則會影響重構效果；此外，BP分解時需要使用所有原子，所占內存和需要的計算時間都會大大增加，在這個仿真實驗中，BP所用時間是OMP所用時間的20倍，效率較低。由分解結果可以看出，BP在SNR和MSE值方面略優(yōu)于MP，但BP和MP都未能在頻域正確恢復出原始信號。而常規(guī)去噪方法小波和EMD則未能濾掉電磁干擾，這是由于它們都認為噪聲在一些轉換域中的能量會較小，將小于一定閾值的信號當做噪聲去除，顯然無法適應于大幅值、高能量的強脈沖電磁干擾信號。而數(shù)學形態(tài)法盡管能一定程度上濾除脈沖電磁干擾，但去噪后的信號發(fā)生了畸變，丟失了細節(jié)。此外，小波、EMD和數(shù)學形態(tài)法去噪后的信號在頻域都丟失了一些重要分量，并包含大量的噪聲，未能正確恢復原信號。

仿真實驗結果說明，常規(guī)去噪方法不適用于強脈沖電磁干擾信號，采用本文所提出的濾波方法能夠有效濾除脈沖類強電磁干擾和白噪聲，并且在時域和頻域都有較好的表現(xiàn)。

圖3　不同濾波方法去噪后的頻譜圖 Fig.3 Spectra of filtered signals by different denoising methods

去噪方法SNRMSEMP4.51.0OMP8.30.6BP5.50.8小波-8.44.2EMD-8.84.4數(shù)學形態(tài)法4.21.0

3在液壓泵壓力信號中的應用

圖4為液壓泵壓力測試平臺，實測壓力信號不僅受到外界的白噪聲干擾，還受到變頻器、電機和線路的強脈沖電磁干擾。測試時液壓泵轉速為900 r/min，液壓泵有7個柱塞，因此轉頻為15 Hz，壓力脈動基頻為105 Hz。

圖4　液壓泵壓力信號采集平臺 Fig.4 The test rig of pressure signal acquisition for hydraulic pump

原始壓力信號和去噪后信號的時域、頻域波形如圖5所示。由圖中可以看出，原始壓力信號已經(jīng)被脈沖電磁干擾和白噪聲嚴重污染，看不到油壓脈動，并且其頻譜圖中也包含很多的噪聲，使得頻譜分析可信度降低。采用本文方法進行濾波后，不僅濾除了強脈沖電磁干擾，精確地恢復了壓力變化過程，還保留了壓力脈動的細節(jié)和突變。此外，去噪后的信號頻譜圖也沒有了噪聲，而且顯示更多的轉頻和壓力脈動的倍頻。這些結果表明，本文所提出的去噪方法能夠有效濾除強脈沖電磁干擾和白噪聲，同時能完整保留有用信號的細節(jié)部分，滿足了實際需要。

圖5　原始壓力信號與去噪結果對比 Fig.5 The comparison between original pressure signal and denoised signal

4結論

針對強脈沖電磁干擾的結構特征構造的單位脈沖原子能夠很好地匹配脈沖干擾，借助于本文設計的信號重構方案，有效地消除了強脈沖電磁干擾和白噪聲，突破了常規(guī)方法的思維定勢，為濾波提供了一個新的思路。仿真實驗對比和實例應用說明，常規(guī)去噪方法抑制強脈沖電磁干擾尚有不足，而本文所提出的基于OMP的多原子的濾波方法能夠對脈沖電磁干擾以及短時強電磁干擾都有良好的抑制作用，同時很好地保留了信號的細節(jié)，克服了常規(guī)去噪算法的不足。當應用本文濾波方法時，可以根據(jù)實際信號的特點改變原子庫D1的組成，以得到更好的去噪效果。

參考文獻

[1]Liang J, Wang Y, Huang Y, et al. Electromagnetic interference shielding of graphene/epoxy composites[J]. Carbon, 2009, 47 (3): 922-925.

[2]Al-Saleh M H, Sundararaj U. Electromagnetic interference shielding mechanisms of CNT/polymer composites[J]. Carbon, 2009, 47 (7): 1738-1746.

[3]Tarateeraseth V, See K Y, Canavero F G, et al. Systematic electromagnetic interference filter design based on information from in-circuit impedance measurements[J]. Electromagnetic Compatibility, IEEE Transactions on, 2010, 52(3): 588-598.

[4]Aggarwal R, Rathore S,Singh J K, et al. Noise reduction of speech signal using wavelet transform with modified universal threshold[J]. International Journal of Computer Applications, 2011, 20(5):14-19.

[5]Inoue T, Saruwatari H, Takahashi Y, et al. Theoretical analysis of musical noise in generalized spectral subtraction based on higher order statistics[J]. Audio, Speech, and Language Processing, IEEE Transactions on, 2011, 19 (6): 1770-1779.

[6]Wu Z, Huang N E. Ensemble empirical mode decomposition: a noise-assisted data analysis method[J]. Advances in adaptive data analysis, 2009,1(1): 1-41.

[7]Feng Z, Chu F. Application of atomic decomposition to gear damage detection[J]. Journal of Sound and Vibration, 2007, 302 (1): 138-151.

[8]Olshausen B A. Emergence of simple-cell receptive field properties by learning a sparse code for natural images[J]. Nature, 1996, 381 (6583): 607-609.

[9]Mallat S G, Zhang Z. Matching pursuits with time-frequency dictionaries[J]. Signal Processing, IEEE Transactions on, 1993, 41(12): 3397-3415.

[10]Tropp J A, Gilbert A C. Signal recovery from random measurements via orthogonal matching pursuit[J]. Information Theory, IEEE Transactions on, 2007, 53(12): 4655-4666.

[11]Chen S S, Donoho D L, Saunders M A. Atomic decomposition by basis pursuit[J]. SIAM Journal on Scientific Computing, 1998, 20(1): 33-61.

[12]Donoho D L, Tsaig Y, Drori I, et al. Sparse solution of underdetermined systems of linear equations by stagewise orthogonal matching pursuit[J]. Information Theory, IEEE Transactions on, 2012, 58(2): 1094-1121.

[13]Needell D, Vershynin R. Uniform uncertainty principle and signal recovery via regularized orthogonal matching pursuit[J]. Foundations of Computational Mathematics, 2007, 9(3): 317-334.

[14]Wright J, Yang A Y, Ganesh A, et al. Robust face recognition via sparse representation[J]. Pattern Analysis and Machine Intelligence, IEEE Transactions on, 2009, 31(2): 210-227.

[15]Flandrin P, Rilling G, Goncalves P. Empirical mode decomposition as a filter bank[J]. Signal Processing Letters, IEEE, 2004, 11(2): 112-114.

[16]Huang N E, Wu M L, Long S R, et al. A confidence limit for the empirical mode decomposition and Hilbert spectral analysis[J]. Proceedings of the Royal Society of London. Series A: Mathematical, Physical and Engineering Sciences, 2003, 459(2037): 2317-2345.

[17]沈長青,謝偉達,朱忠奎, 等. 基于 EEMD 和改進的形態(tài)濾波方法的軸承故障診斷研究[J]. 振動與沖擊, 2013, 32 (2): 39-43.

SHEN Chang-qing, XIE Wei-da,ZHU Zhong-kui, et al.Rolling element bearing fault diagnosis based on EEMD and improved morphological filtering method[J].Journal of Vibration and Shock,2014,33(4):11-16.

[18]鐘先友,趙春華,陳保家,等. 基于形態(tài)自相關和時頻切片分析的軸承故障診斷方法[J]. 振動與沖擊, 2014,33(4):11-16.

ZHONG Xian-you, ZHAO Chun-hua, CHEN Bao-jia, et al.Bearing fault diagnosis method based on morphological filtering,time-delayed autocorrelation and time-frequency slice analysis[J].Journal of Vibration and Shock,2014,33(4):11-16.