强跟踪-容积卡尔曼滤波在弹道式再入目标跟踪中的应用

强跟踪-容积卡尔曼滤波在弹道式再入目标跟踪中的应用

张龙;崔乃刚;王小刚;白俞亮

【摘要】对于具有一定机动能力的弹道式再入目标跟踪问题,稳定性好、鲁棒性强、收敛精度高的估计方法是保证跟踪精度的关键.针对再入运动模型和测量体制的强

非线性以及目标机动引起的滤波精度下降问题,提出一种将强跟踪滤波(STF)和基于三阶球面-向径容积规则的容积卡尔曼滤波(CKF)相结合的强跟踪-容积卡尔曼滤波(STCKF).通过将强跟踪算法中的自适应渐消因子引入到滤波时间更新和测量更新方程中,在线实时调整滤波增益矩阵,能有效避免模型失准造成的滤波性能下降,使该算法兼具CKF滤波精度高和STF鲁棒性强的优点.通过数学仿真表明,改进后的STCKF可以实现对具有机动的弹道式再入目标的高精度跟踪,相对于CKF精度提高50％,并且具有更强的鲁棒性和自适应能力.

【期刊名称】《中国惯性技术学报》

【年(卷),期】2015(023)002

【总页数】8页(P211-218)

【关键词】弹道式再入目标跟踪;容积卡尔曼滤波;自适应渐消因子;非线性系统

【作者】张龙;崔乃刚;王小刚;白俞亮

【作者单位】哈尔滨工业大学航天工程系,哈尔滨150001;哈尔滨工业大学航天工

程系,哈尔滨150001;哈尔滨工业大学航天工程系,哈尔滨150001;哈尔滨工业大学

航天工程系,哈尔滨150001

【正文语种】中文

【中图分类】V557

弹道式再入目标的跟踪是攻防对抗体系中的关键环节，有效的跟踪可以延长防御方的反应时间和提高拦截能力。再入稠密大气的弹道式目标具有速度快、飞行时间短、运动环境复杂的特点，降低了再入目标的跟踪精度。同时目标再入过程中为突防进行的机动未知，无法准确建模，必然会造成在一段时间内跟踪模型失准，进一步增加跟踪的难度。因此，开发鲁棒性强、稳定性好、精度高的估计算法是解决该问题的关键[1]。

理论上，在先验信息（如弹道系数、目标机动大小和时间）充分的前提下，可以对再入运动模型准确建模。但通常情况下，防御方无法获取目标弹道系数的真值，可将弹道系数作为状态量进行联合估计，弹道系数初值可借助一定的先验知识结合已搜集到的情报数据库来确定[2]。但是，目标机动大小、时间和机动形式不能通过

先验信息获取，因此无法建立准确的目标机动模型和应用多模型跟踪方法[3]。

由于目标再入运动方程和雷达对目标的测量方程均是待估状态量的非线性方程，再入跟踪系统呈现出较强的非线性特性。传统的解决非线性实时跟踪滤波问题的方法是扩展卡尔曼滤波（EKF）。该方法采用泰勒级数展开近似非线性函数，引入了线性化误差，对于非线性强的系统精度较低[4]。不敏卡尔曼滤波（UKF）利用UT变换计算采样点，直接对状态的概率分布近似[5]，可得到比EKF高的精度。但对于

高维（维数大于3）非线性系统，该算法在滤波过程中可能出现协方差非正定情况，导致滤波数值不稳定甚至发散，同时易出现非局部效应，严重影响滤波精度[6]。

粒子滤波（PF）基于Monte Carlo采样策略，通过大量的随机粒子逼近概率分布函数，可获得较高的滤波精度。但是由于粒子数较多，造成计算量较大，且易出现粒子退化问题[7]。容积卡尔曼滤波（CKF）采用三阶球面-相径容积规则来近似经非线性函数传递的后验均值和协方差[8]。与UKF相比，CKF算法中各容积点的权值均为正，不会出现协方差非正定的情况，且数值稳定性好，同时计算量远小于

PF算法，适用于再入目标跟踪问题。

