来源:市场资讯

(来源:EW Frontier)

【EW Frontier】——雷达通信AI科研人的一站式技术平台!

✅ 超1000+实战代码:DOA/调制识别/ISAC/抗干扰/无人机等雷达、通信、电子战全方向(MATLAB+Python)

✅ 专属科研辅导:论文专利选题/仿真/写作、项目定制全程答疑

✅ 优质社群资源:985/211硕博同行交流,最新技术干货实时同步

资源获取通道

知识星球(全部资源无限看):https://wx.zsxq.com/group/15554455154582

面包多(单个代码精准购):https://mbd.pub/o/EWFrontier/work

辅导/答疑:

‍客服微信: EWFrontier

面向强对抗环境的空空导弹末制导多源信息融合理论框架

摘要

本文针对现代空空导弹在电子对抗环境下对高机动目标的精确跟踪需求,系统阐述了一套多源信息融合制导理论与模型框架。该框架包含目标机动模型、异构传感器观测模型、干扰机理建模、联邦扩展卡尔曼滤波以及制导参数重构等关键理论环节。通过引入 Singer 机动模型、Van Loan 精确离散化方法以及基于信息分配的联邦滤波结构,显著提升了系统在红外诱饵、雷达箔条和数据链降级等复杂干扰下的状态估计鲁棒性。本文重点推导各模块的数学模型,并通过典型实例说明其物理含义与工作机理。

打开网易新闻 查看精彩图片

1 引言

现代空战环境中,目标飞行器普遍装备有先进的自卫电子对抗系统,可投放红外诱饵弹和雷达箔条,并可能对数据链实施干扰。空空导弹的末制导阶段需在极短时间内融合雷达、红外、相机及数据链等多源异质信息,以精确估计目标的运动状态,并生成满足比例导引及其改进形式所需的视线角速率、接近速度和目标加速度等制导指令参数。为此,需要建立一个能够描述目标复杂机动、模拟传感器物理特性与干扰效应、并实现鲁棒信息融合的完整理论框架。本文将围绕上述功能需求,逐层阐释其背后的数学模型与物理原理。

2 目标运动学与动力学建模

2.1 目标机动指令生成

目标的战术机动由机动调度器描述,依据典型任务剖面(如战斗机 “拦截 - 规避”、预警机 “跑道形” 巡航)输出三个独立控制量:切向加速度、水平转弯角速率和爬升率指令。切向加速度沿速度方向改变速率;水平转弯角速率通过产生法向(向心)加速度实现航向变化;爬升率指令驱动高度变化。

2.2 飞行力学积分模型

将控制指令转化为惯性系下的三维运动,采用恒加速度积分模型。记时刻的状态为位置、速度、加速度,时间步长为,则离散递推公式为:

加速度矢量由指令综合产生:

其中为瞬时偏航角,为垂直通道增益。姿态角由速度矢量和转弯率通过运动学关系反算。该模型为仿真提供 “真值” 数据,不包含在滤波器内部。

2.3 滤波器内部预测模型(CV/CA/Singer)

在目标状态估计器中,常采用不同的运动模型来预测目标行为。为统一描述,考虑单轴向状态向量,即位置、速度、加速度。连续时间状态空间模型为:

其中是功率谱密度为的连续白噪声。三种常用模型如下:

常速度模型(CV):将速度建模为随机游走过程,加速度被白噪声直接驱动速度变化,在工程实现中常附加一个快速衰减的加速度状态以维持状态维数一致。

常加速度模型(CA):加加速度(jerk)为白噪声,状态矩阵为

噪声强度。

Singer 模型:加速度为指数自相关的一阶马尔可夫过程,机动时间常数表示目标机动持续时间的期望,机动强度由稳态加速度方差刻画。其状态方程:

相应的矩阵为

2.4 Van Loan 精确离散化方法

滤波器的离散预测步需要精确的离散状态转移矩阵和过程噪声协方差。对于线性时不变系统,可利用矩阵指数进行一次性求解。构造增广矩阵:

计算其矩阵指数,则,,其中表示分块子矩阵。与直接泰勒展开近似相比,Van Loan 方法可保证数值精度,尤其适用于大时间步长或快变动力学。

数值示例:取s,s,,Singer 模型单轴离散转移矩阵:

可见加速度状态在 0.1 s 内仅衰减约 0.5%,反映了 Singer 模型对持续机动的 “记忆” 特性,优于 CA 模型将加速度视为常数、CV 模型完全忽略加速度的做法。

获得单轴矩阵后,利用各轴向运动解耦的假设,通过克罗内克积形式扩展为三维 9×9 矩阵,并排列为的顺序。

3 多传感器观测模型

3.1 雷达观测方程

雷达提供距离、方位角、俯仰角和多普勒速度四个测量量。在球坐标系下,非线性观测方程为:

测量噪声假设为零均值高斯白噪声,功率谱密度矩阵,典型参数:m,,m/s。

3.2 红外传感器

红外(IR)传感器仅提供高精度的角度测量,其方位和俯仰与雷达角度真值相同,但噪声远小于雷达,。观测矢量简化为。

3.3 光学相机

相机采用针孔透视投影模型。给定相机焦距(以像素为单位)、主点坐标,目标在相机坐标系中的坐标为。投影方程为:

仅在目标位于相机前方且像素坐标在有效幅面内时才输出有效测量。测量噪声标准差通常为 1~2 像素。

3.4 数据链

