绘制卡尔曼滤波运行流程图

绘制卡尔曼滤波运行流程图

flyfish

2024-06-28
2024-07-01 修正错别字
绘图软件是LibreOffice Draw,里面的数学公式用的是TexMaths插件,语法是LaTeX

下面的文章都是为了理解卡尔曼滤波
积分的可视化

详细的介绍匀加速运动的物理方程是如何转化为卡尔曼滤波的状态空间模型的

(目标跟踪算法)中的状态向量与状态转移矩阵

(目标跟踪算法)中卡尔曼滤波器中的更新

(目标跟踪算法)中卡尔曼增益的理解

(目标跟踪算法)中的卡尔曼滤波 - 看了就会的状态转移矩阵

(目标跟踪算法)中的初始化卡尔曼滤波器的状态向量和协方差矩阵

(目标跟踪算法)中卡尔曼滤波器中的预测

(目标跟踪算法)中自由度决定卡方分布的形状

(目标跟踪算法)中的解三角方程计算标准化残差(解线性方程组)

(目标跟踪算法)中的计算观测值与状态估计的马氏距离

(目标跟踪算法)中的马氏距离详解(很详细)

(目标跟踪算法)中 可以设置阈值进行异常检测或目标跟踪的原因

(目标跟踪算法)中的数值表格与调参的关系

从最小二乘法的角度来理解卡尔曼滤波(1)

从最小二乘法的角度来理解卡尔曼滤波(2)

(目标跟踪算法)卡尔曼滤波中的贝叶斯定理

Determinant 一个翻译很不友好的名字 行列式 det(A)

(目标跟踪算法) 卡尔曼滤波 状态向量是如何映射到观测向量(测量向量)的即观测矩阵的构建方式

(目标跟踪算法) 卡尔曼滤波 状态空间模型的理解

包含公式

在这里插入图片描述

图上卡尔曼滤波的五个主要公式的说明

卡尔曼滤波器由两个阶段组成:预测阶段和更新阶段。每个阶段有一组方程。

1. 初始化

在第一个时刻,需要初始化状态向量和误差协方差矩阵:

  • 状态向量(State Vector) x 0 \mathbf{x}_0 x0:对系统初始状态的估计。
  • 误差协方差矩阵(Error Covariance Matrix) P 0 \mathbf{P}_0 P0:初始状态估计的不确定性。

2. 预测步骤(Prediction Step)

预测系统在下一个时刻的状态和误差协方差矩阵。

  • 预测状态(Predicted State) x ^ k ∣ k − 1 \hat{\mathbf{x}}_{k|k-1} x^kk1 x ^ k ∣ k − 1 = A k − 1 x ^ k − 1 ∣ k − 1 + B k − 1 u k − 1 \hat{\mathbf{x}}_{k|k-1} = \mathbf{A}_{k-1} \hat{\mathbf{x}}_{k-1|k-1} + \mathbf{B}_{k-1} \mathbf{u}_{k-1} x^kk1=Ak1x^k1∣k1+Bk1uk1其中, A k − 1 \mathbf{A}_{k-1} Ak1 是状态转移矩阵, B k − 1 \mathbf{B}_{k-1} Bk1 是控制输入矩阵, u k − 1 \mathbf{u}_{k-1} uk1 是控制输入。
  • 预测误差协方差矩阵(Predicted Error Covariance Matrix) P k ∣ k − 1 \mathbf{P}_{k|k-1} Pkk1 P k ∣ k − 1 = A k − 1 P k − 1 ∣ k − 1 A k − 1 T + Q k − 1 \mathbf{P}_{k|k-1} = \mathbf{A}_{k-1} \mathbf{P}_{k-1|k-1} \mathbf{A}_{k-1}^T + \mathbf{Q}_{k-1} Pkk1=Ak1Pk1∣k1Ak1T+Qk1其中, Q k − 1 \mathbf{Q}_{k-1} Qk1 是过程噪声协方差矩阵。

3. 更新步骤(Update Step)

