基于差分模型的近場無源定位算法

趙 研，陶海紅，暢 鑫，董春曦，趙國慶

（1.西安電子科技大學雷達信號處理全國重點實驗室，陜西西安 710071；2.西安電子科技大學電子工程學院，陜西西安 710071）

0 引 言

輻射源定位是目標探測領域的重要應用，在雷達、地震、聲學、海洋和聲吶等專業有著廣泛的需求。無源定位問題按照陣列孔徑、目標距離和工作波長的關系可以分為兩分支，（a）典型天線陣列遠場條件（r≥2D2/λ）和（b）天線陣列近場條件，其中r是目標到天線陣列的距離，D是陣列孔徑，λ是目標信號的波長。在遠場條件下，波達方向（DOA）是分析天線陣列平行波前［1］信號傳播過程,經過多年研究目前已經取得了長足發展和大量成果。當滿足近場條件或者滿足菲涅爾區條件時，電波波前不再滿足平行假設，傳統單要素線性DOA 角度估計方法不再適用，但也具備了解算目標未知的可能（同時估計目標角度和距離）。本文針對典型近場輻射源定位問題，通過對于距離和角度的非線性聯合分析可獲得目前位置信息。

近場目標定位近些年是較為熱門的問題，傳統上采用的方法有：1）經典采用多重用戶分離（MUSIC）的二維拓展算法［2］；2）利用菲涅爾近似波前的二階表達式估計［3］，在此基礎上利用陶布利茲矩陣旋轉不變性，進行ESPRIT 的DOA 估計［4］或高階累積量的DOA 估計［5］；3）通過菲涅爾近似的特殊陣列結構協方差矩陣構造，可以建立角度的單維問題，并由此通過MUSIC 估計距離［6-8］,利用MUSIC 代價函數估計最佳值［9-10］，此外基于子陣遠場估計模型，還提出一種三角函數近場逼近方法［11］。然而以上的方法都是通過二維要素的一維空間化簡求解過程。

不同于傳統近場目標輻射傳播的菲涅爾近似模型，基于近場空間幾何結構和三角函數關系推導得到新的差分迭代模型［12］，本文建立全新高精度目標關系模型并形成定位解析表達式，通過陣列處理估計TOA 定位觀測量，利用新定位算法快速準確地獲取近場多個獨立同分布目標的定位信息。

1 近場菲涅爾近似目標定位

對于均勻線陣（ULA），近場目標二維定位的幾何模型如圖1 所示。其中均勻線陣包含N個天線陣元，相鄰陣元的間距為d，其中原點陣元、第i個陣元和第k個目標，在三角形ΔOIK中，從O到K的距離為rk，從第I個陣元到K的距離為rik，而陣列方向和rk之間夾角為θk，陣列方向與rik之間的夾角為θik，對于三角形ΔONK，依據余弦定理有

圖1 近場目標定位結構示意圖

參見圖1，(rnk,θnk)代表第k個輻射源對第n個節點的距離和到達角，假設λ為電波波長，則傳輸的相位延遲為

式中Rs是目標信號協方差矩陣。求解近場目標數量和位置的問題變成了對采樣數據序列M個樣點的估計。針對式（6）傳統研究提出了多種處理方法，基于WVD 的分析估計［3］,基于子陣遠場估計（MUSIC/ESPRIT）的統計分析［4,6，11］,4 階累積量的子空間分析［5］,通過分析二階協方差矩陣數學模型，構造特殊選擇節點陣列，抵消數學結構中的二次項而只保留一階線性相移項，通過對表達式的算法模型估計［7］,得到最優化目標位置估計逼近值［8-10］。

2 基于差分模型的定位新算法

不同于傳統的菲涅爾近似，文獻［12］中，作者提出一種新的近場電波傳播近似模型，其結構與初值和迭代次數相關，可知［12］

2.1 基于二次迭代的定位分析

2.2 基于三次迭代的定位分析

傳統的定位算法都是基于菲涅爾的二階近似關系式進行推導估計，其性能會繼承近似基礎誤差影響，而按泰勒級數展開三階近似雖可提升近似精度，但所產生的二元三次方程組難以獲得可行解。基于差分迭代新模型，可獲取對于目標定位的三次迭代新算法，有［12］

令x= cosθk，由此建立方程：

令n= 1,2建立方程組

由于表達式（18）存在規則的結構形式，使得方程組（19）消元操作較為便捷，可將二元三次方程轉化為一元三次方程，而后者存在通解公式，并由此獲得(rk,θk)實數形式的非零可行解，由于具體解數學表達式較長，可通過數學工具軟件解得，此處不再贅述。

新算法的特征優勢在于基于較高精度的近場傳播近似模型條件下建立了定位方程，使得模型自身的基本估計誤差優于傳統菲涅爾近似的結果［12］，也使得其數學表達式非線性特征要更為強烈。但是正如上面所推導，其結構的規律性使得二元三次復雜方程組依舊可推導和存在閉式結構的解析可行解。基于傳統菲涅爾近似條件下，以往算法表達式具有了更良好的線性特征和簡潔數學結構，通過經典的最優化估計、子空間分解或累積量分析都可以得到一定準則下的最優定位參數解，但其算法運算量往往較大，對很多實際客觀存在的非線性誤差適應性低，新算法采用了更加精準的近場模型，并且推導定位參數閉式解，在保證精度的同時使得后續實際應用中數值計算存在很大的自由度，為進一步精簡高效運算提供可能。

