NLS次優簡化算法及其在單站無源定位中的應用

梅勇兵,南建設

(中國西南電子技術研究所,成都 610036)

NLS次優簡化算法及其在單站無源定位中的應用

梅勇兵,南建設

(中國西南電子技術研究所,成都 610036)

為了簡化基于外輻射源固定單站無源定位求解算法以用于工程實現,建立了定位原理極坐標模型,給出了極坐標系下非線性最小二乘(NLS)定位算法,根據極坐標系下定位算法的特點,為簡化計算,提出了NLS次優簡化算法。仿真結果表明,NLS次優簡化算法有較好的穩定度和定位精度,為工程實現提供了參考。

單站無源定位;外輻射源;極坐標;非線性最小二乘;次優算法

1 引 言

無源定位對于提供系統在電子戰環境下的生存能力和作戰能力有重要的作用[1-6]。單站無源定位具有設備量少和系統相對獨立等優點。目前,單站無源定位主要利用兩種類型的信號:一種是基于定位目標自身輻射的信號,另一種是基于定位目標之外的外輻射源的信號。利用外輻射源可以對低空目標和采用射頻隱身技術的目標進行定位,本文主要討論基于外輻射源的單站無源定位算法。

基于外輻射源的單站無源定位[1]主要是靠定位站的測向信息以及外輻射源信號的直達與經待定位目標散射的時間差來進行定位,定位解為非線性方程。

目前,單站無源定位算法都是非線性濾波問題,主要采用的是擴展卡爾曼濾波算法[1-2]和非線性最小二乘[2-4]。擴展卡爾曼濾波及其他基于動態系統的濾波算法都需要估計動態方程的轉移矩陣并進行自適應修正,算法相對復雜。最小二乘算法只需要利用觀測方程,通過對非線性的觀測方程進行線性化處理,通常采用Taylor級數展開的方法,取一次項進行近似,近似之后的線性方程可以利用最小二乘進行求解。

基于外輻射源的單站無源定位觀測到的目標相對定位站的方位角、俯仰角和相對時延,極坐標模型可以避免觀測值與坐標值之間的頻繁轉換。特別是利用多組數據進行定位估計的時候,每組數據都需要進行直角坐標與方位角、俯仰角之間的轉換。而采用極坐標模型,僅在有需要時對估計結果進行一次坐標轉換即可。另外,NLS定位算法是一個多維尋優的問題,其計算中包含矩陣求逆等矩陣運算,而矩陣求逆工程實現比較困難。

本文根據基于外輻射源固定單站無源定位觀測量包含方位角和俯仰角的特點,采用極坐標建立了定位原理的模型;推導出了極坐標系下的非線性最小二乘(NLS)算法;根據極坐標系下NLS算法的特點,提出了NLS次優簡化算法,將三維尋優問題簡化為一維RLS求解;最后,通過計算機仿真對NLS次優簡化算法的收斂速度和定位精度進行仿真驗證。

2 基于外輻射源單站定位原理

2.1 定位原理

基于外輻射源單站定位是利用已知外輻射源的信號直接到達定位站與經定位目標散射后達到定位站的時間差信息,以及定位站利用散射信號或目標自身的輻射信號測向來進行定位。

通過測定時間差以及已知外輻射源與定位站的距離,可以確定待定位目標到外輻射源和到定位站的距離和;到兩點距離和一定的點分布在以兩點為焦點的一個橢球上。因此可以確定,待定位目標在該橢球面上。同時,通過定位站測向可以確定待定位目標在由測得俯仰角和方位角決定的一條示向線上。由上述信息決定的橢圓和示向線的交點就是待定位目標所在的位置(示向線起點為橢球內部焦點,有方向性,故只有一個交點)。

2.2 極坐標模型

設外輻射源A與定位站P之間的距離為d,且定位站測得待定位目標的方位角為θ、俯仰角為 φ,三者空間位置關系如圖1所示。

圖1 基于外輻源位單站無源定位位置關系圖Fig.1 Position relationship of single-observer passive location based on outer radiant point

以定位站P為極坐標原點,P到A的方向0°角方向建立極坐標系,則P的極坐標為(d,0,0)。

設待定位目標O的極坐標為(r,θ,φ),信號到達定位站的時間差為τ,則有如下關系:

