基于RCB與鎖相環(huán)的跟蹤鑒別測向方法

任曉松，楊嘉偉，崔 嵬，吳嗣亮，安建波，歐春湘

（1．北京遙感設備研究所，北京100854；2．中國航天科工防御技術研究院，北京100854；3．北京理工大學信息與電子學院，北京100081）

基于RCB與鎖相環(huán)的跟蹤鑒別測向方法

任曉松1，楊嘉偉2，崔 嵬3，吳嗣亮3，安建波1，歐春湘1

（1．北京遙感設備研究所，北京100854；2．中國航天科工防御技術研究院，北京100854；3．北京理工大學信息與電子學院，北京100081）

為了彌補陣列天線導向矢量失配和相位測量噪聲對測向性能的影響，提出一種基于穩(wěn)健capon波束形成技術（robust capon beamforming，RCB）和鎖相環(huán)的矢量最優(yōu)估計與跟蹤鑒別測向方法。首先基于RCB與鎖相環(huán)原理，對目標來波信號導向矢量進行最優(yōu)估計與跟蹤測量；然后在穩(wěn)定、準確跟蹤導向矢量的基礎上，借鑒擴頻接收機偽碼鑒相原理確定來波信號方向。仿真分析表明，該方法能夠彌補陣列流型失配和相位噪聲的影響、突破角度相關間隔的限制，準確測量來波信號的方向。

陣列信號處理；穩(wěn)健波束形成；鎖相環(huán)；角度鑒別；導向矢量；測向

0 引 言

相關干涉測向方法適用于多種陣列天線，具有測向準確度高的特點，在軍事和民用領域中得到了廣泛應用［1-2］。其基本原理是測量所接收信號在各陣元間的相位差，并與樣本相位差進行相關運算，以代價函數(shù)最大方向作為來波信號方向的測量值。導向矢量與目標信號方向之間存在一一對應關系，與相關干涉測向方法中的相位差測量樣本也存在對應關系，因此可以利用導向矢量的相關性進行測向。然而，在實際應用中陣列天線不可避免地存在各種誤差（如陣元位置擾動誤差、通道幅相誤差）和相位噪聲，這些因素會導致實際的陣列形式與理想狀態(tài)存在偏差［35］，因此，如果簡單基于相位差測量值構(gòu)造與陣列相關的導向矢量，勢必會影響到測向性能。

在過去30年里，有很多優(yōu)秀的改善波束形成的穩(wěn)健性方法，主要可分為基于特征空間的自適應波束形成、多點約束和對角加載技術等幾類［6-10］，可以有效克服系統(tǒng)誤差帶來的不利影響，如信號被抑制、輸出信干噪比下降等。文獻［10］提出了基于不確定集的穩(wěn)健波束形成方法（robust capon beamforming，RCB）；該方法將要估計的信號功率作為目標函數(shù)，在給定的誤差范圍內(nèi)對導向矢量進行搜索，選擇與波束形成輸出功率的最大值對應的導向矢量作為真實導向矢量的估計，并用估計的導向矢量進行波束形成。該方法在提高信干噪比抑制系統(tǒng)誤差等干擾的同時，對目標信號真實導向矢量實現(xiàn)了最優(yōu)估計。因此，可以借鑒RCB原理，將最優(yōu)估計得到的導向矢量應用到相關干涉測向方法中。

相關測向方法中角度的測量精度與相關角度間隔有關，角度間隔越小測量精度越高，但是角度間隔不能無限減小，并且角度間隔越小相關搜索帶來的運算量越大。數(shù)字接收機鎖相環(huán)（phase-locked loop，PLL）［11-14］直接對載波相位進行跟蹤，通過鑒相器提取并輸出相位估計誤差，經(jīng)過環(huán)路濾波器得到相位誤差控制量，當環(huán)路穩(wěn)定時具有較高的跟蹤精度。而導向矢量與陣元相位差之間存在對應關系，并且在來波信號方向發(fā)生變化時，導向矢量與相位差同時產(chǎn)生對應變化；因此，可以借鑒數(shù)字接收機PLL原理對來波信號導向矢量從相位角度進行閉環(huán)跟蹤測量。在環(huán)路鎖定精確跟蹤導向矢量之后，借鑒偽碼跟蹤環(huán)路鑒相器原理確定穩(wěn)定跟蹤環(huán)路的最終結(jié)果（導向矢量）所對應的角度，從而突破角度間隔限制、測量精度與相關搜索之間的矛盾。

