KalmanFilter卡尔曼滤波

news2024/11/28 2:34:33

预测(prediction)

状态变量

x k = A x k − 1 + B u k − 1 + w k − 1 x_k = Ax_{k-1} + Bu_{k-1} + w_{k-1} xk=Axk1+Buk1+wk1
其中n维向量 x k x_k xk为k时刻的系统状态变量,n维向量 x k − 1 x_{k−1} xk1是k-1时刻的系统状态变量。

A A A是状态转移矩阵或者过程增益矩阵,是 n × n n×n n×n阶方阵,它将k-1时刻状态和当前的k时刻状态联系起来。
在这里插入图片描述
B B B是可选的控制输入 u ∈ R l u\in R^l uRl的增益,在大多数实际情况下并没有控制增益,所以 B u k − 1 Bu_{k−1} Buk1这一项很愉快的变成零了。

w k − 1 w_{k-1} wk1是n维向量,代表过程激励噪声,它对应了中每个分量的噪声,是期望为0,协方差为Q的高斯白噪声, w ~ N ( 0 , Q ) w~N(0,Q) wN(0Q)
在这里插入图片描述

观测变量

z k = H x k + v k z_k = Hx_k + v_k zk=Hxk+vk
其中观测值 z k z_k zk是m阶向量,状态变量 x k x_k xk是n阶向量。对于目标检测, z k z_k zk就是检测结果转化后的值
在这里插入图片描述
观测矩阵 H H H m × n m×n m×n阶矩阵,代表状态变量 x k x_k xk对测量变量 z k z_k zk的增益
在这里插入图片描述
观测噪声 v k v_k vk是期望为0,协方差为R的高斯白噪声, v k ~ N ( 0 , R ) v_k~N(0,R) vkN(0,R)

误差协方差

P ^ k ˉ = A P ^ k − 1 A T + Q \hat P_{\bar k} = A \hat P_{k-1}A^T + Q P^kˉ=AP^k1AT+Q
在这里插入图片描述
在这里插入图片描述

P ^ k ˉ − 1 \hat P_{\bar k-1} P^kˉ1:表示k时刻的先验估计协方差,这个协方差矩阵只要确定了一开始的 P ^ 0 \hat P_0 P^0,后面都可以递推出来,而且初始协方差矩阵 P ^ 0 \hat P_0 P^0只要不是为0,它的取值对滤波效果影响很小,都能很快收敛。

P ^ k ˉ \hat P_{\bar k} P^kˉ:表示k时刻的后验估计协方差,是滤波结果之一。

Q Q Q:表示过程激励噪声的协方差,它是状态转移矩阵与实际过程之间的误差。这个矩阵是卡尔曼滤波中比较难确定的一个量,一般有两种思路:一是在某些稳定的过程可以假定它是固定的矩阵,通过寻找最优的Q值使滤波器获得更好的性能,这是调整滤波器参数的主要手段,Q一般是对角阵,且对角线上的值很小,便于快速收敛;二是在自适应卡尔曼滤波(AKF)中Q矩阵是随时间变化的。
在这里插入图片描述

更新(update/correction)

卡尔曼增益

K k = P ^ k ˉ H T H P ^ k ˉ H T + R K_k = \frac {\hat P_{\bar k}H^T}{H\hat P_{\bar k}H^T + R} Kk=HP^kˉHT+RP^kˉHT
在这里插入图片描述
在这里插入图片描述

H H H:表示量测矩阵,是 m × n m×n m×n阶矩阵,它把m维测量值转换到n维与状态变量相对应。
在这里插入图片描述
R R R:测量噪声的协方差矩阵,这个值过大过小都会使滤波效果变差,且取值越小,收敛越快。
在这里插入图片描述

由观测变量 z k z_k zk更新估计

x ^ k = x ^ k ˉ + K k ( z k − H x ^ k ˉ ) \hat x_k = \hat x_{\bar k} +K_k(z_k-H \hat x_{\bar k}) x^k=x^kˉ+Kk(zkHx^kˉ)
在这里插入图片描述
在这里插入图片描述

在这里插入图片描述
在这里插入图片描述

