摘要
前导脉冲信号到达时间(Time of arrival,TOA)估计精确度直接影响多点定位精确度。在差分匹配滤波器数学模型基础上,分析理想脉冲信号差分匹配滤波器TOA精确度与采样频率的关系。讨论非理想脉冲信号之间的时间结构关系,提出基于S模式前导四脉冲信号构造差分匹配滤波器校准时间戳的改进算法,实现过零点的准确估计。从信噪比、采样频率和复杂度仿真数据输出的均方根误差对比可见,通过四脉冲滤波获得的TOA精确度要比单脉冲和双脉冲滤波获得TOA精确度高。匹配滤波器和差分匹配滤波器TOA精确度对比分析结果表明,修正的差分匹配滤波器算法有更好的效果。在采样频率 40 MHz、信噪比为40 dB时,四脉冲信号差分匹配滤波器TOA精确度为0.402 7 μs,远远小于国际民航组织规定的时间误差精度最大值25 μs。最后实际空管运行数据验证结果证明了本文修正算法估算TOA的精确度和可用性。
S模式前导脉冲信号到达时间估计精确度是多点定位准确实现机场场面监视的基础和前提。机场区域和航路上的飞机受噪声等因素的影响,针对S模式前导脉冲信号到达时间的精确估计存在诸多难题。Torrieri等提出通过设定自适应门限检测信号到达时间的方
针对上述问题,利用前导脉冲信号之间的时间结构关系,提出用四脉冲信号差分匹配滤波器过零点来标记时间戳的改进算法,通过匹配滤波器输出信号峰值对应的采样点来限制差分匹配滤波器输出信号过零点的范围,实现较低信噪比时TOA的测量。对比分析双脉冲和四脉冲信号使用该算法的TOA精确度,与第1个脉冲的TOA精确度作对比,结果表明四脉冲信号的TOA估计精确度要优于单脉冲和双脉冲TOA估计精确度。
S模式信号有固定的格式,如
(1) |
式中:T为脉冲重复周期;p为S模式脉冲序列,p=[1,0,1,0,0,0,0,1,0,1,0,0,0,0,0,0]。

图1 S模式信号脉冲串
Fig.1 S-mode signal pulse train
利用前导脉冲信号上升沿所用时间和下降沿所用时间的过渡过程,展开对单脉冲、双脉冲和四脉冲信号差分匹配滤波时间戳标记精确度的讨论。机场附近和航路上飞机的S模式信号信噪比一般在15 dB以
首先将S模式前导单脉冲、双脉冲和四脉冲信号分别输入对应的匹配滤波器和差分匹配滤波器。然后提取出匹配滤波器输出函数的峰值,将峰值对应的采样点j保存下来。提取差分匹配滤波器输出函数的第j位,在第j位输出信号so( j)和第j+1位输出信号so( j+1)间循环进行过零点检测,如果so( j)>0,so( j+1)<0,则提取出这两个采样点j,j+1及其对应的信号输出值so( j),so( j+1),然后在点( j,so( j))和点( j+1,so( j+1))间做直线方程。最后令直线方程等于零,求取过零点,标记时间戳。具体的算法流程图如

图2 差分匹配滤波器算法流程图
Fig.2 Flow chart of differential matched filter algorithm
如果输入脉冲信号为s(n),则s(n)可表示为
(2) |
式中:N为信号长度;为信号幅值;表示冲激响应序列。由匹配滤波器可知,在n时刻对应的传输函数可表示
(3) |
在n-1时刻对应的传输函数可表示为
(4) |
则n时刻差分匹配滤波器的传输函数可表示为
(5) |
式中:,。
对应的匹配滤波器的输出信号so(n)为
(6) |
对应的差分匹配滤波器的输出信号sod(n)为
(7) |
如果输入信号为单脉冲信号,则差分匹配滤波器的传输函数中,。
差分匹配滤波器是匹配滤波器的传输函数作差分运算,将输出信噪比最大的点转化为零点。根据国际民航组织标准,在无干扰条件下,滤波器的脉冲响应为理想脉冲信号的导数,此时到达时间为滤波器的过零点对应时
仿真条件如下:无噪声,采样频率为40 MHz,前导单脉冲的上升沿所用时间为0.1 μs,下降沿所用时间为当0.125 μs。差分匹配滤波器的传输函数波形如

