中国科技核心期刊      中国指挥与控制学会会刊     军事装备类重点期刊
理论研究

基于最优数据压缩的雷达和ESM间断点迹融合跟踪

  • 王永安 1 ,
  • 李世忠 2 ,
  • 王国宏 2
展开
  • 1.91880部队, 山东 胶州 266300
  • 2.海军航空大学信息融合研究所, 山东 烟台 264001

王永安(1972-),男,山东即墨人,高级工程师,研究方向为雷达数据处理。

李世忠(1980-),男,博士,工程师。

王国宏(1963-),男,教授,博士生导师。

收稿日期: 2017-09-16

  修回日期: 2017-09-20

  网络出版日期: 2022-05-09

Fusion Precision of Intermittent Radar and ESM Data Based on Optimal Data Compression

  • WANG Yong-an 1 ,
  • LI Shi-zhong 2 ,
  • WANG Guo-hong 2
Expand
  • 1. Unit 91880, Jiaozhou 266300, China
  • 2. Institute of Information Fusion Technology, NAAU, Yantai 264001, China

Received date: 2017-09-16

  Revised date: 2017-09-20

  Online published: 2022-05-09

摘要

雷达和ESM传感器均是飞机重要的信息获取手段,但当对己方雷达的辐射进行控制,且敌方雷达也间歇工作时,可能出现雷达和ESM均测量间断的问题。针对该问题,提出了一种基于最优数据压缩的雷达和ESM间断点迹融合跟踪方法,有效解决了雷达和ESM量测间断的问题。并且针对传感器不同测量精度、传感器不同数据率和雷达不同关机时长对融合精度的影响进行了仿真分析。仿真结果表明,ESM测角精度是否小于雷达测角精度对融合精度有重要影响,而融合误差则受雷达测量误差的影响相对较大。

本文引用格式

王永安 , 李世忠 , 王国宏 . 基于最优数据压缩的雷达和ESM间断点迹融合跟踪[J]. 指挥控制与仿真, 2017 , 39(6) : 38 -42 . DOI: 10.3969/j.issn.1673-3819.2017.06.009

Abstract

Both radar and ESM fusion are important sensors in airborne information processing. When the radar sensors are under the condition of radiation, the radar and the ESM sensors cannot get the measurements continuously. Aiming at this problem, a method of fusion precision of intermittent radar and ESM data based on optimal data compression is proposed in this paper. And the problem of radar and ESM measurement interruption could be effectively solved using this method. In addition, this paper studies the influence of the radar and ESM sensors’ precision and data rate as well as the length of shut down time on the fusion performance. Simulation results show that whether the ESM angle measuring accuracy is smaller than the radar angle measuring accuracy has an important influence on the fusion accuracy, and the fusion error is affected by the radar measurement error largely.

雷达和电子支援措施(ESM)是最主要的两种机载传感器,具有性能互补的特点,通过对二者的数据进行融合,可实现信息互补,取得单传感器无法完成的功能和效果,改善对目标的探测和跟踪性能[1-3]。文献[4-11]研究了雷达与ESM的数据关联和融合问题,其中,文献[4-5]研究了机载雷达与ESM的关联问题,文献[6]研究了雷达与ESM的误差配准问题,文献[7]研究了雷达与ESM的融合跟踪问题,文献[8]研究了雷达与ESM的综合检测、跟踪和识别问题,文献[9]研究了ESM测量间歇情况下的雷达与ESM的协同跟踪及雷达辐射控制问题,文献[10]研究了修正极坐标系下雷达与ESM航迹对准关联的问题,文献[11]研究了数据链航迹与ESM间断航迹的融合跟踪。
在实际应用中,当敌方雷达关机时,ESM就不能连续获取目标的测量,从而得到间断的目标信息量测;当对己方雷达的电磁辐射进行控制时,雷达也面临间断测量的情况,从而面临雷达和ESM均为间断测量情况下的数据融合问题,而关于这方面的研究尚未见到相关报道。因此,本文研究了雷达和ESM间断量测情况下的数据融合,并对影响雷达和ESM融合精度的多种因素进行了仿真分析。