仿真結(jié)果表明，本文提出的基于RCB與PLL的跟蹤鑒別測向方法，不僅彌補了陣元位置擾動、通道幅相誤差和相位噪聲的影響，還突破了角度相關間隔與測量精度之間的矛盾。

1 陣列信號模型

假設K個窄帶遠場信號入射到M元均勻線陣（uniform linear array，ULA）上（K＜M），則t時刻陣列接收的數(shù)據(jù)為

式中，快拍矢量X（t）＝［x1（t），x2（t），…，xM（t）］T、噪聲矢量N（t）＝［n1（t），n2（t），…，nM（t）］T均為M×1維的陣列接收數(shù)據(jù)；信號包絡S（t）＝［s1（t），s2（t），…，sK（t）］T為K×1維，第i個信號為具有幅度α和初始相位φ的非隨機信號，即

基于均勻線陣的導向矢量矩陣為A＝［a1，a2，…，aK］，其中，ai為與第i個來波信號對應的導向矢量，對應的來波方向為θi?；陉囋g距為d的均勻線陣和窄帶來波信號si的陣列模型示意如圖1所示。

圖1 均勻線陣示意

以第一個陣元為參考陣元，則與等距線陣和si對應的導向矢量ai可表示為

理想情況下，假設單目標信號（即K＝1）、入射角為θ情況下，根據(jù)式（2）和式（3）可以得出t時刻陣列接收的數(shù)據(jù)為

由于通道幅相誤差和陣元位置擾動誤差的存在，陣列接收的數(shù)據(jù)變?yōu)?/p>

可以記a～＝Γa，表示由于幅相誤差和陣元位置擾動誤差形成的真實導向矢量，其中

定義R＝E｛X（t）XH（t）｝為陣列接收信號的相關矩陣，其中，E｛·｝表示數(shù)學期望。

2 基于RCB與PLL的跟蹤測向

2.1 算法提出

眾所周知，對于確定的天線陣和目標信號，有式（1）和式（3）所示的接收信號模型和與目標信號方向有關的導向矢量，導向矢量與相關干涉測向方法中使用的相位差樣本具有一一對應關系?？梢越梃b相關干涉儀的測角原理［15－17］，選取導向矢量樣本進行相關匹配處理，以相關序列極大值方向作為目標信號方向。

利用導向矢量進行目標信號方向的相關測量時，需要準確測量導向矢量，可以說來波信號導向矢量測量越準確，測向越精確。然而，由式（5）可知，接收信號受到通道幅相誤差、陣元位置擾動等誤差和噪聲的影響，來波信號的導向矢量與理想導向矢量產(chǎn)生偏差；并且角度的估計精度與相關測向的角度間隔有關。

因此，為了彌補通道幅相誤差和陣元位置擾動誤差的影響，突破角度相關間隔的限制，消除噪聲的影響，提高測向精度，借鑒RCB算法和PLL原理提出基于RCB與數(shù)字接收機PLL的設計方法。利用RCB對來波信號真實導向矢量進行最優(yōu)估計；然后，將估計得到的導向矢量送給導向矢量閉環(huán)跟蹤系統(tǒng)，使本地導向矢量始終跟隨目標信號的導向矢量。在環(huán)路鎖定之后，借鑒偽碼跟蹤環(huán)路鑒相器原理設計角度鑒別器，確定目標信號角度。

本文所提方法由3部分構(gòu)成，分別是：基于RCB的導向矢量估計單元、基于PLL的導向矢量跟蹤單元和角度鑒別單元。算法結(jié)構(gòu)框圖如圖2所示。

