基于无人机跟踪场景的智能传感器

(整期优先)网络出版时间:2022-06-09
/ 3


基于无人机跟踪场景的智能传感器

管理算法

卢开钦 、于冬洋、吕途、 王猛

(郑州航天电子技术有限公司,河南省郑州市高新区,450001)

摘要 针对无人机跟踪场景下,由于目标风险评估的不确定性对无人机造成潜在威胁的问题,提出一种面向风险控制的自适应传感器管理算法。首先,通过改进威胁评估算法,提出利用动态模糊贝叶斯网络的方法对目标威胁风险进行评估;其次,利用目标预测协方差对传感器的跟踪风险进行评估;最后,通过融合跟踪风险以及威胁风险构建效能函数,以风险最小化为优化准则,合理地分配传感器资源。仿真及分析表明,所提算法有效地控制无人机作战风险,提高了对目标的跟踪能力。

关键词 无人机跟踪 风险控制 传感器管理 效能函数

1.引言

在现代化信息战领域,利用无人机实施对目标的精准稳定跟踪,对目标态势信息进行准确估计,有助于完成目标预警、监视、定位、拦截和打击等任务需求[1]。传感器作为无人机的“眼睛”,在任务执行中发挥着无与伦比的作用。但是,由于单一传感器受到探测范围、环境、自身性能等限制,无法满足对机动目标的连续性跟踪,为了提高跟踪性能,需要合理调度多传感器来协同探测。但是整个跟踪系统的传感器资源往往有限,如何合理调度有限的传感器资源满足目标跟踪要求,是传感器管理研究的意义所在[2]

传感器管理技术的研究最早由于军事需要而诞生,随着科技的不断发展,仅仅提高目标跟踪精度已经无法满足任务需求,使得传感器管理涌现出了新的问题:(1)当任务不要求打击目标,仅保证传感器不丢失目标即可,此时过于追求跟踪精度就会大大浪费传感器资源;(2)随着目标武器装备不断升级,在跟踪的过程中,由于传感器辐射功率,己方无人机承担着来自目标的威胁风险也不断增强,给目标跟踪带来了极大的挑战。针对风险控制的优化问题具有良好的应用前景,也是本文重点研究内容。

针对上述问题,本文研究内容将关注重点放在了控制作战风险上,运用贝叶斯理论、模糊推理、状态估计等技术手段,对跟踪过程中目标威胁风险进行预测评估。此外,结合状态估计协方差对目标跟踪风险进行预测,进而融合两者构建自适应的效能函数,对由于传感器决策所造成的潜在风险进行控制,并以风险最小化为优化准则,实现传感器资源的合理分配。

2. 目标威胁风险评估策略

2.1目标威胁评估模型构建

目标威胁风险评估是根据跟踪目标的特征,确定该目标对我方无人机构成的威胁程度。在实际跟踪环境中,由于信息保密性,只有结合传感器探测设备的量测信息,以及已有的目标信息进行估计[3]

62a182bc97166_html_7db2c6932cfd3765.gif

1 目标威胁风险评估模型

本文构建目标威胁评估模型[4]如图1所示,从目标威胁能力、打击意图、攻击机会三个方面对目标威胁不确定性进行评估。威胁能力、打击意图、攻击机会属于一级属性。机载武器、干扰能力、目标类型、航向角、目标高度、速度、距离属于二级属性,可以通过数据预处理和目标探测评估得到。

2.2目标威胁风险评估算法

本文在模糊推理、贝叶斯网络以及多源信息融合等理论基础上,通过改进目标威胁评估策略,提出动态模糊贝叶斯网络的方法评估目标威胁风险。算法的具体过程如框图2所示。

62a182bc97166_html_3f233268c9ab6f21.gif

62a182bc97166_html_6c8389580449565d.gif62a182bc97166_html_ff82d837f3fd52b1.gif62a182bc97166_html_11ce5d379b2b582a.gif62a182bc97166_html_c1085fc7f55ecec.gif62a182bc97166_html_e2a8d17c135ba4a9.gif62a182bc97166_html_17295ba246c8cc77.gif62a182bc97166_html_bcea604db332c3b4.gif62a182bc97166_html_6828efd029563ed0.gif62a182bc97166_html_7b30a62b12501210.gif2 动态模糊贝叶斯网络算法框图