1 基于最优点迹数据压缩融合处理的流程

雷达点迹和ESM点迹数据融合处理主要有两种方式:一种是序贯处理方式,一种是基于点迹数据最优压缩的处理方式。但序贯处理方式计算量较大,计算复杂度较高,因此从实际工程应用方面考虑,本文主要以第二种方式为对象进行研究讨论。
雷达和ESM点迹数据压缩融合处理的流程如图1所示,主要包括以下步骤:
图1 ESM点迹+雷达点迹数据压缩滤波流程图
1)以ESM测量点迹的时间为基准,对雷达点迹进行外推或内插平滑;
2)雷达点迹与ESM点迹进行关联;
3)以外推或内插的雷达点迹为基准,采用最优压缩方法对关联上的ESM和雷达点迹进行融合;
4)对融合后的点迹进行滤波以得到融合航迹,为了避免由极坐标到直角坐标转换带来的误差,可采用无偏不敏变换(UUKF)滤波算法。

2 雷达和ESM点迹的数据关联和最优压缩融合

本文以方位角θ为例给出数据关联算法。假设在雷达的一个测量周期[tk-1,tk]内ESM有n个量测,即
Ze= θ e ( 1 ) θ e ( 2 ) ··· θ e ( n ) T
设雷达相应量测为θr(k-1)、θr(k),求得ESM量测与θr(k-1)的时间差值,记为Δtn,将ESM的n个量测用下式内插到雷达测量周期内:
θr(i)= θ r ( k ) - θ r ( k - 1 ) T RΔtnr(k-1),i=1,···,n
则得到雷达的量测集合:
Zr= θ r ( 1 ) θ r ( 2 ) ··· θ r ( n ) T
于是,可构造如下检验统计量:
d= i = 1 n θ r ( i ) - θ e ( i ) σ r 2
则该统计量服从自由度为nχ2分布。设置信度为1-α,利用χ2分布表可得门限 χ α 2(n),若
d< χ α 2 (n)
则说明该雷达周期内的ESM量测与雷达量测为同一目标的量测。
N个雷达周期进行上述方法的检测,若有M次使得该式成立(如N=4,M=3),则认为ESM量测与雷达量测为同一目标的量测,从而实现了点迹的关联。
首先对雷达和ESM测量的角度数据进行时间对准,然后利用最优压缩方法将对准后的雷达和ESM角度测量数据进行最优压缩,从而得到融合方位角测量值。以目标方位角为例,最优融合方位角为:
$\hat{\theta}_{i}=\left(\frac{\theta_{i}}{\sigma_{R}^{2}}+\frac{\tilde{\theta}_{i}}{\sigma_{E}^{2}}\right) /\left(\frac{1}{\sigma_{R}^{2}}+\frac{1}{\sigma_{E}^{2}}\right)$
同理,可对俯仰角数据进行压缩,从而获得融合后的目标测量点迹( θ ^ i, β ^ i,ri)。

3 航迹起始

为了实现对目标的跟踪,需要对融合测量进行航迹起始。兼顾航迹起始性能和航迹起始时间,这里采用3/4逻辑法对航迹进行起始,在航迹起始中采取如下限制:
1)首先,目标速度要介于最小速度Vmin最大速度Vmax之间;
2)其次,目标加速度绝对值要小于最大加速度amax;当有多个回波的加速度满足约束时,则选择最小加速度绝对值对应的回波。

4 雷达和ESM点迹滤波算法