2.2 基于RCB的導向矢量估計

基于接收到的快拍數(shù)據(jù)建立球形RCB［10］的優(yōu)化模型為

在式（8）約束下，可以避免平凡解as＝0的出現(xiàn)，除非ε＝（即此時as＝0位于約束集合的邊界上）。對于式（8）中等式約束最優(yōu)化問題，利用Lagrange乘數(shù)法進行有效求解有

式中，λL為Lagrange乘數(shù)。對as求導，并令其等于零，可得最優(yōu)解利用矩陣求逆引理可得

圖2 算法框圖

最優(yōu)Lagrange乘數(shù)λL可以利用牛頓迭代法通過求解約束方程g（λL）＝‖（I＋λLR）－1a～s‖2＝ε獲得；然后代入最優(yōu)權(quán)矢量表達式從而獲得對目標信號真實導向矢量的最優(yōu)估計值。

圖3 RCB算法原理

式中，Γa為真實導向矢量；D＝diag［ejδ1，ejδ2，…，ejδM］為殘留相位噪聲引起的導向矢量偏差。

因此，基于RCB算法的導向矢量最優(yōu)求解方法，初步彌補了通道幅相誤差和陣元位置擾動誤差的影響，抑制了相位噪聲的影響，可以實現(xiàn)對目標信號真實導向矢量的估計。

2.3 導向矢量跟蹤環(huán)路

為了進一步消除RCB算法估計所得導向矢量中殘留相位噪聲的影響，實現(xiàn)對來波信號導向矢量的動態(tài)跟蹤，基于PLL原理提出導向矢量的閉環(huán)跟蹤方法。跟蹤環(huán)路主要由導向矢量鑒別器、環(huán)路濾波器和本地導向矢量數(shù)控振蕩器（numerically controlled oscillator，NCO）3部分組成，組成如圖2所示。鑒別器用于鑒別RCB估計導向矢量與本地導向矢量的復指數(shù)誤差信號（與相位樣本對應）；環(huán)路濾波器用于濾除導向矢量誤差信號中的高頻諧波分量，進一步濾除殘留相位噪聲的影響，對輸入導向矢量進行精確估計；本地導向矢量NCO根據(jù)導向矢量初始值以及濾波器輸出值，不斷調(diào)整NCO的頻率字保持對接收信號導向矢量的精確跟蹤。

2．3．1 導向矢量鑒別

根據(jù)式（10）可知，與當前時刻來波方向?qū)膶蚴噶抗烙嬛禐?/p>

本地NCO產(chǎn)生的導向矢量為

式中，~D為本地NCO生成的初始導向矢量偏差，是D經(jīng)環(huán)路濾波和NCO累加消除高頻噪聲所得。對式（12）取共軛，然后與式（11）相乘可得

式中

為來波信號角度變化引起的導向矢量偏差；

為相位噪聲殘差引起的導向矢量偏差。

2.3.2 環(huán)路濾波

與數(shù)字接收機載波PLL一樣，導向矢量跟蹤環(huán)路濾波器［12］的作用是濾除導向矢量誤差信號中的高頻諧波分量，使其在輸出端保持對來波信號導向矢量的精確跟蹤估計。

在角度變化不快的情況下，跟蹤環(huán)路濾波器的階數(shù)不需要太高，可以采用一階環(huán)路濾波器。導向矢量跟蹤環(huán)路濾波器設計如圖4所示。

圖4 導向矢量一階濾波器數(shù)字實現(xiàn)

濾波器傳遞函數(shù)為式中，Kd為鑒頻器增益；Ko為數(shù)控振蕩器增益；K為環(huán)路增益；F（z）為環(huán)路濾波器轉(zhuǎn)移函數(shù)；wn為自然頻率；T為環(huán)路更新時間。

2.3.3 本地跟蹤導向矢量NCO

本地跟蹤導向矢量NCO根據(jù)環(huán)路濾波器的輸出結(jié)果，不斷調(diào)整NCO的頻率字，保持對輸入信號導向矢量的精確跟蹤。導向矢量NCO的基本原理如圖5所示。

