冯 奇 曲长文 周 强
多运动站异步观测条件下的直接定位算法
冯 奇①曲长文*①周 强②
①(海军航空工程学院电子信息工程系 烟台 264001)②(海军航空工程学院科研部 烟台 264001)
针对多运动站无源融合定位系统中的异步观测问题,该文提出一种基于空间谱分析的直接定位算法。首先利用子空间分析理论建立异步观测直接定位模型,构造多站协同定位联合空间谱函数,采用多级搜索方法实现定位解算。然后推导了该算法的定位误差协方差矩阵及克拉美罗下限。最后对该算法在异步观测条件下的定位性能进行蒙特卡洛仿真。仿真结果表明算法定位精度能够逼近克拉美罗下限,定位性能优于基于定位参数的分步定位算法,尤其在低信噪比情况下优势更明显。
无源定位;异步观测;直接定位;阵列信号处理;脉冲到达角
辐射源被动定位一直是战场信息获取的重要研究方向。从观测站的数量划分,无源定位可分为单站无源定位和多站无源定位。单站无源定位系统战术使用灵活,但面对日益复杂的空间电磁环境,不易及时、准确地获取战场信息。而多站无源定位系统利用站间协同,能够快速、准确地获取辐射源状态信息,已经成为国内外研究的热点[1,2]。
随着电子对抗技术的不断发展,在对窄波束辐射源(如相控阵雷达等)无源定位应用中,可能遇到多个观测站无法同时落入辐射源天线主瓣或者副瓣照射波束的情况。此时各观测站接收到的信号发射时间不同,此时时差(Time Difference Of Arrival, TDOA)、频差(Frequency Difference Of Arrival, FDOA)等定位参数无法估计,相应的定位算法也失去作用。同一观测时间内辐射源信号总是不能被所有观测平台都侦收到,称为异步观测。
针对异步观测条件下的多站无源定位问题,文献[6]提出一种相干定位系统。该系统从脉冲列信号的产生原理出发,还原每个搜索网格点期望的信号到达时间(Time Of Arrival, TOA)序列,通过迭代法逼近实际侦测的TOA序列,最终得到目标位置估计。该方法假设辐射源信号存在一个内在不变的脉冲重复间隔(Underlying Constant Pulse Repetition Interval, UCPRI),而这种重复间隔的存在性和估计方法都有待进一步验证和研究。辐射源的脉冲到达角(Angle Of Arrival, AOA)信息的估计不需要接收信号具有相同的发射时间,可用于实现异步观测条件下的辐射源定位。文献[7-9]提出基于异步AOA的运动目标定位跟踪算法,由于AOA观测量中未包含辐射源速度信息,因此需要通过累积多次观测来估计运动辐射源的位置和速度。文献[10]中提出一种利用异步AOA估计辐射源位置和速度的滤波算法。
上述定位方法都将定位过程分为两步:第1步,通过接收的辐射源信号估计定位参数(如脉冲到达角);第2步,利用这些参数进行定位解算得到辐射源状态信息。两步定位法中定位观测量在各观测站独立估计得到,忽略了定位参数与辐射源之间的约束关系,因此这类方法是次优的[11]。新兴的基于信号层次的直接定位算法(Direct Position Determination, DPD),不需要估计定位参数,直接利用接收信号估计辐射源状态。
本文提出一种基于空间谱分析的直接定位算法,用以解决多运动站异步观测条件下的无源定位问题。建立异步观测条件下的直接定位算法模型,构造多站协同定位的空间谱函数,然后通过多级搜索方法实现辐射源位置估计。最后通过数字仿真对比了本文算法与文献[9]中基于异步AOA的分步定位算法的定位性能并分析各参数对定位精度的影响。
图1 多运动站异步观测定位场景示意图
图1中实心表示该观测站能够接收到辐射源信号;空心表示该观测站不能接收到辐射源信号。在该假设场景中由于各观测站接收到信号的发射时间不同,文献[11-14]中提出的直接定位法也无法实现定位。若图中直线表示描述波达方向的空间谱,则谱线相交处即为目标估计位置。
(2)
条件3:为了便于表述,假设观测站只有两种状态:能接收到所有辐射源信号;不能接收到任一辐射源信号。
条件3只是为了方便下文表述,不影响本文算法的一般性。异步观测定位问题就是利用各观测平台侦收到的异步观测信号和观测平台位置信息估计各辐射源位置。
3.1 基于空间谱分析的直接定位算法
(5)
结合式(4)可得
(7)
因此
进而可得
(9)
(11)
(13)
易得空间谱函数是关于波达方向的函数,为了实现目标定位可以利用式(2)将空间谱函数转化为目标位置的函数,则目标位置估计值为
从算法实现原理来看本文算法也可以称为异步观测条件下测向交叉直接定位算法,该算法具有自身的特点。
(1)与文献[9]中基于AOA的异步观测条件下定位算法相比,本文算法不需要估计AOA参数,属于一种直接定位算法。
(2)与文献[11-14]中直接定位法相比不要求观测信号具有相同的发射时间,即对静止或慢速运动目标本文算法能实现异步观测条件下的无源定位。
3.2 定位误差分析
克拉美罗下限(Cramer-Rao Lower Bound, CRLB)是无偏估计器所能达到的估计误差下限。设观测站接收到的信号表示为,辐射源位置矩阵,则关于辐射源位置的似然函数表示为
(16)
(17)
式中
当样本数据足够多时,式(17)可以简化为
(19)
利用式(14)多站协同定位空间谱函数作为代价函数估计辐射源位置坐标,推导定位协方差矩阵过程中会出现不同观测站噪声向量在信号空间的正交投影之间的互相关。利用条件1和条件2可得不同观测站正交投影互相关为,即
(21)
(23)
(24)
由式(24)可得,整理后定位误差方差求解不再需要特征分解运算。同理可得轴方向定位误差方差。
本节通过数字仿真对比本文算法与文献[9]中基于AOA的分步定位算法性能并分析各参数对定位精度的影响。仿真定位场景如图1,观测站初始位置坐标为,和,观测站速度为,且观测时间内不变。两个辐射源真实位置坐标分别为和,辐射源均静止,假设辐射源信号为固定重频的脉冲列信号,脉冲重频分别为和,脉冲幅度为单位1,脉冲占空比为1%,载波频率可测且分别为1 GHz和3 GHz。如图1所示,时刻观测站1和观测站2能接收到辐射源信号,时刻只有观测站3能够接收到辐射源信号,时间间隔。传感器天线阵列阵元个数,阵元间距,数据样本个数。为方便分析假设辐射源信号能量相等,信号之间互不相关(仿真中只给出了对辐射源1的定位性能分析,辐射源2具有相同的性质)。
(1)算法定位性能对比分析: 仿真对比本文提出的基于空间谱分析的直接定位算法与基于定位参数的分步定位算法定位性能,如图2所示。图2中“DPD_RMSE”表示本文算法均方根误差(Root Mean Square Error, RMSE),“Two-Steps”表示分步定位算法均方根误差。蒙特卡洛仿真次数为1000次。
图2 算法定位性能随信噪比的变化
由图2得,直接定位算法定位性能优于基于定位参数的分步定位算法,尤其在低信噪比情况下优势更明显。当信噪比为4 dB时,直接定位算法的定位精度能够达到CRLB,与分步定位算法相比,定位精度提高了约0.2 km;当信噪比为-10 dB时,与分步定位算法相比,直接定位算法定位精度提高了约1 km。
(2)阵元个数和样本数对算法精度的影响: 定义算法有效性参数为均方根误差与克拉美-罗下限的比值。
由图3(a)得到:阵元越多,样本数越大,定位性能越好;由图3(b)得到:与样本数相比,阵元个数对有效性参数的影响比较强。即样本数相同情况下阵元个数越多,RMSE越逼近CRLB,反之,阵元个数相同情况下,增加样本数对有效性参数影响不大。
(3)阵元间距对定位性能的影响: 利用空间谱估计测向时,要求阵元间距应满足“半波长”条件,否则相位差可能大于,而产生方向模糊,即和都可以是信号的波达方向。如果能够确定目标是位于天线阵列的“前面”还是“后面”,则可以取消“半波长”条件的限制。当阵元间距分别取和时,辐射源位置估计的RMSE和CRLB对比图如图4所示。
由图4得,不考虑“半波长”条件限制情况下,阵元间距越大辐射源位置估计精度越高,且定位精度的提升速率随着阵元间距的增加而减小。
本文提出一种基于空间谱分析的直接定位算法,用以解决多运动站异步观测情况下的无源定位问题。建立异步观测条件下的直接定位算法模型,推导了该算法的定位误差协方差矩阵和Cramer- Rao下限。通过数字仿真对比分析了本文算法与分步定位算法的定位性能,分析了不同参数对定位精度的影响。仿真结果表明异步观测条件下对静止或慢速运动目标,直接定位算法性能优于基于AOA的分步定位算法,尤其在低信噪比情况下优势更明显;阵元个数越多、样本数越大定位精度越高;在不产生方向模糊情况下,阵元间距越大,定位精度越高。
图3 阵元个数和样本数对定位性能的影响 图4 阵元间距对定位性能的影响
[1] 徐征. 机载多平台无源定位跟踪算法研究[D]. [博士论文],海军航空工程学院, 2014.
XU Zheng. Research on algorithms for airborne multi- platform passive localization and tracking[D]. [Ph.D. dissertation], Naval Aeronautical and Astronautical University, 2014.
[2] 雷文英, 陈伯孝, 杨明磊, 等. 基于TOA 和TDOA的三维无源目标定位方法[J]. 系统工程与电子技术, 2014, 36(5): 816-823. doi: 10.3969/j.issn.1001-506X.2014.05.02.
LEI Wenying, CHEN Boxiao, YANG Minglei,. Passive 3D target location method based on TOA and TDOA[J]., 2014, 36(5): 816-823. doi: 10.3969/j.issn.1001-506X.2014.05.02.
[3] CAO Yalu, PENG Li, LI Jinzhou,. A new iterative algorithm for geolocating a known altitude target using TDOA and FODA measurements in the presence of satellite location uncertainty[J]., 2015, 28(5): 1510-1518. doi: 10.1016/j.cja.2015.08.015.
[4] LI Jing, ZHAO Yongjun, LI Donghai,. Accurate single-observer passive coherent location estimation based on TDOA and DOA[J]., 2014, 27(4): 913-923. doi: 10.1016/j.cja.2014.06.004.
[5] LI Jinzhou, GUO Fucheng, and JIANG Wenli. Source localization and calibration using TDOA and FDOA measurements in the presence of sensor location uncertainty [J]., 2014, 57: 042315(1)-042315(12). doi: 10.1007/s11432-013-4800-2.
[6] SCHIFFMILIER R, ADLER H, and CARROLL M. Multi- ship coherent geolocation system[P]. USA Patent, 124385, 2005-12-29.
[7] MUSICKI D. Multi-target tracking using multiple passive bearings-only asynchronous sensors[J]., 2008, 44(3): 1151-1160. doi: 10.1109/TAES.2008.4655370.
[8] MUSICKI D. Bearings only multi-sensor maneuvering target tracking[J].&, 2008, 57(3): 216-221. doi: 10.1016/j.sysconle.2007.08.008.
[9] SATHYAN T and SINHA A. A two-stage assignment-based algorithm for asynchronous multisensor bearings-only tracking[J]., 2011, 47(3): 2153-2167. doi: 10.1109/TAES.2011. 5937289.
[10] HANSELMANN T and MORELANDE M. Multiple target tracking with asynchronous bearings only measurements[C]. Proceeding of the International Conference on Information Fusion, Quebec, Canada, 2007. doi: 10.1109/ICIF.2007. 4408056.
[11] AMAR A and WEISS A I. Localization of narrowband radio emitters based on doppler frequency shifts[J]., 2008, 56(11): 5500-5508. doi: 10.1109/TSP.2008.929655.
[12] AMAR A and WEISS A I. Direct geolocation of stationary wideband radio signal based on time delays and Doppler shifts[C]. IEEE/SP Workshop on Statistical Signal Processing, Cardiff, UK, 2009: 101-104. doi: 10.1109/SSP. 2009.5278630.
[13] OISPUU M. Direct detection and position determination of multiple sources with intermittent emission[J]., 2010, 90(12): 3056-3064. doi: 10.1016/j.sigpro. 2010.05.010.
[14] 王云龙, 吴瑛. 联合时延与多普勒频率的直接定位改进算法[J]. 西安交通大学学报, 2015, 49(4): 123-129. doi: 10.7652/ xjtuxb201504020.
WANG Yunlong and WU Ying. An improved direct determination algorithm with combined time delay and doppler[J]., 2015, 49(4): 123-129. doi: 10.7652/ xjtuxb201504020.
[15] 张贤达. 现代信号处理[M]. 北京: 清华大学出版社, 2002: 126-135.
ZHANG Xianda. Modern Signal Processing[M]. Beijing: Tsinghua University Press, 2002: 126-135.
[16] STOICA P and NEHORAI A. Music, maximum likelihood, and Cramer-Rao bound: Further results and comparisons[J].,,, 1990, 38(12): 2140-2150. doi: 10.1109/29.61541.
Direct Position Determination Using Asynchronous Observations of Multiple Moving Sensors
FENG Qi①QU Changwen①ZHOU Qiang②
①(,,264001,)②(,,264001,)
A Direct Position Determination (DPD) algorithm based on spatial spectrum analysis for the issue of asynchronous observation in the multiple moving sensors passive fusion location system is proposed. At first, a direct determination mathematical model of the asynchronous scene is constructed and a multiple sensors cooperative spatial spectrum function is composed. Furthermore, the position of the emitter can be got by a two or three dimensional search, and then the derivations of the estimators’ variance and the Cramer-Rao Lower Bound (CRLB) are given. Finally, the Monte Carlo simulations over the asynchronous scene indicate that the accuracy of the proposed method is close to the Cramer-Rao lower bound and superior to the two-step location method based on the location parameters, especially in the low SNR.
Passive location; Asynchronous observation; Direct Position Determination (DPD); Array signal processing; Angle of arrival
TN958.97
A
1009-5896(2017)02-0417-06
10.11999/JEIT160314
2016-04-01;改回日期:2016-09-01;
2016-10-17
曲长文 qcwwby@sohu.com
冯 奇: 男,1990年生,博士生,研究方向为无源定位跟踪技术.
曲长文: 男,1963年生,教授,博士生导师,研究方向为无源定位跟踪技术、信号与信息处理、SAR目标检测与识别技术.
周 强: 男,1984年生,助理研究员,研究方向为SAR目标检测与识别技术、无源定位跟踪技术.