基于似然比检验的目标机动检测方法

2024-02-22 00:00:00邵雷贺杨超赵锦
航空兵器 2024年6期

摘 要:""""" 针对博弈对抗过程目标意图识别中的机动检测问题, 设计了一种基于似然比检验的目标机动检测方法。 该方法利用目标实时观测数据进行机动检测, 将目标机动检测问题转化为目标观测时间序列数据检验问题。 首先, 基于极大似然原理对目标运动状态进行回归估计; 其次, 利用似然比检验原理设计检测模型, 构建目标机动检测方法。 仿真结果表明, 该方法在噪声环境下具有较好的误检测和检测延迟性能。

关键词:"""" 机动检测; 似然比; 回归估计; 时间序列估计

中图分类号:"""""" TJ760

文献标识码:""" A

文章编号:"""" 1673-5048(2024)06-0064-06

DOI: 10.12132/ISSN.1673-5048.2024.0052

0 引" 言

随着信息化技术的发展, 人工智能技术在空战领域博弈对抗中的应用越来越广泛, 智能化空战成为未来空战的一种重要形式[1-2。 如何根据目标的运动状态快速评估目标的意图, 为己方决策提供支持, 进而制定出针对性对抗策略, 成为智能化空战过程取得胜利的重要性影响因素之一[3-5。 基于此, 国内外众多学者在目标意图提取方面进行了大量的探索。 赵智伟等[6将目标机动识别信息应用于机动意图分析, 为空战态势评估提供支持; Li等[7将基于目标意图评估应用于目标轨迹预测, 可有效提高轨迹预测精度, 为目标拦截方案制定提供支撑; 针对给定运动模式, 乔殿峰等[8采用分段识别方式, 通过引入STI算法对模型差异小的非预分段区航迹进行精细化识别, 提高目标行为意图判断的准确性和快速性。 针对预定的机动模式, 吴家湖等[9采用循环神经网络识别方法实现对目标机动模式的识别; 罗艺等[10在分析目标运动特性的基础上, 将目标的运动过程建模为一个马尔科夫过程, 以分布概率的方式实现了目标的意图轨迹预测; 张凯等[11通过引入意图代价函数, 评估目标机动模式, 提高对目标机动模式的识别准确度, 进而提高目标轨迹预测精度。

上述方法提取意图过程中, 以目标先验机动模式为基础, 通过观测数据与先验目标机动模式的匹配, 实现对目标的机动检测与模式识别。 然而, 目标先验机动模式在对抗中往往难以准确获取, 如何直接利用观测数据进行机动检测与识别成为对抗过程中的一个关键性问题。 在假定目标意图攻击要地的条件下, 张凯等[12利用贝叶斯理论辨识目标机动模式, 借助采样方法实现轨迹预测; 李世杰等[13通过分析要地重要性和目标可达域的关联关系, 对目标轨迹参数进行先验预测, 缩小了目标机动的可能范围, 实现了较高精度的中长期轨迹预测; Hu等[14利用贝叶斯原理, 在对目标机动规律进行辨识的基础上, 借助目标意图得到积分预测轨迹。

在上述基于目标运动状态的意图评估中, 目标机动的检测与快速判断是一个重要前提, 快速的机动检测可有效提高对目标运动轨迹的跟踪预测精度[15-16。 对于目标机动检测问题, 以提高目标跟踪精度为目标, 国内外学者对具有目标机动检测的滤波跟踪算法进行研究[17, 试图在检测到目标机动后通过目标机动补偿提高目标跟踪精度。 这类算法通常将目标机动检测与滤波跟踪算法融合设计[17-18, 以提高跟踪精度为目标, 但机动检测的虚警、 延迟等能力有待进一步加强[19-20

在基于意图识别的博弈对抗中, 目标机动检测的虚警将导致意图识别错误, 影响后续的决策。 为进一步提高目标机动检测性能, 苏建敏等[21对小波分解后的信号进行统计, 并通过与设定阈值的对比实现卫星轨道机动检测, 但没有给出阈值选取方法; 崔红正等[22以观测值残差比为特征, 分析了不同推力下的目标机动特征; 王庆瑞等[23采用概率统计方法, 基于Neyman-Pearson准则, 利用相对速度变化量对发动开机带来的目标机动进行检测。 上述机动目标检测通常以相对速度为特征量, 在高精度相对速度测量条件下, 对于具有阶跃过载机动特性目标, 具有较好的检测性能。 然而, 多数防空雷达通常仅能直接测量目标位置信息, 同时空气动力目标的机动过载通常是一个连续变化量, 往往不会产生瞬时阶跃机动。 针对空气动力目标机动检测的这些不同特征, 本文以目标位置测量信息为特征, 将基于极大似然原理的目标运动状态回归估计与似然比检验原理相结合, 设计目标检测方法, 为目标运动状态意图分析提供支撑。

1 问题描述

利用目标跟踪获取的观测数据(xk, yk)(k=1, 2, …, N)对目标机动与否进行判断实际上是一种时间序列数据的检测, 本质上是一个假设检验问题。 若目标运动状态为没有机动, 则目标运动轨迹呈线性关系。 为简化问题分析, 文中以二维平面为例进行研究, 假定目标运动轨迹为

y=a·x+b+v(1)

式中: x, y分别为目标水平面纵向和横向的运动轨迹坐标; a, b为目标运动状态参数; v为观测噪声。

若目标在此期间进行了机动, 则运动轨迹将相对该状态发生偏移, 可表示如下:

y=a·x+b+f(x, y)+v(2)

式中: f(x, y)为目标机动带来的偏移量。

从上述分析可以看出, 根据目标观测数据对目标进行机动检测, 可描述为如下关于新观测数据相比历史运动状态偏移量是否存在的判断问题, 即

H0:yk-a·xk-b=vkk=1, 2, …, N

H1:yk-a·xk-b≠vkk=1, 2, …, N(3)

式中: H0表示目标未机动; H1表示目标进行了机动。

2 目标机动检测原理

考虑在实际对抗过程中, 对目标数据的观测是一个周期性的动态过程; 因此, 对目标的机动检测, 可描述为新观测数据相对历史运动状态变化判断的动态过程。 基于上述分析对目标进行机动检测, 可分为两个环节: 一是对目标历史运动状态参数的估计; 二是以新获取目标数据为输入, 基于目标历史运动状态信息对目标机动与否进行判断。

2.1 目标运动状态参数估计模型

假定目标在观测窗口期间没有发生机动, 同时假定测量噪声服从均值为0、 方差为σ2的独立正态分布, 则目标的运动状态参数估计可转化为回归问题。

假定目标历史状态没有机动, 定义如下的条件概率密度函数:

p((x, y)|H0)=

1(2πσ2)N/2exp-12σ2∑Nk=1(yk-a·xk-b)2(4)

利用极大似然原理, 可以得到目标运动状态参数估计为[24

a^b^=A-1B(5)

式中: A=∑Nk=1x2k∑Nk=1xk

∑Nk=1xkN; B=∑Nk=1xkyk

∑Nk=1yk。

2.2 目标机动检测模型

基于假设检验原理, 若观测序列(xk, yk)(k=1, 2, …, N)在条件H0和条件H1下的概率密度函数分布为p((x, y)|H0)和p((x, y)|H1), 则似然比函数[25

p((x, y)|H1)p((x, y)|H0)=

1(2πσ2)N/2exp-12σ2∑Nk=1(yk-a·xk-b-μ-)21(2πσ2)N/2exp-12σ2∑Nk=1(yk-a·xk-b-μ)2(6)

式中: μ为目标运动状态测量误差均值, 对于给定假设vk, 该值为0; μ-为目标观察数据误差的最大似然估计均值。

似然比描述了观测序列(xk, yk)(k=1, 2, …, N)属于H0和H1的可能性比值。 根据似然比检验原理, 选取

p(x|H1)p(x|H0)gt;c(7)

式中: c∈(0, 1)为检测阈值。

根据式(7), 若拒绝假设H1, 接受假设H0, 即认为目标未机动; 否则, 认为目标进行了机动。

定义误差序列:

Δk=yk-a·xk-b "k=1, 2, …, N (8)

将式(6), (8)代入式(7), 并对两边取对数, 简化可得

-12σ2∑Nk=1(Δk-μ-)2+12σ2∑Nk=1(Δk-μ)2gt;lnc(9)

利用恒等式:

∑Nk=1(Δk-μ-)2=∑Nk=1(Δk-μ)2+N(μ-μ-)2(10)

则式(9)变为

(μ-μ-)2lt;-2σ2Nlnc (11)

从式(11)可以看出, 目标机动的检测问题实际变为新的观察数据与目标历史运动状态之间的轨迹统计值判断问题。 满足式(11)则可认为目标没有发生机动。

2.3 目标机动检测方法

基于上述模型进行目标机动检测, 可将观测数据(xk, yk)(k=1, 2, …, N)分为两段, 将前段作为目标运动状态参考窗口, 描述目标当前运动状态, 认为目标没有机动, 采用式(5)获取目标运动状态信息; 将包含最新观测数据的后段作为检测窗口, 当目标没有机动时, 目标运动状态不发生变化, 可认为误差统计量不发生变化, 按照式(11)对目标进行机动检测判断。

基于上述思想进行目标机动检测, 基本步骤如下:

(1) 在每一工作周期获取最新观测数据, 构建检测窗口数据集合{(xk, yk)(k=1, 2, …, N)}, 并分割观测数据, 得到检测参考窗口Rk=(xk, yk)(k=1, 2, …, N1)和机动检测窗口Tk=(xk, yk)(k=N-M+1, …, N);

(2) 目标运动状态估计, 利用参考窗口Rk, 按照式(5)估计目标历史状态参数;

(3) 目标运动状态误差统计量估计, 利用参考窗口Rk, 根据式(8)计算误差量, 并采用最大似然原理估计参考窗口中的统计量μ, σ2;

(4) 利用步骤(3)中的误差统计量和机动检测窗口Tk, 计算检测窗口误差统计量μ-, 根据式(11)对目标进行机动检测判断。

检测算法流程如图1所示。

需要说明的是, 本检测算法以目标滤波跟踪数据为输入, 通过运动参数估计与误差统计进行目标检测, 每一次检测与后续检测之间具有相互独立性, 降低了前一次检测过程中目标运动状态估计不准对后续检测的影响; 同时, 尽管在实际运动过程中很难保证目标不机动, 但文中检测过程实质是以窗口中的一段数据作为基准, 对另一段数据进行检测, 当测量数据率较高且窗口不是很长时, 可近似认为参考段数据所反映目标运动没有机动。 故而, 此检测方法具有一定的适用性。

3 仿真验证与分析

采用仿真分析方法验证本文设计目标机动检测方法的有效性。

3.1 仿真场景设置

假定观测设备配置于原点, 测量误差为均值为0、 均方差为20 m的白噪声, 数据率选择为50 Hz, 目标速度为300 m/s, 目标运动过程按照表1中的方式进行机动。 目标运动轨迹与运动状态变化如图2~3所示。

采用经典kalman滤波方法对目标轨迹进行跟踪滤波处理, 跟踪误差如图4所示。

3.2 仿真结果分析

3.2.1 检测性能分析

为验证本文检测方法的检测性能, 选取目标观测数据长度为N=200, 与文献[16]、 文献[20]、 文献[21]、 文献[23]中的方法进行对比仿真分析。 仿真过程中, 选取前半段数据为参考数据, 后半段数据为检测数据, 对目标进行机动检测, 对比检测结果如图5所示。

从图5中的对比仿真可以看出, 本文检测方法与文献中方法均能在目标机动时检测出目标的机动, 但由于受到噪声的影响都存在错误检测和检测延迟的情况。

(1) 错误检测情况分析

为了进一步分析本文所设计检测方法的检测性能, 采用蒙特卡洛仿真方法对错误检测率进行仿真分析。 仿真条件同上, 仿真次数设置为500次, 仿真结果如图6所示。

这里, 错误检测率定义为错误检测累计时间与总时间的比值。 错误检测率仿真统计结果如表2所示。

从图6与表2可以看出, 文献[16]与文献[20]中的方法直接利用滤波跟踪过程新息与其协方差之间的关系进行机动检测, 受到噪声的影响, 错误检测率较高; 文献[21] 中的方法采用小波变换方法进行处理, 可在一定程度上降低错误检测率; 文献[23]与本文方法均基于极大似然原理进行检测, 错误检测率相对较低; 本文在基于极大似然原理基础上引入回归估计可进一步降低错误检测率。

(2) 检测延迟情况分析

检测过程中, 当目标刚开始机动时, 一方面机动带来的偏差相对较小, 另一方面由于噪声的存在, 小的偏差可能会被噪声淹没, 难以快速检测出来, 使得对目标的机动检测总存在一定的延迟, 往往不能在目标刚机动时就检测出目标机动。 为分析本文设计方法对目标机动检测的延迟特性, 采用蒙特卡洛仿真方法进行延迟时间仿真分析。 仿真条件同上, 仿真次数设置为500次, 仿真结果如图7所示。

这里, 延迟时间定义为目标状态发生变化的时间与机动检测出目标状态变化的时间之差。 延迟时间仿真统计结果如表3所示。 文献[16]中的方法直接利用当前滤波新息进行机动检测判断, 若满足条件则认为目标发生机动。 因此, 从检测角度只考虑错误检测的情况, 检测延迟情况未进行分析。

从图7与表3可以看出, 在干扰条件下目标机动检测过程中, 本文方法相比文献中几种方法, 检测延迟时间均有一定提升。 需要说明的是, 检测延迟时间受到窗口宽度的影响, 减小窗口宽度可能进一步缩短延迟时间。

3.2.2 检测性能影响因素分析

从图6与图7可以看出, 由于噪声的存在, 以及从机动开始到位置产生显著偏差需要经历一个过程等因素的影响, 目标机动检测会有一定的延迟, 同时也存在一定的错误检测率。 本文采用蒙特卡洛仿真方法对检测性能影响因素进行分析。

(1) 噪声对检测性能影响分析

采用与3.2.1节中同样的条件, 并选择检测窗口宽度为200, 对噪声影响进行蒙特卡洛仿真分析。 不同噪声环境下的检测性能仿真结果如图8~9所示。

从图8~9的仿真结果可以看出, 随着噪声标准差变大, 其强度变大, 错误检测率呈现上升趋势, 即测量误差越大, 目标越不容易准确检测出来。 通过对比可以看出, 采用基于极大似然原理的本文方法与文献[23]中的方法相比文献[16], [20], [21]在错误检测概率与平均检测延迟时间方面具有更好的性能; 同时, 本文引入回归估计之后可进一步降低噪声对检测的影响, 检测性能相比文献[23]有进一步改善。

(2) 检测窗口宽度对检测性能影响分析

检测窗口宽度描述了检测过程采用的数据的数量。 采用与3.2.1节中同样的条件, 并且选择噪声标准差为30 m, 对检测窗口宽度进行蒙特卡洛仿真分析。 不同检测窗口宽度下的检测性能仿真结果如图10~11所示。

从图10~11的仿真结果可以看出, 随着检测窗口宽度增加, 检测延迟总体变大, 本文检测方法相比文献中方法均有所改善。 当检测窗口过小时, 由于数据量较少, 检测过程受到噪声影响变大, 导致错误检测概率增加, 因此在检测过程中需要合理选择检测窗口宽度。

4 结" 论

针对博弈对抗过程中, 目标意图识别对目标机动提出的高效快速检测需求, 本文设计了一种基于似然比检验的目标机动检测方法。 将目标机动检测问题转化为目标观察时间序列数据检测问题, 通过引入似然比检验原理, 设计检验方法。 该方法利用目标实时观测数据进行机动检测, 仿真结果分析表明, 该方法在噪声环境下的误检测率和检测延迟时间方面具有较好的检测性能。

参考文献:

[1] 孙智孝, 杨晟琦, 朴海音, 等. 未来智能空战发展综述[J]. 航空学报, 2021, 42(8): 525799.

Sun Zhixiao, Yang Shengqi, Piao Haiyin, et al. A Survey of Air Combat Artificial Intelligence[J]. Acta Aeronautica et Astronautica Sinica, 2021, 42(8): 525799.(in Chinese)

[2] 樊会涛, 闫俊. 空战体系的演变及发展趋势[J]. 航空学报, 2022, 43(10): 527397.

Fan Huitao, Yan Jun. Evolution and Development Trend of Air Combat System[J]. Acta Aeronautica et Astronautica Sinica, 2022, 43(10): 527397.(in Chinese)

[3] Yang Q M, Zhang J D, Shi G Q, et al. Maneuver Decision of UAV in Short-Range Air Combat Based on Deep Reinforcement Learning[J]. IEEE Access, 2020, 8: 363-378.

[4] 纪毅, 王伟, 张宏岩, 等. 面向高机动目标拦截任务的空空导弹制导方法综述[J]. 航空兵器, 2022, 29(6): 15-25.

Ji Yi, Wang Wei, Zhang Hongyan, et al. A Survey on Guidance Method of Air-to-Air Missiles Facing High Maneuvering Targets[J]. Aero Weaponry, 2022, 29(6): 15-25.(in Chinese)

[5] Berntorp K. Online Bayesian Inference and Learning of Gaussian-Process State-Space Models[J]. Automatica, 2021, 129: 109613.

[6] 赵智伟, 袁伟伟, 关东海. 基于Multivariate BiLSTM-FCNs的机动实时识别方法[J]. 航空计算技术, 2023, 53(2): 65-69.

Zhao Zhiwei, Yuan Weiwei, Guan Donghai. A Real-Time Identification Method for Maneuvers Based on Multivariate BiLSTM-FCNs[J]. Aeronautical Computing Technique, 2023, 53(2): 65-69.(in Chinese)

[7] Li M J, Zhou C J, Shao L, et al. Intelligent Trajectory Prediction Algorithm for Reentry Glide Target Based on Intention Inference[J]. Applied Sciences, 2022, 12(21): 10796.

[8] 乔殿峰, 梁彦, 张会霞, 等. 具有自动回溯的机动目标航迹精细化分段识别[J]. 航空学报, 2021, 42(4): 524744.

Qiao Dianfeng, Liang Yan, Zhang Huixia, et al. Automatic Backtracking-Based Refined Segment Recognition of Maneuvering Target Track[J]. Acta Aeronautica et Astronautica Sinica, 2021, 42(4): 524744.(in Chinese)

[9] 吴家湖, 熊华, 宗睿, 等. 基于循环神经网络的目标转弯机动类型识别[J]. 广东工业大学学报, 2020, 37(2): 67-73.

Wu Jiahu, Xiong Hua, Zong Rui, et al. Target Turning Maneuver Type Recognition Based on Recurrent Neural Networks[J]. Journal of Guangdong University of Technology, 2020, 37(2): 67-73.(in Chinese)

[10] 罗艺, 谭贤四, 王红, 等. 一种高超声速飞行器攻击意图预测方法[J]. 西安电子科技大学学报, 2019, 46(5): 113-119.

Luo Yi, Tan Xiansi, Wang Hong, et al. Method for Predicting the Attack Intention of Hypersonic Vehicles[J]. Journal of Xidian University, 2019, 46(5): 113-119.(in Chinese)

[11] 张凯, 熊家军. 高超声速滑翔飞行器长期轨迹预测问题探讨[J]. 战术导弹技术, 2018(4): 13-17.

Zhang Kai, Xiong Jiajun. Discussion on Long-Term Trajectory Prediction of Hypersonic Gliding Reentry Vehicle[J]. Tactical Missile Technology, 2018(4): 13-17.(in Chinese)

[12] 张凯, 熊家军, 李凡, 等. 基于意图推断的高超声速滑翔目标贝叶斯轨迹预测[J]. 宇航学报, 2018, 39(11): 1258-1265.

Zhang Kai, Xiong Jiajun, Li Fan, et al. Bayesian Trajectory Prediction for a Hypersonic Gliding Reentry Vehicle Based on Intent Inference[J]. Journal of Astronautics, 2018, 39(11): 1258-1265.(in Chinese)

[13] 李世杰. 临近空间高超声速滑翔目标跟踪及轨迹预测方法研究[D]. 西安: 空军工程大学, 2020.

Li Shijie. Research on Tracking and Trajectory Prediction Technology for Near Space Hypersonic Glide Vehicles[D]. Xi’an: Air Force Engineering University, 2020. (in Chinese)

[14] Hu Y D, Gao C S, Li J L, et al. Novel Trajectory Prediction Algorithms for Hypersonic Gliding Vehicles Based on Maneuver Mode on-Line Identification and Intent Inference[J]. Measurement Science and Technology, 2021, 32(11): 115012.

[15] Li M J, Zhou C J, Shao L, et al. An Intelligent Trajectory Prediction Algorithm for Hypersonic Glide Targets Based on Maneuver Mode Identification[J]. International Journal of Aerospace Engineering, 2022, 2022: 4625001.

[16] 张娜, 王锐, 蔡炯. 基于机动检测的参数自适应跟踪算法[J]. 信号处理, 2022, 38(2): 367-374.

Zhang Na, Wang Rui, Cai Jiong. Parameter Adaptive Tracking Algorithm Based on Maneuver Detection[J]. Journal of Signal Processing, 2022, 38(2): 367-374.(in Chinese)

[17] 范红旗, 王胜, 付强. 目标机动检测算法综述[J]. 系统工程与电子技术, 2009, 31(5): 1064-1070.

Fan Hongqi, Wang Sheng, Fu Qiang. Survey of Algorithms of Target Maneuver Detection[J]. Systems Engineering and Electronics, 2009, 31(5): 1064-1070.(in Chinese)

[18] Karsaz A, Khaloozadeh H, Pariz N, et al. A New Algorithm Based on Generalized Target Maneuver Detection[C]∥ IEEE International Conference on Control and Automation, 2007: 3034-3039.

[19] Xu Y, Fang Y W, Wu Y L, et al. A Weight-Based Multi-Innovation Intelligent Tracking Algorithm with Maneuver Detection[C]∥" 38th Chinese Control Conference, 2019: 3367-3372.

[20] Chan Y T, Couture F. Manoeuvre Detection and Track Correction by Input Estimation[J]. IEE Proceedings F-Radar and Signal Processing, 1993, 140(1): 21-28.

[21] 苏建敏, 董云峰. 利用小波分析在轨识别空间目标轨道机动[J]. 空间科学学报, 2012, 32(3): 412-416.

Su Jianmin, Dong Yunfeng. Detection of Space Target Orbit Maneuver on Board by Wavelet Analysis[J]. Chinese Journal of Space Science, 2012, 32(3): 412-416.(in Chinese)

[22] 崔红正, 刘文玲, 唐歌实, 等. 不同推力下的非合作空间目标轨道机动检测[J]. 宇航学报, 2016, 37(3): 253-261.

Cui Hongzheng, Liu Wenling, Tang Geshi, et al. Different Thrust Maneuvers Detection of Uncooperative Space Objects[J]. Journal of Astronautics, 2016, 37(3): 253-261.(in Chinese)

[23] 王庆瑞, 邹江威, 吴文振, 等. 基于Neyman-Pearson准则的空间目标轨道机动检测方法[J]. 中国空间科学技术, 2021, 41(2): 96-103.

Wang Qingrui, Zou Jiangwei, Wu Wenzhen, et al. Orbital Maneuver Detection Method of Space Target Based on Neyman-Pearson Criterion[J]. Chinese Space Science and Technology, 2021, 41(2): 96-103.(in Chinese)

[24] 王格格, 鲁钰, 吴刘仓. 偏正态数据下众数混合专家回归模型的参数估计[J]. 应用数学, 2021, 34(4): 929-939.

Wang Gege, Lu Yu, Wu Liucang. The Parameter Estimation of the Mixture of Expert Regression for Mode Models with Skew-Normal Data[J]. Mathematica Applicata, 2021, 34(4): 929-939.(in Chinese)

[25] 桑会平. 基于线性模型广义似然比检验的突发信号检测方法[J]. 通信技术, 2023, 56(9): 1051-1055

Sang Huiping. Burst Signal Detection Method Based on Generali-zed Likelihood Ratio Test of Linear Model[J]. Communications Technology, 2023, 56(9): 1051-1055.(in Chinese)

Target Maneuver Detection Method Based on Likelihood Ratio Test

Shao Lei1*, He Yangchao1, Zhao Jin2

(1. Air and Missile Defense College, Air Force Engineering University, Xi’an 710051, China;

2. The Institute of Huanghe Group, Xi’an 710043, China)

Abstract: Aiming at the target maneuver detection in antagonistic process, a detection method is proposed based on likelihood ratio test. Based on the target real time observation data, the target detection problem is translated into the problem of time series data validation. Firstly, the maximum likelihood principle is applied to target regression state estimation. Then, the detection model is designed based on the principle of likelihood ratio detection, and the target maneuver detection method is constructed. Simulation result shows the improved misdetection and detection delay performance of the proposed method in the interference environment.

Key words: maneuver detection; likelihood ratio; regression estimation; time series estimation