如图所示,算法的具体内容如下:

1)目标二级属性值量化

结合实际细化需求对定性因素进行量化处理。本文对定性属性的量化过程如表1所示。

1 定性属性量化过程


机载武器

其它设备

0

常规武器

0.5

导弹

1


干扰能力

1/3

2/3

1


目标类型

地面小型目标

1/3

地面大型目标

2/3

空中目标

1

对于定量属性,结合线性映射函数对其进行具体量化,使归一化后的属性值越大,对威胁值的贡献越大,定量因素量化过程如表2所示。

2 定量属性量化过程

属性因素

归一化属性函数

自变量

航向角(62a182bc97166_html_2754db69ed6903c1.gif

62a182bc97166_html_c6d4ef3f1f8a6155.gif

62a182bc97166_html_29f98167555f5dd0.gif62a182bc97166_html_15489cce6c488b7f.gif

高度(m

62a182bc97166_html_7ebcf4570edb6693.gif

62a182bc97166_html_d2f0129fecb664b2.gif62a182bc97166_html_83df6dc0ba266d19.gif

速度 (m/s)

62a182bc97166_html_2cf3158cb2c370fa.gif

62a182bc97166_html_3e5da154552710b7.gif62a182bc97166_html_e26571a02e12a967.gif

距离(km)

62a182bc97166_html_d54a89bb57116dab.gif

62a182bc97166_html_9488f0463bab07f3.gif62a182bc97166_html_771bcedcd2768fc5.gif

2)利用模糊隶属度函数计算一级属性

根据上述量测信息计算出的威胁度影响因素值62a182bc97166_html_7fb64a3f3c7fe692.gif ,从而计算出隐含层变量节点威胁能力62a182bc97166_html_d98b9451328cdff5.gif 、打击意图62a182bc97166_html_44621dafd80181ef.gif 、攻击机会62a182bc97166_html_c0f26aff3366fc38.gif 。由于影响因素之间相互独立,因此采用以下高斯隶属度函数结合影响因素初始值求解隐含层节点概率[6],即:

威胁能力分为强、中、弱三个等级,由影响因素62a182bc97166_html_c315c767f8759407.gif 融合获得,则威胁能力隶属度函数如下:

62a182bc97166_html_24561b569561b8cb.gif(1)

其中,62a182bc97166_html_8589bf0aad590d54.gif ,因此,62a182bc97166_html_6e6f2e846ee986fb.gif62a182bc97166_html_426ef09449e3271a.gif 为标准差,62a182bc97166_html_6e07d09d7838caac.gif 为各影响因素所对应的权重。

打击意图模糊子集分为大、中、小三个等级。由影响因素62a182bc97166_html_3153a5bc5f74debb.gif 融合而成,同理,则打击意图隶属度函数如下:

62a182bc97166_html_6d3ddf399568bce8.gif(2)

其中,62a182bc97166_html_ed26c2c848ed0897.gif ,因此,62a182bc97166_html_e17a4bc8c03c8719.gif62a182bc97166_html_67fc6202934c6aa1.gif 为标准差,62a182bc97166_html_5d63e961ff080237.gif 为各影响因素所对应的权重。

攻击机会模糊子集分为大、中、小三个等级。由影响因素62a182bc97166_html_87493ccb69ffdbd.gif 融合而成,同理,则攻击机会隶属度函数如下:

62a182bc97166_html_1e30c278cbabdd2b.gif(3)

其中,62a182bc97166_html_938007cda041468a.gif ,因此,62a182bc97166_html_ef9fee923c8f8d57.gif62a182bc97166_html_eeea01c46548ff9e.gif 为标准差,62a182bc97166_html_cf27addd42756224.gif 为各影响因素所对应的权重。

由于需要采用贝叶斯网络推理的方法获得最终的威胁值,因此需要将确定度转换为概率,其具体计算公式如下:

62a182bc97166_html_d6292ddf20ac930f.gif(4)

62a182bc97166_html_22a06717d3a3a0f8.gif(5)

式中,62a182bc97166_html_5af829e96116cd57.gif 代表威胁能力、打击意图、攻击机会等不同的属性,62a182bc97166_html_8685255aeb71f105.gif 代表属性下对应的等级,62a182bc97166_html_e0cdfb1470aef23f.gif 为一致性参数,其范围为62a182bc97166_html_befd1ce98a010f1b.gif ,本文取62a182bc97166_html_34484cc82e09ba61.gif

3)利用贝叶斯推理求解威胁值

