卡尔曼滤波算法的深入理解

2023-05-16

        以前对于这个算法只会调用,并没有深入的理解。现在碰到一个项目,没办法只能硬着头皮去理解了。网上看了很多资料,感觉写的也一般般,现在我讲自己的理解。

        卡尔曼滤波算法的核心是动态调整权值。用过互补滤波的应该知道它的权值是静态的,而卡尔曼是动态的。刚刚接触卡尔曼也不要紧张,我来一步步剖析这个东西。

       以上五条公式是卡尔曼的核心。它的本质就是通过预测结合测量来估计当前系统的状态。举个例子,假如我们要估计一架飞行器的姿态,可以通过IMU来实时测量,但是测量值有一定的风险是不准确的,所以并不能完全依赖传感器。任何一个满足物理规律的系统应当是连续的,所以我们还可以通过上一状态来预测当前状态。Kalman Filter正是结合这两条进行状态估计,到底是相信哪一个多一点,还要根据Kt来决定,我们定义Kt为卡尔曼增益,它是根据 测量和预测的协方差来计算的。

      先解释下每个公式所要表达的含义以及变量的含义:

line 2: 首先通过上一状态最优值和将要施加的控制量来预测当前状态,由假设一可以得到:

这里写图片描述

因为我们只是求均值,而高斯噪声均值为0,所以可省去最后一项。

At指当前时刻的状态转移矩阵(就是指从上一状态转变为下一状态的关系矩阵);

Bt指当前时刻的控制矩阵(就是指影响控制量的控制矩阵);

ut指当前时刻的控制量;

ut-1指上一时刻最优估计值;

指当前时刻估计值;

line 3: 除了预测均值之外,我们还需要预测值的协方差来计算Kalman增益。

这里写图片描述

指上一时刻的预测值协方差矩阵;

Rt指当前时刻测量值噪声矩阵;

line 4:根据预测值的协方差这里写图片描述,测量值和状态的比例系数,测量值的协方差来计算Kalman增益。

这里写图片描述

Kt指卡尔曼增益(指的就是权值);

Ct指的是当前时刻的测量方程;

Qt表示观察量的协方差矩阵;

line 5:这一行可以说是Kalman Filter 的精华了,现在我们有了对状态的预测值和协方差,同时也收集到了对状态的测量值。这时就可以通过kalman增益来计算状态估计值了。

这里写图片描述

增益越大,表明我们越相信测量值。

line 6: 根据 line3 ,预测当前状态需要用到上一状态的协方差,所以我们还需要计算当前状态的协方差用于下一次迭代。它同样要根据Kalman增益来计算:

这里写图片描述

相信到这里,大家应该对kalman Filter的原理有了一个大致的了解,算法中,从初始状态开始,不断计算当前状态的均值和方差来迭代,直至系统结束。

        上面解释了各个公式以及各个变量的含义。其实扩展卡尔曼的主要作用还是在不断的迭代中求出最接近真实值的那个值。卡尔曼滤波的作用有以下两种。第一,如果卡尔曼用作单种数据滤波(或者多种数据分开),那么将数据作为测量量传入模型中,卡尔曼模型会通过上一次的值估计出下一次的值,然后将此次的估计值和测量值分别取一定的权值(模型自己所计算的权值(卡尔曼增益)),求出这次的最优值。第二,是多数据的融合。可将一种数据作为测量量,另一种数据作为估计值进行融合。一般来说观测的数据要准确一些,估计的数据要不准确一些。

附上视频:https://www.bilibili.com/video/av4356232?from=search&seid=12641787542715420206

卡尔曼数据融合matalab代码:https://download.csdn.net/my/uploads/1/2

本文内容由网友自发贡献,版权归原作者所有,本站不承担相应法律责任。如您发现有涉嫌抄袭侵权的内容,请联系:hwhale#tublm.com(使用前将#替换为@)

