空間角度相近多干擾源DOA估計*

王強，張永順，司文濤

(空軍工程大學 防空反導學院，陜西 西安　710051)

?

空間角度相近多干擾源DOA估計*

王強，張永順，司文濤

(空軍工程大學 防空反導學院，陜西 西安710051)

摘要：針對空間角度相近多個干擾源DOA估計時，經(jīng)典MUSIC算法估計性能與分辨效果不佳的問題，提出了一種改進的MUSIC算法。該算法利用干擾與目標回波信號的特征，劃分出干擾子空間、信號與噪聲子空間，重新構(gòu)造了空間譜函數(shù)。通過對經(jīng)典MUSIC算法與改進的MUSIC算法DOA估計的仿真分析，結(jié)果表明，在小快拍數(shù)、多干擾源空間角度相近、低信干比的情況下，改進的MUSIC算法估計性能與分辨精度更佳。

關(guān)鍵詞：多干擾源；DOA估計；MUSIC算法

0引言

隨著科技的迅猛發(fā)展，現(xiàn)代戰(zhàn)場電磁干擾變得日益復雜，對雷達天線實施主瓣干擾已變得可能。主瓣干擾的一大特征是干擾機同目標空間方位分離度小，倘若進入接收主波束內(nèi)的干擾源方向未知，將使得雷達從信號處理角度實施主瓣干擾抑制，變得更加困難。對于主瓣內(nèi)存在單個干擾源，即使方向未知，因其強度遠大于目標回波信號與噪聲的強度，利用簡單的空間譜估計方法(最小方差法)就能獲取方向信息。然而對于主瓣內(nèi)存在多個干擾源，尤其是天線波束為針狀波束的情況，干擾源之間空間方位相差甚微。就經(jīng)典的DOA(direction of arrival)估計算法——MUSIC(multiple signal classification)算法而言，估計性能與分辨率均顯著變差，影響到主瓣內(nèi)干擾源個數(shù)的確定[1-2]。結(jié)合文獻[3]中提到的多基地雷達聯(lián)合信息處理抗主瓣干擾方法，該方法能夠不抑制目標回波而實現(xiàn)主瓣干擾。但前提條件是接收站的個數(shù)至少比干擾源的個數(shù)多1個。如果主瓣內(nèi)干擾源個數(shù)不確定，會影響多基地雷達的布站配置，限制了該方法優(yōu)勢的發(fā)揮。因此主瓣內(nèi)多干擾源的DOA估計，是一個亟待解決的問題。

針對以上分析，本文提出了一種改進MUSIC算法，該算法在小快拍數(shù)、多干擾源空間角度相近、低信干比的情況下，對干擾源DOA的估計與分辨性能更佳。

1陣列天線信號模型

