劉春靜

(中國電子科技集團公司第三十八研究所 合肥 230088)

一種改善和差測角盲區(qū)的快速算法

劉春靜

(中國電子科技集團公司第三十八研究所 合肥 230088)

傳統(tǒng)的和差測角方法在雷達、對抗、通信等多個領(lǐng)域得到了廣泛的應用，但受限于信噪比的影響，角度測量盲區(qū)較大。超分辨算法(如MUSIC等)能夠達到非常高的角度測量精度，但計算復雜度隨著陣元或波束個數(shù)的增大而增加，為了降低實時處理的難度，同時考慮到實際環(huán)境中都是針對單目標的角度測量，提出了一種基于兩波束的超分辨MUSIC算法，該算法涉及到的二維矩陣特征分解在角度估計過程中具有閉式形式，不僅運算復雜度很低，而且適用于目前大多數(shù)的陣列處理系統(tǒng)。

和差測角;超分辨;快速算法;MUSIC

0 引言

和差測角是一種廣泛使用的角度測量技術(shù)，由于測量精度高、計算簡便易用等優(yōu)點，應用于各種工程設(shè)計領(lǐng)域。和差測角的基本原理是利用和差波束之間幅度比較估計目標偏離天線軸中心的角度，根據(jù)和差波束之間的相位比較估計目標偏離天線軸中心的方向。當和波束接收到的目標信噪比不夠高時，差波束接收到的信號將被噪聲淹沒，達不到正確測量的要求，從而導致固有的角度測量盲區(qū)。測量盲區(qū)的大小與被檢測目標的信噪比有關(guān)，信噪比越高，測量盲區(qū)越小，而常用的和差測角精度與信噪比的表達式并沒有考慮到測量盲區(qū)的影響。

文獻［1，6－8，10，12，15，17］研究了和差測角性能與目標信噪比、天線幅相不一致性、阻抗匹配誤差等因素的關(guān)系;文獻［3，4，14］研究了二維陣列的和差角度測量及工程化設(shè)計方法，提出了快捷簡便的查表法測角;文獻［2，9，16］針對高精度測量雷達中的和差跟蹤測角進行了應用分析，提出了基于跟蹤濾波器的角度測量;文獻［5，11，13］討論了干擾環(huán)境下的和差測角問題，給出了可行的解決方案。雖然研究和差測角的文獻資料很多，但尚未有關(guān)于測角盲區(qū)的分析，本文從和差測角的基本原理出發(fā)，分析了角度測量盲區(qū)與目標信噪比的關(guān)系，提出采用兩波束的超分辨角度測量算法(多重信號分類 MUSIC，MUtiple SIgnal Classification)，不僅能夠大大提高角度測量精度，徹底消除和差測量盲區(qū)，而且推導了兩波束MUSIC算法的快速求解閉式，具有運算復雜度低、實時處理容易實現(xiàn)的優(yōu)點。此外，該方法不限于一維陣列的使用，能夠適用于任意二維或三維陣列。

1 和差測角盲區(qū)分析

設(shè)一維N元均勻直線陣列的陣元間距為d，工作頻率為 f0，波長為 λ0，波束指向為 θ0，a(θ) 為天線陣列的N維導向矢量，s為接收信號的包絡(luò)波形，n為天線陣列的接收機噪聲向量，則和差波束的輸出信號分別為:

其中wΣ、wΔ分別為和差波束的N維加權(quán)向量，x為天線陣列接收的N維信號向量，即:

(·)H表示矩陣的共軛轉(zhuǎn)置。

為了使得加權(quán)后的和差波束輸出保持相同的噪聲能量，若第i個天線接收通道的噪聲為零均值、方差為σ2

i的高斯白噪聲，則和差加權(quán)系數(shù)應當滿足以下條件:

當N個接收通道的噪聲方差近似相等時，則上述條件簡化為

從和差波束的接收信號表達式可以看出，接收信號的能量由加權(quán)信號和加權(quán)噪聲兩部分組成，當和差波束接收到的加權(quán)信號能量都很強時，和差比幅測角能夠獲得較好的測量精度;然而在實際系統(tǒng)的目標檢測中，如果以12dB信噪比為和波束的檢測門限，則對于略大于檢測門限的目標信號，差波束在零軸附近區(qū)域的接收信號幅度可能非常小，此時差波束中的加權(quán)噪聲占據(jù)了主導地位，不僅不能夠正確測量出和差波束的幅度比值，也不能夠正確判斷目標偏離零軸的方向。

