网刊加载中。。。

使用Chrome浏览器效果最佳,继续浏览,你可能不会看到最佳的展示效果,

确定继续浏览么?

复制成功,请在其他浏览器进行阅读

S模式前导脉冲过零点提取与TOA时间戳精确度研究  PDF

  • 宫峰勋
  • 曹雅茹
中国民航大学电子信息与自动化学院,天津,300300

中图分类号: V351.37TN959.6

最近更新:2020-11-06

DOI:10.16356/j.1005-2615.2020.05.004

  • 全文
  • 图表
  • 参考文献
  • 作者
  • 出版信息
目录contents

摘要

前导脉冲信号到达时间(Time of arrival,TOA)估计精确度直接影响多点定位精确度。在差分匹配滤波器数学模型基础上,分析理想脉冲信号差分匹配滤波器TOA精确度与采样频率的关系。讨论非理想脉冲信号之间的时间结构关系,提出基于S模式前导四脉冲信号构造差分匹配滤波器校准时间戳的改进算法,实现过零点的准确估计。从信噪比、采样频率和复杂度仿真数据输出的均方根误差对比可见,通过四脉冲滤波获得的TOA精确度要比单脉冲和双脉冲滤波获得TOA精确度高。匹配滤波器和差分匹配滤波器TOA精确度对比分析结果表明,修正的差分匹配滤波器算法有更好的效果。在采样频率 40 MHz、信噪比为40 dB时,四脉冲信号差分匹配滤波器TOA精确度为0.402 7 μs,远远小于国际民航组织规定的时间误差精度最大值25 μs。最后实际空管运行数据验证结果证明了本文修正算法估算TOA的精确度和可用性。

S模式前导脉冲信号到达时间估计精确度是多点定位准确实现机场场面监视的基础和前提。机场区域和航路上的飞机受噪声等因素的影响,针对S模式前导脉冲信号到达时间的精确估计存在诸多难题。Torrieri等提出通过设定自适应门限检测信号到达时间的方[

1⁃2],该方法可以克服脉冲幅度起伏对到达时间的影响,简单、易于硬件实现,但是在低信噪比情况下,信号容易埋没在噪声中,造成较大的测时误差。Chan等提出了自卷积算[3⁃4],通过脉冲信号的自动卷积确定卷积峰值位置,然后使用最小二乘法进一步细化,估计脉冲信号的到达时间(Time of arrival,TOA),但是该算法估计精度较低,对信噪比的要求较高。席轶敏等提出自相关检测算[5],采用双门限判决TOA的方法实时性好,但是涉及复数计算,运算过程比较复杂。胡国兵等将自相关检测算法和Haar小波变换相结合来估计信号的TOA[6],无任何先验条件下,信噪比较低时,实现对TOA的精度估计,但是小波变换尺度较小时,抑制噪声的能力减弱,易造成TOA的误判,形成较大的测时误差。文献[7⁃10]提出在匹配滤波后加差分的TOA估计方法,它的阈值为固定的零点,避免由门限带来的测量误差,但其仅讨论第1个脉冲在高信噪比时TOA精确度的测量。

针对上述问题,利用前导脉冲信号之间的时间结构关系,提出用四脉冲信号差分匹配滤波器过零点来标记时间戳的改进算法,通过匹配滤波器输出信号峰值对应的采样点来限制差分匹配滤波器输出信号过零点的范围,实现较低信噪比时TOA的测量。对比分析双脉冲和四脉冲信号使用该算法的TOA精确度,与第1个脉冲的TOA精确度作对比,结果表明四脉冲信号的TOA估计精确度要优于单脉冲和双脉冲TOA估计精确度。

1 模型分析

S模式信号有固定的格式,如图1所示。S模式信号的数学模型[

11]