图3 单脉冲差分匹配滤波器传输函数波形
Fig.3 Transfer function waveform of monopulse differential matched filter
假设理想脉冲信号的宽度为,采样频率为,采样间隔,采样点数。
匹配滤波器的输出信号为
(9) |
式中Rs为自相关函数。
由于自相关偶函数特性
(10) |
差分匹配滤波的输出为
(11) |
根据自相关偶函数特性
(12) |
当n=N时,式中的前者取最大值,此时的差分匹配滤波输出为
(13) |
当n=N+1时,
(14) |
所以说差分匹配滤波器过零点的值在N和N+1之间。求(N, Rs(0)-Rs(1))与(N+1, Rs(1)-Rs(0))间的直线方程,该直线过零点时刻即为所求差分匹配滤波的到达时间。
直线方程为
(15) |
令y=0,过零点x
(16) |
则差分匹配滤波的输出均方根误差为
(17) |
如
仿真条件如下:信噪比15 dB,不同频率,理想前导脉冲信号通过差分匹配滤波器的输出均方根误差如

图4 理想前导脉冲差分匹配滤波输出均方根误差
Fig.4 RMSE of ideal preamble pulse differential matched filter
仿真条件如下:信噪比15 dB,采样频率40 MHz和50 MHz,脉冲的上升沿所用时间为t1,脉冲的下降沿所用时间为t2,t1=t2的脉冲信号通过各自对应的差分匹配滤波器均方根误差如

图5 t1=t2的脉冲信号差分匹配滤波的均方根误差
Fig.5 RMSE of pulse signal differential matched filtering at t1=t2
仿真条件如下:脉冲宽度为T,采样频率为fs1,离散化后脉冲信号的采样个数为Nn=T×fs1,上升沿采样个数为N1=t1×fs1,下降沿采样个数为N2=t2×fs1。无噪声,采样频率为40 MHz和50 MHz,t1=0.05 μs、t2不定,t1=0.1 μs、t2不定的脉冲信号,通过各自匹配滤波器和差分匹配滤波器的输出均方根误差如

图6 不同频率不同时间脉冲信号TOA测量精确度对比图
Fig.6 Comparison of TOA measurement accuracy of pulse signals at different frequencies and different time
由
首先对模型性能进行仿真分析,本节将从TOA测量均方根误差的角度分析单脉冲、双脉冲和四脉冲信号的时间戳测量精确度。
脉冲积累是信号处理中有效的能量积累方法,随着脉冲个数的增加,相较于单脉冲,四脉冲信号积累的能量更多,又因脉冲信号的噪声之间互不相关,所以四脉冲信号的信噪比比单脉冲的信噪比大,四脉冲匹配滤波输出的峰值点和差分匹配滤波输出的过零点比单脉冲的精确,进而可以降低输出均方根误差,提高TOA精确度。
仿真条件如下:信噪比15 dB,采样频率40 MHz,t1= 0.1 μs,t2= 0.125 μs的脉冲信号,采用本文设计的差分匹配滤波和匹配滤波模型,经过1 000次仿真,得到单、双和四脉冲的匹配滤波和差分匹配滤波均方根误差对比关系,如
仿真条件如下:采样频率40 MHz,t1=0.1 μs,t2=0.125 μs的脉冲信号,采用本文设计的差分匹配滤波和匹配滤波模型,经过1 000次仿真,得到单脉冲、双脉冲和四脉冲信号的均方根误差随信噪比变化的曲线,如

图7 差分匹配滤波器与匹配滤波器输出误差随信噪比变化的曲线
Fig.7 Error curves of differential matched filter and matched filter with SNR
仿真条件如下:信噪比15 dB,t1= 0.1 μs,t2= 0.125 μs的脉冲信号,采用本文设计的差分匹配滤波器和匹配滤波器模型,经过1 000次仿真,得到单脉冲、双脉冲和四脉冲信号均方根误差随采样频率变化的曲线,如

图8 差分匹配滤波与匹配滤波输出误差随频率变化曲线
Fig.8 Output error curves of differential matched filter and matched filter with frequency
仿真条件如下:信噪比20 dB,采样频率40 MHz,t1=0.1 μs,t2= 0.125 μs的单脉冲、双脉冲和四脉冲信号,采用本文设计的差分匹配滤波器,1 000次运行后仿真所需时间如
在天津机场终端区设置1个广播式自动相关监视接收机,实时接收的真实数据如