2 改善和差測角盲區(qū)的快速算法

為了減小和差測角盲區(qū)，在信噪比一定的條件下，采用波束域的超分辨MUSIC算法提高測角精度;由于只針對一個目標測量角度，避免矩陣特征分解導致的巨大運算量，將對稱的天線陣面(一般情況下都能夠滿足該條件)劃分為兩個相同的陣面，分別形成兩個指向相同、具有恒定相位差的空間接收波束。超分辨測量的角度范圍很小，不會超過一個波束寬度，因此在陣面允許的條件下，兩個波束的相位中心可以選擇較大的間距差，而不會受到角度模糊的影響。

設(shè)兩個空間接收波束的采樣數(shù)據(jù)向量為:

其中a為2×1的空間信號導向矢量，s(t)為入射信號，n(t)為高斯白噪聲向量，采樣數(shù)據(jù)向量的協(xié)方差矩陣估計為

其中M為選擇的數(shù)據(jù)樣本數(shù)，數(shù)據(jù)樣本的選擇必須剔除掉其它疑似目標，保證所有的數(shù)據(jù)樣本中只有一批期望的目標，協(xié)方差矩陣Rx的特征多項式為:

則特征值λ的兩個根分別為

其中E為單位矩陣，則協(xié)方差矩陣最小特征值(即噪聲特征值)對應的特征向量為:

根據(jù)超分辨MUSIC算法的原理，空間譜計算表達式為:

其中a(θ)為方向θ的陣列導向矢量，空間譜掃描的峰值位置則為目標角度的估計值。

根據(jù)以上算法步驟估計兩波束超分辨MUSIC的運算量，在M個樣本的采樣時間內(nèi)，即MTs(Ts為采樣時間間隔)，需要完成協(xié)方差矩陣估計(4M次復乘法、M次復加法、1次復除法)、特征值計算(6次復乘法、6次復加法、1次求根和1次復除法)、噪聲特征向量計算(1次復加法、1次復除法)、空間譜計算(若空間掃描L個角度，則為2L次復乘法、L次復除法和復加法);整個運算量為(4M+6+2L)次復乘法、(M+7+L)次復加法、(3+L)次復除法、1次求根運算。為了進一步降低運算量，可以將空間譜表達式的峰值搜索簡化為求根 MUSIC，通過對PMUSIC(θ)求導直接估計出目標的最優(yōu)角度。

3 仿真實驗

假設(shè)均勻直線陣列的陣元個數(shù)為64，工作頻率為300MHz，單元間距為半波長，和波束采用30dB的泰勒加權(quán)，差波束采用30dB的Bayliss加權(quán)，則方位法線方向上角度測量盲區(qū)與目標信噪比的關(guān)系如圖1所示。從圖中可以看出，隨著目標信噪比的增加，角度測量盲區(qū)不斷減小;當信噪比足夠大時(例如30dB)，角度測量盲區(qū)已經(jīng)能夠達到0.05度以下。

將64元均勻直線陣劃分為兩個相同的對稱陣面，從而形成方位向的左波束和右波束，兩個波束之間在某相同角度上存在恒定相位差，兩波束超分辨MUSIC算法的協(xié)方差矩陣估計樣本數(shù)為30，角度掃描間隔為0.001°，Monte－Carlo仿真次數(shù)為500次，圖2為兩波束超分辨MUSIC算法的角度估計均方根誤差(RMSE)隨目標信噪比變化的關(guān)系曲線，從圖中可以看出，即使在只有30個樣本的條件下，目標信噪比為10dB時已經(jīng)能夠達到0.035°的角度測量精度，遠遠優(yōu)于傳統(tǒng)的和差測角方法。

4 結(jié)論