兼顾性能和计算量,同时也为了避免坐标转换引起的固定偏差,这里采用无偏转换测量[12]Kalman滤波方法。设k时刻雷达的量测向量为
Z(k)= r ( k ) θ ( k ) ϕ ( k )
其协方差为
R= σ r 2 0 0 0 σ θ 2 0 0 0 σ ϕ 2
由于Z(k)为三维向量,因此选取7个sigma点,即
Zpi(k)= r i ( k ) θ i ( k ) ϕ i ( k ) T=Z(k)+ R·Ai(i=1,2,···,7)
其中
A= 0 η - η 0 0 0 0 0 0 0 η - η 0 0 0 0 0 0 0 η - η
η= 3 + λ
其中,λ为自由度,此处取λ=0。
Zmi(k)= c 1 r i ( k ) c o s ( ϕ i ( k ) ) c o s ( θ i ( k ) ) c 1 r i ( k ) c o s ( ϕ i ( k ) ) s i n ( θ i ( k ) ) c 2 r i ( k ) s i n ( ϕ i ( k ) )
其中
c 1 = ε 1 - 1 e x p ( σ θ 2 + σ ϕ 2 ) 2 c 1 = ε 2 - 1 e x p σ ϕ 2 2 ε 1 = ( 1 + c o s ( 3 σ θ ) + c o s ( 3 σ ϕ ) ) 3 ε 2 = ( 2 + c o s ( 3 σ ϕ ) ) 3
k时刻雷达新的量测向量为
Zm(k)= 1 6 i = 2 7Zmi(k)
相应的协方差阵为
Rm(k)= 1 6 i = 2 7[Zmi(k)-Zm(k)]·[Zmi(k)-Zm(k)]T
由此可得到量测方程,然后可按照卡尔曼滤波方法进行滤波处理,即可得到目标航迹。

5 仿真结果分析

5.1 仿真场景及初始参数

假设我方飞机编队有2架飞机,坐标系选为NED坐标系,目标与我方飞机均做匀加速直线飞行,仿真时间100s,等间隔间断4次,导航数据间隔为0.1s,融合数据输出间隔为0.5s,雷达测量周期为0.5s;ESM传感器等间隔间断5次,测量周期为0.1s,ESM和雷达关机时间比例均为0.2;目标初始位置、速度与文献[7]相同。

5.2 传感器误差对跟踪精度的影响分析

设目标与我方飞机均做匀加速直线运动;雷达测距均方根误差σd=100m,雷达方位和俯仰角测量均方根误差均为σR;ESM的方位和俯仰角测量均方根误差均为σE;Monte Carlo仿真100次,分4种情况进行仿真, 4种场景的设置如表1所示,对4种场景进行仿真后的精度比较如表2所示。
表1 不同测量精度的设置
场景 1 2 3 4
σR/(°) 0.3 0.6 0.6 0.8
σE/(°) 0.3 0.3 0.6 1.5
表2 不同场景下雷达误差与融合跟踪误差比较
精度场景 距离误差 方位误差 俯仰误差 位置误差
均值 极大值 极小值 均值 极大值 极小值 均值 极大值 极小值 均值 极大值 极小值
1 雷达误差 61.49 205.3 0.466 0.036 0.208 0.001 0.040 0.134 0.001 131.7 392.5 44.78
融合误差 57.83 201.1 0.918 0.022 0.111 0.001 0.056 0.166 0.001 138.4 315.0 49.30
2 雷达误差 59.16 266.9 1.814 0.080 0.174 0.011 0.130 0.270 0.006 300.9 475.6 25.41
融合误差 62.28 308.5 4.616 0.035 0.132 0.002 0.050 0.161 0.002 137.6 385.5 33.16
3 雷达误差 46.80 199.1 0.100 0.120 0.305 0.017 0.123 0.282 0.004 326.0 583.1 123.2
融合误差 53.33 197.1 0.763 0.116 0.213 0.021 0.060 0.273 0.007 252.7 481.9 79.01
4 雷达误差 53.38 215.6 0.557 0.173 0.256 0.109 0.171 0.520 0.001 447.9 991.3 239.5
融合误差 55.71 286.8 11.98 0.209 0.444 0.080 0.391 0.934 0.111 791.5 0.001 385.7
由仿真结果可得如下结论:
1)不管雷达与ESM的测角精度如何变化,融合后距离误差基本不变,说明雷达与ESM点迹融合对测距精度影响较小,这是因为距离精度主要受雷达测量距离测量精度的影响所致。
2)在方位、俯仰和总体上的融合精度随着ESM精度的下降稍有下降,随着雷达精度的下降,在方位、俯仰和总体上融合后的精度下降明显,由此可见,雷达测量精度对融合精度的影响要明显大于ESM测量精度的影响。
3)当ESM测角精度与雷达测角精度相等时,融合后精度略有提高;当ESM测角精度好于雷达的测角精度时,融合精度高于雷达测量精度,当ESM测角精度低于雷达测角精度时,融合精度低于雷达测量精度,说明当σER≤1时,即ESM精度高于雷达的测角精度时,雷达点迹与ESM点迹融合可以提高跟踪精度。