圖5 數(shù)控振蕩器結(jié)構(gòu)圖

2.4 角度鑒別

在導向矢量跟蹤環(huán)路鎖定之后，借鑒偽碼跟蹤環(huán)鑒相器原理設計角度鑒別器，確定目標信號角度。以上一時刻本地導向矢量為中心，生成左偏和右偏角度θd／2的導向矢量；然后，分別計算左偏導向矢量與跟蹤環(huán)路實時得到導向矢量之間的距離幅值pls，右偏導向矢量與跟蹤環(huán)路實時得到導向矢量之間的距離幅值prs?；趯蚴噶?范數(shù)進行歸一化并取反加一操作之后的角度鑒別原理如圖6所示。

圖6 角度鑒別原理

基于本地導向矢量構(gòu)造左偏角導向矢量和右偏角導向矢量為

左偏幅值與右偏幅值為

與幅值偏差對應的角度誤差鑒別為

從而確定當前時刻測向結(jié)果為

式中，θ（n）和θ（n－1）分別代表當前時刻角度鑒別算法確定的來波信號方向和上一時刻的來波信號方向。

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

設置仿真條件：天線為10陣元的均勻直線陣、每個陣元均為全向天線、理想狀態(tài)下相鄰陣元間距為半個波長。陣元通道間幅度誤差服從均值為0，方差為0.1的高斯分布；相位誤差服從均值為0，方差為20°的高斯分布。陣元位置擾動引起的d／λ服從均值為0．5，方差為0．03的高斯分布。

3.1 算法性能分析

在上述仿真條件下假設來波信號信噪比為15 dB，采樣快拍數(shù)為40，目標來波信號方位向為－10°，對本文算法進行仿真分析如圖7所示?；赗CB估計得到的導向矢量進行相關測向的結(jié)果如圖7中虛線所示；利用鎖相環(huán)路進一步消除導向矢量中高頻干擾的測向結(jié)果如圖7中實線所示，也即本文算法得到的最終測量值?？梢园l(fā)現(xiàn)，經(jīng)過環(huán)路跟蹤處理后，測量值更趨近于真實來波信號方向。

圖7 算法性能對比示意

這說明本文提出的基于PLL的導向矢量閉環(huán)跟蹤方法進一步消除了噪聲干擾的影響，從而提高了測向精度。

在不同信噪比和采樣快拍數(shù)條件下，對靜態(tài)來波信號進行80次蒙特卡羅仿真，統(tǒng)計測角均方根誤差如圖8所示?？梢园l(fā)現(xiàn)，在相同信噪比情況下，隨著快拍數(shù)增加，測角均方根誤差具有逐漸減小的趨勢；在相同采樣快拍數(shù)情況下，隨著信噪比增加，測角均方根誤差具有逐漸減小的趨勢；并且當信噪比達到20 dB時，采樣快拍數(shù)的影響減小，測角均方根誤差趨于一致。

圖8 影響測向因素分析

因此，利用本算法測量來波信號方向時，可以結(jié)合信噪比大小對采樣快拍數(shù)進行調(diào)整，從而實現(xiàn)最優(yōu)的參數(shù)配置。

3.2 靜態(tài)仿真分析

在上述仿真條件下，設置來波信號信噪比為15 dB，采樣快拍數(shù)為40；在［－10°，10°］區(qū)間內(nèi)，以2°為間隔設置靜態(tài)來波信號。利用本文方法、基于相位重構(gòu)導向矢量的傳統(tǒng)相關測向方法和基于RCB算法獲得導向矢量進行相關測向的方法，對各個來波方向分別進行80次蒙特卡羅仿真，統(tǒng)計測向結(jié)果如圖9所示。

可以發(fā)現(xiàn)，基于相位重構(gòu)導向矢量的傳統(tǒng)相關測向均方根誤差數(shù)量級最大；利用RCB算法獲得的導向矢量進行相關測向，其相關測向結(jié)果具有10－2數(shù)量級的測向精度；而本文算法則突破了角度相關間隔的限制，并具有較高的測向精度，達到了10－3數(shù)量級的測向精度。

