章耀文， 陳少昌

(海軍工程大學(xué) 電子工程學(xué)院，湖北 武漢 430000)

基于壓縮感知的虛擬暗室測試方法

章耀文， 陳少昌

(海軍工程大學(xué) 電子工程學(xué)院，湖北 武漢 430000)

摘要:為了在現(xiàn)場電磁發(fā)射測試中有效濾除背景噪聲，針對存在相關(guān)信號的環(huán)境中多重信號分類(MUSIC)算法失效的問題，提出基于壓縮感知(CS)的虛擬暗室測試方法。通過CS理論獲得受試設(shè)備的輻射信號和環(huán)境中其他干擾信號的來波方向，并采取最小方差無畸變波束形成算法進(jìn)行波束形成，從而對干擾信號實現(xiàn)有效抑制。該方法不僅適用于存在相干信號的情況下，而且不受陣列結(jié)構(gòu)限制，應(yīng)用更加靈活。通過仿真驗證:方法在保證受試設(shè)備特性不失真的情況下能夠有效抑制干擾噪聲。

關(guān)鍵詞：虛擬暗室； 壓縮感知； 波束形成

0引言

對于大型系統(tǒng)或已固定設(shè)備的電磁干擾(electromagneticinterference，EMI)問題，因為屏蔽室容量有限或設(shè)備已經(jīng)固定，EMI測試無法在屏蔽室內(nèi)完成，需要采用現(xiàn)場測試來進(jìn)行EMI測試。但是，普通環(huán)境中的背景噪聲會超出電磁兼容(electromagneticcompatibility，EMC)的極限值，并且被測設(shè)備(equipmentundertest，EUT)產(chǎn)生的電磁信號可能會被環(huán)境中的背景噪聲所調(diào)制。所以，想在普通環(huán)境下進(jìn)行測試就必須設(shè)法消除背景噪聲的影響。針對以上難題，MarinoJr和MichaelA于2002年提出虛擬暗室測試?yán)碚?，并?005年申請和獲得了虛擬暗室的專利[1]。該理論自適應(yīng)噪聲抵消技術(shù)，使用雙通道分布測試克服單通道測試帶來的時間差問題。基于虛擬暗室理論設(shè)計的CASSPER系統(tǒng)被美國軍方采購和使用。FrenchA和BraunS等人還提出了基于時域自適應(yīng)噪聲抵消的電磁輻射發(fā)射現(xiàn)場測試方法[2]。但這些基于自適應(yīng)噪聲抵消的虛擬暗室技術(shù)要求測試環(huán)境中不能有與EUT同頻的信號分量以及參考通道中不能存在EUT信號等。DmochowskiJP提出基于波束形成的干擾抵消技術(shù)[3]，不需要參考通道并且具有很好的干擾抑制效果。盧中昊和林銘團(tuán)等人提出了基于陣列信號的虛擬暗室技術(shù)[4]，通過實驗證明方法的可行性。但文中采用的算法只能用于獨立信號并且受到陣列結(jié)構(gòu)限制，使得應(yīng)用條件受到限制。

當(dāng)EUT放置在室內(nèi)環(huán)境，如艙室環(huán)境中，EUT信號和背景電磁噪聲在空間中不斷地發(fā)生反射、散射。由于空間相對狹小，反射后的強(qiáng)度并沒有減弱多少，會產(chǎn)生相干信號。同時，由于實際的測量環(huán)境中，陣列的排布方式往往是受實際尺寸限制的，而絕大多數(shù)實際條件不能夠滿足均勻陣布陣要求[5,6]。針對以上情況，本文根據(jù)空域濾波原理，引入壓縮感知(compressivesensing,CS)理論和波束成形技術(shù)，實現(xiàn)環(huán)境干擾的濾除，既可以有效處理相干信號，并且可以使用于任意陣列。

1基本原理

如圖1所示，在EUT周圍設(shè)置天線陣列，每個天線單元均可以接收到由EUT輻射信號和背景干擾信號所組成的混合信號，這個混合過程與各輻射源到天線的電磁傳播路徑有關(guān)，利用基于CS理論的空間譜估計技術(shù)獲取各個信號的波達(dá)方向(directionofarrival，DOA)。根據(jù)信號的空間譜估計結(jié)果對各陣列接收信號做最優(yōu)加權(quán)，在保證EUT輻射信號不失真的前提下，使陣列波束在干擾信號的來向形成“零陷”，實現(xiàn)對信號的空域濾波，進(jìn)而來準(zhǔn)確反映出EUT的真實輻射信號特性。