2.3 相關處理與算法性能分析

在近場陣列輻射源定位的整體信號處理過程中，各通道輸入信號間一般都滿足相互獨立且隨機同分布，噪聲為加性高斯的前提假設，因此對于實際陣列條件下式（4）可對等表示為

式中△tnl==rnk/C。

對于式（5）有

由此可得所提算法輸入參數△Rnk= △tnkC。對于TOA 的精確估計還可以參考其他算法文獻［13-15］，本文主要重點在后面的基于TOA 定位觀測量的定位算法研究，這里不再贅述。

從文獻［12］的性質分析，可以看到新的算法模型（二次迭代）與傳統菲涅爾近似模型（二階近似）算法仿真分析，兩者在象限區間中各有一般的近似優勢；而如果采用優化初值的三次迭代新算法則近似程度明顯提高，在0°～82°和97°～180°近場范圍內近似精度都顯著提高，參見圖2。由此可知在高信噪比條件下（數值精度可分辨情況下），新算法的定位精度應優于菲涅爾近似的定位精度。

圖2 所提算法與TSECL2rd定位誤差對比圖

3 對所提算法仿真和性質分析

針對均勻線陣的近場區域范圍內目標定位性質進行仿真分析，對比分析實驗：1）針對近場陣列的子陣遠場分解及多個參數回歸分析估計算法［11］；2）基于菲涅爾近似建立的定位方程二階條件求解［3］，參見式（16）；3）基于高階累計量的陣列統計估計定位方法［5］；4）文獻［12］提出基于二次差分迭代計算定位方程新算法，參見式（14）；5）本文所提出基于三次差分迭代計算定位方程新算法，由于其目標定位參數閉式解結構較為復雜，具體參見附錄。

由圖3 可知，隨著信噪比條件（噪聲比S/N=Psignal/Pnoise，其中Psignal為信號能量，Pnoise為與信號時、頻等域相匹配噪聲能量）從-20 dB 到30 dB 以1 dB 的步長逐漸變換，高階累計量算法的測距精度最高誤差在1%左右；所提出三次迭代比二次迭代算法精度略高，也逼近1%的誤差水平，指標性能分別為第2 和第3；菲涅爾近似算法與遠場子陣定位算法由于具有噪聲分解與原始數據信息測量優勢，在低信噪比條件下效果相對較好，但整體定位性能相對較差。

圖3 各種算法測距誤差隨信噪比變化性能圖

由圖4可知，隨著信噪比條件從-20 dB到30 dB以1 dB 的步長逐漸變換，從曲線結果可知所提出三次迭代定位算法具有優良的近場模型非線性特征適應能力，在大于0 dB 信噪比具有良好的測角性質，測角誤差約1.1%；而高階累計量算法的測角精度也較為優秀，誤差也逼近1%，兩者性能較優且逼近性能漸近線；文獻［12］提出二次迭代算法精度要好于遠場子陣定位算法和菲涅爾近似定位算法，而這兩種算法在低信噪比條件下的性能要優于前面的算法。

圖4 各種算法測角誤差隨信噪比變化性能圖

通過對于近場目標距離和角度的綜合測量評估，5種典型方法的定位誤差對比如圖5所示。

圖5 各種算法定位誤差隨信噪比變化性能圖

由圖5 可知，5 種近場定位方法整體定位誤差隨信噪比提升而改善。而其中新提出的三次迭代估計誤差與高階累積量估計方法的定位指標最好，均在1%左右，且兩者誤差值十分接近，且在負信噪比條件下性能較差，并隨著信噪比提升誤差快速減小；新提出二次迭代方法的誤差估計性能弱于前兩種算法，但在高信噪比條件時，性能要優于遠場子陣估計和菲涅爾近似定位方法；而這兩種定位算法在低信噪比條件性能相對穩定，誤差惡化程度減小。

4 結束語

本文所提算法具有3個優點：1）具備代數差分迭代結構形式，可以實現近似精度更高的閉式解析定位解表達式，其三次迭代解結構是同階次菲涅爾近似難以推導得到的高精度解；2）通過高精度逼近的傳播模型迭代推導得到的近場定位算法，可以得到正信噪比條件下和高階累積量定位算法誤差基本相同的定位結果能力；3）在相同陣列近場正信噪比工作條件下，新三次迭代算法與高階累計量定位誤差性能較好，相比其他方法定位性能改善0.3%～3%，而本方法是代數計算與結果統計，相比于最優化分析和最佳點估計過程，其計算分析過程與數值計算量大為降低，處理復雜度和時間代價更低。本文所提算法優勢和特征較為突出，對于近場DOA問題的解決提供新思路。