5.3 传感器数据率对跟踪精度的影响分析

设雷达测距均方根误差σd=0.1km,雷达方位角、俯仰角均方根误差为0.3°,雷达的周期用TR表示; ESM方位角、俯仰角均方根误差为0.5°,ESM的周期用TE表示;目标与我方飞机均做匀加速直线飞行, Monte Carlo仿真100次,3种场景的设置如表3所示,对3种场景进行仿真后的精度比较如表4所示。
表3 不同传感器数据率情况下的场景设置
场景 1 2 3
TR/s 0.1 0.5 0.5
TE/s 0.05 0.05 0.1
表4 不同场景下雷达误差与融合跟踪误差比较
精度场景 距离误差 方位误差 俯仰误差 位置误差
均值 极大值 极小值 均值 极大值 极小值 均值 极大值 极小值 均值 极大值 极小值
1 雷达误差 45.33 188.6 0.026 0.017 0.063 0.001 0.022 0.053 0.001 73.44 202.5 6.269
融合误差 45.08 197.3 0.409 0.029 0.091 0.002 0.032 0.123 0.001 100.9 298.1 35.24
2 雷达误差 54.47 496.4 0.314 0.039 0.229 0.001 0.033 0.172 0.001 126.2 636.5 6.514
融合误差 54.50 474.8 2.085 0.093 0.260 0.005 0.103 0.307 0.006 264.3 699.5 26.32
3 雷达误差 48.60 176.9 0.578 0.033 0.166 0.001 0.055 0.134 0.001 135.1 303.3 24.86
融合误差 52.75 156.8 1.497 0.063 0.252 0.002 0.049 0.225 0.006 277.7 440.5 59.91
由仿真结果可得如下结论:
1)当ESM数据率不变时,提高雷达精度可以提高雷达及融合的位置精度;
2)当雷达数据率不变时,系统跟踪精度主要与雷达数据率有关,与ESM数据率关系相对较小;
3)如果ESM的测角精度差于雷达的测角精度,则位置融合精度均略差于雷达位置精度。

5.4 传感器关机时长对跟踪精度的影响分析