本文通過詳細分析和差測角盲區(qū)產(chǎn)生的原理，闡述了傳統(tǒng)測角方法依賴于信噪比的缺陷，由于實際環(huán)境中總能夠針對單個目標的角度進行測量，因此通過采用兩波束超分辨MUSIC算法實現(xiàn)一個波束寬度內(nèi)的目標角度估計，不僅可以達到非常高的角度測量精度，而且具有簡單可行的閉式表達式，能夠應用于各種形式的天線陣列系統(tǒng)中。

［1］馬振球，崔嵬.相位和差單脈沖雷達測角性能分析［J］.北京理工大學學報，2009，29(8):726－731.

［2］員琳紅.和差多波束跟蹤算法及工程實現(xiàn)［D］.西安:西安電子科技大學，2012.

［3］保的.相控陣二維和差測角方法及其實現(xiàn)研究［D］.西安:西安電子科技大學，2010.

［4］楊蓓蓓.一種二維數(shù)字陣列雷達的和差波束測角方法［J］.雷達與對抗，2014，34(3):6－10.

［5］楊玲.基于自適應干擾抑制與和差測角性能的陣形評估［D］.西安:西安電子科技大學，2011.

［6］王玉濤.單脈沖測角技術(shù)及工程實現(xiàn)［D］.西安:西安電子科技大學，2007.

［7］周明宇.毫米波導引頭信號處理若干關(guān)鍵技術(shù)研究［D］.南京:南京理工大學，2015.

［8］孔永飛.相控陣自適應和差波束測角研究［D］.西安:西安電子科技大學，2012.

［9］孔曉丹.單脈沖跟蹤雷達跟蹤測量技術(shù)應用［D］.西安:西安電子科技大學，2011.

［10］陳知明，魏述剛.相控陣雷達測角誤差自適應計算方法［J］. 現(xiàn)代雷達，2013，35(7):13－15.

［11］王峰.雷達主瓣干擾環(huán)境下的測角方法研究［J］.現(xiàn)代雷達，2015，37(1):21－23.

［12］李文臣，李青山，馬飛.相位和差單脈沖相控陣天線方向圖仿真與性能分析［J］.中國電子科學研究院學報，2011，6(4):336－339.

［13］張忠傳，邰新軍，孟慶昌，孫東.主/旁瓣干擾對單脈沖測角得影響分析［J］.雷達科學與技術(shù)，2012，10(2):207 －211.

［14］毛祺，安紅，周先敏.二維相位和差單脈沖雷達的測角性能分析［J］.電子信息對抗技術(shù)，2007，22(6):15 －18.

［15］Wei Zhou et al.Angle measurement accuracy analysis of sum-difference amplitude-comparison monopulse in onshore or shipborne ISAR［J］.Proceedings of the 19thInternational Conference on Digital Signal Processing，2014，8:20－23.

［16］Evgeny Markin.Method of automatic target angle tracking by sum-and-difference monopulse radar invariant against the polarization jamming［J］.Proceedings of the 7thEuropean Radar Conference，2010，10:499 －502.

［17］David A.Shnidman.An optimization technique for sum and difference array factors［J］.IEEE Transactions on Aerospace and Electronic Systems，2004，40(1):360－370.

A Fast Algorithm for Improving Blind Zone of Sum-Difference Angle Estimation

Liu Chunjing
(The No.38 Research Institute of CETC，Hefei 230088)

The traditional sum-difference angle estimation method is widely applicable in many fields such as radar，electronic countermeasures and communication;however，blind zone of angle estimation is too wide due to impact of signal-to-noise ratio.Although super-resolution algorithm(such as MUSIC etc.)can achieve very high angle measurement accuracy，computational complexity increases with the number of array elements or beams.In order to reduce difficulty of real-time processing(generally angular measurement is to single target under real environment)，a super-resolution MUSIC algorithm based on two beams is proposed.Eigen-decomposition of two-dimensional matrix involved in the proposed algorithm can be expressed as a closed form to achieve angle estimation，therefore，the proposed algorithm can be implemented with very low complexity，and applicable for most array processing systems.

sum-difference angle estimation;super-resolution;fast algorithm;MUSIC

TN911.7

A

1008-8652(2017)01-005-03

2016-10-18

劉春靜(1985－)，男，高級工程師，主要研究方向為陣列信號處理和自適應信號處理。