公开/公告号CN112653410A
专利类型发明专利
公开/公告日2021-04-13
原文格式PDF
申请/专利权人 中国地质大学(武汉);
申请/专利号CN202011478686.7
申请日2020-12-15
分类号H03H7/01(20060101);H03H7/075(20060101);
代理机构42238 武汉知产时代知识产权代理有限公司;
代理人郝明琴
地址 430000 湖北省武汉市洪山区鲁磨路388号
入库时间 2023-06-19 10:35:20
技术领域
本发明涉及图像处理领域,尤其涉及一种噪声相关和量测一步随机延时系统高斯滤波器的设计方法。
背景技术
滤波器作为系统状态估计的实现方法,广泛应用于信号处理、组合导航、目标监测与跟踪等领域。现有的滤波算法基本思想是在递推贝叶斯估计框架下,利用当前量测信息构造状态后验概率分布函数,进而在最小方差准则下获得系统状态的最优估计。在线性离散系统状态估计问题中,很容易实现状态后验概率分布函数的计算,从而获得后验概率密度函数的最优解。然而,对于一般的非线性离散系统,由于其自身的非线性特性,满足最优解的条件不再成立,通常难以得到后验概率密度函数的闭环解,故求解过程必须采用近似策略,因此,非线性高斯滤波算法通常只能达到对系统状态的次优估计。然而,非线性高斯滤波算法,对系统模型中的噪声特性和量测实时性都做了理想假设,即在量测数据实时获取、过程噪声和量测噪声相互独立等理想条件下,所设计的算法滤波性能才能得到保证;但在实际系统中,量测实时性、噪声特性通常难以满足上述要求。实际情况下,如在网络控制系统中,由于传感器的老化、灵敏度不够等原因,使得系统中量测数据存在随机延时问题;与此同时,由于系统受同一背景噪声源的干扰和连续系统离散化等原因,均会致使过程噪声和量测噪声不再独立而具有相关性。由于实际系统中存在上述非理想条件,故而采用常规高斯滤波器将会出现精度下降,甚至滤波结果严重发散的情况。
发明内容
针对噪声相关和量测一步随机延时下非线性离散系统的状态估计问题,基于理想量测条件的关于先验概率密度函数的标准高斯滤波假设不能适用,同时,一般形式的高斯滤波框架中,要求其统计特性满足系统的状态噪声和量测噪声互相独立。现有的标准高斯滤波器,具有量测随机延时的高斯滤波器,具有噪声相关的高斯滤波器都不适用于实现具有量测随机一步延时和相关噪声的非线性系统状态估计。为了解决上述问题,本发明的实施例提供了一种噪声相关和量测一步随机延时系统高斯滤波器的设计方法。基于贝叶斯滤波理论,通过构造正交变换矩阵实现噪声解相关,递推得到高斯滤波框架形式的最优估计算法,采用球径容积法则近似计算所提框架中的高斯加权积分,构建考虑量测一步随机延时和噪声相关的改进容积卡尔曼滤波器框架,利用该改进容积卡尔曼滤波器框架,优化了滤波器的滤波处理。
本发明提供的技术方案带来的有益效果是:对比于传统的只考虑一步延迟的滤波器或者只考虑噪声相关的滤波器,或者已出现的同时考虑一步延迟及相关噪声的滤波器,具有更高的精确性和实用性。
附图说明
下面将结合附图及实施例对本发明作进一步说明,附图中:
图1是本发明实施例中建立的高斯滤波器框架示意图;
图2是本发明实施例中一种噪声相关和量测一步随机延时系统高斯滤波器的设计方法的流程图;
图3是本发明实施例中p
图4是本发明实施例中p
图5是本发明实施例中p
图6是本发明实施例中p
图7是本发明实施例中p
图8是本发明实施例中p
图9是本发明实施例中p
图10是本发明实施例中p
图11是本发明实施例中p
具体实施方式
为了对本发明的技术特征、目的和效果有更加清楚的理解,现对照附图详细说明本发明的具体实施方式。
针对噪声相关和量测一步随机延时下非线性离散系统的状态估计问题,基于理想量测条件的关于先验概率密度函数的标准高斯滤波假设不能适用,同时,一般形式的高斯滤波框架中,要求其统计特性满足系统的状态噪声和量测噪声互相独立。现有的标准高斯滤波器,具有量测随机延时的高斯滤波器,具有噪声相关的高斯滤波器都不适用于实现具有量测随机一步延时和相关噪声的非线性系统状态估计。为了解决上述问题,本发明的实施例提供了一种噪声相关和量测一步随机延时系统高斯滤波器的设计方法。基于贝叶斯滤波理论,通过构造正交变换矩阵实现噪声解相关,递推得到高斯滤波框架形式的最优估计算法,采用球径容积法则近似计算所提框架中的高斯加权积分,构建考虑量测一步随机延时和噪声相关的改进容积卡尔曼滤波器框架。
对于非线性离散系统:
x
z
其中x
由前所述,系统在某些情况下产生量测延迟,因此实际量测y
{γ
为了实现上述系统状态的估计,将状态噪声、量测噪声增广为状态量,与系统原状态量一起共同基于实际量测y
如图2所示,本实施例公开了一种噪声相关和量测一步随机延时系统高斯滤波器的设计方法,该设计方法包括:
S1:将非线性离散系统的状态噪声和量测噪声作为状态增量,基于贝叶斯滤波理论,通过构造正交变换矩阵实现噪声解相关,递推得到高斯滤波框架形式的最优估计算法;
S2:根据所述最优估计算法,采用三阶球径容积法则计算出所述高斯滤波器的框架中的高斯加权积分,得到高斯滤波框架下的CKF算法;
S3:基于所述CKF算法,构建量测一步随机延时和噪声相关的改进容积卡尔曼滤波器框架,通过使用该改进容积卡尔曼滤波器框架,优化了滤波器的滤波处理。
具体如下:
1.高斯滤波器框架
所提噪声相关和量测一步随机延时下的高斯滤波框架的示意图如图1所示,高斯滤波器的推导将分为三部分构成:(1)状态和噪声的一步预测,即新状态量的一步预测;(2)延时量测的一步预测;(3)新状态量的滤波更新。
(1)状态和噪声的一步预测,即新状态量的一步预测
对于非线性离散系统,将状态噪声、量测噪声增广为状态量,经过一步步推导得到得到k时刻系统状态a
其中,x
其中
进而
其中,状态x
其中f()为非线性函数,
(2)延时量测的一步预测
其中
(3)新状态量的滤波更新
上述非线性离散系统中,已知k-1时刻的状态后验概率密度p(a
其中,
其中,
2.数值积分近似实现
对上述高斯滤波框架,采用三阶球径容积法则近似计算所提框架中的高斯加权积分,得到高斯滤波框架下的CKF算法。算法主要包括时间更新和量测更新两部分。具体实现过程如下:
(1)时间更新
Step 1:对k-1时刻得到的状态量a
Step 2:容积点计算:
其中,ξ
Step 3:容积点扩散,计算容积点经非线性函数f(·)和h(·)的传播结果:
Step 4:计算式中状态一步预测
其中,ξ
(2)量测更新
Step 1:对系统状态量x
Step 2:容积点计算:
Step 3:容积点扩散,计算容积点χ
Z
Step 4:计算式量测估计中中的高斯加权积分部分:
其中,U
将上式代入所述高斯滤波框架形式的最优估计算法中,可得噪声相关和一步随机延时下系统状态a
仿真实验:
以含有三角函数和平方函数的模型的为例,该模型为普遍非线性非平稳增长模型[],通过分别比较本实施例提出的CKF-RDMCN,和已出现的一步延迟噪声相关CKF算法、和一步延迟CKF算法以及噪声相关CKF算法,可显示本实施例所提出的CKF-RDMCN的优势。该模型如下:
其中w
表1多种情况下,4种算法的RMSE的均值
当p
1)由图3-11和表一可知,本实施例所提的改进的CKF算法在三种不同的情况下:RMSE的结果基本上与已存在一步延迟及噪声相关的CKF算法持平,RMSE的均值甚至还要比其小点;跟一步延迟CKF算法和噪声相关CKF算法相比,RMSE的结果偏低,说明本文的所提的CKF算法精度更高。且随着噪声相关系数的增大,其RMSE结果变化不大,说明在噪声相关条件下,所提出CKF算法具有更好的鲁棒性和更优的滤波精度。
本发明的有益效果是:主要运用于信号处理、组合导航、目标监测与跟踪等领域,对比于传统的只考虑一步延迟的滤波器或者只考虑噪声相关的滤波器,或者已出现的同时考虑一步延迟及相关噪声的滤波器,具有更高的精确性和实用性。
以上所述仅为本发明的较佳实施例,并不用以限制本发明,凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
机译: 用于建立高斯或排列类型的随机数值模型的迭代方法,例如在油气勘探中,涉及将与伪数据的公差相关的第一高斯白噪声与第二高斯白噪声相结合
机译: 用于建立高斯或排列类型的随机数值模型的迭代方法,例如在油气勘探中,涉及将与伪数据的公差相关的第一高斯白噪声与第二高斯白噪声相结合
机译: 利用局部统计量的噪声检测方法及通过局部统计量的改进广义高斯滤波器去除噪声的方法和系统