设雷达测距均方根误差σd=0.1km,雷达方位角、俯仰角均方根误差为0.3°,雷达的周期TR=0.5s,等间隔间断4次,ESM方位角、俯仰角均方根误差为0.5°,ESM的周期用TE=0.1s,等间隔间断5次,ESM关机时间占总工作时间的20%;。目标与我方飞机均做匀加速直线飞行,Monte Carlo仿真100次。雷达关机时间占总工作时间的百分比用TClose表示,分3种情况进行分析,3种场景的设置如表5所示,对3种场景进行仿真后的精度比较如表6所示。
表5 不同传感器关机时长的场景设置
场景 1 2 3
TClose/s 20% 40% 60%
表6 不同场景下雷达误差与融合误差比较
精度场景 距离误差 方位误差 俯仰误差 位置误差
均值 极大值 极小值 均值 极大值 极小值 均值 极大值 极小值 均值 极大值 极小值
1 雷达误差 58.78 218.4 0.970 0.022 0.126 0.001 0.029 0.157 0.001 101.9 289.9 16.27
融合误差 53.83 226.1 2.586 0.061 0.183 0.002 0.061 0.149 0.005 176.0 357.2 49.09
2 雷达误差 79.64 409.9 8.044 0.054 0.190 0.001 0.085 0.150 0.003 204.8 536.8 54.53
融合误差 71.56 351.5 1.144 0.067 0.258 0.005 0.055 0.317 0.006 184.6 684.6 45.17
3 雷达误差 56.59 207.7 0.346 0.097 0.321 0.009 0.078 0.241 0.001 239.3 841.1 39.03
融合误差 90.27 479.0 1.016 0.114 0.831 0.011 0.098 0.502 0.005 315.1 1438 69.75
由仿真结果可得如下结论:
1)当雷达的关机时间占总时间的比例提高时,雷达及融合后的距离、方位、俯仰和位置精度整体上在逐步变差,位置精度的均值下降。
2)如果ESM的测角精度差于雷达的测角精度,雷达关机时间扩大时,融合后位置精度仍然低于融合前雷达的位置精度。

6 结束语

本文对雷达和ESM数据不连续情况下基于最优压缩的融合滤波进行了仿真分析,分别在不同测量精度、不同传感器数据链和不同关机时长情况下对雷达测量精度与融合精度进行了比较。结果表明:ESM测角精度是否小于雷达测角精度对融合精度有重要影响,当ESM的测角误差小于雷达测角误差或者大致相当时,融合航迹误差要小于雷达航迹误差,反之则前者略大于后者;相对而言,融合误差受雷达测量误差的影响相对较大;当雷达和ESM数据率降低时,融合误差均增加;关机时长增加时,融合误差也显著增大。
[1]
何友, 王国宏, 陆大金. 多传感器信息融合及应用[M]. 第2版. 北京: 电子工业出版社, 2009.

[2]
关欣, 彭彬彬, 衣晓. 基于区间重合度的雷达与ESM航迹关联算法[J]. 雷达科学与技术, 2017, 15(1): 61-67.

[3]
肖林, 王华. 改进的ESM目标与雷达航迹关联方法[J]. 信息化研究, 2016, 42(1): 35-39.

[4]
董红波, 王刚, 刘进忙, 等. 机载雷达与ESM航迹关联技术研究[J]. 传感器与微系统, 2010, 29(1): 12-14.

[5]
陈中华, 王国宏, 刘德浩, 等. 基于几何法的雷达与ESM航迹关联算法[J]. 电光与控制, 2012, 19(4): 10-12,22.

[6]
刘德浩, 王国宏, 陈中华. 基于UKF的雷达与ESM误差配准技术研究[J]. 电子测量技术, 2011, 34(6): 26-29.

[7]
张翔宇, 王国宏, 王娜, 等. 基于序贯滤波的机载雷达和ESM融合跟踪[J]. 电子测量技术, 2011, 34(9): 44-47.

[8]
石绍应, 杜鹏飞, 张靖, 等. 雷达与ESM综合多目标检测、跟踪与识别[J]. 系统工程与电子技术, 2016, 38(7): 1524-1531.

[9]
吴巍, 王国宏, 李世忠, 等. ESM量测间歇下雷达/ESM协同跟踪与辐射控制[J]. 现代防御技术, 2011, 39(3): 132-138.

[10]
关欣, 彭彬彬, 衣晓. 修正极坐标系下雷达与ESM航迹对准关联[J]. 航空学报, 2017, 38(5): 221-232.

[11]
王永安, 李世忠, 关成斌. 数据链航迹与ESM间断航迹融合跟踪精度研究[J]. 海军航空工程学院学报, 2012, 27(4): 401-408.

[12]
S. J. Julier, J. K. Uhlmann, H. F. Durrant-Whyte. A new method for the nonlinear transformation of means and covariances in filters and estimators[J]. IEEE Trans on Automatic Control, 2000, 45(3): 477-482.

DOI

文章导航

/