几种非线性滤波算法的研究-内附程序

  • 格式:doc
  • 大小:343.50 KB
  • 文档页数:9

下载文档原格式

  / 9
  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。

2017 年秋季学期研究生课程考核

(读书报告、研究报告)

考核科目:雷达系统导论

学生所在(系):电子与信息工程学院

学生所在学科:电子与同学工程

学生姓名:

学号:

学生类别:

考核结果阅卷人

第 1 页(共页)

几种非线性滤波算法的介绍与性能分析

作者姓名:学号:

专业院系:电信学院电子工程系

电子邮件:

摘要—非线性滤波算法在雷达目标跟踪中有着重要的应用,对雷达的跟踪性能有着至关重要的影响。好的滤波算法有利于目标航迹的建立及保持,能够得到较精确的目标位置,为发现目标后的后续工作提供可靠的数据依据。本文重点介绍了雷达数据处理中的几种非线性滤波算法:扩展卡尔曼滤波(EKF)、不敏卡尔曼滤波(UKF)、粒子滤波(PF),并且给出了一个利用这三种算法进行数据处理的一个实例,通过这个实例对比分析了这三种算法的性能以及优劣。

关键字—非线性滤波算法;扩展卡尔曼滤波;不敏卡尔曼滤波;粒子滤波;

I.概述(一级表题格式)

在雷达对目标进行跟踪前要先对目标进行检测。对于满足检测条件的目标就需要进行跟踪,在跟踪的过程中可以利用新获得的数据完成对目标的进一步检测比如去除虚假目标等,同时利用跟踪获得数据可以进一步完成对目标动态特性的检测和识别。因此对目标进行准确的跟踪是雷达性能的一个重要指标。在检测到满足条件的目标后,根据目标运动状态建立目标运动模型,然后对目标跟踪算法进行设计,这是雷达目标跟踪中的核心部分。

目前主要的跟踪算法包括线性自回归滤波,两点外推滤波、维纳滤波、-

αβ滤波、加权最小二乘滤波、维纳滤波和卡尔曼滤波[1]。对于线性系统而言最优滤波的方法就是卡尔曼滤波,卡尔曼滤波是线性高斯模型下的最优状态估计算法。但是实际问题中目标的运动模型往往不是线性的,因此卡尔曼滤波具有很大的局限性。目前主要用的非线性滤波算法可以分为高斯滤波和粒子滤波[2]。不敏卡尔曼滤波和扩展卡尔曼滤波就是高斯滤波中的典型代表,也是应用相对较为广泛的。粒子滤波的应用范围比高斯滤波的适用范围要广,对于系统状态非线性,观测模型非高斯等问题都有很好的适用性。本文具体分析阐述了扩展卡尔曼滤波算法,不敏卡尔曼滤波算法,粒子滤波算法,并且通过一个实例利用仿真的方法分析了这三种算法在滤波性能上的优劣,最后对这三种算法做了一定的总结。

我本科毕业设计题目为《基于历史数据的路径生成算法研究》,由于我是跨专业保研到电信学院,该课题所研究内容不属于雷达系统研究范围,是一种城市路网最快路径生成算法。

II.几种非线性滤波算法

A.扩展卡尔曼滤波

扩展卡尔曼滤波是将非线性系统转换为近似的线性系统的一种方法,其核心思想是围绕滤波值将非线性函数展开成泰勒级数并略去二阶及以上的项,得到一个近似的线性化模型,然后应用卡尔曼滤波完成状态估计。

扩展卡尔曼滤波状态空间模型:

k

k

k

w

x

f+

=

+

)

(

x

1

状态方程

k

k

k

v

x

h+

=)

(

z观测方程

其中(.)

f和(.)

h为非线性函数

在扩展卡尔曼滤波中,状态的预测以及观测值的预测由非线性函数计算得出,线性卡尔曼滤波中的状态转移矩阵A阵和观测矩阵H阵由f和h函数的雅克比矩阵代替。

(.)

f和(.)

h Taylor展开,只保留一次项有:

)

ˆ

(

)

ˆ(

)

(

k

k

k

k

k

x

x

A

x

f

x

f-

+

)

ˆ

(

)

ˆ(

)

(

k

k

k

k

k

x

x

H

x

h

x

h-

+

其中:

k

k

x

x

k

k dx

df

A

ˆ=

=为f对

1-

k

x求导的雅克比矩阵

k

k

x

x

k

k dx

dh

H

ˆ=

=为h对

1-

k

x求导的雅克比矩阵

)

ˆ(

ˆ

1-k

k

x

f

x=,于是可以得出:

k

k

k

k

k

k

k

w

x

A

x

f

x

A

x+

-

+

+

)

ˆ

)

ˆ(

(

1

k

k

k

k

k

k

k

v

x

H

x

h

x

H

z+

-

+

+

)

ˆ

)

ˆ(

(

1

通过以上变换,将非线性问题线性化。接下来EKF 滤波过程同线性卡尔曼滤波相同,公式如下:

))

|

(ˆ(

)

|1

(

Xˆk

k

X

f

k

k=

+

)

(

)

(

)

|

(

)

(

)

|1

(P k

Q

k

k

k

P

k

k

k+

Φ'

Φ

=

+

)1

(

)1

(

)

|1

(

)1

(

)1

(S+

+

+

'

+

+

=

+k

R

k

H

k

k

P

k

H

k

)1

(

)1

(

)

|1

(

)1

(

K1+

+

'

+

=

+-k

S

k

H

k

k

P

k