DeepSORT(目标跟踪算法) 卡尔曼滤波的完整流程

2024-06-14 06:04

本文主要是介绍DeepSORT(目标跟踪算法) 卡尔曼滤波的完整流程,希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

DeepSORT(目标跟踪算法) 卡尔曼滤波的完整流程

flyfish

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

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

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

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

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

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

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

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

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

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

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

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

卡尔曼滤波的完整流程

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. 循环迭代

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

各个字母的含义

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:预测的误差协方差矩阵。

这篇关于DeepSORT(目标跟踪算法) 卡尔曼滤波的完整流程的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



http://www.chinasem.cn/article/1059574

相关文章

工作流Activiti初体验—流程撤回【二】

已经玩工作流了,打算还是研究一下撤回的功能。但是流程图里面并不带撤回的组件,所以需要自己动态改造一下,还是延续上一个流程继续试验撤回功能。《工作流Activiti初体验【一】》 完整流程图 我们研究一下分发任务撤回到发起任务,其他环节的撤回类似 撤回的原理大概如下: 将分发任务后面的方向清空,把发起任务拼接到原来的判断网关,然后结束分发任务,这样流程就到发起任务了 此时的流程如上图,

ROS话题通信流程自定义数据格式

ROS话题通信流程自定义数据格式 需求流程实现步骤定义msg文件编辑配置文件编译 在 ROS 通信协议中,数据载体是一个较为重要组成部分,ROS 中通过 std_msgs 封装了一些原生的数据类型,比如:String、Int32、Int64、Char、Bool、Empty… 但是,这些数据一般只包含一个 data 字段,结构的单一意味着功能上的局限性,当传输一些复杂的数据,比如:

代码随想录算法训练营:12/60

非科班学习算法day12 | LeetCode150:逆波兰表达式 ,Leetcode239: 滑动窗口最大值  目录 介绍 一、基础概念补充: 1.c++字符串转为数字 1. std::stoi, std::stol, std::stoll, std::stoul, std::stoull(最常用) 2. std::stringstream 3. std::atoi, std

Android我的二维码扫描功能发展史(完整)

最近在研究下二维码扫描功能,跟据从网上查阅的资料到自己勉强已实现扫描功能来一一介绍我的二维码扫描功能实现的发展历程: 首页通过网络搜索发现做android二维码扫描功能看去都是基于google的ZXing项目开发。 2、搜索怎么使用ZXing实现自己的二维码扫描:从网上下载ZXing-2.2.zip以及core-2.2-source.jar文件,分别解压两个文件。然后把.jar解压出来的整个c

人工智能机器学习算法总结神经网络算法(前向及反向传播)

1.定义,意义和优缺点 定义: 神经网络算法是一种模仿人类大脑神经元之间连接方式的机器学习算法。通过多层神经元的组合和激活函数的非线性转换,神经网络能够学习数据的特征和模式,实现对复杂数据的建模和预测。(我们可以借助人类的神经元模型来更好的帮助我们理解该算法的本质,不过这里需要说明的是,虽然名字是神经网络,并且结构等等也是借鉴了神经网络,但其原型以及算法本质上还和生物层面的神经网络运行原理存在

大林 PID 算法

Dahlin PID算法是一种用于控制和调节系统的比例积分延迟算法。以下是一个简单的C语言实现示例: #include <stdio.h>// DALIN PID 结构体定义typedef struct {float SetPoint; // 设定点float Proportion; // 比例float Integral; // 积分float Derivative; // 微分flo

3月份目标——刷完乙级真题

https://www.patest.cn/contests/pat-b-practisePAT (Basic Level) Practice (中文) 标号标题通过提交通过率1001害死人不偿命的(3n+1)猜想 (15)31858792260.41002写出这个数 (20)21702664840.331003我要通过!(20)11071447060.251004成绩排名 (20)159644

LeetCode 算法:二叉树的中序遍历 c++

原题链接🔗:二叉树的中序遍历 难度:简单⭐️ 题目 给定一个二叉树的根节点 root ,返回 它的 中序 遍历 。 示例 1: 输入:root = [1,null,2,3] 输出:[1,3,2] 示例 2: 输入:root = [] 输出:[] 示例 3: 输入:root = [1] 输出:[1] 提示: 树中节点数目在范围 [0, 100] 内 -100 <= Node.

【Java算法】滑动窗口 下

​ ​    🔥个人主页: 中草药 🔥专栏:【算法工作坊】算法实战揭秘 🦌一.水果成篮 题目链接:904.水果成篮 ​ 算法原理 算法原理是使用“滑动窗口”(Sliding Window)策略,结合哈希表(Map)来高效地统计窗口内不同水果的种类数量。以下是详细分析: 初始化:创建一个空的哈希表 map 用来存储每种水果的数量,初始化左右指针 left

ROS2从入门到精通4-4:局部控制插件开发案例(以PID算法为例)

目录 0 专栏介绍1 控制插件编写模板1.1 构造控制插件类1.2 注册并导出插件1.3 编译与使用插件 2 基于PID的路径跟踪原理3 控制插件开发案例(PID算法)常见问题 0 专栏介绍 本专栏旨在通过对ROS2的系统学习,掌握ROS2底层基本分布式原理,并具有机器人建模和应用ROS2进行实际项目的开发和调试的工程能力。 🚀详情:《ROS2从入门到精通》 1 控制插