第 27 卷 第 1 期
Vol. 27 No. 1
控 制 与 决 策
Control and Decision
2012 年 1 月
Jan. 2012
基于均差滤波与高斯和的非线性非高斯系统滤波算法
文章编号: 1001-0920 (2012) 01-0129-06
李振华
1
, 宁 磊
1
, 徐胜男
2
(1. 山东大学 控制科学与工程学院,济南 250061;2. 济南广播电视大学 信息技术学院,济南 250001)
摘 要: 针对一类非线性非高斯系统的滤波问题, 在分析均差滤波算法和高斯和滤波算法的基础上, 提出一种基于
均差滤波的高斯和滤波算法, 适于处理非线性非高斯系统的滤波问题. 对于似然密度位于条件转移概率密度拖尾处
的情况, 与传统的粒子滤波算法相比, 所提算法能提高滤波的精度和实时性. 仿真实验验证了新算法的有效性.
关键词: 非线性非高斯滤波;贝叶斯统计;均差滤波;高斯和滤波
中图分类号: TP242 文献标识码: A
Nonlinear non-Gaussian system filtering based on Gaussian sum and
divided difference filter
LI Zhen-hua
1
, NING Lei
1
, XU Sheng-nan
2
(1. School of Control Science and Engineering,Shandong University,Ji’nan 250061,China;2. School of
Information Technologh,Ji’nan Radio and TV University,Ji’nan 250001,China. Correspondent:LI Zhen-hua,
E-mail:zhenhua.li@sdu.edu.cn)
Abstract: Based on analyzing divided difference filter(DDF) and Gaussian sum filter(GSF), a GSF-based DDF algorithm
is developed for nonlinear dynamic state space(DSS) models with non-Gaussian noise, which is suitable for the filtering
problem of nonlinear/non-Gaussian systems. When the likelihood function appeares at the tail of the transfer probability
density, the proposed algorithm can improve the precision of nonlinear/non-Gaussian filtering compared with the traditional
particle filter(PF). Experiments show that the proposed method works well in the filtering for DSS models with non-Gaussian
noise.
Key words: nonlinear non-Gaussian niltering;Bayesian estimation;divided difference filter;Gaussian sum filter
1 引引引 言言言
滤波问题广泛存在于自动控制、信号处理、金融
统计、导航、目标跟踪、信息融合、故障检测等领域,
在过去的半个世纪中得到了较为深入的研究. 其基本
方法是通过含有噪声的观测值来递推估计系统的未
知状态, 这种方法中的系统模型被称为动态状态空间
模型 (DSSM). 根据贝叶斯统计理论, 后验密度包含了
状态变量的全部统计信息, 因此, 滤波问题可转变为
在每次获得新的量测值后, 如何根据现有的全部量测
信息递推地估计状态变量的后验密度.
当 DSSM 是线性的且与之相关的噪声满足高斯
分布时, 卡尔曼滤波算法 (KF)
[1]
是贝叶斯递推估计问
题的最优解. 当 DSSM 变为非线性非高斯时, 最优滤
波便失去了吸引力, 因此需要寻找新的次优滤波算法
来替代 KF. 非线性非高斯滤波算法大致可分为两类:
一类是将后验密度视为高斯分布的局部逼近方法.
例如文献 [1] 中提到的利用一阶或二阶 Taylor 级数展
开的扩展卡尔曼滤波算法 (EKF); [2-3] 提到的应用无
味变换 (UT) 的无味卡尔曼滤波算法 (UKF); [4-6] 提
到 的 利 用 Stiring 内 插 公 式 的 均 差 滤 波 算 法 (DDF);
[7] 提 到 的 利 用 Gauss-Hermite 积 分 的 Gauss-Hermite
滤波算法 (GHQ-KF) 等. 由于这些算法都假定后验密
度为高斯分布, 它们又被称为高斯滤波器. 另一类是
直接对后验密度进行分析的全局逼近方法. 例如 [8-
10] 提到的综合 Monte Carlo 数值积分理论和序贯重
要性采样 (SIS) 重采样思想的粒子滤波算法 (PF); [11]
提到的用一簇高斯分布之和来近似逼近后验密度的
高斯和滤波算法 (GSF) 等. 因为概率密度函数 (PDF)
收稿日期: 2010-09-07;修回日期: 2010-12-12.
基金项目: 国家自然科学基金项目(60904097);教育部留学回国人员科研启动基金项目;国防基础科研计划项目
(B1420080209-08).
作者简介: 李振华(1976−), 男, 副教授, 博士, 从事多传感器信息融合、图像处理等研究;宁磊(1986−), 男, 硕士生, 从
事非线性滤波估计的研究.