扩展卡尔曼滤波

扩展卡尔曼滤波(EKF)是一种处理非线性系统的滤波方法,通过线性化非线性函数来近似求解。本文详细介绍了EKF的系统描述、线性化过程以及滤波和观测更新的步骤,适用于高精度的动态系统预测和跟踪。在高非线性情况下,EKF可能性能下降,此时可考虑采用迭代EKF以提高滤波效果。
摘要由CSDN通过智能技术生成

EKF

目录

前言

一、系统描述

二、扩展卡尔曼滤波

1.线性化

2.算法流程

1)初始化

2)滤波、观测更新 

3)预测更新 


前言

针对非线性系统,卡尔曼滤波算法的效果比较差甚至不可使用。于是就有了扩展卡尔曼滤波,针对非线性系统也可以做到比较好的预测。


一、系统描述

对一个离散的非线性随机系统的状态空间的描述如下:

x_{k+1}=f_{k}\left ( x_{k} \right )+w_{k}                                k=0,1,2,.......

z_{k}=h_{k}\left ( x_{k} \right )+v_{k}                                     k=0,1,2,......

其中x_{k}是一个n维的列向量,z_{k}m维的列向量,分别是在时间t的状态和观测值。函数f_{k}n维的列向量映射到n维的列向量空间中去。函数h_{k}n维的状态列向量映射到观测空间的m维的列向量。w_{k}v_{k}分别是在状态空间和观测空间的相互独立的白噪声。白噪声的概率密度函数符合高斯分布,均值是0,方差分别是\Sigma _{k}^{w}\Sigma _{k}^{v}

二、扩展卡尔曼滤波

1.线性化

首先是对非线性函数的线性化,利用一阶泰勒函数的展开把非线性函数进行线性化。

h_{k}\left ( x_{k} \right )\approx h_{k}\left ( x\widehat{}_{k|k-1} \right )+H_{k}\left ( x_{k}-x\widehat{}_{k|k-1} \right )

f_{k}\left ( x_{k} \right )\approx f_{k}\left ( x\widehat{}_{k|k} \right )+F_{k}\left ( x_{k}-x\widehat{}_{k|k} \right )

其中H_{k}=\frac{\partial h\left ( x_{k} \right )}{\partial x_{k}}|x_{k=x\widehat{}_{k|k-1}}h_{k}\left ( \cdot \right )x_{k}=x\widehat{}_{k|k-1}处的估计值。 

其中F_{k}=\frac{\partial f\left ( x_{k} \right )}{\partial x_{k}}|x_{k=x\widehat{}_{k|k}}f_{k}\left ( \cdot \right )x_{k}=x\widehat{}_{k|k}处的估计值。

2.算法流程

1)初始化

x\widehat{}_{0|-1}=x\widehat{}_{0}

P_{0|-1}=P_{0}

2)滤波、观测更新 

K_{k}=P_{k|k-1}H_{k}^{T}\left ( H_{k}P_{k|k-1}H_{k}^{T}+\Sigma _{k}^{v} \right )^{-1}

x\widehat{}_{k|k}=x\widehat{}_{k|k-1}+K_{k}\left ( z_{k}-h_{k}\left (x\widehat{}_{k|k-1}\right ) \right )

P_{k|k}=P_{k|k-1}-K_{k}H_{k}P_{k|k-1}

3)预测更新 

x\hat{}_{k+1|k}=f_{k}\left (x\hat{}_{k|k}\right )

P_{k+1|k}=F_{k}P_{k|k}F_{k}^{T}+\Sigma _{k}^{w}

观测方程如果是强非线性的则滤波器的性能会恶化,此时迭代扩展卡尔曼滤波的效果会更好。


  • 0
    点赞
  • 3
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

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

余额充值