圖1　基于空域濾波的電磁輻射發(fā)射現(xiàn)場測試體系Fig 1　Emission field of electromagnetic radiation testing system based on spatial filtering

2基于CS的DOA估計

CS[7]理論表明:當(dāng)信號具有稀疏性或可壓縮性時，通過采集少量的信號投影值就可以實現(xiàn)信號的準(zhǔn)確或近似重構(gòu)[8]。

(1)

式中S僅含有K個(K?N)個非零值，此時稱信號X是可壓縮的。利用M×N(M?N)維觀測矩陣Φ對可壓縮信號進(jìn)行觀測，其表達(dá)式為

Y=ΦX=ΦΨZ=ΘS

(2)

式中Φ∈RM×N為測量矩陣，Ψ∈RN×N為信號的稀疏表示基矩陣。S∈RN×1為稀疏表示系數(shù)，y∈RM×1為測量向量，Θ∈RM×N為由測量矩陣和稀疏表示基矩陣構(gòu)成的字典，也稱之為信息算子。由CS理論可知，如果Θ滿足有限等距性(restricted isometry property，RIP)等稀疏重構(gòu)條件，則可以求解以下L0范數(shù)為題和極高的概率精準(zhǔn)重構(gòu)系數(shù)信號S。

假設(shè)有K個遠(yuǎn)場窄帶信號入射到有M個的任意平面天線的陣列上，第k個信號的入射角度為θk。入射信號的波長為λ。t時刻陣列接收的快拍數(shù)據(jù)矢量可以表示為

X(t)=A(θ)S(t)+N(t),t∈{1,2,…,T}

(3)

式中X=[x(1),…,x(T)]T,N=[n(1),…,n(T)]T，A為M×K維的陣列導(dǎo)向矢量矩陣，表達(dá)式為

A=[a(θ1)a(θ2) …a(θK)]

(4)

a(θi)=[e-jωτ1i,e-jωτ2i,…,e-jωτMi]T

(5)

式中ω為入射信號的中心頻率，τli為第i個信號在第l個陣列位置相對參考陣元的延時(波程差)，且ωτli=2π(dxlcosθl+dylsinθi)/λ，(dxl,dyl)為第l個陣元的位置坐標(biāo)。

輻射信號的潛在方位按照等角度空間劃分方法，分成均勻的網(wǎng)格,抽樣網(wǎng)格數(shù)量為n，設(shè)第i個潛在信號信號的入射角度為θi，且θi∈[-π,π]，則網(wǎng)格索引為{θ1,θ2,…,θn}，假定每個方位都對應(yīng)潛在的信號si(i=1,2,…,p)，構(gòu)造出一個n×1維的信號S=[s1,s2,…,sn]，顯然S中只有k個位置有非零的元素，其余n-k的位置均為零，由于潛在信源個數(shù)遠(yuǎn)遠(yuǎn)大于實際信源個數(shù)，因此，信號S是稀疏的[9]。

針對公式的求解可以采用MOMP算法、MORMP算法、MFOCUSS算法和RMFOCUSS算法等，文獻(xiàn)[10]對上述集中算法做了分析比較,在同等條件下RMFOCUSS算法比較其他算法恢復(fù)成功率高且誤差小。

本文采用RMFOUSS算法實現(xiàn)稀疏信號的恢復(fù)和重構(gòu)[11]，其算法如下

Wk+1=diag(ck[i]1-p/2)

(6)

其中

(7)

(8)

其中

(9)

其中,p的選取與信號稀疏度有關(guān),p越小表示信號越稀疏，λ為反映噪聲功率大小的變量，通過調(diào)節(jié)參數(shù)λ可以除去部分噪聲。當(dāng)噪聲服從標(biāo)準(zhǔn)正太分布時，則有λ=σ2，其中，σ2為噪聲的方差。

通過RMFOCUSS算法可以求出譜估計公式為

