1. 项目背景与核心问题
航天器末端追逃博弈是空间对抗领域的关键课题,其本质是追踪方与逃逸方在有限时间内的动态策略对抗。传统研究通常假设双方完全掌握对方的动力学参数和控制策略,但在实际太空任务中,这种理想假设往往难以成立。逃逸方可能通过主动干扰、参数伪装或通信屏蔽等手段隐藏真实控制特性,导致博弈进入不完全信息状态。
当追踪方基于错误参数制定策略时,会出现两个典型问题:
- 拦截轨迹偏离实际目标
- 燃料消耗超出预期
最终导致任务失败率上升。本研究的创新点在于将扩展卡尔曼滤波(EKF)与博弈论相结合,通过实时参数估计和策略调整,使系统在信息不对称条件下仍能保持接近最优的拦截性能。
需要模型API调用? 免费领10W Token,多模型网关一键接入 Claude、DeepSeek 等主流模型。
2. 关键技术解析
2.1 航天器相对运动建模
采用Clohessy-Wiltshire(C-W)方程描述近地轨道上的相对运动:
code复制dx'' - 2ωdy' - 3ω²dx = ax
dy'' + 2ωdx' = ay
dz'' + ω²dz = az
其中ω为轨道角速度,(dx,dy,dz)表示相对位置分量,(ax,ay,az)为控制加速度。该模型的优势在于:
- 线性化处理简化了计算复杂度
- 保持足够的精度(适用于相对距离<100km)
- 状态方程形式适合EKF处理
2.2 EKF参数估计实现
将逃逸方的未知控制矩阵B扩展为系统状态变量:
code复制X_hat = [dx, dy, dz, dx', dy', dz', B11, B12,..., B33]^T
EKF预测与更新步骤如下:
-
状态预测:
matlab复制
X_pred = f(X_hat_prev) + w P_pred = F*P_prev*F' + Q -
卡尔曼增益计算:
matlab复制
K = P_pred*H'/(H*P_pred*H' + R) -
状态更新:
matlab复制
X_hat = X_pred + K*(Z - h(X_pred)) P = (I - K*H)*P_pred
关键参数设置经验:
- 过程噪声Q取10^-6量级
- 测量噪声R取10^-2量级
