小结

我们已经推导出了卡尔曼滤波全部的五个方程的矩阵形式,我们现在把它们汇总在一起。

卡尔曼滤波执行的是一个“预测-修正”循环,如下图所示。

Predict-Update Diagram

一旦初始化完成,卡尔曼滤波首先预测下一时刻的系统状态,并且同时给出这个预测的不确定性。

随后拿到测量结果,卡尔曼滤波会更新(或修正)这个预测值及其不确定性,并且同时再进行下一时刻的预测,以此类推。

下图给出了卡尔曼滤波运算的全貌。

The Kalman Filter Diagram

下表给出了卡尔曼滤波全部的方程。

方程 方程名 别名
预测 x^n+1,n=Fx^n,n+Gun 状态外插 预测器方程
转移方程
预测方程
动态模型
状态空间模型
Pn+1,n=FPn,nFT+Q 协方差外插 预测器协方差方程
更新
(修正)
x^n,n=x^n,n1+Kn(znHx^n,n1) 状态更新 滤波方程
Pn,n=(IKnH)Pn,n1(IKnH)T+KnRnKnT 协方差更新 修正方程
Kn=Pn,n1HT(HPn,n1HT+Rn)1 卡尔曼增益 权重方程
辅助方程 zn=Hxn 测量方程
Rn=E(vnvnT) 测量协方差 测量误差
Qn=E(wnwnT) 过程噪声协方差 过程噪声误差
Pn,n=E(enenT)=E((xnx^n,n)(xnx^n,n)T) 估计协方差 估计误差

下表汇总了符号(包括其他文献中的不同记法)和对应维度。

名称 其他记法 维度
x 状态向量 nx×1
z 输出向量 y nz×1
F 状态转移矩阵 Φ,A nx×nx
u 输入向量 nu×1
G 控制矩阵 B nx×nu
P 估计协方差 Σ nx×nx
Q 过程噪声协方差 nx×nx
R 测量协方差 nz×nz
w 过程噪声向量 y nx×1
v 测量噪声向量 nz×1
H 观测矩阵 C nz×nx
K 卡尔曼增益 nx×nz
n 离散时间变量 k

维度注解:

  • nx 是状态向量中的状态个数
  • nz 是测量到的状态个数
  • nu 是输入向量中的元素个数
上一章 下一章