在基于贝叶斯推理的目标威胁评估过程中,首先根据目标属性间的因果关系,构建贝叶斯网络拓扑结构[7],如图3所示。

62a182bc97166_html_bf62c2e43a5db656.gif

3 贝叶斯推理网络

在图3中,节点62a182bc97166_html_84cd3eee1a93a3ee.gif 表示威胁度变量,隐含层节点62a182bc97166_html_3d462aa4c28946.gif62a182bc97166_html_4dda476be47f25e2.gif62a182bc97166_html_a2dccb4dbd11a8ef.gif 表示一级属性状态变量,节点62a182bc97166_html_28e39043c165dca1.gif 表示二级属性(观测变量)。隐含层节点变量和观测节点变量体现为威胁评估相关的所有因素。

假设把目标威胁分为62a182bc97166_html_faf1e9e7b287d3c1.gif 个等级:62a182bc97166_html_8bbb826421e8273.gif ,经过贝叶斯网络推理,获得推理结果:

62a182bc97166_html_61ecc5b6c42fc731.gif(6)

其中,62a182bc97166_html_e904086027cc4ab2.gif 表示威胁等级属于62a182bc97166_html_b1d5c66cf466827d.gif 的概率,计算方法如下:

62a182bc97166_html_f6b05946f8aa4638.gif(7)

62a182bc97166_html_580b4aa4ef524c13.gif (8)

其中,62a182bc97166_html_a529870737d6f159.gif 表示隐含层节点变量个数,62a182bc97166_html_229a740b6c529b77.gif 表示隐含层节点到威胁值的条件概率。

根据上述推论,设目标各威胁等级上的威胁权重值为62a182bc97166_html_705694b4969ec392.gif ,则最终的威胁值为:

62a182bc97166_html_cf2d6119af66540a.gif (9)

4)传感器功率截获概率计算

所谓功率截获概率就是电磁探测系统在给定虚警概率62a182bc97166_html_975cd371d89bba7d.gif 下,对传感器辐射信号的探测概率[8]。计算方法如下:

  1. 计算目标探测系统接收机接收到的传感器信号辐射功率,即

62a182bc97166_html_f73745114322970a.gif(10)

式中,62a182bc97166_html_548b4c5e38274327.gif 为目标电磁探测系统接收机接收的辐射功率,62a182bc97166_html_caeb2d90e5c994a1.gif 为无人机传感器平台辐射脉冲的峰值功率,62a182bc97166_html_c188e72d51941ee5.gif 为传感器工作波长,62a182bc97166_html_3420fcfcdd3551b.gif 为无人机与目标之间的距离,62a182bc97166_html_ee9f6758dbb3c192.gif 为目标传感器发射天线增益,62a182bc97166_html_859364c2e7771c5c.gif 为电磁无源探测系统接收增益,62a182bc97166_html_c8648ba672f23446.gif 为电磁无源探测系统接收机处理器净增益。

  1. 计算目标电磁探测系统接收机输出端的信噪比,即

62a182bc97166_html_dd882bcb7d7c3108.gif(11)

式中,62a182bc97166_html_d6885a5c290fc745.gif 为目标电磁探测系统接收机输出端噪声,62a182bc97166_html_b4a12357d29d2e7c.gif 为噪声系数,62a182bc97166_html_699a9a661e327f2e.gif 为玻尔兹曼常数,其值为62a182bc97166_html_6cc8ed4b956fe57b.gif ,62a182bc97166_html_b25150b6c826c8e0.gif 为噪声温度,62a182bc97166_html_1a15193c3e11db00.gif 为电磁无源探测系统接收机带宽。

