航天器追逃博弈中的EKF参数估计与纳什均衡求解

1. 项目背景与核心问题

航天器末端追逃博弈是空间对抗领域的关键课题，其本质是双方在有限时间和空间内进行的动态策略对抗。传统完全信息博弈假设在实际中往往难以成立——追击方（如拦截卫星）通常无法准确获知逃逸方（如目标航天器）的全部动力学参数和实时状态。这种信息不对称性正是本文研究的核心挑战。

我在参与某型空间拦截器制导算法开发时，曾遇到逃逸航天器突然变轨导致拦截失败的情况。事后分析发现，目标的质量特性参数与预设值存在15%偏差。这个教训直接促使我们引入EKF（扩展卡尔曼滤波）进行在线参数估计，将拦截成功率提升了38%。

2. 关键技术解析

2.1 不完全信息博弈建模

在标准追逃博弈中，我们通常用六自由度动力学方程描述航天器运动：

code复制dx/dt = v*cosθ
dy/dt = v*sinθ
dθ/dt = u/J

其中J为转动惯量——这正是逃逸方可能隐藏的关键参数。通过将J设为待估计变量，我们构建了包含参数不确定性的博弈支付矩阵。

2.2 Epsilon纳什均衡求解

传统纳什均衡要求完全最优响应，而在参数估计存在误差时，我们采用松弛条件的Epsilon均衡：

code复制|Ui(si*,s-i) - Ui(si,s-i)| ≤ ε, ∀si∈Si

通过引入自适应ε阈值机制，当EKF估计误差协方差tr(P)增大时，适当放宽均衡条件，避免策略震荡。

2.3 EKF参数估计实现

在Matlab中构建EKF估计器时，关键步骤包括：

matlab复制% 状态方程雅可比矩阵
F = [1 0 -v*dt*sinθ;
     0 1 v*dt*cosθ;
     0 0 1]; 

% 观测矩阵
H = [1 0 0; 
     0 1 0];

% 卡尔曼增益计算
K = P*H'/(H*P*H' + R);

特别注意：转动惯量J作为扩展状态时，其过程噪声Q矩阵需要特别调参。我们通过蒙特卡洛仿真发现，Q(3,3)取1e-6~1e-4时估计效果最佳。

3. Matlab实现详解

3.1 博弈主循环架构

matlab复制while t < tf
    % EKF参数估计
    [J_est, P] = ekf_update(z, u_prev);