噪声相关和量测一步随机延时系统高斯滤波器的设计方法技术方案

本发明专利技术提供了一种噪声相关和量测一步随机延时系统高斯滤波器的设计方法，将状态噪声和量测噪声作为状态增量,并基于贝叶斯滤波理论，通过构造正交变换矩阵实现噪声解相关，递推得到高斯滤波框架形式的最优估计算法，采用球径容积法则近似计算所提框架中的高斯加权积分，构建考虑量测一步随机延时和噪声相关的改进容积卡尔曼滤波器框架。本发明专利技术的有益效果是：主要运用于信号处理、组合导航、目标监测与跟踪等领域，对比于传统的只考虑一步延迟的滤波器或者只考虑噪声相关的滤波器，或者已出现的同时考虑一步延迟及相关噪声的滤波器，具有更高的精确性和实用性。更高的精确性和实用性。更高的精确性和实用性。

【技术实现步骤摘要】
噪声相关和量测一步随机延时系统高斯滤波器的设计方法


[0001]本专利技术涉及图像处理领域，尤其涉及一种噪声相关和量测一步随机延时系统高斯滤波器的设计方法。

技术介绍

[0002]滤波器作为系统状态估计的实现方法，广泛应用于信号处理、组合导航、目标监测与跟踪等领域。现有的滤波算法基本思想是在递推贝叶斯估计框架下，利用当前量测信息构造状态后验概率分布函数，进而在最小方差准则下获得系统状态的最优估计。在线性离散系统状态估计问题中，很容易实现状态后验概率分布函数的计算，从而获得后验概率密度函数的最优解。然而，对于一般的非线性离散系统，由于其自身的非线性特性，满足最优解的条件不再成立，通常难以得到后验概率密度函数的闭环解，故求解过程必须采用近似策略，因此，非线性高斯滤波算法通常只能达到对系统状态的次优估计。然而，非线性高斯滤波算法，对系统模型中的噪声特性和量测实时性都做了理想假设，即在量测数据实时获取、过程噪声和量测噪声相互独立等理想条件下，所设计的算法滤波性能才能得到保证；但在实际系统中，量测实时性、噪声特性通常难以满足上述要求。实际情况下，如在网络控制系统中，由于传感器的老化、灵敏度不够等原因，使得系统中量测数据存在随机延时问题；与此同时，由于系统受同一背景噪声源的干扰和连续系统离散化等原因，均会致使过程噪声和量测噪声不再独立而具有相关性。由于实际系统中存在上述非理想条件，故而采用常规高斯滤波器将会出现精度下降，甚至滤波结果严重发散的情况。

技术实现思路