(c)计算最终的功率截获概率。

在给定虚警概率62a182bc97166_html_351725b08e08b82e.gif ,可计算出功率截获概率为:

62a182bc97166_html_2035b45ec4b0c723.gif (12)

62a182bc97166_html_1dc8f3d4ea05b213.gif(13)

根据传感器功率截获概率,结合目标威胁值,可得传感器平台受到目标最终的威胁风险为:

62a182bc97166_html_fd8ce3aced07789e.gif (14)

62a182bc97166_html_5f55e94cebc84eb6.gif (15)


3 目标跟踪风险评估策略

目标跟踪风险,即在任务执行过程中,由于传感器决策造成的目标失跟、误跟的风险。对于目标跟踪风险评估的关键是对目标状态估计误差的评估,计算方式如下:

62a182bc97166_html_fc40f7fb061056b4.gif(16)

其中,62a182bc97166_html_f6ed0d0b9680d709.gif62a182bc97166_html_2600089216e318fb.gif 时刻传感器控制序列为62a182bc97166_html_5d7d50adcc5c4c2c.gif 时的状态预测误差协方差,由滤波算法求得,在此不做赘述。62a182bc97166_html_3d074c91e6cfdac0.gif62a182bc97166_html_f7c5b13a65fb5131.gif 时刻传感器的量测噪声方差,62a182bc97166_html_704205eef2fccfac.gif 表示求矩阵的迹,62a182bc97166_html_4cbf85c7d7ad5935.gif 为量测噪声系数。此式的物理含义为在当前时刻的量测噪声下,对应传感器对目标的状态估计误差,误差越大,目标跟踪风险越高。

4效能函数的构建

本文结合上述跟踪风险以及威胁风险的评估策略,以风险控制为原则,构建效能函数。算法的核心在于对未来时刻由于任务决策带来的潜在风险进行预测评估,以风险最小化对传感器资源进行合理地分配。风险的一步预测值计算方法如下:

62a182bc97166_html_e6332a3304b38be6.gif(17)

其中,62a182bc97166_html_c0ddaf673209fe6b.gif62a182bc97166_html_38fa168eaa78bc8c.gif62a182bc97166_html_d7c3b37573e16f0a.gif 表示在执行传感器控制策略62a182bc97166_html_9b7bc2632aef508.gif 后,预测在62a182bc97166_html_2a35455172640c4.gif 时刻传感器平台所面临的总风险,跟踪风险以及威胁风险,62a182bc97166_html_3916855f4ecef04c.gif 为效能平衡权重,用于平衡两种风险对总风险的贡献。

进一步,可以得到多步预测总风险计算方法如下:

62a182bc97166_html_ac1a23de26f6d253.gif(18)

式中62a182bc97166_html_7a493a42d3f9f851.gif 为任务决策步长。

本文采用长期管理方式,利用多步预测总风险,构建效能函数如下:

62a182bc97166_html_6b3edada6a8d35a2.gif(19)

该效能函数的最优解62a182bc97166_html_33cab5e34675dc97.gif 表示使得未来时刻62a182bc97166_html_c0bb21ed668af1f0.gif62a182bc97166_html_35aa81c7d134d1a4.gif 总风险最小的传感器控制序列,由于该时间段处于同一传感器控制周期内,因此62a182bc97166_html_a394e84c920e4abd.gif

5 仿真结果与分析

5.1 仿真场景

为了验证本文所提算法的有效性,本次仿真在单机对抗的跟踪场景下,无人机携带三种不同的主动传感器对目标进行跟踪定位,控制周期为10s。无人机探测系统参数如表4所示,初始状态为

62a182bc97166_html_2201631ac0a0f0dd.gif ,单位(km)。目标初始状态为62a182bc97166_html_7aa4401927a671f0.gif ,单位(km)。目标的属性参数值为:常规设备、干扰能力强、空中目标、航向角60°。

仿真时长为1200s,量测噪声系数62a182bc97166_html_6be704b9d1355bd3.gif 。仿真实验采用公式(19)所示的效能函数对传感器资源进行合理地分配,取步长62a182bc97166_html_82d3c9b7a8f52090.gif