假設(shè)在一個平面天線陣列上有M個陣元，天線主波束內(nèi)存在K(K

X′=S(t)+A(θ)G(t)+N(t),

(1)

式中：S(t)=(S1(t),S2(t),…,SM(t))T為M個陣元輸入端的目標回波信號矢量；G(t)=(G1(t),G2(t),…,GK(t))T為K個干擾源形成的干擾信號矢量；A(θ)=(a(θ1),a(θ2),…,a(θk))為干擾導向矩陣，a(θk)=(1,e-j2πdsinθk/λ,…,e-j(M-1)2πdsinθk/λ)T，θk為第k個干擾源的波達方向;N(t)為M個陣元內(nèi)部噪聲。

由式(1)得到陣列輸出總信號的協(xié)方差矩陣為[5]

R′=E[XXH]=Rs+ARgAH+σ2I，

(2)

式中：Rs，Rg分別為目標回波信號、干擾信號的協(xié)方差矩陣；σ2為噪聲功率；I為M維單位陣。

(3)

式中：L為譜估計的快拍數(shù)。

2經(jīng)典MUSIC算法原理

經(jīng)典MUSIC算法的數(shù)據(jù)協(xié)方差R由信號、噪聲各自協(xié)方差矩陣2部分組成，首先對R進行特征分解求出其所有特征值，然后將特征值進行大小依次排序，選取較大的特征值構(gòu)成信號對角陣Λs，其對應的特征向量構(gòu)成信號子空間Ωs，將較小的特征值構(gòu)成噪聲對角陣Λn，對應的特征向量構(gòu)成噪聲子空間Ωn[7-8]。則特征分解后，R可表示為[9]

(4)

結(jié)合信號的導向矢量ρ(θ)與噪聲子空間Ωn的正交特性，MUSIC空間譜函數(shù)可寫為[10]

(5)

選取一定的搜索補償Δθ，則θi=(i-1)Δθ。通過計算PMUSIC(θi)的所有峰值，其對應的θi便是待估計信號的波達方向。

3改進的MUSIC算法

為了有效估計空間內(nèi)多干擾源的空間方位，本文對MUSIC算法的步驟流程作出了改進?；谇懊嫣岬降年嚵刑炀€信號模型，改進的MUSIC算法首先對數(shù)據(jù)協(xié)方差矩陣R′進行特征分解，由于干擾功率大于目標回波信號以及接收機內(nèi)部噪聲功率， 因此在估計干擾源方位時，可將目標回波信號同噪聲作為一個整體考慮。協(xié)方差矩陣R′可分為干擾信號協(xié)方差、目標回波信號同噪聲協(xié)方差2部分，然后依照大小順序排列所有特征值，選取較大的特征值構(gòu)成干擾對角陣Λg，其對應的特征向量構(gòu)成干擾子空間Ωg，將較小的特征值構(gòu)成信號與噪聲對角陣Λsn，對應的特征向量構(gòu)成信號與噪聲子空間Ωsn。則特征分解后，R′可表示為

(6)

結(jié)合文獻[11-13]對MUSIC空間譜函數(shù)作出改進，其表示式具體為

(7)

式中：a(θ)為干擾的導向矢量；ν為調(diào)節(jié)因子，表達式為

(8)

4仿真結(jié)果與分析

仿真采用8×1等距直線陣，陣元間距設(shè)為半波長，搜索補償Δθ=0.02π，仿真分2種情況(圖1～5)：

(1) 接收機內(nèi)部噪聲為高斯白噪聲，有2個獨立干擾源掩護目標，對地面雷達天線實施干擾，圖1中2個干擾源的來襲方向為-10°，10°。

(2) 圖2中2個獨立干擾源的來襲方向為5°，7°，圖3中2個獨立干擾源的來襲方向為5°，6°，其他條件設(shè)定同上。

圖1～3中設(shè)定信干比為-2 dB，快拍數(shù)為500，圖4中快拍數(shù)為500，圖5中信干比為-2 dB。

圖1　兩干擾源空間角為20°時，不同算法的空間譜比較Fig.1　Space spectrum comparison of different algorithms　　　in the two interference sources space angel of 20

圖2　兩干擾源空間角為2°時，不同算法的空間譜比較Fig.2　Space spectrum comparison of different algorithms　　　　in the two interference sources space angel of 2

圖3　兩干擾源空間角為1°時，不同算法的空間譜比較Fig.3　Space spectrum comparison of different algorithms　　　in the two interference sources space angel of 1

圖4　兩干擾源空間角為1°時，信干比對算法測向誤差的影響Fig.4　Algorithm direction error on SIR in the two　　　interference sources space angel of 1

圖5　兩干擾源空間角為1°時，快拍數(shù)對算法　　　測向誤差的影響Fig.5　Algorithm direction error on snapshots in the　　　two interference sources space angel of 1

前3幅圖主要是算法分辨性能的仿真對比，通過觀察可以看出：①當2干擾源空間角度為20°時，經(jīng)典MUSIC算法與改進MUSIC算法測向誤差均不大，且改進的MUSIC算法的估計性能更佳。②當2干擾源空間角度由2°下降為1°時，經(jīng)典的MUSIC算法估計性能嚴重變差，且隨著2個干擾源空間方位的接近，分辨性能明顯下降，而改進的MUSIC算法有效估計出2個干擾源的來襲方向。

后2幅圖主要是算法測角誤差的仿真對比，通過觀察可以看出：①當設(shè)定2干擾源空間角度為1°，且快拍數(shù)一定時，2種算法的測向誤差隨信干比的增加而逐漸減小，改進的算法性能更優(yōu)。②當設(shè)定2干擾源空間角度為1°，且信干比一定時，2種算法的測向誤差隨快拍數(shù)的增加而逐漸減小，改進的算法性能更優(yōu)。

5結(jié)束語

基于空間角度相近多個未知方向干擾源的DOA估計，經(jīng)典的MUSIC算法的估計性能與分辨效果不佳。對此，本文根據(jù)干擾與目標回波信號的特征，劃分了干擾子空間、信號與噪聲子空間，提出了一種改進的MUSIC算法。通過仿真對比分析，結(jié)果表明，改進的MUSIC算法估計性能與分辨精度更佳，對干擾源DOA的估計與分辨性能更佳。

參考文獻：

[1]HE Jin, SWAMY M N S, AHMAD M O. Efficient Application of MUSIC Algorithm Under the Coexistence of Far-Field and Near-Field Sources[J]. IEEE Transa-ctions on Signal Processing, 2012, 60(4): 2066-2067.

[2]Yoke Leen Sit, Christian Sturm, Johannes Baier, et al. Direction of Arrival Estimati-on Using the MUSIC Algorithm for a MIMO OFDM Radar[J]. 2012 IEEE, 2012, 4(12): 26-227.

[3]Victor S Chernyak. Fundamentals of Multisite Radar Systems:Multistatic Radars and Multiradar Systems[M]. Beijing: Publishing House of Elsctronics Industry,2011: 70-222.

[4]梁國龍, 張鍇, 范展, 等. 單矢量傳感器MUSIC算法的DOA估計及性能評價[J].哈爾濱工程大學學報, 2012, 33(1):31-32.

LIANG Guo-long, ZHANG Kai, FAN Zhan, et al. Performance Evaluation of DOA Estimation Using a Single Acoustic Vector-Sensor Based on an Improved MUSIC Algorithm[J]. Journal of Harbin Engineering University, 2012, 33(1): 31-32.

[5]趙謙, 董民, 梁文娟. DOA估計算法的一種修正MUSIC算法的研究[J]. 計算機工程與應用, 2012, 48(10): 103-104.

ZHAO Qian, DONG Min, LIANG Wen-juan. Research on Modified MUSIC Algorithm of DOA Estimation[J]. Computer Engineering and Applications, 2012, 48(10): 103-104.

[6]王陽, 魯祖坤, 石宇, 等. MIMO雷達中基于Cyclic-Music算法的DOA估計[J]. 中國雷達, 2012, 30(4): 30-31.

WANG Yang, LU Zu-kun, SHI Yu, et al. DOA Estimation Based on Cyclic-Music Algorithm in MIMO Radar[J]. China Radar, 2012, 30(4): 30-31.

[7]司偉建, 藍曉宇, 劉學. 提高二維DOA估計分辨率的改進MUSIC算法[J]. 西安電子科技大學學報:自然科學版, 2012, 39(2): 139-140.

SI Wei-jian, LAN Xiao-yu, LIU Xue. Modified MUSIC Algorithm for Improving the Resolution of 2-D Direction of Arrival Estimation[J]. Journal of Xidian University：Nature Science ed, 2012, 39(2): 139-140.

[8]AKKAR S, GHARSALLAH A. Reactance Domains Unitary MUSIC Algorithms Based on Real-Valued Orthogonal Decompositi-on for Electronically Steerable Parasitic Array Radiator Antennas[J]. IET Micro-Waves, Antennas & Propagation, 2012, 6(2): 226-227.

[9]楊永建, 王晟達, 馬健, 等. 基于MUSIC 算法的圓柱共形陣DOA估[J]. 空軍工程大學學報:自然科學版, 2012, 13(5): 66-68.

YANG Yong-jian, WANG Sheng-da, MA Jian, et al. DOA Estimation for Cylindrical Coonformal Array Antenna Based on MUSIC Algorithm[J]. Journal of Air Force Engineering University:Nature Science ed, 2012, 13(5): 66-68.

[10]劉芬, 胡文彪. 基于MUSIC 算法的變頻電機轉(zhuǎn)子故障診斷研究[J]. 船電技術(shù), 2013, 33(1): 40-41.

LIU Fen, HU Wen-biao. Diagnosis Research of Rotor Fault in Variable Frequency Motor Based on MUSIC Algorithm[J]. Marine Electric & Electronic Engineering, 2013, 33(1): 40-41.

[11]楊志偉, 賀順, 廖桂生. 加權(quán)偽噪聲子空間投影的修正MUSIC算法[J]. 信號處理, 2011, 27(1): 1-5.

YANG Zhi-wei, HE Shun, LIAO Gui-sheng. Modified MUSIC Approach with Weighted Pseudo-Noise Subspace Projection[J]. Signal Processing, 2011, 27(1): 1-5.

[12]梁國龍, 張鍇, 張光普，等. 基于二元矢量陣MUSIC算法的相干源分辨研究[J]. 高技術(shù)通訊, 2012, 22(4): 335-336.

LIANG Guo-long, ZHANG Kai, ZHANG Guang-pu, et al. Research on the Resolution in Coherent Sources Using a MUSIC Algorithm Based on a Dual Vector-Sensor Array[J]. High Technology Letters, 2012, 22(4): 335-336.

[13]劉志剛, 汪晉寬, 王福利. 基于實值分解技術(shù)的循環(huán)root-MUSIC算法[J]. 系統(tǒng)仿真學報, 2006, 18(9): 2438-2442.

LIU Zhi-gang, WANG Jin-kuan, WANG Fu-li. Cyclic root-MUSIC Algorithm Based on Real-Valued Decomposition[J]. Journal of System Simulation, 2006, 18(9): 2438-2442.

DOA Estimation of Space Angle Close Multiple Interference Sources

WANG Qiang, ZHANG Yong-shun, SI Wen-tao

(AFEU，Air and Missile Defense School, Shaanxi Xi’an 710051,China)

Abstract:Focusing on the problem that classical MUSIC(multiple signal classification)algorithm has poor estimation performance and bad resolution effect on DOA (direction of arrival) estimation of close space angle with multiple interference sources, a kind of improved MUSIC algorithm is proposed. It classifies interference subspace, signal with noise subspace, and gives the new structure of space spectrum function by using the features of interference and target echo signal. Through the simulation analysis on classical MUSIC algorithm and the improved MUSIC algorithm, the result states the latter has the better estimation performance and resolution effect, especially, for instance, a small number of snapshots, the low SIR(signal interference ratio) and closer space angel

Key words:multiple interference sources; direction of arrival(DOA) estimation; multiple signal classification(MUSIC) algorithm

中圖分類號：TN973

文獻標志碼：A

文章編號：1009-086X(2015)-01-0120-04

doi:10.3969/j.issn.1009-086x.2015.01.020

通信地址：710051陜西省西安市長樂東路甲字1號防空反導學院研2隊E-mail：1019611183@qq.com

作者簡介：王強(1989-)，男，陜西三原人。碩士生，主要從事多基地雷達信號處理研究。

基金項目：陜西省自然科學基金資助項目(2010JQ8007)

收稿日期：2014-01-14；
修回日期：2014-03-29