利用新测量数据更新预测的状态和误差协方差矩阵。

  • 卡尔曼增益(Kalman Gain) K k \mathbf{K}_k Kk K k = P k ∣ k − 1 H k T ( H k P k ∣ k − 1 H k T + R k ) − 1 \mathbf{K}_k = \mathbf{P}_{k|k-1} \mathbf{H}_k^T (\mathbf{H}_k \mathbf{P}_{k|k-1} \mathbf{H}_k^T + \mathbf{R}_k)^{-1} Kk=Pkk1HkT(HkPkk1HkT+Rk)1其中, H k \mathbf{H}_k Hk 是观测矩阵, R k \mathbf{R}_k Rk 是观测噪声协方差矩阵。
  • 更新状态(Updated State) x ^ k ∣ k \hat{\mathbf{x}}_{k|k} x^kk x ^ k ∣ k = x ^ k ∣ k − 1 + K k ( z k − H k x ^ k ∣ k − 1 ) \hat{\mathbf{x}}_{k|k} = \hat{\mathbf{x}}_{k|k-1} + \mathbf{K}_k (\mathbf{z}_k - \mathbf{H}_k \hat{\mathbf{x}}_{k|k-1}) x^kk=x^kk1+Kk(zkHkx^kk1)其中, z k \mathbf{z}_k zk 是测量向量。
  • 更新误差协方差矩阵(Updated Error Covariance Matrix) P k ∣ k \mathbf{P}_{k|k} Pkk P k ∣ k = ( I − K k H k ) P k ∣ k − 1 \mathbf{P}_{k|k} = (\mathbf{I} - \mathbf{K}_k \mathbf{H}_k) \mathbf{P}_{k|k-1} Pkk=(IKkHk)Pkk1其中, I \mathbf{I} I 是单位矩阵。

4. 循环迭代

重复预测步骤和更新步骤,处理每个新的测量数据。

各个字母的含义

1. 预测步骤(Prediction Step)

预测步骤包括状态预测和误差协方差预测。

状态预测(State Prediction)

x ^ k ∣ k − 1 = A k − 1 x k − 1 + B k − 1 u k − 1 \hat{\mathbf{x}}_{k|k-1} = \mathbf{A}_{k-1} \mathbf{x}_{k-1} + \mathbf{B}_{k-1} \mathbf{u}_{k-1} x^kk1=Ak1xk1+Bk1uk1

  • x ^ k ∣ k − 1 \hat{\mathbf{x}}_{k|k-1} x^kk1:预测的状态向量(在时间 k k k 时的先验估计)。
  • A k − 1 \mathbf{A}_{k-1} Ak1:状态转移矩阵,描述系统从时间 k − 1 k-1 k1 到时间 k k k 的动态。
  • x k − 1 \mathbf{x}_{k-1} xk1:在时间 k − 1 k-1 k1 时的状态向量(后验估计)。
  • B k − 1 \mathbf{B}_{k-1} Bk1:控制输入矩阵,描述控制输入对系统状态的影响。
  • u k − 1 \mathbf{u}_{k-1} uk1:在时间 k − 1 k-1 k1 时的控制输入向量。
误差协方差预测(Error Covariance Prediction)

P k ∣ k − 1 = A k − 1 P k − 1 ∣ k − 1 A k − 1 T + Q k − 1 \mathbf{P}_{k|k-1} = \mathbf{A}_{k-1} \mathbf{P}_{k-1|k-1} \mathbf{A}_{k-1}^T + \mathbf{Q}_{k-1} Pkk1=Ak1Pk1∣k1Ak1T+Qk1

  • P k ∣ k − 1 \mathbf{P}_{k|k-1} Pkk1:预测的误差协方差矩阵(在时间 k k k 时的先验估计的不确定性)。
  • P k − 1 ∣ k − 1 \mathbf{P}_{k-1|k-1} Pk1∣k1:在时间 k − 1 k-1 k1 时的误差协方差矩阵(后验估计的不确定性)。
  • A k − 1 \mathbf{A}_{k-1} Ak1:状态转移矩阵。
  • Q k − 1 \mathbf{Q}_{k-1} Qk1:过程噪声协方差矩阵,描述系统过程噪声的不确定性。

2. 更新步骤(Update Step)

更新步骤包括计算卡尔曼增益、更新状态估计和更新误差协方差矩阵。

计算卡尔曼增益(Kalman Gain Calculation)