圖9 靜態(tài)測向均方根誤差對比

3.3 動態(tài)仿真分析

設置目標信號入射角度范圍為［－3°，3°］，以某項目為背景設置目標信號以5（°）／s的角速度勻速運動；來波信號信噪比為15 dB，采樣快拍數(shù)為40。利用本文算法可以實現(xiàn)對來波信號的穩(wěn)定跟蹤，環(huán)路穩(wěn)定跟蹤動態(tài)來波信號角度仿真結(jié)果如圖10所示。

圖10 動態(tài)測向示意

利用本文方法在［－3°，3°］區(qū)間范圍內(nèi)對動態(tài)來波信號進行80次蒙特卡羅仿真，統(tǒng)計測向結(jié)果，并與基于相位重構(gòu)導向矢量的傳統(tǒng)相關測向方法和基于RCB算法獲得導向矢量進行相關測向的結(jié)果進行比較，如圖11所示。動態(tài)仿真表明，本文算法能夠?qū)聿ㄐ盘柕慕嵌茸兓M行穩(wěn)定跟蹤測量，具有較高的測量精度。

圖11 動態(tài)測向均方根誤差對比

4 結(jié) 論

本文所提方法，首先利用RCB算法對來波信號導向矢量進行最優(yōu)估計；為了進一步消除系統(tǒng)噪聲的影響，基于PLL原理實現(xiàn)了對來波信號導向矢量的準確跟蹤測量?；诮邮諜C偽碼跟蹤環(huán)路鑒相器原理設計的角度鑒別器，有效突破了角度相關間隔的限制。通過仿真分析表明，該方法彌補了陣列天線系統(tǒng)誤差和測量相位隨機噪聲的影響，可以實現(xiàn)來波信號方向的高精度測量。

［1］Kebeli M．Extended symmetrical aperture direction finding using correlative interferometer method［C］∥Proc.of the 7th IEEE International Conference on Electrical and Electronics Engineering，2011：209- 213．

［2］Wei H W，Shi Y G．Performance analysis and comparison of correlative interferometers for direction finding［C］∥Proc.of the 10th IEEE International Conference on Signal Processing，2010：393- 396.

［3］Zhu L Y，Meng H E，A robust minimum variance beamformer with new constraint on uncertainty of steering vector［J］.Signal Processing，2006，86（9）：2243- 2254.

［4］Jiang X，Zeng W J，Yasotharan A，et al．Robust beamforming by linear programming［J］.IEEE Trans.on Signal Processing，2014，62（7）：1834- 1849.

［5］Robert G L，Stephen P B．Robust minimum variance beamforming［J］.IEEETrans.on Signal Processing，2005，53（5）：1684- 1696.

［6］Ruan H，De L，Rodrigo C．Robust adaptive beamforming using a low-complexity shrinkage-based mismatch estimation algorithm［J］. IEEESignal Processing Letters，2014，21（1）：60- 64.

［7］Vorbyov S A，Gershman A B，Luo Z Q．Robust adaptive beamforming using worst case performance optmization a solution to the signal mismatch problem［J］.IEEE Trans.on Signal Processing，2003，51（2）：313- 324.

［8］Elnashar A，Elnoubi S M，El-Mikati H A．Further study on robust adaptive beam forming with optimum diagonal loading［J］. IEEE Trans.on Antennas Propagation，2006，54（12）：3647- 3658.

［9］Gong C，Huang L，Xu D，et al．Knowledge-aided robust adaptive beamforming with small snapshots［J］.Electronics Letters，2013，49（20）：1258- 1259.

［10］Li J，Stoica P，Wang Z S．On robust Capon beamforming and diagonal loading［J］.IEEE Trans.on Signal Processing，2003，51（7）：1702- 1715.

