一种基于扩展核递归最大相关熵准则的系统状态自适应估计方法技术方案

本发明专利技术属于自适应信号处理技术领域，提供一种基于扩展核递归最大相关熵准则的自适应系统状态估计方法，该方法能够在非高斯噪声条件下，对缓慢时变的非线性系统的状态进行估计和跟踪。该方法包括：1)提出适合于刻画非线性、缓慢时变的系统的数学模型与基于最大相关熵准则的代价函数；2)在再生核希尔伯特空间中，对系统状态进行迭代；3)利用核技巧，在输入空间中对系统状态进行迭代求解，实现系统状态的自适应估计。实验证明本发明专利技术算法性能良好，在真实的工程应用中，具有较好的应用前景。

【技术实现步骤摘要】

本专利技术属于自适应信号处理
，涉及到非高斯噪声条件下，非线性、缓慢时变系统的状态自适应估计方法，特别涉及到一种基于扩展核递归最大相关熵准则的自适应系统状态估计方法。
技术介绍
自适应滤波技术是信号处理领域的理论之一，自从自适应滤波算法被提出以来，就被各种实际工程广泛地应用。例如，自适应时延估计、自适应噪声抵消、自适应系统辨识、自适应信道均衡、自适应混响抵消等等。可见，在日常生活中，自适应滤波方法是极其重要的。在理论研究和应用方法层面，针对高斯分布噪声和非高斯分布噪声的情况，针对线性系统与非线性系统的情况，乃至针对系统状态时变与时不变等情况，都有不同的技术成果，但是针对非高斯分布噪声条件下，状态缓慢时变的非线性系统这种非常复杂的情况，仍存在理论研究和技术方法的空白，这也是本专利技术要解决的问题。
技术实现思路
针对现有技术的不足，本专利技术提供一种基于扩展核递归最大相关熵准则的系统状态自适应估计方法,该方法是一种对非高斯噪声具有较强抑制能力的自适应方法，并用该方法对非线性的缓慢时变系统状态进行估计。本专利技术提出适合于刻画非线性、缓慢时变的系统的数学模型与基于最大相关熵准则的代价函数，然后在再生核希尔伯特空间中，对系统状态进行迭代，最后利用核技巧，在输入空间中对系统状态进行迭代求解，实现系统状态的自适应估计。为了达到上述目的,本专利技术的技术方案为:一种基于扩展核递归最大相关熵准则的系统状态自适应估计方法,包括以下步骤：第一步,提出适合于刻画非线性、缓慢时变的系统的数学模型与基于最大相关熵准则的代价函数。(1)提出适合于刻画非线性、缓慢时变的数学模型；(2)针对具体的应用，例如慢衰落瑞利信道，选择适当的状态转移矩阵，将数学模型具体化；(3)针对该模型，基于最大相关熵原则，提出代价函数。第二步，在再生核希尔伯特空间中，对系统状态进行迭代。(1)设定初始条件；(2)计算当前步骤下输出信号与期望信号的误差；(3)计算当前步骤下的增益；(4)对系统状态进行更新；(5)对状态误差相关矩阵进行更新，并重新进行迭代。第三步，利用核技巧，在输入空间中对系统状态进行迭代求解，实现系统状态的自适应估计。(1)设定初始条件；(2)计算当前步骤下输出信号与期望信号的误差；(3)对系统状态进行更新；(4)对其他变量、矩阵进行更新。本专利技术的有益效果为：该方法能够在非高斯噪声条件下，对缓慢时变的非线性系统的状态进行估计和跟踪。实验证明本专利技术算法性能良好，在真实的工
程应用中，例如：慢衰落信道参数估计问题，具有较好的应用前景。附图说明图1是该方法的流程图；图2是洛伦兹时间序列的三维图；图3是洛伦兹时间序列的X轴坐标状态变化图；图4是本专利技术应用于洛伦兹时间序列估计学习曲线；图5是本专利技术应用于缓慢衰落瑞利信道的状态估计学习曲线。具体实施方式为使本专利技术实施例的目的、技术方案及其优点更加清楚，下面结合本专利技术实施例中的附图，对本专利技术实施例中的技术方案进行清楚完整的描述，整体算法流程图如图1所示：第一步，得到系统的数学模型和代价函数1)用公式(1)状态方程和公式(2)观测方程表示系统的数学模型，数学模型适合于刻画非线性、缓慢时变的系统； x ~ ( n ) = A x ~ ( n - 1 ) + n ~ 1 ( n ) - - - ( 1 ) ]]> d ~ ( n ) = h ( u ~ ( n ) , x ~ ( n ) ) + n ~ 2 ( n ) - - - ( 2 ) ]]>其中，是系统的状态，A是状态变化矩阵，系统的输入，是系统的输出，是状态噪声，是观测噪声，h表示与输入和系统状态有关的非线性变换；2)针对具体应用，选择状态转移矩阵，将数学模型具体化，具体化后的数学模型为：x(n)＝ηx(n-1)+n1(n) (3)其中，η取值0.9999或者更接近1，表示系统转移参数；表示对输入信号的非线性变换，该非线性变换是未知和非显式的；在上述具体化后的模型中，公式(3)、公式(4)中没有“～”号的各变量和向量表示的物理含义与公式(1)、(2)中有“～”号的物理含义相同；3)针对具体化后的数模模型，基于最大相关熵原则，得到公式(5)所示的代价函数；其中，κ表示满足Mercer条件的核函数，通常选用高斯核函数表示，高斯核函数的定义为： κ σ ( e ( m ) ) = exp ( - e 2 ( m ) / 2 σ 2 ) / 2 π σ - - - ( 6 ) ]]>其中，表示观测误差，σ表示核长；将公式(5)转化为下述等价关系： m i n Σ m = 1 n κ σ ( e ( m ) ) | e本文档来自技高网...