5.2 仿真结果

根据公式(19),本文构建多步预测效能函数对未来时刻的潜在风险进行预测,总风险的预测结果如图5所示。

62a182bc97166_html_ec0ca47452ab70ff.jpg

图5预测总风险曲线图

图5为跟踪过程中总风险变化曲线图,图中,风险预测值为通过对目标状态估计预测未来时刻的总风险变化,风险真实值是在执行传感器控制方案后通过真实量测值所计算的总风险变化。由图可以看出,风险预测值与真实值大致相同,验证了本文对目标风险评估的有效性,证明了本文所提算法的可行性。

为了详细说明本文所提管理算法(proposed management approach,PMA)对控制跟踪过程总风险的有效性,结合以下两种常见的管理算法进行对比分析。

(1) 随机管理方法 (random management approach,RMA),随机分配传感器资源执行对目标的跟踪探测,主要用于态势紧急,来不及解算目标函数的情况。

(2) 短期管理方法 (myopic management approach,MMA),根据公式(17),以一步风险预测值为准则,对传感器资源进行管理。

根据上述三种管理方案,对目标跟踪过程中的风险进行预测评估,其结果如图7、图8所示。

62a182bc97166_html_a511aa44232473a0.jpg

7 传感器总风险变化

62a182bc97166_html_cca3ad16c99fb634.gif

8 传感器归一化风险对比

7所示为传感器总风险随时间变化曲线对比图,图8所示为整个跟踪过程中累积总风险、跟踪风险、威胁风险对比图。从图中可以看出,RMA由于传感器分配的随机性,使得作战风险无法得到有效控制,因此三种风险均为最高。MMA根据一步预测效能函数对传感器资源进行分配,在一定程度上可以对风险进行控制,但是相对于整个控制时刻,一步预测的结果具有片面性,而且风险评估准确度不高,无法得到最佳的风险控制效果。本文所提算法,由于通过对一定时间的累积风险进行预测,完成对传感器的管理,使得在整个跟踪时间内的累积风险最低,每一时刻的风险也为最低。说明了本文所提算法能较好地评估跟踪风险和威胁风险,进而控制跟踪过程中的总风险,提高了无人机的跟踪性能。

6 结论

本文以风险最小化为原则,提出一种自适应的传感器管理算法,所提算法改进了威胁评估策略,提出动态模糊贝叶斯网络的方法对目标的威胁风险进行评估,同时融合目标跟踪风险构建多步预测效能函数对传感器资源进行合理地分配,仿真结果表明该算法有效地降低了目标跟踪过程中无人机所面临的任务风险,提高了跟踪性能。

参考文献

  1. 王学佼. 机载平台多机动目标协同跟踪中传感器管理算法研究[D].哈尔滨工程大学,2017.

  2. 高晓光, 李飞, 万开方. 数据丢包环境下的多传感器协同跟踪策略研究[J]. 系统工程与电子技术, 2018, 40(11):63-71.

  3. 孔尚萍. 基于多源信息融合的目标航迹估计与威胁评估[D].中国航天科技集团公司第一研究院,2017.

  4. 朱建益. 空战中的威胁估计与态势评估研究[D].西安电子科技大学,2013.

  5. 庞策,单甘霖,段修生.多传感器协同识别跟踪多目标管理方法研究[J].北京航空航天大学学报:1-8[2019-05-29].

  6. 刘文博,梁敏,张欧亚,代海峰.基于贝叶斯网络的侦察威胁等级评估研究[J].火力与指挥控制,2019,44(03):72-76.

  7. 韩其松,余敏建,高阳阳,宋帅,陈双艳.云模型和距离熵的TOPSIS法空战多目标威胁评估[J].火力与指挥控制,2019,44(04):136-141.

  8. Anton Milan, Konrad Schindler, Stefan Roth. Multi-Target Tracking by Discrete-Continuous Energy Minimization[J]. IEEE Transactions on Pattern Analysis & Machine Intelligence, 2016, 38(10):2054-2068.