卡尔曼滤波算法的深入理解 的相关文章

  • 嵌入式面试刷题

    1 表示一年有多少秒 define SECONDS PER YEAR 606024 365 UL 2 写一个标准宏 MIN define MIN A B A lt 61 B A B 3 指针数组 int a 10 优先级高所以是a 10 数
  • ubuntu18配置ftp

    安装ftp 修改文件 sudo gedit etc vsftpd conf span class token comment Example config file etc vsftpd conf span span class token
  • win11切换win10资源管理器

    HKEY LOCAL MACHINE SOFTWARE Microsoft Windows CurrentVersion Shell Extensions 右键 Blocked 选择 新建 字符串值 名称为 e2bf9676 5f8f 43
  • 联想小新pro13笔记本外接显示屏没信号

    step1 关机 step2 拔下电源 step3 安住 fn 43 s 43 v键 xff0c 开机 xff08 开不了机 xff0c 我重复了几次 xff09 step4 插电源开机 xff0c 扩展屏幕正常
  • C51内存类型

    bdata bdata内存类型只能用于声明变量 您不能声明bdata函数 该存储器使用8位地址直接访问 xff0c 是8051的片内位可寻址RAM 用bdata类型声明的变量是位可寻址的 xff0c 可以用位指令读写 code 代码存储器类
  • MQTT问题

    是否存在c gt ping outstanding 61 1 的后一秒就触发TimerIsExpired amp c gt last received span class token keyword int span span class
  • Failed to start apt-news.service Failed to start esm-cache.service

    luozw 64 luozw vpc etc apt apt conf d span class token function sudo span span class token function apt get span update
  • stm中断优先级理解+抢占优先级和相应优先级

    一 抢占优先级比子优先级的优先权更高 xff0c 这意味抢占优先级更高的中断会先执行 xff0c 而不管子优先级的优先权 xff0c 数值越低优先级越高 二 同理 xff0c 如果抢占优先级相同 xff0c 那么就会比较子优先级 xff0c
  • Realsense D435基于ROS跑通ORBSLAM2

    Realsense D435基于ROS跑通ORBSLAM2 系统ubuntu16 04 ROS Kinetic 相机RealsenseD435 SLAM系统 xff1a ORBSLAM2 一 安装Realsense的SDK 官方链接 htt
  • Qt学习:综合案例应用-上(翻金币小游戏)

    本案例是对Qt的基本控件 xff0c 事件处理 xff0c 资源文件的使用等知识的综合应用 以及一些开发思想和逻辑控制 首先了解下案例的文件构成 头文件 xff1a mainwindow h chooselevelscene h plays
  • 在TX2上运行realsenseD435摄像头

    在TX2上运行realsenseD435 先给出相关的链接在TX2上安装realsense SDK库在TX2上安装realsense SDK库 先给出相关的链接 github 上的一些链接 realsense SDK库 xff1a http
  • docker build 后面为什么要跟个 .

    我们在构建镜像文件时无非是使用 xff1a docker build t test ubuntu v1 或者 docker build f docker test Dockerfile 来进行构建镜像 xff0c 用第一个命令时任务 指代的
  • 微机原理中地址总线、数据总线与内存容量之间的关系

    今天在复习微机原理的时候 xff0c 看到一个概念 xff1a 存储总量 61 存储单元个数 存储字长 xff0c 然后存储单元个数 61 2 地址总线位数 xff0c 存储字长和数据总线位数有关 xff0c 如果是这样 xff0c 那么
  • HDLC——高级数据链路控制(HDLC,High-level Data Link Control)

    一 HDLC概述 1 1 HDLC的发展历史 高级数据链路控制 xff08 High Level Data Link Control或简称HDLC xff09 xff0c 是一个在同步网上传输数据 面向比特的数据链路层协议 xff0c 它是
  • 差分技术:LVDS(低压差分信号)、MLVDS(多点低压差分信号)的区别与应用场景

    差分传输是一种信号传输的技术 xff0c 区别于传统的一根信号线一根地线的做法 xff0c 差分传输在这两根线上都传输信号 xff0c 这两个信号的振幅相同 xff0c 相位相反 在这两根线上的传输的信号就是差分信号 信号接收端比较这两个电
  • 小白能理解的奈奎斯特采样及延伸出的理论

    一 取样定理 其实奈奎斯特采样有两种方式 xff0c 一种是矩形脉冲采样 xff0c 一种是冲激采样 xff0c 采样方式如下图 我们在不计算数学公式的情况下来讲解 xff0c 只是让大家明白是这么回事 xff0c 具体为什么是这样 xff
  • 单边谱和双边谱

    实际中 xff0c 只会有单边谱 xff0c 并不会有负频率的信号 在引入欧拉公式后 xff0c 出现了双边谱 单边谱转换为双边谱 xff0c 幅度会降低一半 xff0c 其他不变
  • 小白也能搞通UDP通信(88E1111 RGMII 接口)

    一 网络协议 下表描述了整个从上到下的网络协议层 xff1a 这些网络协议在FPGA实际开发的过程中用到的就是传输层 网络层 数据链路层和物理层 xff0c 在我们的举例中用到UDP IP ARP协议 xff0c 物理层就用88E1111
  • 485通讯和modbus通讯协议

    485通信 xff1a 采用差分信号 xff1a A比B电压高是1 xff0c A比B电压低是0 xff0c 电压高低值在0 2V 6V之间 硬件连接上 xff1a 所有A接到一起 xff0c 所有B接到一起AB之间要加匹配电阻100欧到1
  • MODBUS RTU

    Modbus xff1a 是一种单主 从通信协议 MODBUS网络上只有一个主站 xff0c 主站在MODBUS网络没有地址 xff0c 从站的地址范围为0 247 xff0c 其中0为广播地址 xff0c 从站的实际地址为1 247 MO