图9 真实数据
Fig.9 Real data
如
本文研究了S模式前导脉冲信号的过渡过程与时间戳测量精确度问题。在匹配滤波和差分匹配滤波基本原理的基础上,利用脉冲间的时间结构关系,提出用四脉冲信号差分匹配滤波过零点来标记信号TOA时间戳的改进算法,实验结果表明四脉冲信号的时间戳标记精确度优于单脉冲信号。与匹配滤波的时间戳标记精确度作对比,差分匹配滤波的时间戳标记精确度高于匹配滤波。在采样频率40 MHz、信噪比40 dB时,四脉冲信号差分匹配滤波器的时间戳标记精确度达到0.402 7 μs。在天津机场终端区接收实时信号的采样频率为80 MHz,信噪比约为15.12 dB,四脉冲信号差分匹配滤波器的时间戳标记精确度达到了8.948 75 μs,远小于国际民航组织规定的时间误差精度最大值25 μ
参考文献
Torrieri D J. Adaptive thresholding systems[J]. IEEE Transactions on Aerospace and Electronic Systems, 1977, AES-13(3): 273-280. [百度学术]
Zehavi E. Estimation of time of arrival for rectangular pulses[J]. IEEE Transactions on Aerospace and Electronic Systems, 1984, AES-20(6): 742-747. [百度学术]
Chan Y T, Lee B H, Inkol R, et al. Estimation of pulse parameters by convolution[C]/ [百度学术]
/Proceedings of Electrical and Computer Engineering, CCECE’06 Canadian Conference on.[S.l.]: IEEE, 2006: 17-20. [百度学术]
Chan Y T, Lee B H, Inkol R, et al. Estimation of pulse parameters by autoconvolution and least squares[J]. IEEE Transactions on Aerospace and Electronic Systems, 2010, 46(1): 363-374. [百度学术]
席轶敏,刘渝,靖晟.电子侦察信号实时检测算法及性能分析[J]. 南京航空航天大学学报,2001, 33(3): 277-281. [百度学术]
XI Yimin, LIU Yu, JING Sheng. Analysis on probability of false alarm and detection[J]. Journal of Nanjing University of Aeronautics & Astronautics, 2001, 33(3): 277-281. [百度学术]
胡国兵,刘渝,邓振淼.基于Haar小波变换的信号到达时间估计[J].系统工程与电子技术,2009,31(7): 1615-1619. [百度学术]
HU Guobing, LIU Yu, DENG Zhenmiao. Arrival time estimation of signals based on Haar wavelets transform[J].Systems Engineering and Electronics,2009,31(7): 1615-1619. [百度学术]
GASPARE G, MAURO L, Patrizio D M, et al. New time of arrival estimation method for multilateration target location[C]//Proceedings of Joint International Symposium on Sensors and Systems for Airport Surveillance. Paris, France:[s.n.], 2005: 116-119. [百度学术]
Lin W C, Liu K C, Wang C K. Differential matched filter architecture for spread spectrum communication systems[J]. Electronics Letters, 1996, 32(17): 1539. [百度学术]
Petrochilos N, Galati G, Piracci E. Separation of SSR signals by array processing in multilateration systems[J]. IEEE Transactions on Aerospace and Electronic Systems, 2009, 45(3): 965-982. [百度学术]
Galati G, Leonardi M, Tosti M. Multilateration (local and wide area) as a distributed sensor system: Lower bounds of accuracy[C]//Proceedings of European Radar Conference.[S.l.]:IEEE, 2008. [百度学术]
王洪,金尔文,刘昌忠,等. 多点定位TOA精确估计及同步误差校正算法[J]. 系统工程与电子技术, 2013, 35(4): 835-839. [百度学术]
WANG Hong, JIN Erwen, LIU Changzhong, et al, Accurate estimation of TOA and calibration of synchronization error for multilateration[J]. Systems Engineering and Electronics, 2013, 35(4): 835-839. [百度学术]
徐江, 李艳, 纪晶华.统计学基础[M]. 北京:机械工业出版社, 2015. [百度学术]
XU Jiang, LI Yan, JI Jinghua. Statistical basis[M]. Beijing: Mechanical Industry Press,2015. [百度学术]
Eurocae. Minimum operational performance specification for mode S multilateration system for use in advanced surface movement guidance and control system(A-SMCGS): ED-117[S]. [S.l.]: Eurocae, 2003: 36-46. [百度学术]
Galati G, Leonardi M, Magaro P, et al. Wide area surveillance using SSR mode S multilateration: Advantages and limitations[C]//Proceedings of European Radar Conference, EURAD 2005. Paris: [s.n.], 2005: 225-229. [百度学术]