式中,c為電磁傳播速度。

在極坐標下,O點坐標(r,θ,φ)中的 θ和φ已通過測向獲得,故確定目標位置僅需確定PO的距離r。

由式(1)知:

由可觀測值(τ,θ,φ)決定。

3 定位NLS求解

3.1 NLS算法

根據第2節內容,且 c常數,為記法簡單,定義觀測量 y=(cτ,θ′,φ′),其與目標的極坐標 z=(r,θ,φ)有如下關系:

將非線性方程f(z)進行線性化處理,在z0附近泰勒展開,并保留線性部分得到

偏導數函數如下:

根據文獻[3],式(6)的非線性最小二乘解的迭代形式如下:

式中,N為觀測誤差協方差矩陣,表達式為

3.2 簡化算法

式(10)中,f′(z(k))是矩陣形式,一般情況下不能像遞

推

最小二乘一樣避免矩陣求逆,用遞推方式進行求解。因此,提出一種降維的簡化方法,將三維尋優簡化為一維交替迭代尋優,降維后的最小二乘可以采用遞推最小二乘(RLS)進行計算,進一步簡化計算量。

在N=diag(λ1,λ2,λ3)為對角矩陣的時候,即觀測量是相互獨立的,極坐標模型的NLS算法的目標函數為

(1)對方位角θ進行最小二乘求解;

(2)對俯仰角 φ進行最小二乘求解;

(3)根據 θ、φ的優化結果,對r進行非線性最小二乘求解。

對方位角θ和俯仰角φ的搜索求解可采用RLS算法思想,此時,其第i步新息的權因子為。在完成θ和φ的更新后,作為新的估計值,代入對r進行搜索,可以采用RLS進行計算。

極坐標系下,基于外輻射源的單站無源定位NLS簡化算法如下。

(1)初始化

3.3 小結

新提出的NLS簡化算法是根據極坐標系下,方位角、俯仰角的觀測量與估計量是系數為1的線性關系的特點,將NLS三維尋優的問題簡化為一維尋優交互迭代的問題。通過簡化,可以避免NLS中矩陣求逆的問題;同時,簡化后可以采用RLS方法進行求解,提高穩定度和減少計算量。

簡化的NLS方法僅需要將第一次觀測作為初始值就可以迭代求解,遞推求解過程是采用已經成熟應用于工程的RLS算法實現的,初始值的計算只需要用到完全能工程化的乘加算法,故簡化的NLS算法能滿足工程實現的要求。

4 仿真分析

設外輻射源與定位站距離150 km,待定位目標初始位置距離定位站81.2 km,方位角1.446 rad,俯仰角0.123 rad。觀測時間間隔1 s,假設到達時差測量誤差均方差為20 ns,角度測量誤差均方差為1°。分別對靜態和慢時變目標進行定位仿真。

(1)靜態定位仿真

假設目標真實位置在觀測時間內部未發生變化或變化可以忽略不計,此時通過公式(4)的計算、采用NLS簡化遞推算法得到的結果如圖2所示。

圖2 靜態定位仿真結果Fig.2 Location results of static target

從圖2可知,NLS簡化算法有較快的收斂速度和穩定度,10次迭代之后,距離估計值偏差就在1 km以內,方位角頻差在0.005 rad以內,俯仰角偏差在0.008 rad以內。

(2)慢時變目標定位仿真

假設待定位目標在直角坐標系內x軸4 m/s、y軸10m/s、z軸6m/s進行運動,此時通過公式(4)的計算、采用NLS簡化遞推算法得到的結果如圖3所示。

圖3 慢時變目標定位仿真結果Fig.3 Location results of slow moving target

從圖3可知,對于慢時變目標,NLS簡化算法也有較好的收斂效果,但是隨著迭代次數的增加,估計結果慢慢偏離真實值,這主要是由于無法跟蹤上目標真實位置的變化。

(3)基于初始更新的定位仿真

為了校正偏離,可以采取周期性對NLS簡化算法的初始值進行更新的方式,進行定位計算。仿真結果如圖4所示。

圖4 周期初始化的NLS簡化算法結果Fig.4 Location results of low-complexity periodic initialization algorithm

從圖4可知,通過周期性初始化NLS簡化算法的初始值,能在一定程度上進行偏離校正。

