卡尔曼滤波 - 状态空间模型中的状态方程

这篇具有很好参考价值的文章主要介绍了卡尔曼滤波 - 状态空间模型中的状态方程。希望对大家有所帮助。如果存在错误或未考虑完全的地方,请大家不吝赐教,您也可以点击"举报违法"按钮提交疑问。

卡尔曼滤波 - 状态空间模型中的状态方程

flyfish

状态方程和观测方程统称为状态空间模型

位移

陀螺仪 加速度计 卡尔曼滤波 状态空间方程,多目标跟踪,算法,概率论,线性代数,卡尔曼滤波
 位移 = Δ x = x f − x 0 \text { 位移}=\Delta x=x_f-x_0  位移=Δx=xfx0
x 0 x_0 x0 是起始位置
x f x_f xf 是终止位置
在坐标轴里,右边是正,左边是负

面积等于物体的位移

陀螺仪 加速度计 卡尔曼滤波 状态空间方程,多目标跟踪,算法,概率论,线性代数,卡尔曼滤波

绿色矩形的高度为 v 0 v_0 v0
宽度为 t t t
所以面积等于 v 0 v_0 v0 t t t

黄色三角形的底是 t t t
高度为 v − v 0 v-v_0 vv0
黄色三角形的面积为
1 2 t ( v − v 0 ) \large \frac {1}{2} t (v - v_0) 21t(vv0)

两者求和时,我们得到位移公式
Δ x = v 0 t + 1 2 t ( v − v 0 ) \large \Delta x = v_0 t + \frac {1}{2} t (v - v_0) Δx=v0t+21t(vv0)

加速度
( 1 ) a = Δ v Δ t (1)\Large a = \frac {\Delta v}{\Delta t} 1a=ΔtΔv
速度差,也就是三角形的高
( 2 ) Δ v = v − v 0 (2)\Large \Delta v = v - v\normalsize{_0} 2Δv=vv0
(2)式代入(1)式
a = v − v 0 Δ t \Large a = \frac {v - v_0}{\Delta t} a=Δtvv0

v = v 0 + a Δ t \Large v = v_0 + \Large a \Delta t v=v0+aΔt
也就是我们常见的
v = v 0 + a t \LARGE v = v_0 + at v=v0+at

  1. 简化位移公式
    Δ x = v 0 t + 1 2 a t 2 \large \Delta x = v_0 t + \frac {1}{2} at^2 Δx=v0t+21at2

黄色和绿色左右分布也一样

陀螺仪 加速度计 卡尔曼滤波 状态空间方程,多目标跟踪,算法,概率论,线性代数,卡尔曼滤波

定义系统的状态变量

x ( t ) x(t) x(t)是位置
x ˙ ( t ) \dot{x}(t) x˙(t)是速度
a = w ( t ) a=w(t) a=w(t)是加速度

x ( t ) = [ x ( t ) x ˙ ( t ) ] = [  位置   速度  ] \boldsymbol{x}(t)=\left[\begin{array}{l} x(t) \\ \dot{x}(t) \end{array}\right]=\left[\begin{array}{l} \text { 位置 } \\ \text { 速度 } \end{array}\right] x(t)=[x(t)x˙(t)]=[ 位置  速度 ]

x ( t ) = x 0 + v 0 t + 1 2 a t 2 {x}(t) = x_0 + v_0 t + \frac {1}{2} at^2 x(t)=x0+v0t+21at2 (位置)

x ˙ ( t ) = v ( t ) = v 0 + a t \dot{x}(t)=v(t) =v_0+a t x˙(t)=v(t)=v0+at (速度)

x ¨ ( t ) = v ˙ ( t ) = a ( t ) = a \ddot{x}(t)=\dot{v}(t)=a(t) =a x¨(t)=v˙(t)=a(t)=a (加速度)