K k = P k ∣ k − 1 H k T ( H k P k ∣ k − 1 H k T + R k ) − 1 \mathbf{K}_k = \mathbf{P}_{k|k-1} \mathbf{H}_k^T (\mathbf{H}_k \mathbf{P}_{k|k-1} \mathbf{H}_k^T + \mathbf{R}_k)^{-1} Kk=Pkk1HkT(HkPkk1HkT+Rk)1

  • K k \mathbf{K}_k Kk:卡尔曼增益矩阵,描述测量值对状态估计的影响。
  • P k ∣ k − 1 \mathbf{P}_{k|k-1} Pkk1:预测的误差协方差矩阵。
  • H k \mathbf{H}_k Hk:观测矩阵,描述状态向量到测量向量的映射。
  • R k \mathbf{R}_k Rk:测量噪声协方差矩阵,描述测量噪声的不确定性。
更新状态估计(State Update)

x k ∣ k = x ^ k ∣ k − 1 + K k ( z k − H k x ^ k ∣ k − 1 ) \mathbf{x}_{k|k} = \hat{\mathbf{x}}_{k|k-1} + \mathbf{K}_k (\mathbf{z}_k - \mathbf{H}_k \hat{\mathbf{x}}_{k|k-1}) xkk=x^kk1+Kk(zkHkx^kk1)

  • x k ∣ k \mathbf{x}_{k|k} xkk:更新后的状态向量(在时间 k k k 时的后验估计)。
  • x ^ k ∣ k − 1 \hat{\mathbf{x}}_{k|k-1} x^kk1:预测的状态向量。
  • K k \mathbf{K}_k Kk:卡尔曼增益矩阵。
  • z k \mathbf{z}_k zk:在时间 k k k 时的实际测量向量。
  • H k \mathbf{H}_k Hk:观测矩阵。
  • H k x ^ k ∣ k − 1 \mathbf{H}_k \hat{\mathbf{x}}_{k|k-1} Hkx^kk1:预测测量值。
更新误差协方差矩阵(Error Covariance Update)

P k ∣ k = ( I − K k H k ) P k ∣ k − 1 \mathbf{P}_{k|k} = (\mathbf{I} - \mathbf{K}_k \mathbf{H}_k) \mathbf{P}_{k|k-1} Pkk=(IKkHk)Pkk1

  • P k ∣ k \mathbf{P}_{k|k} Pkk:更新后的误差协方差矩阵(在时间 k k k 时的后验估计的不确定性)。
  • I \mathbf{I} I:单位矩阵。
  • K k \mathbf{K}_k Kk:卡尔曼增益矩阵。
  • H k \mathbf{H}_k Hk:观测矩阵。
  • P k ∣ k − 1 \mathbf{P}_{k|k-1} Pkk1:预测的误差协方差矩阵。

不包含公式的版本
在这里插入图片描述

  • 51
    点赞
  • 27
    收藏
    觉得还不错? 一键收藏
  • 打赏
    打赏
  • 0
    评论
卡尔曼滤波流程图可以使用图解来更加清晰地展示。 卡尔曼滤波流程图由一系列的步骤组成,如下所示: 1. 初始化:设置初始状态和初始误差协方差矩阵。 2. 预测步骤:根据系统模型和控制输入,预测下一个状态和误差协方差矩阵。 3. 测量步骤:获取传感器测量值。 4. 更新步骤:使用测量值和预测值之间的差异,计算卡尔曼增益和更新后的状态和误差协方差矩阵。 5. 返回到第2步,重复预测、测量和更新步骤,以实现连续的状态估计和测量更新。 卡尔曼滤波的核心思想是通过使用当前的状态和测量值的信息,结合系统模型和控制输入,来估计出最优的状态估计值,并最小化估计误差。这使得卡尔曼滤波在估计动态系统状态时非常有效和准确。 请注意,卡尔曼滤波流程图可以根据具体的应用和系统模型进行调整和扩展,以满足具体需求。<span class="em">1</span><span class="em">2</span><span class="em">3</span> #### 引用[.reference_title] - *1* *2* *3* [一文图解卡尔曼滤波(Kalman Filter)](https://blog.csdn.net/qq_42722197/article/details/121586553)[target="_blank" data-report-click={"spm":"1018.2226.3001.9630","extra":{"utm_source":"vip_chatgpt_common_search_pc_result","utm_medium":"distribute.pc_search_result.none-task-cask-2~all~insert_cask~default-1-null.142^v93^chatsearchT3_2"}}] [.reference_item style="max-width: 100%"] [ .reference_list ]
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

西笑生

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付:¥1
获取中
扫码支付

您的余额不足,请更换扫码支付或充值

打赏作者

实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值