文献[9-11]对几种弹道式再入目标的跟踪滤波方法作了对比分析。文献[9]在建立

二维非线性目标再入模型中，将弹道系数作为一个已知的先验信息来考虑，并对比了EKF、UKF、PF、CADET四种滤波方法。Ming Xin[12]针对高超声速再入飞行器的状态估计问题，在Bayesian滤波框架的基础上，采用任意自由度容积规则计算采样点积分，推导了五阶球面-相径容积规则的ICKF模型，并与EKF、UKF、三阶CKF进行对比，获得了更高的估计精度，但是由于采用更高自由度的积分规则，相对于三阶CKF计算量增大了几乎n（状态维数）倍。文献[13]基于三阶球面-相

径积分规则，结合强跟踪滤波STF算法，推导出强跟踪容积滤波算法，研究了初

值不确定并且噪声方差不准确的GPS/INS组合导航问题，有效提高了跟踪精度和系统的鲁棒性。文献[14]为了解决无源传感器机动目标跟踪系统非线性较强、传统的跟踪滤波方法不稳定容易发散的缺陷，在Quadrature 卡尔曼滤波（QKF）基

础上提出了一种带渐消因子的 QKF（ FQKF）算法。该算法通过引入时变渐消因

子来实时调整误差协方差阵，对状态传播积分点和量测传播积分点进行渐消，达到实时调整滤波器增益矩阵的目的，具有STF的优良性能，能够克服QKF算法的缺陷，对于无源传感器机动目标跟踪中系统的突变状态具有较强的跟踪能力，但是基于Gauss-Hermite积分规则对于非线性系统的积分点个数呈指数增长，易造成维数灾难。

本文针对弹道系数未知、带有一定机动的弹道式再入目标跟踪问题，建立再入目标运动模型和雷达测量模型，将弹道系统作为状态量进行实时估计。建立非线性系统的STF算法，结合三阶球面-相径积分的CKF算法，设计了带有自适应渐消因子的强跟踪容积卡尔曼滤波算法STCKF。仿真结果表明，针对弹道式再入目标跟踪问题，本文提出的改进CKF算法的滤波性能优于EKF、CKF算法。

1.1 再入目标运动模型

弹道式目标再入时，受地球引力和空气动力作用。考虑到弹道式目标一般保持零度攻角再入，因此受到的空气动力表现为大气阻力，大气阻力加速度方向与再入速度方向相反。为便于描述再入目标相对地面雷达站的运动关系，在地面雷达站坐标系下建立再入目标的运动方程。定义雷达站坐标系o-xyz的原点o位于雷达站，ox

轴指向北方，oy轴垂直地面指向上方，oz轴指向东方，且与ox、oy轴构成右手坐标系。由于再入时间较短，可忽略地球自转角速度的影响。对再入目标建模如图1所示，图中ag、ad分别是引力加速度和气动力加速度。

假设地球为标准球体，目标的位置和速度分量为[x,y,z,vx,vy,vz]T，则弹道式再入

目标的运动方程如下：

式中：μ = 3.986×1014m3/s2为地球重力常数，Re=为目标的地心距，

Re=6371.11 km为地球平均半径，为目标的再入速度，ρ为大气密度，CD为气

动阻力系数，A为有效面积，m为目标质量。定义弹道系数β，使

由于式(2)中三个参数均事先无法获取，在滤波过程中需要对弹道系数建模。文献[15]对弹道系数模型进行了系统的总结，不同的建模方法对非线性次优滤波算法的性能有直接影响。为了保证弹道系数的非负性，防止滤波器发散，本文选用指数模型，即

式中：0β为弹道系数初值，具体数值可根据经验选定。将α的变化率用零均值高斯白噪声αω表示，即

综上，选取目标的位置、速度和弹道系数参数α作为状态量x=[x,y,z,vx,vy,vz,α]T，建立离散的状态方程：

式中：ΔT为仿真步长；wk为系统噪声向量，满足N(wk;0,Qk)的统计特性。

注意到，式(1)和(5)中均没有考虑目标机动加速度aj，因此当目标发生机动时，易出现模型失准。

1.2 测量模型