[0003]针对噪声相关和量测一步随机延时下非线性离散系统的状态估计问题，基于理想量测条件的关于先验概率密度函数的标准高斯滤波假设不能适用，同时，一般形式的高斯滤波框架中，要求其统计特性满足系统的状态噪声和量测噪声互相独立。现有的标准高斯滤波器，具有量测随机延时的高斯滤波器，具有噪声相关的高斯滤波器都不适用于实现具有量测随机一步延时和相关噪声的非线性系统状态估计。为了解决上述问题，本专利技术的实施例提供了一种噪声相关和量测一步随机延时系统高斯滤波器的设计方法。基于贝叶斯滤波理论，通过构造正交变换矩阵实现噪声解相关，递推得到高斯滤波框架形式的最优估计算法，采用球径容积法则近似计算所提框架中的高斯加权积分，构建考虑量测一步随机延时和噪声相关的改进容积卡尔曼滤波器框架，利用该改进容积卡尔曼滤波器框架，优化了滤波器的滤波处理。
[0004]本专利技术提供的技术方案带来的有益效果是：对比于传统的只考虑一步延迟的滤波器或者只考虑噪声相关的滤波器，或者已出现的同时考虑一步延迟及相关噪声的滤波器，具有更高的精确性和实用性。
附图说明
[0005]下面将结合附图及实施例对本专利技术作进一步说明，附图中：
[0006]图1是本专利技术实施例中建立的高斯滤波器框架示意图；
[0007]图2是本专利技术实施例中一种噪声相关和量测一步随机延时系统高斯滤波器的设计方法的流程图；
[0008]图3是本专利技术实施例中p
k
＝0.5，S
k
＝0.3时，本申请中CKF和已存在的同类型一步延迟噪声相关CKF状态x的RSME的示意图；
[0009]图4是本专利技术实施例中p
k
＝0.5，S
k
＝0.3时，本申请中CKF和一步延迟CKF状态x的RSME的示意图；
[0010]图5是本专利技术实施例中p
k
＝0.5，S
k
＝0.3时，本申请中CKF和噪声相关CKF状态x的RSME的示意图；
[0011]图6是本专利技术实施例中p
k
＝0.8，S
k
＝0.3时，本申请中CKF和已存在的同类型一步延迟噪声相关CKF状态x的RSME的示意图；
[0012]图7是本专利技术实施例中p
k
＝0.8，S
k
＝0.3时，本申请中CKF和一步延迟CKF状态x的RSM的示意图；
[0013]图8是本专利技术实施例中p
k
＝0.8，S
k
＝0.3时，本申请中CKF和噪声相关CKF状态x的RSME的示意图；
[0014]图9是本专利技术实施例中p
k
＝0.2，S
k
＝0.7时，本申请中CKF和已存在的同类型一步延迟噪声相关CKF状态x的RSME的示意图；
[0015]图10是本专利技术实施例中p
k
＝0.2，S
k
＝0.7时，本申请中CKF和一步延迟CKF状态x的RSM的示意图；
[0016]图11是本专利技术实施例中p
k
＝0.2，S
k
＝0.7时，本申请中CKF和噪声相关CKF状态x的RSME的示意图。
具体实施方式
[0017]为了对本专利技术的技术特征、目的和效果有更加清楚的理解，现对照附图详细说明本专利技术的具体实施方式。
[0018]针对噪声相关和量测一步随机延时下非线性离散系统的状态估计问题，基于理想量测条件的关于先验概率密度函数的标准高斯滤波假设不能适用，同时，一般形式的高斯滤波框架中，要求其统计特性满足系统的状态噪声和量测噪声互相独立。现有的标准高斯滤波器，具有量测随机延时的高斯滤波器，具有噪声相关的高斯滤波器都不适用于实现具有量测随机一步延时和相关噪声的非线性系统状态估计。为了解决上述问题，本专利技术的实施例提供了一种噪声相关和量测一步随机延时系统高斯滤波器的设计方法。基于贝叶斯滤波理论，通过构造正交变换矩阵实现噪声解相关，递推得到高斯滤波框架形式的最优估计算法，采用球径容积法则近似计算所提框架中的高斯加权积分，构建考虑量测一步随机延时和噪声相关的改进容积卡尔曼滤波器框架。
[0019]对于非线性离散系统：
[0020]x
k
＝f(x
k
‑1)+w
k
‑1[0021]z
k
＝h(x
k
)+v
k
[0022]其中x
k
∈R
n
为系统状态量，z
k
∈R
m
为无延迟的量测向量。f(
·
)和h(
·
)为已知的非线性函数，状态噪声w
k
和量测噪声v
k
为相关的零均值高斯白噪声。
[0023]由前所述，系统在某些情况下产生量测延迟，因此实际量测y
k
：
[0024][0025]{γ
k
；k＞1}表示独立同分布的伯努利随机变量(取值为0或1)。
[0026]为了实现上述系统状态的估计，将状态噪声、量测噪声增广为状态量，与系统原状态量一起共同基于实际量测y
k
进行更新。
[0027]如图2所示，本实施例公开了一种噪声相关和量测一步随机延时系统高斯滤波器的设计方法，该设计方法包括：
[0028]S1：将非线性离本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种噪声相关和量测一步随机延时系统高斯滤波器的设计方法，其特征在于：该设计方法包括：S1：将非线性离散系统的状态噪声和量测噪声作为状态增量，基于贝叶斯滤波理论，通过构造正交变换矩阵实现噪声解相关，递推得到高斯滤波框架形式的最优估计算法；S2：根据所述最优估计算法，采用三阶球径容积法则计算出所述高斯滤波器的框架中的高斯加权积分，得到高斯滤波框架下的CKF算法；S3：基于所述CKF算法，构建量测一步随机延时和噪声相关的改进容积卡尔曼滤波器框架，通过使用该改进容积卡尔曼滤波器框架，优化了滤波器的滤波处理。2.如权利要求1所述的一种噪声相关和量测一步随机延时系统高斯滤波器的设计方法，其特征在于：步骤S1中，高斯滤波器框架的建立包括状态和噪声的一步预测、延时量测的一步预测及新状态量的滤波更新3个部分，具体过程如下：(1)状态和噪声的一步预测，即新状态量的一步预测对于非线性离散系统，将状态噪声、量测噪声增广为状态量，经过一步步推导得到k时刻状态x
k
的一步预测和一步预测误差协方差矩阵P
k|k
‑1：：其中f()为非线性函数，为k
‑
1时刻得到的状态噪声w
k
‑1的预测值，为k
‑
1时刻得到的系统状态量x
k
‑1的预测值；P
k
‑
1|k
‑1、和分别为k
‑
1时刻得到的系统状态量x
k
‑1的预测误差协方差矩阵，系统状态量x
k
‑1与状态噪声w
k
‑1的预测误差协方差矩阵和状态噪声w
k
‑1的预测误协差方差矩阵；N()表示高斯分布，d[]表示积分；(2)延时量测的一步预测(2)延时量测的一步预测(2)延时量测的一步预测(2)延时量测的一步预测(2)延时量测的一步预测(2)延时量测的一步预测
其中和分别为实际量测y
k
的一步预测和一步预测误差方差矩阵，p
k
为量测一步随机延时概率；和分别为量测z
k
的一步预测和一步预测误差协方差矩阵；和分别为k
‑
1时刻得到的量测z
k
‑1的预测值和预测误差协方差矩阵；h()为非线性函数；x
k
和x
k
‑1分别表示k和k
‑
1时刻的系统状态量，v
k
‑1表示k
‑
1时刻的量测噪声；为k
‑
1时刻得到的系统状态量x
k
‑1的预测值；表示k
‑
1时刻的量测噪声的一步预测，和分别表示系统状态量x
k
‑1与量测噪声v
k
‑1的一步预测误差协方差矩阵和量测噪声v
k
‑1的一步预测误差协方差矩阵；P
k
‑
1|k
‑1为k
‑
1时刻得到的系统状态量x
k
‑1的预测误差协方差矩阵；P
k|k
‑1表示系统状态量x
k
的一步预测误差协方差矩阵；表示k
‑
1时刻得到的系统状态量x
k
的预测值；R
k
表示量测噪声v
k
的方差系数；(3)新状态量的滤波更新已知k
‑
1时刻的状态后验概率密度p(a
k
‑1|Y
k
‑1)及k时刻的一步预测概率密度p(a
k
|Y
k
‑1)，则k时刻状态量a
k
的后验估计为：的后验估计为：的后验估计为：的后验估计为：其中，表示状态量a
k
的一步预测，K
k
为滤波增益，为状态量a
k
与实际量测y
k
的一步预测误差协方差矩阵，表示实际量测y
k
的一步预测，R
k
和S
k
分别表示k时刻量测噪声v
k
的方差系数和状态噪声w<...

【专利技术属性】
技术研发人员：潘世晔，王新梅，刘珍珠，司洋荣，
申请(专利权)人：中国地质大学武汉，
类型：发明
国别省市：