t t t t − 1 t-1 t1 时刻之间的时间差为 Δ t \Delta t Δt
{ x t = x t − 1 + x ˙ t − 1 Δ t + 1 2 a Δ t 2 x ˙ t = x ˙ t − 1 + a Δ t \left\{\begin{array}{l} x_t=x_{t-1}+\dot{x}_{t-1} \Delta t+\frac{1}{2} a \Delta t^2 \\ \dot{x}_t=\dot{x}_{t-1}+a \Delta t \end{array}\right. {xt=xt1+x˙t1Δt+21aΔt2x˙t=x˙t1+aΔt

矩阵形式
[ x t x ˙ t ] = [ 1 Δ t 0 1 ] [ x t − 1 x ˙ t − 1 ] + [ Δ t 2 2 Δ t ] a \left[\begin{array}{l} x_t \\ \dot{x}_t \end{array}\right]=\left[\begin{array}{cc} 1 & \Delta t \\ 0 & 1 \end{array}\right]\left[\begin{array}{l} x_{t-1} \\ \dot{x}_{t-1} \end{array}\right]+\left[\begin{array}{c} \frac{\Delta t^2}{2} \\ \Delta t \end{array}\right] a [xtx˙t]=[10Δt1][xt1x˙t1]+[2Δt2Δt]a

A = [ 1 Δ t 0 1 ] , B = [ Δ t 2 2 Δ t ] A=\left[\begin{array}{cc} 1 & \Delta t \\ 0 & 1 \end{array}\right], B=\left[\begin{array}{c} \frac{\Delta t^2}{2} \\ \Delta t \end{array}\right] A=[10Δt1],B=[2Δt2Δt]
x t = A x t − 1 + B u t − 1 (状态方程 ) {x_t} = A{x_{t - 1}} + B{u_{t - 1}}(状态方程) xt=Axt1+But1(状态方程)
t t t 换成 k k k
x k = A x k − 1 + B u k − 1 (状态方程 ) {x_k} = A{x_{k - 1}} + B{u_{k - 1}}(状态方程) xk=Axk1+Buk1(状态方程)

x k = A x k − 1 + B u k − 1 + w k − 1 (加入随机项状态方程 ) {x_k} = A{x_{k - 1}} + B{u_{k - 1}} + {w_{k - 1}}(加入随机项状态方程) xk=Axk1+Buk1+wk1(加入随机项状态方程)

其他式子
先验估计
− 代表先验,ˆ代表估计
x ^ k − = A x ^ k − 1 + B u k − 1 \hat{x}_{k}^{-}=A \hat{x}_{k-1}+B u_{k-1} x^k=Ax^k1+Buk1文章来源地址https://www.toymoban.com/news/detail-784758.html

到了这里,关于卡尔曼滤波 - 状态空间模型中的状态方程的文章就介绍完了。如果您还想了解更多内容,请在右上角搜索TOY模板网以前的文章或继续浏览下面的相关文章,希望大家以后多多支持TOY模板网!

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处: 如若内容造成侵权/违法违规/事实不符,请点击违法举报进行投诉反馈,一经查实,立即删除!

领支付宝红包 赞助服务器费用

相关文章

  • 卡尔曼滤波模型及Matlab模型建立

    目录 一、卡尔曼滤波 1.概念解析: 2.卡尔曼滤波的最优估计模型 3.实例——小车匀加速直线运动 4.Matlab建模  二、扩展卡尔曼滤波(EKF:Extended KAlman Filter) 1.非线性系统的局部线性化 2.扩展卡尔曼滤波模型         卡尔曼滤波的目的:由于人的主观认识(数学模型的建立而产

    2023年04月10日
    浏览(29)
  • 【状态估计】变分贝叶斯近似的递归噪声自适应卡尔曼滤波(Matlab代码实现)

    💥💥💞💞 欢迎来到本博客 ❤️❤️💥💥 🏆博主优势: 🌞🌞🌞 博客内容尽量做到思维缜密,逻辑清晰,为了方便读者。 ⛳️ 座右铭: 行百里者,半于九十。 📋📋📋 本文目录如下: 🎁🎁🎁 目录 💥1 概述 📚2 运行结果 🎉3 参考文献 🌈4 Matlab代码及文献 文献来

    2024年02月09日
    浏览(32)
  • 了解卡尔曼滤波器4--非线性状态估算器(EKF,UKF,PF)

            一般来说,我们希望我们的生活是线性的,就像这条线,这可能表示成功、收入或者幸福。但实际上,生活并不是线性的,它充满了起伏,有时甚至更复杂。         如果您是工程师,您经常会需要处理非线性系统,为了帮助您,我们将讨论非线性状态估算器

    2023年04月20日
    浏览(42)
  • 【状态估计】粒子滤波器、Σ点滤波器和扩展/线性卡尔曼滤波器研究(Matlab代码实现)

    💥💥💞💞 欢迎来到本博客 ❤️❤️💥💥 🏆博主优势: 🌞🌞🌞 博客内容尽量做到思维缜密,逻辑清晰,为了方便读者。 ⛳️ 座右铭: 行百里者,半于九十。 📋📋📋 本文目录如下: 🎁🎁🎁 目录 💥1 概述 📚2 运行结果 2.1 扩展卡尔曼滤波 2.2 线性卡尔曼滤波 

    2024年02月09日
    浏览(29)
  • 【状态估计】基于线性卡尔曼滤波器和粒子滤波器无人机估计地形高度(Matlab代码实现)

      💥💥💞💞 欢迎来到本博客 ❤️❤️💥💥 🏆博主优势: 🌞🌞🌞 博客内容尽量做到思维缜密,逻辑清晰,为了方便读者。 ⛳️ 座右铭: 行百里者,半于九十。 📋📋📋 本文目录如下: 🎁🎁🎁 目录 💥1 概述 📚2 运行结果 🎉3 参考文献 🌈4 Matlab代码实现 本文模

    2024年02月16日
    浏览(33)
  • 扩展卡尔曼滤波在目标跟踪中的应用(1)

    前面几节相继介绍的不同维度的 卡尔曼滤波算法 , 在线性高斯模型的条件下,具有很好的跟踪效果 ,当时在实际生活中, 不可能 所有的运动都呈现一种线性关系,也存在着 非线性 关系,比如说:飞机的飞行状态,导弹的制导问题等等。 为了保证结果的准确性和精确度

    2024年02月10日
    浏览(37)
  • 无迹卡尔曼滤波在目标跟踪中的作用(一)

    在前一节中,我们介绍了 扩展卡尔曼滤波算法EKF 在目标跟踪中的应用,其 原理是 将非线性函数局部线性化,舍弃高阶泰勒项,只保留一次项 ,这就不可避免地会影响结果的准确性,除此以外,实际中 要计算雅各比矩阵 也 不是特别的容易 ,因此有必要研究其他的滤波算法

    2024年02月11日
    浏览(29)
  • 无迹卡尔曼滤波在目标跟踪中的作用(二)

    在上一节的内容中,我们介绍了 UKF 中最重要的内容— 无迹变换UT ,今天我们将具体介绍 UKF 是如何实现的。 好了,话不多说,开整!!! 我们知道,我们可以使用 状态方程 和 观测方 程来 对系统进行描述 ,那么一个 非线性系统 可以用以下的方程进行描述: { X ( k + 1 )

    2024年02月12日
    浏览(37)
  • 【鲁棒、状态估计】用于电力系统动态状态估计的鲁棒迭代扩展卡尔曼滤波器研究(Matlab代码实现)

    💥💥💞💞 欢迎来到本博客 ❤️❤️💥💥 🏆博主优势: 🌞🌞🌞 博客内容尽量做到思维缜密,逻辑清晰,为了方便读者。 ⛳️ 座右铭: 行百里者,半于九十。 📋📋📋 本文目录如下: 🎁🎁🎁 目录 💥1 概述 📚2 运行结果 2.1 IEEE14节点  ​ 2.2 IEEE30节点  ​  2.3 

    2024年02月07日
    浏览(29)
  • 基于卡尔曼滤波检测GNSS 数据中的智能手机位移附matlab代码

     ✅作者简介:热爱科研的Matlab仿真开发者,修心和技术同步精进,代码获取、论文复现及科研仿真合作可私信。 🍎个人主页:Matlab科研工作室 🍊个人信条:格物致知。 更多Matlab完整代码及仿真定制内容点击👇 智能优化算法       神经网络预测       雷达通信       无

    2024年04月25日
    浏览(27)

觉得文章有用就打赏一下文章作者

支付宝扫一扫打赏

博客赞助

微信扫一扫打赏

请作者喝杯咖啡吧~博客赞助

支付宝扫一扫领取红包,优惠每天领

二维码1

领取红包

二维码2

领红包