P(θi)=‖S(i,:)‖2,i=1,2,…,L

(10)

3陣列波束控制與零點形成

適應(yīng)波束成形在某種最優(yōu)準(zhǔn)則下通過自適應(yīng)算法來實現(xiàn)權(quán)集尋優(yōu)，它能適應(yīng)各種環(huán)境的變換，實時地將權(quán)集調(diào)整到最佳位置附近。

波束形成算法是在一定準(zhǔn)則下綜合各輸入信息來計算最優(yōu)權(quán)值的數(shù)學(xué)方法。如表1，這些準(zhǔn)則主要包括最小均誤差(MMSE)準(zhǔn)則、最大信噪比(MSNR)準(zhǔn)則和線性約束最小方差(LCMV)準(zhǔn)則等[12]。

在電磁兼容現(xiàn)場測量中,由于需要提取EUT輻射信號的準(zhǔn)確數(shù)據(jù),所以,必須在抑制噪聲干擾的情況下保持期望信號不失真。上述LCMV準(zhǔn)則在滿足線性約束條件wHα(θ)=1時,可以保證提取的EUT輻射信號不失真，這種方法也就是最小方法無畸變響應(yīng)(MVDR)波束形成器。

為了有效抑制干擾，在獲得空間信號DOA后，采用MVDR準(zhǔn)則確定各陣元的最佳加權(quán)系數(shù)，對陣列方向圖進(jìn)行綜合，使陣列主波束指向被測信號，并在干擾來向形成零點，盡可能保證被測信號不失真，最大化抑制干擾信號功率。MVDR在期望信號來向上增益為1。代價函數(shù)為J(w)=wHRw,約束條件為wHa(θ)=1。a(θ)為期望信號的方向相量。令

(11)

求式(11)的梯度并令其為零，得到

(12)

解得其權(quán)重向量為w=λR-1a(θ)

(13)

將式(13)帶入約束條件中求得

(14)

將式(14)帶入式(13)中可以求出最優(yōu)權(quán)重向量為

(15)

4仿真實驗結(jié)果

設(shè)接收陣列為8陣元的均勻線陣，間距為200 m。同時設(shè)空中有3個不同的信號，數(shù)學(xué)形式如下

(16)

信號的來向分別為40°，20°，-20°，信噪比為10 dB。

4.1EUT和干擾信號來向仿真

通過經(jīng)典MUSIC和CS算法對EUT輻射信號的DOA進(jìn)行Matlab仿真得到實驗結(jié)如圖2。

圖2　MUSIC和CS的DOA估計結(jié)果Fig 2　DOA estimation results of MUSIC and CS

由圖2可知：MUSIC在存在相干信號的情況下對EUT輻射信號和干擾信號的測向結(jié)果出現(xiàn)偏差，而CS-DOA算法可以準(zhǔn)確獲得EUT輻射信號和干擾信號的來向。同時CS-DOA估計算法對陣列結(jié)構(gòu)沒有限制，可以運用于非均勻限制情況，更加貼近實際應(yīng)用。

4.2波束成形仿真

采用MVDR準(zhǔn)則通過Matlab進(jìn)行波束形成仿真得到結(jié)果如圖3。

圖3　MVDR波束形成結(jié)果Fig 3　Result of MVDR beamforming

如圖3所示，陣列空間譜估計的主波束指向被測信號，不僅準(zhǔn)確獲取EUT輻射信號和干擾信號的DOA，而且在干擾來向處形成兩個極深的“零陷”。

從圖4、圖5和圖6是輸入輸出頻譜對比，在保證期望信號失真的情況下，在300 Hz有效抑制噪聲100 dB，在400 MHz的頻譜峰值處與EUT原輻射信號抑制，該頻率點的相干信號噪聲基本被抑制，同時對其他頻率的高斯噪聲有0～25 dB的有效抑制。

圖4　EUT輻射信號的頻譜Fig 4　Spectrum of EUT radiated signal

圖5　陣列信號輸入的頻譜Fig 5　Spectrum of array signal input

圖6　陣列信號輸出的頻譜Fig 6　Spectrum of array signal output

5結(jié)論