5 結束語

單站無源定位的處理算法比較多,本文從工程實現簡單化的需求出發,給出了基于極坐標的非線性最小二乘的簡化遞推算法,將傳統的非線性最小二乘多維尋優簡化為一維交替迭代尋優,并利用成熟應用于工程的RLS算法來解決一維尋優的問題。仿真結果表明,極坐標系下NLS簡化算法對于慢時變系統有較好的收斂速度與穩定性,可用于基于外輻射源單站無源定位求解,但由于沒有目標的狀態轉移信息,多次迭代后會偏離真實值,本文提出的周期初始化的方法有一定校正能力。簡化算法避免了矩陣求逆、向量運算等較復雜算法,對慢時變系統有較好的收斂性和精度,能滿足工程實現的要求。

[1]張芹.固定單站無源定位系統的可行性研究[D].成都:電子科技大學,2007.

ZHANG Qin.The Possibility of Single Fixed Observer Passive Location System[D].Chengdu:University of Electronic Science and Technology of China,2007.(in Chinese)

[2]李宗華,馮道旺,周一宇,等.固定單站對三維運動輻射源無源定位算法及可觀測性分析[J].信號處理,2004,20(2):117-121.

LI Zong-hua,FENG Dao-wang,ZHOU Yi-yu,et al.Passive Location and Tracking Algorithm for 3-D Moving Emitter By Single Non-moving Observer and Its observability Analysis[J].Signal Processing,2004,20(2):117-121.(in Chinese)

[3]付忠.基于單站的無源定位技術研究[D].成都:電子科技大學,2007.

FU Zhong.The Studies on Single-observer Passive Location[D].Chengdu:University of Electronic Science and Technology of China,2007.(in Chinese)

[4]李悅,柳征,伍渝江.基于非線性最小二乘的三星座時差定位綜合算法[J].電訊技術,2010,50(3):27-32.

LI Yue,LIU Zheng,WU Yu-jian.Integration Algorithm of Three-Satellite Constellation TDOA Location Based on NLS[J].Telecommunication Engineering,2010,50(3):27-32.(in Chinese)

[5]廖平,付忠,劉剛.采樣濾波算法在單站無源定位中的應用[J].電訊技術,2006,46(4):30-33.

LIAO Ping,FU Zhong,LIU Gang.Application of Sample Filtering Algorithm in Single Observer Passive Location[J].Telecommunication Engineering,2006,46(4):30-33.(in Chinese)

[6]廖平,楊中海,姜道安.基于概率的單站多目標無源定位算法[J].電訊技術,2006,46(1):47-51.

LIAO Ping,YANG Zhong-hai,JIANG Dao-an.A Single Observer Passive Location Algorithm Based on Probability in Multi-target Environment[J].Telecommunication Engineering,2006,46(1):47-51.(in Chinese)

MEI Yong-bing was born in Xinzhou,Hubei Province,in 1979.He received the M.S.degree from Sichuan University in 2005.He is now an engineer.His research concerns signal processing and electronic warfare.

Email:meiybmail@sina.com

南建設(1954—),男,山東沂南人,研究員,主要從事系統總體設計。

NAN Jian-she was born in Yinan,Shandong Province,in 1954.He is now a senior engineer of professor.His research concerns system design.

A Low-complexity Suboptimal NLS Algorithm and its Application in Single-observer Passive Location

MEI Yong-bing,NAN Jian-she
(Southwest China Institute of Electronic Technology,Chengdu 610036,China)

To simplify the algorithm of single-observer passive location based on outer radiant point,the polar coordinate model is built.The algorithm based on nonlinear least square(NLS)algorithm is presented in polar coordinate.To reduce the computation complex of NLS algorithm,a low-complexity suboptimal algorithm is proposed.The simulation results demonstrate the location precision and stabilization of the suboptimal algorithm,which provides reference for engineering realization.

single-observer passive location;outer radiant;polar coordinate;nonlinear least square(NLS);suboptimal algorithm

TN971

A

10.3969/j.issn.1001-893x.2012.06.016

1001-893X(2012)06-0913-05

2012-04-26;

2012-05-11

梅勇兵(1979—),男,湖北新洲人,2005年于四川大學獲碩士學位,現為工程師,主要研究方向為信號處理與電子對抗;