【技术保护点】
一种基于扩展核递归最大相关熵准则的系统状态自适应估计方法，其特征在于，包括以下步骤：第一步，得到系统的数学模型和代价函数1)用公式(1)状态方程和公式(2)观测方程表示系统的数学模型，数学模型适合于刻画非线性、缓慢时变的系统；x~(n)=Ax~(n-1)+n~1(n)---(1)]]>d~(n)=h(u~(n),x~(n))+n~2(n)---(2)]]>其中，是系统的状态，A是状态变化矩阵，系统的输入，是系统的输出，是状态噪声，是观测噪声，h表示与输入和系统状态有关的非线性变换；2)针对具体应用，选择状态转移矩阵，将数学模型具体化，具体化后的数学模型为：x(n)＝ηx(n‑1)+n1(n)   (3)其中，η取值0.9999或者更接近1，表示系统转移参数；表示对输入信号的非线性变换，该非线性变换是未知和非显式的；在上述具体化后的模型中，公式(3)、公式(4)中没有“～”号的各变量和向量表示的物理含义与公式(1)、(2)中有“～”号的物理含义相同；3)针对具体化后的数模模型，基于最大相关熵原则，得到公式(5)所示的代价函数；s.t.x(m)＝ηx(m‑1)+n1(m)   (5)其中，κ表示满足Mercer条件的核函数，通常选用高斯核函数表示，高斯核函数的定义为：κσ(e(m))=exp(-e2(m)/2σ2)/2πσ---(6)]]>其中，x(m)表示观测误差，σ表示核长；将公式(5)转化为下述等价关系：s.t.x(m)＝ηx(m‑1)+n1(m)   (7)为提升方法的鲁棒性，将公式(7)改写为：s.t.x(m)＝ηx(m‑1)+n1(m)   (8)其中，λ表示遗忘因子，是正则化因子，q是状态噪声和观测噪声的能量之比；第二步，在由高斯核函数所引入的再生核希尔伯特空间中，对系统状态进行迭代；1)设定初始条件；在m＝0时，设x(0)＝0，P(0)＝I；2)当迭代开始后，即m≥1时，按照公式计算当前步骤下输出信号与期望信号间的观测误差e(m)；3)根据公式计算得到当前步骤下的增益k(m)；其中：r^(m)=r(m)κσ-1(e(m))---(10)]]>4)根据公式x(m)＝ηx(m‑1)+k(m)e(m)更新当前的系统状态；5)根据公式更新状态误差的相关矩阵，并按照步骤2)、步骤3)、步骤4)、步骤5)的重新开始下一次迭代；第三步，利用核技巧，在u(n)所在的输入空间中对系统状态x(n)进行迭代求解，实现系统状态的自适应估计；1)设定初始条件；在m＝1时，设：a(1)=ηd(1)λλ~+κσ(u(1),u(1))---(11)]]>ρ(1)=η2λ~-1+λq---(12)]]>Q(1)=η2λ~-1λλ~+κσ(u(1),u(1))---(13)]]>其中，a(m)，ρ(m)，Q(m)表示中间变量，向量或矩阵；2)按照如下公式计算当前步骤下输出信号与期望信号的误差：h(m)＝[κσ(u(m),u(1)),…,κσ(u(m),u(m‑1))]T   (14)z(m)＝Q(m‑1)h(m)   (15)e(m)＝d(m)‑h(m)Ta(m‑1)   (16)其中，h(m)表示由高斯核函数组成的向量，z(m)是中间向量；3)根据如下公式对系统当前状态进行更新：r(m)＝λn+κσ(e(m))ρ(m‑1)κσ(u(m),u(m‑1))‑hT(m)z(m)   (17)r^(m)=r(m)κσ-1(e(m))---(10)]]>a(m)=ηa(m-1)-z(m)r^-1(m)e(m)ρ(m-1)r^-1(m)e(m)---(18)]]>x(m)＝H(m)a(m)   (19)其中，是由映射后的输入组成的向量；4)对其他变量、矩阵按照公式(20)、公式(21)进行更新，并继续进行迭代直到结束。ρ(m)＝η2ρ(m‑1)+λmq   (20)Q(m)=η2r~-1(m)Q(m-1)r~(m)+z(m)zT(m)-ρ(m-1)z(m)-ρ(m-1)zT(m)ρ2(m-1)---(21)]]>...