数据链(Data Link)提供了来自预警机等平台的超视距目标指示,输出目标在惯性系下的三维位置和速度。由于信号传播和网络延迟,该测量具备不确定的时延(s,s)以及独立的伯努利丢包过程,丢包率可达 10%。接收端基于报文发送时间戳进行时延补偿,线性外推得到当前时刻近似位置:

测量噪声为标准差m、m/s 的高斯噪声,远劣于本机雷达和光电器件。

4 电子对抗环境建模

4.1 红外诱饵弹(Flare)

诱饵弹从载机尾部弹射,初始速度附加一个向下的弹射分量,随后在重力和空气阻力作用下运动。其动力学方程:

式中为阻力系数(典型值 0.5)。采用线性积分可得到诱饵位置,诱饵因其高辐射强度,在红外 / 相机视场中形成非常明亮的干扰源,可能导致跟踪波门偏移。

4.2 雷达箔条(Chaff)

箔条释放后初始速度与载机相同,之后在大气阻滞下迅速减速并弥散。采用指数衰减模型描述其速度变化:

时间常数通常为几秒量级。箔条云形成大范围的雷达强散射体,其回波功率和多普勒谱展宽远大于典型目标,易使雷达自动检测和跟踪算法发生误锁。

4.3 干扰对传感器数据的影响

在干扰条件下,传感器可能输出虚假测量值以替代真实目标信息。例如,红外 / 相机系统可能锁定最亮的诱饵,将诱饵的角位置作为目标测量值;雷达距离 - 速度门可能被箔条云捕获,输出箔条云的平均距离和展宽后的多普勒。因此,滤波器实际接收到的测量量可能被污染,需要融合算法具备鉴别或抑制虚假信息的能力。

5 联邦扩展卡尔曼滤波融合理论

5.1 联邦滤波总体架构

为处理多源异步、不同精度的传感器信息,并增强对局部传感器故障或欺骗的鲁棒性,采用 Carlson 提出的联邦滤波结构。系统包含四个并行的子滤波器:雷达子滤波器、红外子滤波器、相机子滤波器和数据链子滤波器。每个子滤波器均为扩展卡尔曼滤波器(EKF),独立维护目标的全状态(9 维)和协方差矩阵。

在每个滤波周期开始时,融合中心将全局估计信息按信息分配因子分配给各子滤波器():

表征对各子滤波器的置信度,通常依据传感器精度和抗干扰能力预先设定。例如,近距离高精度雷达和相机可赋予较大权重(),红外次之(),数据链因噪声大、时延不确定而权重最小()。

5.2 子滤波器时间更新(预测)

各子滤波器采用相同的目标运动模型(如 Singer 模型)和对应的离散状态矩阵进行一步预测:

其中和由 Van Loan 方法从连续模型生成。

5.3 子滤波器量测更新

对于每个子滤波器,根据其对应的传感器模型进行 EKF 更新。

雷达更新:利用非线性函数和解析雅可比矩阵,计算新息。角度新息约束在内。随后计算新息协方差,卡尔曼增益,最后修正状态与协方差。

红外与相机更新:红外仅利用子集,相机则利用与对应的雅可比,过程类似。

数据链更新:数据链提供线性观测,观测矩阵为。为抑制数据链大噪声与潜在污染,人为设置极大的测量噪声协方差 ,其中 m, m/s。此举本质是在融合时降低数据链的统计权重,使其只在其他传感器不可用时作为保底参考。

5.4 全局融合

各子滤波器独立完成更新后,融合中心通过逆协方差加权实现全局最优估计:

该融合准则是最大似然意义下的最优组合,且能容纳任意子滤波器失效(只需将其信息矩阵置零)。为数值稳健,实践中对协方差求逆时添加微小正则化项(如)。融合后的全局状态与协方差将作为下一周期的分配起点,形成闭环重置。

抗干扰机理:当雷达遭受箔条欺骗而导致测量误差增大时,其新息增大,估计误差协方差也随之增大;相反的,若相机仍捕获真实目标,其保持较小值。逆协方差融合时,雷达的贡献自动衰减,全局估计将主要依从相机和红外,从而抑制了欺骗干扰的影响,实现了 “软切换” 而非硬判决。

6 制导参数重构

滤波器输出的全局状态包含目标位置、速度和加速度,结合导弹自身状态(由导航系统提供),可直接导出比例导引及其增强形式所需的所有参数。

计算视线矢量及相对运动量:

径向æŽè¿‘速率

视线角速率矢量为核心量:

其模长用于经典比例导引指令:,其中,为导航比。

目标加速度估计直接取自,在增强比例导引(APN)中,垂直于视线的加速度分量被用于补偿目标机动:

剩余飞行时间由近似估计(要求m/s)。

7 结论

本文针对强对抗场景下的空空导弹末制导跟踪,完整论述了从目标运动建模、异构传感器测量方程、对抗物干扰机理,到联邦卡尔曼滤波融合与制导参数重构的理论框架。通过 Van Loan 精确离散化的多模型预测、基于信息分配的分布式融合以及针对数据链的统计降权策略,系统在多个层面提升了估计的准确性与鲁棒性。所建立的数学模型与参数配置可为相关仿真实验和算法设计提供理论基础。

本框架可进一步扩展至交互多模型(IMM)以自适应切换运动模式、引入概率数据关联(PDAF)处理多假目标、以及集成被动测角(ESM)等辅助信息,以适应更为复杂的战场环境。