本文提出一種基于CS的虛擬暗室技術(shù)用于解決現(xiàn)場電磁發(fā)射測量中相關(guān)信號干擾問題。該方法將CS理論、MVDR自適應(yīng)波束形成技術(shù)引入到現(xiàn)場電磁兼容測試領(lǐng)域。通過仿真驗證,該方法能夠有效抑制現(xiàn)場輻射發(fā)射測試環(huán)境中的干擾信號，恢復(fù)出的EUT信號具有較好的保真度。該方法不僅可以應(yīng)用在環(huán)境中存在相干信號的情況下，亦不受陣列結(jié)構(gòu)限制，測試陣列可以根據(jù)測量場地的實際環(huán)境調(diào)整，在實際應(yīng)用中更加靈活。

參考文獻(xiàn):

[1]Marino J,Michael A.System and method for measuring RF radiated emissions in the presence of strong ambient signals:US,6980611[P].2000—02—03.

[2]French A,Brauns S,Russer P.Time-domain EMI measurements

in the presence of ambient[C]∥Proc of the International Symposium on Electromagnetic Compatibility,2009:139-142.

[3]Dmochowski J P.Decoupled beamforming and noise cancella-tion[J].IEEE Trans on Instrumentation and Measurement,2007,56(1):80-88.

[4]Lu Z H,Liu P G,Liu J B.A novel method of ambient interference suppressing for in-situ electromagnetic radiated emission test[J].IEEE Trans on Electromagnetic Compatibility,2012,54(6):1205-1215.

[5]蔣伯峰,王文杰,殷勤業(yè).適用于任意陣列的多徑信道二維方向角與相對時延的聯(lián)合估計方法[J].電子學(xué)報,2000(12):1-4.

[6]楊鵬,楊峰,聶在平,等.基于共形天線陣的免搜索來波方向估計算法研究[J].電波科學(xué)學(xué)報,2012,27(2):241-245,325.

[7]Donoho D L.Compressed sensing[J].IEEE Transactions on Information Theory,2006,52(4):1289-1306.

[8]賀亞鵬,李洪濤,王克讓,等.基于壓縮感知的高分辨DOA估計[J].宇航學(xué)報,2011(6):1344-1349.

[9]Malioutov D,Cetin M,Willsky A S.A spare signal reconstruction perspective for source localization with sensor arrays[J].IEEE Transactions on Signal Processing,2005,53(8):3010-3022.

[10] Cotter S F,Rao B D,Engan K,et al.Spare solutions to linear inverse problem with multiple measurement vectors [J].IEEE Transaction on Signal Processing ,2005,53(7):2477-2488.

[11] DaoB,Engan K,Cotter S F,et al.Subset selection in noise based on diversity measure minimization[J].IEEE Transaction on Signal Processing,2003,51(3):760-770.

[12] 張小江，汪飛，陳偉華.陣列信號處理的理論和應(yīng)用[M].2版.北京：國防工業(yè)出版社，2013.

VirtualdarkchambertestmethodbasedonCS

ZHANGYao-wen，CHENShao-chang

(SchoolofElectronicEngineering,NavalUniversityofEngineering,Wuhan430000,China)

Abstract:In order to filtering out background noise effectively in on-site electromagnetic emission test,virtual dark chamber test method based on compressive sensing(CS)is proposed to solve the problem of failure of algorithm for multiple signal classification(MUSIC)when correlation signals exist in environment.Using theory based on CS to achieve directions of arrival(DOA) of radiation signal of equipment under test(EUT)and other interference signal in environment,and using minimum variance distortionless response(MVDR)bream-forming algorithm for beamforming,so as to effectively suppress interference signal.This method can not only be applied in the presence of coherent signals,but also has more flexible applications without restrictions in array structure. Simulation results prove that the method can suppress interference noise effectively without distortion of character of EUT.

Key words:virtual dark chamber; compressive sensing(CS); beamforming

DOI:10.13873/J.1000—9787(2016)04—0059—04

收稿日期：2015—07—08

中圖分類號：TM 931

文獻(xiàn)標(biāo)識碼：A

文章編號：1000—9787(2016)04—0059—04

作者簡介:

章耀文(1992-)，男，江蘇蘇州人，碩士研究生，主要研究方向為電磁兼容。