b(n)=i=1128/240b(i)p(n-iT) (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以[

11],所以文中设置的信噪比为15 dB,噪声类型为高斯白噪声。

1.1 时间戳标记流程图

首先将S模式前导单脉冲、双脉冲和四脉冲信号分别输入对应的匹配滤波器和差分匹配滤波器。然后提取出匹配滤波器输出函数的峰值,将峰值对应的采样点j保存下来。提取差分匹配滤波器输出函数的第j位,在第j位输出信号so( j)和第j+1位输出信号so( j+1)间循环进行过零点检测,如果so( j)>0,so( j+1)<0,则提取出这两个采样点jj+1及其对应的信号输出值so( j),so( j+1),然后在点( j,so( j))和点( j+1,so( j+1))间做直线方程。最后令直线方程等于零,求取过零点,标记时间戳。具体的算法流程图如图2所示。

图2 差分匹配滤波器算法流程图

Fig.2 Flow chart of differential matched filter algorithm

1.2 差分匹配滤波器模型推导

如果输入脉冲信号为s(n),则s(n)可表示为

s(n)=m=0N-1amδ(n-m) (2)

式中:N为信号长度;am为信号幅值;δ(n)表示冲激响应序列。由匹配滤波器可知,在n时刻对应的传输函数可表示

h(n)=m=0N-1amδ(N-n-m) (3)

n-1时刻对应的传输函数可表示为

h(n-1)=m=0N-1amδ(N+1-n-m) (4)

n时刻差分匹配滤波器的传输函数可表示为

hd(n)=h(n)-h(n-1)=m=0Ncmδ(N-n-m) (5)

式中:cN=aN-1, c0=a0, cp=ap-1-ap1pN-1

对应的匹配滤波器的输出信号so(n)为

so(n)=s(n)×h(n) (6)

对应的差分匹配滤波器的输出信号sod(n)为

sod(n)=s(n)×hd(n) (7)

如果输入信号为单脉冲信号,则差分匹配滤波器的传输函数中,cp=ap-1-ap,2pN+1

差分匹配滤波器是匹配滤波器的传输函数作差分运算,将输出信噪比最大的点转化为零点。根据国际民航组织标准,在无干扰条件下,滤波器的脉冲响应为理想脉冲信号的导数,此时到达时间为滤波器的过零点对应时 [

7],即在过零点时刻标记时间戳。

仿真条件如下:无噪声,采样频率为40 MHz,前导单脉冲的上升沿所用时间为0.1 μs,下降沿所用时间为当0.125 μs。差分匹配滤波器的传输函数波形如图3所示。

图3 单脉冲差分匹配滤波器传输函数波形

Fig.3 Transfer function waveform of monopulse differential matched filter

1.3 时间戳精确度衡量公式

设含噪声的前导脉冲信号的采样间隔为T,采样点数为N,无噪脉冲信号经过差分匹配滤波器(匹配滤波器)的输出信号时间戳标记在第N0个采样点,含噪声的脉冲信号经过差分匹配滤波器(匹配滤波器)的输出信号时间戳标记在第Nq位采样点,则时间戳标记的均方根误差(Root mean square error, RMSE)表示[

12]

RMSE=1nq=1nσTq2 (8)

式中:σTq=|Nq-N0|×T,为第q次的TOA测量误差。

2 脉冲信号过零点与时间戳精确度

2.1 理想脉冲信号

假设理想脉冲信号的宽度为τ,采样频率为fs,采样间隔Ts=1fs,采样点数N=τTs=τ×fs

匹配滤波器的输出信号为

so(n)=s(n)×h(n)=Rs(N-n) (9)

式中Rs为自相关函数。

由于自相关偶函数特性

so(n)=Rs(n-N) (10)

差分匹配滤波的输出为

sod(n)=s(n)×hd(n)=Rs(N-n)-Rs(N+1-n) (11)

根据自相关偶函数特性

sod(n)=Rs(n-N)-Rs(n-(N+1)) (12)

n=N时,式中的前者取最大值,此时的差分匹配滤波输出为

sod(N)=Rs(0)-Rs(1) (13)

n=N+1时,式(12)中的后者取最大值,此时的差分匹配滤波输出为

sod(N+1)=Rs(1)-Rs(0)=-sod(N) (14)

所以说差分匹配滤波器过零点的值在NN+1之间。求(N, Rs(0)-Rs(1))与(N+1, Rs(1)-Rs(0))间的直线方程,该直线过零点时刻即为所求差分匹配滤波的到达时间。

直线方程为

y=2[Rs(1)-Rs(0)]x+[Rs(0)-Rs(1)](1+2N)=[Rs(1)-Rs(0)][2x-(1+2N)] (15)

y=0,过零点x

x=1+2N2=12+τfs (16)

则差分匹配滤波的输出均方根误差为

Δt=(x-N)×Ts=x-Nfs=12fs (17)

式(16)所示,理想前导脉冲信号差分匹配滤波器的过零点与脉宽和采样频率有关。如式(17)所示,差分匹配滤波器的时间戳标记精确度只与信号的采样频率有关,成1/2反比关系。

仿真条件如下:信噪比15 dB,不同频率,理想前导脉冲信号通过差分匹配滤波器的输出均方根误差如图4所示,理想脉冲信号通过差分匹配滤波器的输出均方根误差只与采样频率有关,成1/2反比关系。

图4 理想前导脉冲差分匹配滤波输出均方根误差

Fig.4 RMSE of ideal preamble pulse differential matched filter

2.2 非理想脉冲信号过渡过程分析

2.2.1 过渡时间相等

仿真条件如下:信噪比15 dB,采样频率40 MHz和50 MHz,脉冲的上升沿所用时间为t1,脉冲的下降沿所用时间为t2t1=t2的脉冲信号通过各自对应的差分匹配滤波器均方根误差如图5(a)所示, t1=t2的脉冲信号在不同频率下的差分匹配滤波器均方根误差如图5(b)所示。过渡时间相等的脉冲信号通过差分匹配滤波器的输出误差只与采样频率有关,且成反比关系。随着采样频率的增加,输出均方根误差逐渐减小,TOA测量精确度逐渐提高。

图5 t1=t2的脉冲信号差分匹配滤波的均方根误差

Fig.5 RMSE of pulse signal differential matched filtering at t1=t2

2.2.2 过渡时间不等

仿真条件如下:脉冲宽度为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所示。

图6 不同频率不同时间脉冲信号TOA测量精确度对比图

Fig.6 Comparison of TOA measurement accuracy of pulse signals at different frequencies and different time

图6中,Da代表采样频率50 MHz,上升沿所用时间0.05 μs的单脉冲差分匹配滤波均方根误差;Ma代表采样频率50 MHz,上升沿所用时间0.05 μs的单脉冲匹配滤波均方根误差。Db代表采样频率40 MHz,上升沿所用时间0.05 μs的单脉冲差分匹配滤波均方根误差;Mb代表采样频率40 MHz,上升沿所用时间0.05 μs的单脉冲匹配滤波均方根误差。Dc代表采样频率50 MHz,上升沿所用时间0.1 μs的单脉冲差分匹配滤波均方根误差;Mc代表采样频率50 MHz,上升沿所用时间0.1 μs的单脉冲匹配滤波均方根误差。Dd代表采样频率40 MHz,上升沿所用时间0.1 μs的单脉冲差分匹配滤波均方根误差;Md代表采样频率40 MHz,上升沿所用时间0.1 μs的单脉冲匹配滤波均方根误差。

图6可看出,在不含噪声,过渡时间不等的情况下,当t2=t1+1fs1,即N2=N1+1时,过零点恰好为Nn,差分匹配滤波的输出均方根误差为0。以采样频率为50 MHz,t1=0.1 μs的信号为例,当t2=0.12 μs时,滤波器的输出均方根误差为0。

3 过零点提取性能分析

首先对模型性能进行仿真分析,本节将从TOA测量均方根误差的角度分析单脉冲、双脉冲和四脉冲信号的时间戳测量精确度。

脉冲积累是信号处理中有效的能量积累方法,随着脉冲个数的增加,相较于单脉冲,四脉冲信号积累的能量更多,又因脉冲信号的噪声之间互不相关,所以四脉冲信号的信噪比比单脉冲的信噪比大,四脉冲匹配滤波输出的峰值点和差分匹配滤波输出的过零点比单脉冲的精确,进而可以降低输出均方根误差,提高TOA精确度。

仿真条件如下:信噪比15 dB,采样频率40 MHz,t1= 0.1 μs,t2= 0.125 μs的脉冲信号,采用本文设计的差分匹配滤波和匹配滤波模型,经过1 000次仿真,得到单、双和四脉冲的匹配滤波和差分匹配滤波均方根误差对比关系,如表1所示。由单脉冲、双脉冲和四脉冲的顺序来看,信号的能量逐渐增强,滤波器的输出误差逐渐减小,四脉冲信号的输出均方根误差最小,时间戳的标记精确度最高。不管是单脉冲、双脉冲还是四脉冲信号,差分匹配滤波输出误差相较于匹配滤波更小,时间戳标记精确度更高。

表 1 匹配滤波器和差分匹配滤波器均方根误差
Table 1 RMSE of matched filter and differential matched filter

脉冲

类型

差分匹配滤波RMSE/μs匹配滤波RMSE/μs
SNR=10 dBSNR=15 dBSNR=10 dBSNR=15 dB
单脉冲 19.672 7 14.598 9 37.864 9 23.887 8
双脉冲 19.078 3 12.240 4 30.669 6 20.811 7
四脉冲 18.673 3 10.002 3 25.775 5 19.397 2

3.1 信噪比

仿真条件如下:采样频率40 MHz,t1=0.1 μs,t2=0.125 μs的脉冲信号,采用本文设计的差分匹配滤波和匹配滤波模型,经过1 000次仿真,得到单脉冲、双脉冲和四脉冲信号的均方根误差随信噪比变化的曲线,如图7所示。图7中,Done代表单脉冲差分匹配滤波均方根误差随信噪比变化曲线,Mone代表单脉冲匹配滤波均方根误差随信噪比变化曲线。Dtwo代表双脉冲差分匹配滤波均方根误差随信噪比变化曲线,Mtwo代表双脉冲匹配滤波均方根误差随信噪比变化曲线。Dfour代表四脉冲差分匹配滤波均方根误差随信噪比变化曲线,Mfour代表四脉冲匹配滤波均方根误差随信噪比变化曲线。由图7可见,输入信噪比越高,滤波器输出信号的均方根误差越小,时间戳的标记精确度就越高;但是在输入信噪比大于20 dB时,匹配滤波器输出均方根误差基本不发生变化,而差分匹配滤波器输出误差仍一直呈减小趋势,时间戳标记精确度也更高。同一信噪比时,与单脉冲、双脉冲信号相比,四脉冲信号滤波器的输出均方根误差最小,时间戳标记精确度最高,如表2所示。

图7 差分匹配滤波器与匹配滤波器输出误差随信噪比变化的曲线

Fig.7 Error curves of differential matched filter and matched filter with SNR

表 2 信噪比40 dB时输出均方根误差
Table 2 Output RMSE at SNR of 40 dB
脉冲类型差分匹配滤波RMSE/μs匹配滤波RMSE/μs
单脉冲 0.813 1 17.958 3
双脉冲 0.539 4 17.871 1
四脉冲 0.402 7 17.871 0

3.2 采样频率

仿真条件如下:信噪比15 dB,t1= 0.1 μs,t2= 0.125 μs的脉冲信号,采用本文设计的差分匹配滤波器和匹配滤波器模型,经过1 000次仿真,得到单脉冲、双脉冲和四脉冲信号均方根误差随采样频率变化的曲线,如图8所示。图8中,Do代表单脉冲差分匹配滤波输出误差随频率变化曲线,Mo代表单脉冲匹配滤波输出误差随频率变化曲线。Dt代表双脉冲差分匹配滤波输出误差随频率变化曲线,Mt代表双脉冲匹配滤波输出误差随频率变化曲线。Df代表四脉冲差分匹配滤波输出误差随频率变化曲线,Mf代表四脉冲匹配滤波输出误差随频率变化曲线。由图8可见,脉冲信号采样频率越高,滤波器输出的均方根误差越小,时间戳标记精确度就越高。同一采样频率时,与单脉冲、双脉冲信号相比,四脉冲信号滤波输出的均方根误差最小,时间戳标记精确度最高。同一采样频率时,差分匹配滤波的均方根误差相较于匹配滤波更小,时间戳标记精确度更高。虽然匹配滤波的输出均方根误差在减小,但是减小的速度远低于差分匹配滤波,在采样频率大于125 MHz时,单、双和四脉冲差分匹配滤波的输出均方根误差变化已经不太明显,但是明显优于匹配滤波,如表3所示。

图8 差分匹配滤波与匹配滤波输出误差随频率变化曲线

Fig.8 Output error curves of differential matched filter and matched filter with frequency

表 3 采样频率125 MHz时输出均方根误差
Table 3 Output RMSE at sampling frequency of 125 MHz
脉冲类型差分匹配滤波RMSE/μs匹配滤波RMSE/μs
单脉冲 6.906 5 12.529 6
双脉冲 6.713 3 10.670 3
四脉冲 6.595 6 9.840 3

3.3 复杂度

仿真条件如下:信噪比20 dB,采样频率40 MHz,t1=0.1 μs,t2= 0.125 μs的单脉冲、双脉冲和四脉冲信号,采用本文设计的差分匹配滤波器,1 000次运行后仿真所需时间如表4所示。由单脉冲、双脉冲和四脉冲的顺序来看,差分匹配滤波输出的均方根误差逐渐减小,时间戳的标记精确度逐渐提高,但是运行的时间却逐渐增加,加大了运算的复杂度,如四脉冲信号差分匹配滤波的精确度约是单脉冲的1.78倍,但四脉冲运行千次的时间却是单脉冲的8.52倍,增加了运算量。

表 4 程序运行时间和精确度对比
Table 4 Comparison of program running time and accuracy
脉冲类型运行千次的时间/s输出RMSE/μs
单脉冲 95.143 441 8.152 1
双脉冲 270.855 339 5.910 4
四脉冲 810.148 857 4.571 2

3.4 真实数据测试

在天津机场终端区设置1个广播式自动相关监视接收机,实时接收的真实数据如图9所示。接收机自身信号采样频率为80 MHz,信噪比约为15.12 dB,为对比分析方便,仿真信号的采样频率为80 MHz,信噪比为15.12 dB,其单、双和四脉冲信号差分匹配滤波器TOA估算精确度结果如表5所示。

图9 真实数据

Fig.9 Real data

表 5 仿真数据与真实数据均方根误差对比及变化率
Table 5 RMSE comparison between simulation data and real data and their change rates
脉冲类型仿真数据真实数据
RMSE/μs相对单脉冲变化率/%RMSE/μs相对单脉冲变化率/%
单脉冲 9.396 1 11.793 45
双脉冲 8.921 9 5.05 10.602 50 10.10
四脉冲 8.224 7 12.5 8.948 75 24.12

表5所示,在采样频率80 MHz、信噪比15.12 dB时,仿真数据的双脉冲较单脉冲TOA估算准确度提升了5.05%,四脉冲较单脉冲TOA估算准确度提升了12.5%。真实数据的双脉冲较单脉冲TOA估算准确度提升了10.1%,四脉冲较单脉冲TOA估算准确度提升了24.12%,说明四脉冲信号差分匹配滤波器的TOA估算精确度更高,验证了算法的实用性。由于各种因素的限制,真实数据与仿真数据不可能完全一样,真实数据的差分匹配滤波均方根误差比仿真数据的大,但真实数据TOA估算精确度提升速度更快。

4 结 论

本文研究了S模式前导脉冲信号的过渡过程与时间戳测量精确度问题。在匹配滤波和差分匹配滤波基本原理的基础上,利用脉冲间的时间结构关系,提出用四脉冲信号差分匹配滤波过零点来标记信号TOA时间戳的改进算法,实验结果表明四脉冲信号的时间戳标记精确度优于单脉冲信号。与匹配滤波的时间戳标记精确度作对比,差分匹配滤波的时间戳标记精确度高于匹配滤波。在采样频率40 MHz、信噪比40 dB时,四脉冲信号差分匹配滤波器的时间戳标记精确度达到0.402 7 μs。在天津机场终端区接收实时信号的采样频率为80 MHz,信噪比约为15.12 dB,四脉冲信号差分匹配滤波器的时间戳标记精确度达到了8.948 75 μs,远小于国际民航组织规定的时间误差精度最大值25 μs[

13⁃14],仿真数据和真实接收数据均验证算法的实用性。

参考文献

1

Torrieri D J. Adaptive thresholding systems[J]. IEEE Transactions on Aerospace and Electronic Systems1977AES-133): 273-280. [百度学术

2

Zehavi E. Estimation of time of arrival for rectangular pulses[J]. IEEE Transactions on Aerospace and Electronic Systems1984AES-206): 742-747. [百度学术

3

Chan Y TLee B HInkol Ret al. Estimation of pulse parameters by convolution[C]/ [百度学术

/Proceedings of Electrical and Computer EngineeringCCECE’06 Canadian Conference on.[S.l.]: IEEE200617-20. [百度学术

4

Chan Y TLee B HInkol Ret al. Estimation of pulse parameters by autoconvolution and least squares[J]. IEEE Transactions on Aerospace and Electronic Systems2010461): 363-374. [百度学术

5

席轶敏刘渝靖晟.电子侦察信号实时检测算法及性能分析[J]. 南京航空航天大学学报2001333): 277-281. [百度学术

XI YiminLIU YuJING Sheng. Analysis on probability of false alarm and detection[J]. Journal of Nanjing University of Aeronautics & Astronautics2001333): 277-281. [百度学术

6

胡国兵刘渝邓振淼.基于Haar小波变换的信号到达时间估计[J].系统工程与电子技术2009317): 1615-1619. [百度学术

HU GuobingLIU YuDENG Zhenmiao. Arrival time estimation of signals based on Haar wavelets transform[J].Systems Engineering and Electronics2009317): 1615-1619. [百度学术

7

GASPARE GMAURO LPatrizio D Met 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. [百度学术

8

Lin W CLiu K CWang C K. Differential matched filter architecture for spread spectrum communication systems[J]. Electronics Letters19963217): 1539. [百度学术

9

Petrochilos NGalati GPiracci E. Separation of SSR signals by array processing in multilateration systems[J]. IEEE Transactions on Aerospace and Electronic Systems2009453): 965-982. [百度学术

10

Galati GLeonardi MTosti M. Multilateration (local and wide area) as a distributed sensor system: Lower bounds of accuracy[C]//Proceedings of European Radar Conference.[S.l.]:IEEE2008. [百度学术

11

王洪金尔文刘昌忠. 多点定位TOA精确估计及同步误差校正算法[J]. 系统工程与电子技术2013354): 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 Electronics2013354): 835-839. [百度学术

12

徐江李艳纪晶华.统计学基础[M]. 北京机械工业出版社2015. [百度学术

XU JiangLI YanJI Jinghua. Statistical basis[M]. BeijingMechanical Industry Press2015. [百度学术

13

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. [百度学术

14

Galati GLeonardi MMagaro Pet al. Wide area surveillance using SSR mode S multilateration: Advantages and limitations[C]//Proceedings of European Radar ConferenceEURAD 2005. Paris: [s.n.], 2005: 225-229. [百度学术

您是第位访问者
网站版权 © 南京航空航天大学学报
技术支持:北京勤云科技发展有限公司
请使用 Firefox、Chrome、IE10、IE11、360极速模式、搜狗极速模式、QQ极速模式等浏览器,其他浏览器不建议使用!