x ^ k ˉ \hat x_{\bar k} x^kˉ:表示 k k k时刻先验状态估计值,这是算法根据前次迭代结果(就是上一次循环的后验估计值)做出的不可靠估计。

x ^ k \hat x_k x^k:表示k时刻后验状态估计值,也就是要输出的该时刻最优估计值,这个值是卡尔曼滤波的结果。

更新误差协方差

P ^ k = ( I − K k H ) P ^ k ˉ \hat P_k = (I - K_kH)\hat P_{\bar k} P^k=(IKkH)P^kˉ
在这里插入图片描述
在这里插入图片描述

新的最佳估计是基于原最佳估计已知外部影响校正后得到的预测。

新的不确定性是基于原不确定性外部环境的不确定性得到的预测。

在这里插入图片描述

在这里插入图片描述

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.coloradmin.cn/o/940268.html

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈,一经查实,立即删除!

相关文章

Docker容器学习:搭建私有镜像仓库Harbor操作

目录 系统环境: 安装Docker-Compose 下载并安装Harber 启动Harbor 使用Harbor 上传下载镜像到harbor仓库 系统环境: Centos7.9Docker-ce:24 安装Docker-Compose curl -L "https://github.com/docker/compose/releases/download/v2.20.3/do…

【JUC系列-03】熟练掌握Atomic原子系列基本使用

JUC系列整体栏目 内容链接地址【一】深入理解JMM内存模型的底层实现原理https://zhenghuisheng.blog.csdn.net/article/details/132400429【二】深入理解CAS底层原理和基本使用https://blog.csdn.net/zhenghuishengq/article/details/132478786【三】熟练掌握Atomic原子系列基本…

Gossip协议

Gossip协议 一、Gossip协议1.1 工作原理1.2 Gossip优点1.3 Gossip传播方式1.3.1 Anti-Entropy(反熵)1.3.2 Rumor-Mongering(谣言传播)1.3.3 结合 1.4 Gossip协议的通信方式1.4.1 Push1.4.2 Pull1.4.3 Push&Pull 二、手撸简易版…

媒介盒子:医疗软文怎么写才能实现营销效果?

随着互联网的快速发展,医疗行业也逐渐意识到了网络营销的重要性。而作为网络营销的一种形式,医疗软文在传播医疗知识、宣传医疗品牌方面具有独特的优势。本文将从选题、内容、形式等多个方面进行探讨,如何写一篇有效的医疗营销软文? 1、选题非常关键 首…

Python“牵手”天猫商品列表数据,关键词搜索天猫API接口数据,天猫API接口申请指南

天猫平台API接口是为开发电商类应用程序而设计的一套完整的、跨浏览器、跨平台的接口规范,天猫API接口是指通过编程的方式,让开发者能够通过HTTP协议直接访问天猫平台的数据,包括商品信息、店铺信息、物流信息等,从而实现天猫平台…

MySql014——分组的GROUP BY子句排序ORDER BYSELECT子句顺序

前提:使用《MySql006——检索数据:基础select语句(使用products表、查询单列、多列、所有列、DISTINCT去除重复行、LIMIT限制返回结果的行数、了解完全限定)》中创建的products表 一、GROUP BY子句基础用法 SELECT vend_id, COU…

【Debug】解决RecursionError: maximum recursion depth exceeded in comparison报错

🚀Debug专栏 目录 🚀Debug专栏 ❓❓问题: 🔧🔧分析: 🎯🎯解决方案: ❓❓问题: 循环中报错RecursionError: maximum recursion depth exceeded in compari…

IC设计各岗位收入水平对比,看看哪个更适合你?

根据人才招聘平台对于2023年已有数据统计,芯片工程师岗位均薪为26012元,位列全行业第一。 这里需要说明一下,这里的“芯片工程师”涵盖了设计、制造、封测等多环节岗位。并非只有芯片设计岗。 从行业招聘薪酬同比增速来看,电子技…

Bigemap在路桥行业是怎么应用的?

选择Bigemap的原因: 奥维下架了,后来了解到的bigemap,于是测试了这款软件 使用场景: 下载影像、矢量路网做前期策划,下载完数据后导出cad ,做一些标注,最终出图下载等高线,作为前期选址依据 …

