CN104794101A - 一种分数阶非线性系统状态估计方法 - Google Patents
一种分数阶非线性系统状态估计方法 Download PDFInfo
- Publication number
- CN104794101A CN104794101A CN201510164901.9A CN201510164901A CN104794101A CN 104794101 A CN104794101 A CN 104794101A CN 201510164901 A CN201510164901 A CN 201510164901A CN 104794101 A CN104794101 A CN 104794101A
- Authority
- CN
- China
- Prior art keywords
- overbar
- gamma
- represent
- covariance
- value
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
Landscapes
- Feedback Control In General (AREA)
Abstract
本发明公开了一种新的分数阶非线性系统状态估计方法。首先,给出状态预测值和预测误差协方差初始值。接着,利用近似处理方法求解得到非高斯Lévy噪声下状态近似值和量测近似值,并由此计算系统噪声协方差和量测噪声协方差。然后,利用量测噪声协方差计算最优滤波增益,结合已得到的最优滤波增益计算状态估计值和估计误差协方差。最后,利用当前时刻状态估计值更新下一时刻状态预测值,利用当前时刻估计误差协方差和系统噪声协方差更新下一时刻状态预测协方差。本发明由于解决了带有非高斯Lévy噪声的分数阶非线性系统状态估计问题,从而拓展了分数阶理论的应用范围,而且本发明易于与已有的状态估计软件相结合。
Description
技术领域
本发明涉及一种新的分数阶非线性系统状态估计方法,属于系统分析与处理技术领域。
背景技术
系统分析与处理,旨在研究特定系统结构中各部分(各子系统)的相互作用,系统的对外接口与界面,以及系统整体的行为、功能和局限,从而为系统未来的变迁与有关决策提供参考和依据,其目标之一在于改善决策过程及系统性能,以达到系统的整体最优。在系统分析与处理领域,状态估计起着至关重要的作用。状态估计是根据可获取的量测数据估算动态系统内部状态的方法,由于对系统的输入和输出进行量测而得到的数据只能反映系统的外部特性,而系统的动态规律需要用内部状态变量来描述,因此状态估计对于了解和控制一个系统具有重要意义。
卡尔曼滤波作为一种状态估计的方法被广泛应用于各类系统中。传统的卡尔曼滤波只限于整数阶系统,直到分数阶微积分有了很大的发展,分数阶卡尔曼滤波才慢慢地被学者们研究并逐步得到广泛应用。对非线性系统来说,应用最多的是扩展卡尔曼滤波,不管是整数阶系统还是分数阶系统,扩展卡尔曼滤波都得到广泛应用,但是传统的扩展卡尔曼滤波要求系统噪声和量测噪声均为高斯白噪声,这些过于理想化的要求在现实在很难得到满足。因此,研究非高斯白噪声下的扩展卡尔曼滤波问题具有很重要的理论与现实意义。
发明内容
发明目的:基于以上分析,本发明采用分数阶系统理论,提出一种新的分数阶非线性系统状态估计方法,以期提高状态估计的估计精度。
由于非高斯噪声普遍存在于实际非线性系统中,那么利用传统的扩展卡尔曼滤波方法对这些系统进行状态估计时不能得出理想结果。本发明中所提到的一种新的非线性系统状态估计方法通过近似处理得到状态量和量测量的近似值,并由此推导得出系统噪声协方差和量测噪声协方差,再根据量测噪声协方差和已给定的状态预测误差协方差计算最优滤波增益,从而得到状态估计值和估计误差协方差。该方法与常规动态状态估计程序相结合,能很好地处理非高斯噪声情况下非线性系统的状态估计问题。
技术方案:一种分数阶非线性系统状态估计方法,所述方法是在计算机中依次按以下步骤实现的:
(1)、初始化,包括:设定状态预测量的初值和预测误差协方差的初值;
(2)、对非高斯Lévy噪声进行近似处理,推导得出状态近似值和量测输出近似值,计算步骤为:
式中
表示含有Lévy噪声的状态向量,下标k+1表示第k+1时刻,δ1表示所取的第一个阈值,sign(x)表示符号函数,如果x>0,返回1;如果x<0,返回-1。|x|表示取x的绝对值,表示量测输出向量,表示非线性函数,表示第k时刻状态预测向量,δ2表示所取的第二个阈值,在Ω1中,表示非线性函数,uk表示第k时刻控制输入向量,γj表示
αN表示第N个分数阶阶次值,j=1,2,…,k+1。
(3)、利用已得到的状态近似值和量测输出近似值,计算当前时刻系统噪声协方差和量测噪声协方差,计算步骤为:
式中,表示非线性函数,表示第k时刻状态估计向量,上标T表示转置,表示非线性函数在处的雅克比矩阵,表示第k时刻状态估计误差协方差,表示非线性函数处的雅克比矩阵,表示第k时刻状态预测误差协方差。
(4)、利用当前时刻量测噪声协方差和给定的当前时刻预测误差协方差计算当前时刻最优滤波增益,计算步骤为:
(5)、结合已得到的当前时刻最优滤波增益和给定的当前时刻状态预测值计算当前时刻状态估计值,计算步骤为:
(6)、利用当前时刻最优滤波增益和给定的当前时刻预测误差协方差计算当前时刻估计误差协方差,计算步骤为:
(7)、利用当前时刻状态估计值更新下一时刻状态预测值,计算步骤为:
式中Δ表示分数阶算子,上标α表示分数阶阶次向量。
(8)、利用当前时刻估计误差协方差和系统噪声协方差更新下一时刻状态预测协方差,计算步骤为:
(9)、判断k+1是否大于等于步长L,如果是,结束计算,否则返回步骤(2)进行下一次估计。
传统的扩展卡尔曼滤波算法处理的是整数阶系统,并且要求系统噪声和量测噪声均为高斯白噪声,这意味着系统的阶次必须是整数,噪声的概率分布是正态分布,而且满足它的二阶矩不相关。这些理想化的定义在实际系统中很难得到满足,比如一些粘弹性结构,有耗网络和扩散波等,这些本身就有着分数阶的性质,用整数阶模型很难对这些系统进行建模。另一方面,实际系统中的噪声大多为非高斯白噪声,这样的系统处理起来具有一定的复杂性和挑战性。
本发明提出一种新的分数阶非线性系统状态估计方法,在传统卡尔曼滤波的基础上运用分数阶系统理论,通过非高斯Lévy噪声替代传统的高斯白噪声,利用近似代替方法分解得到Lévy噪声的分解值,重新计算系统噪声和量测噪声的协方差,进而得到最优滤波增益,最后计算状态估计值和估计误差协方差。本发明由于结合分数阶系统理论和非高斯Lévy噪声,能有效地解决整数阶系统的局限性问题和非高斯噪声下的状态估计问题。
附图说明
图1为本发明实施例的方法流程图;
图2为本发明提出的状态估计方法得到的真实值与估计值;
图3为真实值与估计值之间的误差。
具体实施方式
下面结合具体实施例,进一步阐明本发明,应理解这些实施例仅用于说明本发明而不用于限制本发明的范围,在阅读了本发明之后,本领域技术人员对本发明的各种等价形式的修改均落于本申请所附权利要求所限定的范围。
如图1所示,分数阶非线性系统状态估计方法,包括如下步骤:
(1)、初始化。包括:设定状态预测量的初值和预测误差协方差的初值。
(2)、对非高斯Lévy噪声进行近似处理,推导得出状态近似值和量测输出近似值。
(3)利用已得到的状态近似值和量测输出近似值,计算当前时刻系统噪声协方差和量测噪声协方差。
(4)、利用当前时刻量测噪声协方差和给定的当前时刻预测误差协方差计算当前时刻最优滤波增益。
(5)、结合已得到的当前时刻最优滤波增益和给定的当前时刻状态预测值计算当前时刻状态估计值。
(6)、利用当前时刻最优滤波增益和给定的当前时刻预测误差协方差计算当前时刻估计误差协方差。
(7)、利用当前时刻状态估计值更新下一时刻状态预测值。
(8)、利用当前时刻估计误差协方差和系统噪声协方差更新下一时刻状态预测协方差。
(9)、判断k+1是否大于等于步长L,如果是,结束计算,否则返回步骤(2)进行下一次估计。
目前分数阶非线性系统状态估计方法主要采用的是20世纪70年代初由R.E.Larson和Debs提出的卡尔曼滤波算法扩展而来的分数阶卡扩展尔曼滤波算法。
卡尔曼滤波是一种递推回归方法,它的基本思想是根据新量测的数据和由前面一步量测数据计算而得的估计值,再来推算出新的估计值,即
新估计值=旧估计值十修正值
分数阶扩展卡尔曼滤波的本质是求出k时刻状态量真值xk的最优估计值估计的准则是以状态量的估计误差协方差阵最小为目标函数,即:
式中E为求期望函数,xk为k时刻状态量的真值,为k时刻状态量的估计值,上标T表示转置。
由于实际中的一些非线性系统具有分数阶性质,并且存在非高斯白噪声,因此传统扩展卡尔曼滤波算法在应用上受到了很大的限制。分数阶扩展卡尔曼滤波算法具有全局相关性,能较好地体现系统函数发展的历史依赖过程,因而能处理具有分数阶性质的非线性系统,非高斯Lévy噪声通过近似分解可以得到具有高斯白噪声的性质。对于具有非高斯Lévy噪声的分数阶非线性系统,模型描述可以表示为:
式中Δ表示分数阶算子,上标α表示分数阶阶次向量,和分别表示非线性函数,表示含有Lévy噪声的状态向量,下标k表示为第k时刻,uk表示控制输入向量,为量测输出向量,和均为非高斯Lévy噪声,两者协方差矩阵分别表示为和γj表示
式中αN表示第N个分数阶阶次值,j=1,2,…,k+1。
在以上所示分数阶非线性系统模型中,系统噪声和量测噪声均为非高斯Lévy噪声,由于其具有无限方差特性,因此不能用来直接计算协方差矩阵,必须用近似处理的方法得出具有高斯白噪声性质的近似值才可以用来计算协方差矩阵。近似处理后状态向量和量测输出向量可以表示为:
式中
表示含有Lévy噪声的状态向量,下标k+1表示第k+1时刻,δ1表示所取的第一个阈值,sign(x)表示符号函数,如果x>0,返回1;如果x<0,返回-1。|x|表示取x的绝对值,表示量测输出向量,表示非线性函数,表示第k时刻状态预测向量,δ2表示所取的第二个阈值,在Ω1中,表示非线性函数,uk表示第k时刻控制输入向量。
第k时刻系统噪声协方差和量测噪声协方差可分别表示为:
利用给定模型中的表达式进行代换,可得到:
将上述两式分别代入系统和量测噪声协方差阵可得:
式中,表示非线性函数,表示第k时刻状态估计向量,上标T表示转置,表示非线性函数在处的雅克比矩阵,表示第k时刻状态估计误差协方差,表示非线性函数在处的雅克比矩阵,表示第k时刻状态预测误差协方差。
由已得到的k时刻的量测噪声协方差,可得k时刻最优滤波增益为:
系统在k时刻的估计状态量和估计误差协方差可分别表示为:
其中,表示k-1时刻之前所有控制量和量测量的累计值,cov表示取协方差函数。
由k时候变到k+1时刻时,k+1时刻状态预测量和预测误差协方差分别为:
状态量估计步骤总结如下:
步骤1:在当前控制量和量测累计值基础上,计算当前时刻的状态估计量;
步骤2:在当前控制量和量测累计值基础上,计算当前时刻估计误差协方差;
步骤3:在当前控制量和量测累计值基础上,利用下一时刻的状态估计量更新状态预测量;
步骤4:在当前控制量和量测累计值基础上,利用下一时刻的估计误差协方差更新预测误差协方差。
步骤5:在步长范围内,利用下一时刻状态预测量更新状态估计量,利用预测误差协方差更新估计误差协方差。
根据以上步骤,k时刻状态量的估计值为:
式中,表示第k时刻状态的预测量,表示第k时刻最优滤波增益。
k时刻的状态估计误差协方差为:
式中Hk表示非线性函数在处的雅克比矩阵,表示k时刻状态预测误差协方差。
接下来利用状态估计量更新状态预测量,计算步骤为:
其中,表示第k时刻状态估计量,uk表示第k时刻控制输入向量。
最后,利用估计误差协方差更新预测误差协方差,计算步骤为:
其中
Fk表示非线性函数在处的雅克比矩阵,表示第k时刻估计误差协方差。
新的分数阶非线性系统状态估计方法计算公式如下:
a.预测步:
b.估计步:
式中:表示系统噪声协方差,表示量测噪声协方差,并且有
下面介绍本发明的一个实施例:
考虑非高斯Lévy噪声的分数阶非线性系统模型:
式中
ck=[0.2 0.2 0.2]T
uk=1
α=[0.3 1.2 0.8]T
运用本发明所提到的状态更新方法进行仿真,结果如图2所示;真实值与估计值之间的比较如图3所示。
Claims (1)
1.一种分数阶非线性系统状态估计方法,其特征在于,包括以下步骤:
(1)、初始化,包括:设定状态预测量的初值和预测误差协方差的初值;
(2)、对非高斯Lévy噪声进行近似处理,推导得出状态近似值和量测输出近似值,计算步骤为:
式中
表示含有Lévy噪声的状态向量,下标k+1表示第k+1时刻,δ1表示所取的第一个阈值,sign(x)表示符号函数,如果x>0,返回1;如果x<0,返回-1;|x|表示取x的绝对值,表示量测输出向量,表示非线性函数,表示第k时刻状态预测向量,δ2表示所取的第二个阈值,在Ω1中,表示非线性函数,uk表示第k时刻控制输入向量,γj表示
αN表示第N个分数阶阶次值,j=1,2,…,k+1。
(3)、利用已得到的状态近似值和量测输出近似值,计算当前时刻系统噪声协方差和量测噪声协方差,计算步骤为:
式中,表示非线性函数,表示第k时刻状态估计向量,上标T表示转置,表示非线性函数在处的雅克比矩阵,表示第k时刻状态估计误差协方差,表示非线性函数在处的雅克比矩阵,表示第k时刻状态预测误差协方差;
(4)、利用当前时刻量测噪声协方差和给定的当前时刻预测误差协方差计算当前时刻最优滤波增益,计算步骤为:
(5)、结合已得到的当前时刻最优滤波增益和给定的当前时刻状态预测值计算当前时刻状态估计值,计算步骤为:
(6)、利用当前时刻最优滤波增益和给定的当前时刻预测误差协方差计算当前时刻估计误差协方差,计算步骤为:
(7)、利用当前时刻状态估计值更新下一时刻状态预测值,计算步骤为:
式中Δ表示分数阶算子,上标α表示分数阶阶次向量。
(8)、利用当前时刻估计误差协方差和系统噪声协方差更新下一时刻状态预测协方差,计算步骤为:
(9)、判断k+1是否大于等于步长L,如果是,结束计算,否则返回步骤(2)进行下一次估计。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510164901.9A CN104794101A (zh) | 2015-04-08 | 2015-04-08 | 一种分数阶非线性系统状态估计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510164901.9A CN104794101A (zh) | 2015-04-08 | 2015-04-08 | 一种分数阶非线性系统状态估计方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN104794101A true CN104794101A (zh) | 2015-07-22 |
Family
ID=53558902
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510164901.9A Pending CN104794101A (zh) | 2015-04-08 | 2015-04-08 | 一种分数阶非线性系统状态估计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104794101A (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105573355A (zh) * | 2016-01-15 | 2016-05-11 | 杭州电子科技大学 | 分数阶状态空间预测函数控制的储液罐液位控制方法 |
CN106878076A (zh) * | 2017-02-20 | 2017-06-20 | 河海大学 | 计及数据丢包和增益扰动的分数阶网络系统状态估计方法 |
CN106936628A (zh) * | 2017-02-16 | 2017-07-07 | 河海大学 | 一种计及传感器故障的分数阶网络系统状态估计方法 |
CN117713750A (zh) * | 2023-12-14 | 2024-03-15 | 河海大学 | 一种基于分数次幂的一致性卡尔曼滤波状态估计方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20050240347A1 (en) * | 2004-04-23 | 2005-10-27 | Yun-Chun Yang | Method and apparatus for adaptive filter based attitude updating |
CN103927436A (zh) * | 2014-04-04 | 2014-07-16 | 郑州牧业工程高等专科学校 | 一种自适应高阶容积卡尔曼滤波方法 |
CN104112079A (zh) * | 2014-07-29 | 2014-10-22 | 洛阳理工学院 | 一种模糊自适应变分贝叶斯无迹卡尔曼滤波方法 |
CN104283529A (zh) * | 2014-09-29 | 2015-01-14 | 宁波工程学院 | 未知测量噪声方差的平方根高阶容积卡尔曼滤波方法 |
CN104462015A (zh) * | 2014-11-26 | 2015-03-25 | 河海大学 | 处理非高斯Lévy噪声的分数阶线性离散系统状态更新方法 |
-
2015
- 2015-04-08 CN CN201510164901.9A patent/CN104794101A/zh active Pending
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20050240347A1 (en) * | 2004-04-23 | 2005-10-27 | Yun-Chun Yang | Method and apparatus for adaptive filter based attitude updating |
CN103927436A (zh) * | 2014-04-04 | 2014-07-16 | 郑州牧业工程高等专科学校 | 一种自适应高阶容积卡尔曼滤波方法 |
CN104112079A (zh) * | 2014-07-29 | 2014-10-22 | 洛阳理工学院 | 一种模糊自适应变分贝叶斯无迹卡尔曼滤波方法 |
CN104283529A (zh) * | 2014-09-29 | 2015-01-14 | 宁波工程学院 | 未知测量噪声方差的平方根高阶容积卡尔曼滤波方法 |
CN104462015A (zh) * | 2014-11-26 | 2015-03-25 | 河海大学 | 处理非高斯Lévy噪声的分数阶线性离散系统状态更新方法 |
Non-Patent Citations (4)
Title |
---|
DOMINIK SIEROCIUK等: "FRACTIONAL KALMAN FILTER ALGORITHM FOR THE STATES, PARAMETERS AND ORDER OF FRACTIONAL SYSTEM ESTIMATION", 《INTERNATIONAL JOURNAL OF APPLIED MATHEMATICS & COMPUTER SCIENCE》 * |
HODA SADEGHIAN等: "On the general Kalman filter for discrete time stochastic fractional systems", 《MECHATRONICS》 * |
SOREN ASMUSSEN等: "Approximations of small jumps of Levy processes with a view towards simulation", 《JOURNAL OF APPLIED PROBABILITY》 * |
刘济等: "基于UKF和神经网络的一类非线性系统状态估计", 《控制与决策》 * |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105573355A (zh) * | 2016-01-15 | 2016-05-11 | 杭州电子科技大学 | 分数阶状态空间预测函数控制的储液罐液位控制方法 |
CN105573355B (zh) * | 2016-01-15 | 2018-06-19 | 杭州电子科技大学 | 分数阶状态空间预测函数控制的储液罐液位控制方法 |
CN106936628A (zh) * | 2017-02-16 | 2017-07-07 | 河海大学 | 一种计及传感器故障的分数阶网络系统状态估计方法 |
CN106936628B (zh) * | 2017-02-16 | 2019-10-18 | 河海大学 | 一种计及传感器故障的分数阶网络系统状态估计方法 |
CN106878076A (zh) * | 2017-02-20 | 2017-06-20 | 河海大学 | 计及数据丢包和增益扰动的分数阶网络系统状态估计方法 |
CN117713750A (zh) * | 2023-12-14 | 2024-03-15 | 河海大学 | 一种基于分数次幂的一致性卡尔曼滤波状态估计方法 |
CN117713750B (zh) * | 2023-12-14 | 2024-05-17 | 河海大学 | 一种基于分数次幂的一致性卡尔曼滤波状态估计方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104462015A (zh) | 处理非高斯Lévy噪声的分数阶线性离散系统状态更新方法 | |
Ferreira et al. | On the block maxima method in extreme value theory: PWM estimators | |
CN105930640A (zh) | 一种处理Lévy噪声的分数阶卡尔曼滤波方法 | |
CN104794101A (zh) | 一种分数阶非线性系统状态估计方法 | |
CN105184027B (zh) | 一种基于交互式多模型算法的电力负荷建模方法 | |
CN104376581A (zh) | 一种采用自适应重采样的高斯混合无迹粒子滤波算法 | |
Bunch et al. | Improved particle approximations to the joint smoothing distribution using Markov chain Monte Carlo | |
CN108170029B (zh) | Mimo全格式无模型控制器基于偏导信息的参数自整定方法 | |
CN105354860A (zh) | 基于箱粒子滤波的扩展目标CBMeMBer跟踪方法 | |
CN109388778A (zh) | 一种迭代容积点无迹卡尔曼滤波方法 | |
Zhou et al. | A New Adaptive Square‐Root Unscented Kalman Filter for Nonlinear Systems with Additive Noise | |
CN107704426A (zh) | 基于扩展小波神经网络模型的水位预测方法 | |
CN115577776A (zh) | 基态能量的确定方法、装置、设备及存储介质 | |
KR20230008005A (ko) | 정보를 확정하기 위한 방법 및 장치 | |
CN108566178A (zh) | 一种非稳态随机机会网络特征值滤波方法 | |
CN104820788A (zh) | 计及Lévy噪声的分数阶扩展卡尔曼滤波方法 | |
CN110610140A (zh) | 人脸识别模型的训练方法、装置、设备及可读存储介质 | |
CN108181809B (zh) | Miso紧格式无模型控制器基于系统误差的参数自整定方法 | |
Hong et al. | Modelling and control of Hammerstein system using B-spline approximation and the inverse of De Boor algorithm | |
CN101872021B (zh) | Gps双频实时星载数据处理方法 | |
CN115577787A (zh) | 量子振幅估计方法、装置、设备以及存储介质 | |
Jia et al. | A Novel Student’st-based Kalman Filter with Colored Measurement Noise | |
CN111241749A (zh) | 一种基于储备池计算的永磁同步电动机混沌预测方法 | |
CN115577783A (zh) | 量子数据处理方法、装置、设备以及存储介质 | |
CN108681621A (zh) | 基于Chebyshev正交多项式扩展RTS Kalman平滑方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
EXSB | Decision made by sipo to initiate substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
RJ01 | Rejection of invention patent application after publication | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20150722 |