随机推荐

  • TM4C123-HWREG()及外设寄存器地址说明

    参考文件 xff1a ti TivaWare C Series 2 1 4 178 inc hw types hti TivaWare C Series 2 1 4 178 inc hw memmap htm4c123gh6pz datas
  • I2C协议

    物理层 xff1a 1 一个I2C总线中可连接多个I2C通信设备 xff0c 支持多个主机及多个从机 2 两线制 xff1a 一条双向串行通信的数据线 xff08 SDA xff09 xff0c 一条串行时钟线 xff08 SCL xff0
  • SPI协议

    物理层 xff1a 1 四线制或三线制 xff1a 四线制时3条总线分别为SCK MOSI MISO xff0c 片选线为NSS xff08 CS xff09 三线制与其不同的是MOSI和MISO合并为一条线 xff0c 端口为双向端口 x
  • 数据处理:修改rosbag topic的frame_id

    1 安装srv tools工具 git clone https github com srv srv tools git 将其作为功能包放到ros工作空间下编译即可 2 使用示例 rosrun bag tools change frame
  • STM32 Cubemax(十二) ——利用状态机实现按键的长短按和双击

    STM32 Cubemax 十二 利用状态机实现按键的长短按和双击 文章目录 STM32 Cubemax 十二 利用状态机实现按键的长短按和双击前言一 状态图二 Cubemax配置1 IO口配置2 定时器配置 三 代码1 编写有关按键的结构
  • python中的闭包函数

    闭包函数初探 通常我们定义函数都是这样定义的 def foo pass 其实在函数式编程中 xff0c 函数里面还可以嵌套函数 xff0c 如下面这样 def foo print 34 hello world in foo 34 def b
  • 树莓派+Ubuntu16.04 MATE+ROS(kinetic)的安装

    树莓派 43 Ubuntu16 04 MATE 43 ROS xff08 kinetic xff09 的安装 零 准备工作 1 树莓派板子 2 电源 xff1a 3A以上 xff0c 供电要充足 xff0c 否则无法启动或影响使用 3 显示
  • canal

    目录 搭建canal测试canal 监控MySQL的binlog的工具 搭建canal 1 开启mysql binlog cp usr share mysql my medium cnf etc my cnf 修改my cnf vim et
  • docker 镜像注册【图文教程】

    docker镜像官网 xff1a https hub docker com 进行登录注册 账号 只能使用4到30个字母和数字 那就名字拼音吧 邮箱 正确的邮箱即可 后续会验证 xff0c 务必真实 密码 密码至少为9个字符 第一个勾选框 偶
  • 使用gazebo实现turtlebot入门级开发

    实验室准备新进两台Turtlebot2 xff0c 为了更快上手 xff0c 便提前开始熟悉一下Turtlebot xff0c 通果查阅相关资料 xff0c 我写了一个demo程序 xff0c 并在gazebo模拟环境下进行了测试 该dem
  • 手写智能指针(类)

    基础知识 xff1a 智能指针的设计与实现 1 智能指针类将一个计数器与类指向的对象相关联 xff0c 引用计数跟踪该类有多少个对象共享同一指针 2 每次创建类的新对象时 xff0c 初始化指针并将引用计数置为1 xff1b 3 当对象作为
  • 求职历程

    前言 xff1a 个人求职有大部分的时间花在了城市与企业性质的纠结上 首先是个人打算 xff0c 因研究生课题偏嵌入式相关 xff0c 所以给自己定下的目标是找嵌入式相关的工作 xff0c 而弱化了互联网求职的刷题关和计算机基础关 而在求职
  • 剑指 Offer 57 - II. 和为s的连续正数序列

    剑指 Offer 57 II 和为s的连续正数序列 难度 简单 输入一个正整数 target xff0c 输出所有和为 target 的连续正整数序列 xff08 至少含有两个数 xff09 序列内的数字由小到大排列 xff0c 不同序列按
  • 使用层次遍历重建二叉树并遍历

    58同城的一道题 xff0c 蛮有意思的 利用层次遍历后的数组 xff0c 进行二叉树的重建 数值 1代表NULL 解题思路 xff1a 那么我们依然使用队列 xff0c 进行层次遍历 xff0c 进行重建 xff0c 这边有的问题是当示例
  • C++秋招知识点总结

    先介绍下本人的情况 xff0c 本人研究生双非院校 xff0c 本科也是普通二本 xff0c 也无实习无项目非科班出生 刷了将近半年的题和自学了相关知识点 xff0c 最终拿到华为 海康 大华 同花顺 深信服 CVTE 紫光展锐offer
  • MPU6050软件姿态解算

    我是勤劳的搬运工 转自 http blog sina com cn s blog c5a00db10102wd7d html http mini eastday com mobile 180306210610472 html 使用MPU60
  • AD布线的线宽

    线宽关系图 xff1a 覆铜厚度与电流的关系 xff1a
  • AD16 Unknown Pin 和Failed to add class member 问题

    转载 xff1a http www 51hei com bbs dpj 39227 1 html 相信很多人在使用DXP 和AD 的时候 xff0c 都遇到过Unknown Pin 和Failed to add classmember 的问
  • 加速度计求二次积分

    我是勤劳的搬运工 xff1a https blog csdn net u011006622 article details 56286833 摘要 此文档描述并使用MMA7260QT三轴加速计和低功耗的9S08QG8八位单片机实现求解位置的
  • 卡尔曼滤波算法的深入理解

    以前对于这个算法只会调用 xff0c 并没有深入的理解 现在碰到一个项目 xff0c 没办法只能硬着头皮去理解了 网上看了很多资料 xff0c 感觉写的也一般般 xff0c 现在我讲自己的理解 卡尔曼滤波算法的核心是动态调整权值 用过互补滤