［11］Hu H，Sun H Z，Ji，Z Y．Study on algorithm and control strategies of GPS carrier tracking loop under high dynamic condition［J］.Journal of Astronautics，2011，32（8）：1805- 1812．（胡輝，孫函子，紀兆云．高動態(tài)GPS載波跟蹤算法和環(huán)路控制策略研究［J］．宇航學報，2011，32（8）：1805- 1812．）

［12］Kaplan E D．Understanding GPS principles and applications［M］．Boston：Artech House，1996．

［13］Ravindra B，Jinling W．Dynamics performance of carrier and code tracking loops in ultra-tight GPS／INS／PL integration［C］∥Proc.of the IEEE India Conference，2005：233- 236.

［14］Roncagliolo P A，Cristian E D B，Carlos H M．GPS digital tracking loops design for high dynamic launching vehicles［C］∥Proc.of the 9th IEEE International Symposium on Spread Spectrum Techniques and Applications，2006：41- 45．

［15］An X J．Research on DF based on improved phase interferometer［J］.Radio Engineering of China，2009，39（3）：59- 61．（安效君．改進的干涉儀測向方法研究［J］．無線電工程，2009，39（3）：59- 61．）

［16］Wu Y W，Rhodes S，Satorius E H．Direction of arrival estimation via extended phase interferometry［J］.IEEE Trans.on Aerospace and Electronic Systems，1995，31（1）：375- 381．

［17］Xiao X L．The theory of interferometer in DOA［J］.China Radio，2006，（5）：43- 49．（肖秀麗．干涉儀測向原理［J］．中國無線電，2006，（5）：43- 49．）

Tracking and discriminating direction finding method based on robust capon beamforming and phase-locked loop

REN Xiao-song1，YANG Jia-wei2，CUI Wei3，WU Si-liang3，AN Jian-bo1，OU Chun-xiang1
（1.Beijing Institute of Remote Sensing Equipment，Beijing 100854，China；2.China Aerospace Defense Technology Institute，Beijing 100854，China；3.School of Information and Electronics，Beijing Institute of Technology，Beijing 100081，China）

In order to compensate the impact of the antenna array steering vector mismatch and the phase measurement noise，the optimal measurement method is proposed for steering vector，and the tracking identification method of direction finding is designed，basing on the robust capon beamforming（RCB）and phase-locked loop．Firstly，using the RCB algorithm estimate the optimal signal steering vector，and then basing on the principle of phase-locked loop to closed-loop track target signal steering vector．After steady and accurate tracking of the steering vectors，referencing the principle of spread spectrum receiver pseudo-code phase detector，the direction of the target signal is determined．Simulation results show that the proposed method can compensate the impact of the steering vector mismatch and phase measurement noise，break through the restrictions of the angle correlation interval，and accurate measure the target signal direction．

array signal processing；robust beamforming；phase-locked loop（PLL）；angle identification；steering vector；direction finding

TN911

A

10．3969／j．issn．1001-506X．2015．04．04

任曉松（1981-），男，工程師，博士研究生，主要研究方向為陣列信號處理、雷達信號處理。E-mail：lengyunrxs＠163．com

楊嘉偉（1963-），男，研究員，博士，主要研究方向為雷達系統(tǒng)。E-mail：yjw＠163．com

崔 嵬（1974-），男，研究員，博士，主要研究方向為遙測遙控。E-mail：cuiwei＠bit．edu．cn

吳嗣亮（1964-），男，研究員，博士，主要研究方向為雷達技術。E-mail：siliangwu＠bit．edu．cn

安建波（1976-），男，高級工程師，碩士，主要研究方向為制導控制。E-mail：anjianbo＠163．com

歐春湘（1988 ），女，工程師，博士研究生，主要研究方向為衛(wèi)星導航。E-mail：ouchunxiang＠163．com

1001-506X（2015）04-0746-06

2014- 05- 28；

2014- 10- 26；網(wǎng)絡優(yōu)先出版日期：2014- 11- 26。

網(wǎng)絡優(yōu)先出版地址：http：／／w ww．cnki．net／kcms／detail／11．2422．TN．20141126．1439．001．html