Linux环境下SVN服务器的搭建与公网访问:使用cpolar端口映射的实现方法

文章目录 前言1. Ubuntu安装SVN服务2. 修改配置文件2.1 修改svnserve.conf文件2.2 修改passwd文件2.3 修改authz文件 3. 启动svn服务4. 内网穿透4.1 安装cpolar内网穿透4.2 创建隧道映射本地端口 5. 测试公网访问6. 配置固定公网TCP端口地址6.1 保留一个固定的公网TCP端口地址6…

在日本做程序员能攒到钱吗?

如果你就是无欲无求,和人合租,自己做饭,不买高级食材,没有业余爱好,那我可以肯定告诉你一定能攒下钱,问题你是吗?能不能攒下钱丰俭由人,拿的少也有人能攒下钱,拿的多的也…

安防视频监控平台EasyNVR视频监控汇聚平台页面无法上传授权文件的问题解决方案

TSINGSEE青犀视频安防监控平台EasyNVR可支持设备通过RTSP/Onvif协议接入,并能对接入的视频流进行处理与多端分发,包括RTSP、RTMP、HTTP-FLV、WS-FLV、HLS、WebRTC等多种格式。在智慧安防等视频监控场景中,EasyNVR可提供视频实时监控直播、云端…

Viobot硬件控制

一.使用上位机控制 TOF版本设备点击TOF ON即可开启TOF,开启后按键会变成TOF OFF,点击TOF OFF即可关闭TOF 补光灯版本设备点击LED ON即可开启LED ,开启后按键会变成LED OFF,点击LED OFF即可关闭LED 设置页面的viobot栏&#xff0c…

13.redis集群、主从复制、哨兵

1.redis主从复制 主从复制是指将一台redis服务器(主节点-master)的数据复制到其他的redis服务器(从节点-slave),默认每台redis服务器都是主节点,每个主节点可以有多个或没有从节点,但一个从节点…

代码随想录算法训练营第五十二天|LeetCode 84.柱状图中最大的矩形

LeetCode 84.柱状图中最大的矩形 代码如下(Java):暴力解法 class Solution {public int largestRectangleArea(int[] heights) {int length heights.length;int[] minLeftIndex new int[length];int[] minRightIndex new int[length];min…

Java实现根据关键词搜索1688商品新品数据方法,1688API节课申请指南

要通过1688的API获取商品新品数据,您可以使用1688开放平台提供的接口来实现。以下是一种使用Java编程语言实现的示例,展示如何通过1688开放平台API获取商品新品数据: 首先,确保您已注册成为1688开放平台的开发者,并创…

【猿灰灰赠书活动 - 03期】- 【RHCSA/RHCE 红帽Linux认证学习指南(第7版) EX200 EX300】

说明:博文为大家争取福利,与清华大学出版社合作进行送书活动 图书:《RHCSA/RHCE 红帽Linux认证学习指南(第7版) EX200 & EX300》 一、好书推荐 图书介绍 《RHCSA/RHCE 红帽Linux认证学习指南(第7版) EX200 & E…

亚马逊云科技CEO分享企业内决策流程、领导力原则与重组下的思考

亚马逊云科技首席执行官Adam Selipsky几乎从一开始就在那里:他于2005年加入,在效力亚马逊11年后于2016年离开,转而经营Tableau,并于2021年成为亚马逊云科技首席执行官。当时亚马逊云科技前首席执行官安迪贾西(Andy Jassy)接替杰夫…

从数据孤岛到企业xPA的演化

“数据孤岛”一直以来是企业在信息化进程中面临的比较头疼的问题,由于数据独立存在于不同部门之中,无法进行相互联动,致使数据库无法兼容,这无形中加大了跨部门合作的沟通成本。在此背景下,一种新兴的规划方法——扩展…

知虾shopee数据分析工具:shopee出单的商机利器

当今数字化时代,数据已经成为商业成功的关键要素之一。而Shopee作为东南亚最大的电商平台之一,其强大的数据分析工具正为商家提供了宝贵的市场洞察和决策支持。本文将深入探讨Shopee数据分析工具如何帮助商家抓住商机并取得成功。 洞察消费者需求&#x…