【技术特征摘要】
1.一种基于扩展核递归最大相关熵准则的系统状态自适应估计方法，其特征在于，包括以下步骤：第一步，得到系统的数学模型和代价函数1)用公式(1)状态方程和公式(2)观测方程表示系统的数学模型，数学模型适合于刻画非线性、缓慢时变的系统； x ~ ( n ) = A x ~ ( n - 1 ) + n ~ 1 ( n ) - - - ( 1 ) ]]> d ~ ( n ) = h ( u ~ ( n ) , x ~ ( n ) ) + n ~ 2 ( n ) - - - ( 2 ) ]]>其中，是系统的状态，A是状态变化矩阵，系统的输入，是系统的输出，是状态噪声，是观测噪声，h表示与输入和系统状态有关的非线性变换；2)针对具体应用，选择状态转移矩阵，将数学模型具体化，具体化后的数学模型为：x(n)＝ηx(n-1)+n1(n) (3)其中，η取值0.9999或者更接近1，表示系统转移参数；表示对输入信号的非线性变换，该非线性变换是未知和非显式的；在上述具体化后的模型中，公式(3)、公式(4)中没有“～”号的各变量和向量表示的物理含义与公式(1)、(2)中有“～”号的物理含义相同；3)针对具体化后的数模模型，基于最大相关熵原则，得到公式(5)所示的代价函数；s.t.x(m)＝ηx(m-1)+n1(m) (5)其中，κ表示满足Mercer条件的核函数，通常选用高斯核函数表示，高斯核函数的定义为： κ σ ( e ( m ) ) = exp ( - e 2 ( m ) / 2 σ 2 ) / 2 π σ - - - ( 6 ) ]]>其中，x(m)表示观测误差，σ表示核长；将公式(5)转化为下述等价关系：s.t.x(m)＝ηx(m-1)+n1(m) (7)为提升方法的鲁棒性，将公式(7)改写为：s.t.x(m)＝ηx(m-1)+n1(m) (8)其中，λ表示遗忘因子，是正则化因子，q是状态噪声和观测噪声的能量之比；第二步，在由高斯核函数所引入的再生核希尔伯特空间中，对系统状态进行迭代；1)设定初始条件；在m＝0时，设x(0)＝0，P(0)＝I；2)当迭代开始后，即m≥1时，按照公式计算当前步骤下输出信号与期望信号间的观测误差e(m)；3)根据公式计算得到当前步骤下的增益k(m)；其中： r ^ ( m ) = r ( m ) κ σ - 1 ( e ( m ) ) - - - ( 10 ) ]]>4)根据公式x(m)＝ηx(m-1)+k(m)e(m)更新当前的系统状态；5)根据公式更新状态误差的相关矩阵，并按照步骤2)、步骤3)、步骤4)、步骤5)的重新开始下一次迭代；第三步，利用核技巧，在u(n)所在的输入空间中对系统状态x(n)进行迭代求解，实现系统状态的自适应估计；1)设定初始条件；在m＝1时，设： a ( 1 ) = η d ( 1 ) λ λ ~ + κ σ ( u ( 1 ) , u ( 1 ) ) - - - ( 11 ) ]]> ρ ( 1 ) = η 2 λ ~ - 1 + λ q - - - ( 12 ) ]]> Q ( 1 ) = η 2 λ ~ - 1 λ λ ~ + κ σ ( u ( 1 ) , u ( 1 ) ) - - - ( 13 ) ]]>其中，a(m)，ρ(m)，Q(m)表示中间变量，向量或矩阵；2)按照如下公式计算当前步骤下输出信号与期望信号的误差：h(m)＝[κσ(u(m),u(1)),…,κσ(u(m),u(m-1))]T (14)z(m)＝Q(m-1)h(m) (15)e(m)＝d(m)-h(m)Ta(m-1...

【专利技术属性】
技术研发人员：邱天爽，栾声扬，朱永杰，于玲，张金凤，刘涛，马济通，宋爱民，
申请(专利权)人：大连理工大学，
类型：发明
国别省市：辽宁;21