扩展卡尔曼滤波EKF_第1页
扩展卡尔曼滤波EKF_第2页
扩展卡尔曼滤波EKF_第3页
扩展卡尔曼滤波EKF_第4页
全文预览已结束

下载本文档

版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领

文档简介

1、第三章 扩展卡尔曼滤波EKF3.1 扩展Kalman滤波原理Kalman滤波能够在线性高斯模型的条件下,可以对目标的状态做出最优的估计,得到较好的跟踪效果。对非线性滤波问题常用的处理方法是利用线性化技巧将其转化为一个近似的线性滤波问题。因此,可以利用非线性函数的局部性特性,将非线性模型局部化,再利用Kalman滤波算法完成滤波跟踪。扩展Kalman滤波就是基于这样的思想,将系统的非线性函数做一阶Taylor展开,得到线性化的系统方程从而完成对目标的滤波估计等处理。非线性系统离散动态方程可以表示为 (3-1-1) (3-1-2)这里为了便于数学处理,假定没有控制量的输入,并假定过程噪声是均值为零

2、的高斯白噪声,且噪声分布矩阵是已知的。其中,观测噪声也是加性均值为零的高斯白噪声。假定过程噪声和观测噪声序列是彼此独立的,并且有初始状态估计和协方差矩阵。和线性系统的情况一样,我们可以得到扩展Kalman滤波算法如下 (3-1-3) (3-1-4) (3-1-5) (3-1-6) (3-1-7)这里需要重要说明的是,状态转移和量测矩阵 是由和的雅克比矩阵代替的。其雅克比矩阵的求法如下: 假如状态变量有n维,即,则对状态方程对各维求偏导, (3-1-8) (3-1-9)3.2 扩展卡尔曼在一维非线性系统中的应用所谓的非线性方程,就是因变量和自变量的关系不是线性的,这类方程很多,例如平方关系,对数

3、关系,指数关系,三角函数关系等等。这类方程可分为两类,一类是多项式方程,一种是非多项式方程。为了便于说明非线性卡尔曼滤波扩展Kalman滤波的原理,我们选用一下系统,系统状态为,它仅包含一维变量,即,系统状态方程为 (3-2-1)观测方程为 (3-2-2)其中,式(3-1-1)是包含分式,平方,三角函数在内的严重非线性的方程,为过程噪声,其均值为0,方差为,观测方程中,观测信号与状态的关系也是非线性的,也是均值为0,方差为的高斯白噪声。因此关于(3-1-1)和(3-2-2)是一个状态和观测都为非线性的一维系统。以此为通用的非线性方程的代表,接下来讲述如何用扩展Kalman滤波来处理噪声问题。第

4、一步:初始化状态,协方差矩阵。第二步:状态预测 (3-2-3)第三步:观测预测 (3-2-4)第四步:一阶线性化状态方程,求解状态转移矩阵 (3-2-5)第五步:一阶线性化观测方程,求解观测矩阵 (3-2-6)第六步:求协方差矩阵预测 (3-2-7)这里需要说明的是,当噪声驱动矩阵不存在的时候,或系统状态方程中,在前没有任何驱动矩阵,这时候,必然和状态的维数一样的方阵,可将式(3-2-7)直接写为。第七步:求Kalman增益 (3-2-8)第八步:求状态更新 (3-2-9)第九步:协方差更新 (3-2-10)以上九步为扩展卡尔曼年滤波的一个计算周期,如此循环下去就是各个时刻EKF对非线性系统的处理过程。

温馨提示

  • 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
  • 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
  • 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
  • 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
  • 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
  • 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
  • 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。

最新文档

评论

0/150

提交评论