2. 西北工业大学 航天学院,陕西 西安 710072
2. School of Astronautics, Northwestern Polytechnical University, Xi′an 710072, China
在现代军事作战中,被攻击目标对来袭武器的探测识别能力和自身动力性能在不断提升,其在避免被我方武器锁定和攻击的情况下需要进行加速逃逸机动,因此对机动目标跟踪的准确性已成为我方武器能否命中目标的关键影响因素之一[1]。描述目标运动的模型选取是否合适直接影响了目标运动数据处理的质量。然而,目前应用于目标跟踪滤波领域的各种目标运动模型,如运用于直线运动目标的CV、CA模型和运用于转弯机动目标的CT模型均不能准确描述目标的运动状态,尤其是目标的加速度信息等。Singer模型作为一种全局统计模型可以较好弥补上述不足,同时兼顾了非机动模型和机动模型的精度,可以很好适用于上述场景下的机动目标逃逸策略模型构建的过程[2]。近年来,为了提高目标跟踪的准确性,魏喜龙等[3]针对高超声速飞行器的周期跳跃运动特点,提出结合Singer模型的扩展卡尔曼滤波方法对目标运动状态进行估计,但是由于高超声速飞行器的运动特点,需要雷达对其连续跟踪;张燕等[4]针对由于Singer模型不准确导致卡尔曼滤波发散的问题,引入渐消因子抑制滤波器长度,从而避免发散,但其对量测方程的数据获取并没有进行深入研究;周琳等[5]为提高空间目标跟踪效果,采用无迹卡尔曼滤波器估计空间目标运行轨迹,但该方法计算难度大,且容易受到量测数据误差等影响,降低目标跟踪精度;贾舒宜等[6]研究了利用Radon-Ambiguity变换估计出目标的径向加速度并将其加入量测向量中,相比传统不带加速度信息的无迹卡尔曼滤波算法精度有所提高,但对于无径向加速度的机动目标提升效果不大;陈立军等[7]提出一种噪声自适应调整的改进交互式多模型算法,仿真结果表明该方法比传统交互式模型算法准确性高,但在仿真中将主动声呐探测周期设计为1 s,不符合现实应用场景;Zhang等[8]引入支持向量机来估计目标运动模式,降低了目标机动时的跟踪误差,但需要选择合适的核函数,由于算法的复杂度随着数据量的增加而增加,在处理大规模数据集时计算和存储的成本较高。
此外,目标跟踪算法的信息来源如果依赖于主动探测获取,实际应用场景下由于长探测周期间隙中获取目标信息的缺失会造成较大的跟踪误差,但若仅依赖于单独的被动探测信息,同样无法对目标进行较为准确的跟踪。Dong等[9]提出一种基于卡尔曼滤波框架的多传感器协同无源定位方案,并通过IMM-SRCKF算法实现对目标的状态估计;刘妹琴等[10]基于水下传感器网络技术,从数据融合的角度论述了水下目标跟踪技术的关键思想和发展现状,为多传感器联合工作发展提供借鉴。
因此,本文基于Singer模型提出将被动探测信息与主动探测信息相融合的方法,在主动观测模式探测信息获取的间隙插入被动模式获取的信息,通过信息的共享和互补提高了对目标的跟踪精度,同时避免对量测值进行非线性的处理导致算法复杂度提升,构建了虚拟量测值,解决目标对抗逃逸机动场景中由于主动信息缺失导致解算误差较大的问题。本文提出的方法在对抗逃逸目标机动模型中相比主动跟踪模式的平均解算精度更高,提升了跟踪算法性能。
1 运动模型与滤波算法 1.1 Singer模型Singer模型是一种基于全局统计的具有普遍意义的机动运动模型,同时兼顾了非机动和机动模型的精度。将目标加速度的分布视为随机过程建模,目标加速度
r(τ)=Ea(t)a(t+τ)=σ2me−α|τ|,α⩾0。 | (1) |
针对水面机动目标,将其状态矩阵
X = [x˙x¨xy˙y¨y]T。 | (2) |
式中:
X(k+1)=FX(k)+w(k)。 | (3) |
式中:
F=diag(Φ11,Φ22), | (4) |
Φ11 = Φ22 = [1T1α2[−1+αT+e−αT]011α[1−e−αT]00e−aT]。 | (5) |
式中:
Q = 2ασ2m[q11q12q13000q12q22q23000q13q23q33000000q11q12q13000q12q22q23000q13q23q33]。 | (6) |
协方差矩阵中
由于量测值的非线性特点,结合离散卡尔曼滤波算法原理,基于信息融合的Singer模型机动目标跟踪滤波算法流程如图1所示。
![]() |
图 1 基于Singer模型的滤波算法流程图 Fig. 1 Flowchart of filtering algorithm based on Singer model |
滤波算法是一种递推算法,随着滤波过程的进行,历史数据会不断累积,若模型不准确,会存在观测数据与模型不匹配的问题,导致滤波结果发散。为了防止可能出现系统状态模型或过程噪声模型建立不准确,影响模型输入,最终导致结果发散的问题,在滤波的基础上引入自适应渐消滤波的思想,降低旧数据的作用,增强跟踪算法的鲁棒性。具体实现为引入渐消因子对预测协方差进行调整:
P(k+1|k)=FP(k|k)FT+λ(k)Q。 | (7) |
式中,渐消因子
N(k)=V0(k)−HQ(k−1)HT−R(k), | (8) |
M(k)=HFP(k−1|k−1)FTHT, | (9) |
c(k)=tr[N(k)]tr[M(k)], | (10) |
λ(k)={c(k),c(k)>1,1,c(k)⩽1。 | (11) |
式中:
通过上述算法流程即可实现基于信息融合的Singer模型的滤波过程。
2 目标机动策略模型及观测数据融合算法 2.1 目标机动策略模型目标在避免被我方武器锁定和攻击的情况下需要进行加速逃逸机动,对典型目标机动对抗逃逸行为进行研究,定义2种目标机动策略用于算法仿真验证,如图2和图3所示。
![]() |
图 2 第1种目标机动策略示意图 Fig. 2 Schematic of the first targeting maneuvering strategy |
![]() |
图 3 第2种目标机动策略示意图 Fig. 3 Schematic of the second targeting maneuvering strategy |
图2为第1种目标机动策略示意图,目标以18 kn的速度沿南偏东60°的方向航行,同时制导站导引水下武器主动运动靠近目标。当水下武器距离目标5 km时,目标发现来袭武器,此时目标开始机动规避,
图3为第2种目标机动策略示意图,当武器到达距离目标5 km的位置时目标报警,A为报警时刻武器的位置,B为报警时刻目标的位置。报警后目标在B点开始机动规避并释放诱饵,诱饵不需要跟踪,只需要通过诱饵计算机动航向,BD为诱饵方向(
为了对目标运动要素进行解算,建立观测方程如下:
Z(K)=HX(K)+W(K)。 | (12) |
式中:
制导站通过发射探测信号对目标进行观测,其初始所在位置为
Rz(k)=√(xt(k)−x0)2+(yt(k)−y0)2+vR(k), | (13) |
βz(k)=arctan(yt(k)−y0xt(k)−x0)+vβ(k)。 | (14) |
式中:
δb(k)=arctan(yt(k)−y0xt(k)−x0)+vδ(k)。 | (15) |
式中:
由于主动探测装置需要发射信号并接收目标反射的回波信号,具有一定时段的信息空档期,在这期间由被动探测装置持续接收来自目标的辐射信号,通过将两者的探测信息进行共享与融合,可以提高目标跟踪的准确性和可靠性。
假定主动探测间隔周期为
R(k)={Rz(k),k=k0,k0+N,k0+2N...,preRz(k),k=k0+1,k0+2...k0+N−1..., | (16) |
β(k)={βz(k),k=k0,k0+N,k0+2N...,δb(k),k=k0+1,k0+2...k0+N−1...。 | (17) |
式中:
传统的纯主动或纯被动目标跟踪算法需要观测方程的线性化,将观测矩阵用雅克比矩阵代替,即式(15)和式(16)的观测矩阵为:
{H=∂Rz(k)∂X(k),H=∂βz(k)∂X(k)。 | (18) |
为了降低算法复杂度,避免对量测信息进行线性化的处理,融合后的观测信息通过构建虚拟量测值
Z=[xy]=[RcosβRsinβ]。 | (19) |
观测矩阵简化为:
H=[100000000100]。 | (20) |
综上,基于信息融合的Singer模型机动目标跟踪,包含了机动目标策略模型、Singer模型、观测数据融合和滤波等4个部分,算法原理框图如图4所示。
![]() |
图 4 算法原理框图 Fig. 4 Schematic diagram of the algorithm |
在仿真过程中,假设主动模式下具有较长的探测周期即10 s,被动探测周期为1 s,主动距离、主被动方位探测精度分别为100 m、3°和3°。制导站位置在[0,0] m处,目标初始位置在[
图5为第一种机动策略轨迹估计情况。在没有被动信息的情况下,信息空档期时的方位量测值为最近一次方位探测结果,从图5(a)中可以看出在纯主动模式下的信息空挡期间估计轨迹与真实轨迹存在明显误差,并且在机动阶段误差更大。图5(b)为基于信息融合的第1种机动策略轨迹估计情况,加入被动信息的Singer模型可以更为准确估计速度大小和方向变化的机动目标的航迹。
![]() |
图 5 第1种机动策略轨迹估计情况 Fig. 5 The trajectory estimation scenario for the first maneuvering strategy |
图6为仿真实验50次的情况下,2种模式下x方向距离估计和y方向距离估计的均方差对比示意图,x方向距离估计的均方差,主动观测模式下区间基本在
![]() |
图 6 第1种机动策略均方差对比图 Fig. 6 Comparison of mean squared error for the first maneuver strategy |
图7(a)为基于主动观测信息的第2种机动策略轨迹估计情况,图7(b)为基于信息融合的第2种机动策略轨迹估计情况,加入被动信息的Singer模型依旧能有效估计速度大小和方向变化的机动目标航迹。
![]() |
图 7 第2种机动策略轨迹估计情况 Fig. 7 The trajectory estimation scenario for the second maneuvering strategy |
图8将2种模式下x方向距离估计和y方向距离估计的均方差进行了对比,x方向距离估计的均方差,主动观测模式下区间基本在
![]() |
图 8 第2种机动策略均方差对比图 Fig. 8 Comparison of mean squared error for the second maneuver strategy |
Singer模型在描述目标运动时,可以很好兼顾机动状态与非机动状态的精度,可以准确地描述典型对抗逃逸机动模型,影响跟踪精度的主要因素为观测信息的完整性与准确性,在面对主动观测信息缺失的情况下,在主动观测模式探测信息获取的间隙插入被动模式获取的信息,通过信息的共享和互补提高了观测信息的完整性,进而提升了平均解算精度。
3.3 渐消因子有效性在模型准确的情况下加入渐消因子解算精度略有提升。为验证渐消因子的有效性,假设系统状态模型建立或过程噪声模型不准确,导致错误输入,反映在系统矩阵发生变化:
Φ′11 = Φ′22 = [101α2[−1+αT+e−αT]001α[1−e−αT]000]。 | (21) |
以2种机动策略为例,跟踪效果对比如图9所示。可知,由于模型状态输入错误,估计轨迹与真实轨迹都存在较大偏差,但在加入渐消因子的情况下增大了新量测值的作用,跟踪效果相比于未加入渐消因子效果有明显提升,使目标不易摆脱跟踪,对滤波发散的情况有抑制作用。
![]() |
图 9 模型错误情况下是否加入渐消因子的估计轨迹对比 Fig. 9 Comparison of trajectory estimation with and without the inclusion of asymptotic factors in case of model error |
针对在目标机动规避场景中的主动观测模式下,由于长周期间隙中探测信息的缺失导致机动目标跟踪误差较大的问题,本文基于Singer模型提出将被动探测信息与主动探测信息相融合的方法,在主动模式探测信息获取的间隙插入被动模式获取的信息,通过信息的共享和互补提高了对目标的跟踪精度,同时避免对量测值进行非线性的处理导致算法复杂度提升,构建了虚拟量测值。同时为了避免边界条件导致算法发散,在滤波中引入渐消因子,提升了算法的鲁棒性。仿真结果表明,在2种典型对抗逃逸目标机动模型下该算法相比于主动跟踪模式的解算结果在平均解算精度上分别提升了39.85%和45.71%,验证了算法的有效性。
本文将主动模式信息缺失阶段的距离探测信息保持为最近一次探测结果,会与实际有偏差,在未来的研究中可以考虑对其进行进一步处理,以更符合实际应用场景。
[1] |
LI Dongdong, LIN Yang, ZHANG Yao. A track initiation method for the underwater target tracking environment[J]. China Ocean Engineering, 2018, 32(2): 206-215. DOI:10.1007/s13344-018-0022-0 |
[2] |
JIA Shuyi, ZHANG Yun, WANG Guohong. Highly maneuvering target tracking using multi-parameter fusion Singer model[J]. Journal of Systems Engineering and Electronics, 2017, 28(5): 841-850. DOI:10.21629/JSEE.2017.05.03 |
[3] |
魏喜庆, 顾龙飞, 李瑞康, 等. 基于Singer模型的高超声速飞行器轨迹跟踪与预测[J]. 航天控制, 2017, 35(4): 62-66+72. |
[4] |
张燕, 柳超, 李云鹏. 基于改进Singer模型的机动目标跟踪方法[J]. 火控雷达技术, 2015, 44(3): 37-40+50. DOI:10.3969/j.issn.1008-8652.2015.03.009 |
[5] |
周琳, 刁伟峰, 王祎. 基于可观性分析的高精度空间目标跟踪方法[J]. 雷达科学与技术, 2021, 19(5): 485-490. DOI:10.3969/j.issn.1672-2337.2021.05.003 |
[6] |
贾舒宜, 张赟. 一种基于径向加速度的Singer-EKF机动目标跟踪算法[J]. 海军航空工程学院学报, 2014, 29(2): 141-145+182. |
[7] |
陈立军, 郝保雷, 郑剑. AUV平台的水下机动目标跟踪[J]. 网络新媒体技术, 2020, 9(5): 48-54. DOI:10.3969/j.issn.2095-347X.2020.05.007 |
[8] |
ZHANG Jiaoxin, LIU Meigin, FAN Zhen. Classify motion model via SVM to track underwater maneuvering target[C]//2019 IEEE International Conference on Signal Processing, Communications and Computing (ICSPCC). IEEE, 2019: 1−6.
|
[9] |
DONG C, GUO Z, CHEN X, et al. Research on passive location of maneuvering target based on IMM-SRCKF algorithm[C]//Signal and Information Processing, Networking and Computers: Proceedings of the 8th International Conference on Signal and Information Processing, Networking and Computers (ICSINC). Singapore: Springer Nature Singapore, 2022: 696−703.
|
[10] |
刘妹琴, 韩学艳, 张森林, 等. 基于水下传感器网络的目标跟踪技术研究现状与展望[J]. 自动化学报, 2021, 47(2): 235-251. |
[11] |
SINGER R A. Estimating optimal tracking filter performance for manned maneuvering targets[J]. IEEE Transactions on Aerospace and electronic systems, 1970(4): 473-483. |