一种具有随机滤波增益变化的有限域滤波器设计方法与流程

本发明涉及一种滤波器设计方法，特别是一种具有随机滤波增益变化的有限域滤波器设计方法。

背景技术：

滤波或状态估计是控制工程和信号处理领域的基础性问题，在航空航天、工业过程控制、自动控制系统中得到了广泛的应用。在实际过程中，通常使用估计误差小于某一上界作为滤波系统的性能指标。比如，在高机动目标的跟踪中，只需要估计误差小于某一上界，而并不需要最小，因此，对方差约束滤波问题的研究具有重要的意义。

网络技术在带来便利的同时，也从以下两个方面对滤波器的设计提出了挑战：一方面，网络技术实现了被控对象和滤波器在地理空间上的分离，滤波参数通过网络传输时，可能会发生微小的随机变化；另一方面，由于数字计算机存储系统的字长是有限的，滤波参数连续信号在进入计算机系统之前必须进行量化，所以会产生截断误差。

现有的滤波器设计研究大都默认滤波器参数能够准确实现，而实际情况中，由于环境变化、仪器精度、未知干扰等因素的影响，滤波器结构参数会发生摄动。keer等证明，对于通过h2，h∞，l1及μ等方法得到的控制器，其参数极其微小的摄动会破坏系统的稳定性。

技术实现要素：

本发明提供一种具有随机滤波增益变化的有限域滤波器设计方法，该方法对具有随机滤波增益变化和量化效应的一类离散时变随机非线性系统，提出基于方差约束下的h∞有限域滤波器设计方法。

一种具有随机滤波增益变化的有限域滤波器设计方法，包括：

步骤1，建立具有随机滤波增益变化和量化效应的离散时变随机非线性系统数学模型；

步骤2，设计方差约束和h∞有限域滤波器；

步骤3，验证方差约束和h∞有限域滤波器设计方法的有效性。

本发明使用方差约束和h∞技术来设计一类离散时变随机非线性系统的有限域滤波器。所设计的滤波器考虑了随机发生的网络诱导滤波增益变化，并利用扇形有界不确定性技术处理量化效应，降低系统运算的复杂度。综合运用schurcomplement和s-procedure引理得到滤波器的lmi表达形式，并给出了有限域滤波器参数求解的迭代算法。同时本发明能够处理网络诱导参数增益变化和量化效应对滤波性能的影响，在滤波器lmis表达形式有解的情况下，保证系统满足h∞性能指标和协方差性能指标。仿真结果验证了算法的有效性，说明达到了预期的设计目标。

下面结合说明书附图对本发明作进一步描述。

附图说明

图1是本发明的方法流程图。

图2是状态x1(k)及其估计示意图。

图3是状态x2(k)及其估计示意图。

图4是输出z(k)及其估计示意图。

具体实施方式

一种具有随机滤波增益变化的有限域滤波器设计方法，按以下步骤实现：

步骤一、建立具有随机滤波增益变化和量化效应的离散时变随机非线性系统数学模型；

步骤二、设计方差约束和h∞有限域滤波器；

步骤三、验证方差约束和h∞有限域滤波器设计方法的有效性。

步骤一中所述建立具有随机滤波增益变化和量化效应的离散时变随机非线性系统数学模型具体为：

考虑定义在k∈[0，n]上的离散时变随机非线性系统：

其中是状态向量，是过程输出，是待估信号。w(k)是定义在概率空间(ω，f，prob)上的一维零均值高斯白噪声序列并且满足a(k)，a1(k)，b(k)，c(k)，d1(k)，d2(k)是维数适当的已知实时变矩阵。

r(k)是取值为1或0且服从如下bernoulli分布的随机变量：

其中是已知常数。

f(k，x(k))是满足下面条件的非线性函数

||f(k，x(k))|²≤θ(k)||g(k)x(k)||²(3)

其中k∈[0，n]，θ(k)＞0为已知正实数，g(k)是已知矩阵。

和是l2[0，n]中的外部扰动且满足

其中w是一已知正定矩阵。

考虑测量信号的量化效应，定义量化器h(·)＝[h1(·)h2(·)…hr(·)]^t，则量化过程的对应关系为：h(y(k))＝[h1(y⁽¹⁾(k))h2(y⁽²⁾(k))…hr(y^(r)(k))]^t

量化器为对数量化器，并满足对称关系，即hj(-y)＝-hj(y)(j＝1，2，…，r)，对每一个hj(·)(1≤j≤r)，量化水平集具有如下形式：

其中χj(j＝1，2，…，r)为量化密度。每个量化水平对应一个区间，则每个量化水平集正好覆盖整个区间。选取如下的量化函数：

其中

由上式可知：hj(y^(j)(k))＝(1+δ(j)(k))y^(j)(k)|δ^(j)(k)|≤δj。所以，可以将量化效应转化成扇形有界的不确定性。

定义δ(k)＝diag{δ⁽¹⁾(k)，δ⁽²⁾(k)…，δ^(r)(k)}，δ＝diag{δ1，δ2，…，δr}，则未知实值时变矩阵满足f(k)f^t(k)≤i。具有量化效应的测量输出可以表示如下：

h(y(k))＝(i+δ(k))y(k)＝(i+δ(k))(b(k)x(k)+d2(k)v(k))(5)

考虑到随机发生的滤波增益变化，采用如下的时变滤波器结构：

其中是状态估计，是估计输出，af(k)，bf(k)，cf(k)是待求的适维滤波器参数矩阵。参数不确定性δaf(k)，δbf(k)，δcf(k)定义为：

δaf(k)=ha(k)δa(k)ea(k)

δbf(k)＝hb(k)δb(k)eb(k)

δcf(k)＝hc(k)δc(k)ec(k)

矩阵ho(k)，eo(k)已知，不确定矩阵δo(k)满足其中o＝a，b，c。

随机变量α(k)，β(k)，γ(k)互不相关且服从bernoulli分布，它们满足

其中是已知实数。

令得到如下增广系统

其中

增广系统(8)的状态协方差矩阵定义为

滤波器设计的目标是使得下列两个条件同时成立：

对于给定的实数γ＞0，矩阵s＞0及初始状态η(0)，系统的h∞性能指标：

其中

对于给定的正定矩阵序列{ψ(k)}0＜k≤n，采样时刻k，估计误差协方差指标满足：

步骤二中所述设计方差约束和h∞有限域滤波器包括以下五个部分

在滤波器设计之前，先给出下面将要用到的引理：

引理1:(schurcomplement)给定常数矩阵s1，s2和s3，其中那么当且仅当

引理2:(s-procedure)n＝n^t，h和e是适当维数的实矩阵，且f^t(t)f(t)≤i。

则不等式n+hfe+(hfe)^t＜0，当且仅当存在一个正实数ε使得n+εhh^t+ε^-1e^te＜0，或者，等价地，

引理3：对于任意向量a,b∈rⁿ，有

ab^t+ba^t≤aa^t+bb^t(14)

引理4：对于任意向量a∈rⁿ，总有

aa^t≤trace(aa^t)i(15)

(一)h∞性能分析，为了方便讨论，做如下假设：

其中，γ为正实数，s为正定矩阵，{τ1(k)}0≤k≤n-1为实数序列，{q(k)}1≤k≤n为正定矩阵序列，且满足满足q(0)≤γ2[i-i]^ts[i-i]，

定义

代入(8)，得到

其中

添加零项到得到

其中

根据(3)，容易得到

对上式两边令k从0到n-1求和，得到

根据上面的不等式可以得到

注意到λ＜0,q(n)＞0及初始条件q(0)≤γ²[i-i]^ts[i-i]，所以j＜0，那么系统的h∞性能指标得到满足。

(二)方差性能分析，为了讨论的方便，做如下假设：

p(k+1)≥φ(p(k))(23)

其中{p(k+1)}0≤k≤n为正定矩阵序列，且满足

由(9)，知

根据引理4，得到

及引理3，有

以，由(24)，得到

运用归纳法，很显然成立，令那么

有如下不等式成立

那么系统的方差约束性能指标得到满足。

(三)在统一的框架下考虑系统方差约束和h∞性能指标，运用schurcomplement引理对(一)(二)的假设条件进行处理，得到假设条件如下等价表述：

其中，(29)对应于(16)，(30)对应于(23)。

(四)方差约束和h∞有限域非脆弱滤波器设计，在(三)工作的基础上，综合运用s-procedure和schurcomplement引理，消除矩阵不等式(29)(30)中的非线性项，从而得到一组线性矩阵不等式(lmis)，如下所示：

其中，

ξ5＝[a1(k)00000]，ξ6＝[c(k)-cf(k)0000]

γ为正实数，s为正定矩阵，{ψ(k)}0≤k≤n+1为方差上界矩阵序列，

及为正定矩阵序列，

{τ1(k)}0≤k≤n，{∈1(k)}0≤k≤n，{∈2(k)}0≤k≤n，{∈3(k)}0≤k≤n，{∈a(k)}0≤k≤n，

{∈b(k)}0≤k≤n，{∈c(k)}0≤k≤n，为实数序列，{p3(k)}1≤k≤n+1，{af(k)}0≤k≤n，{bf(k)}0≤k≤n，{cf(k)}0≤k≤n为实值矩阵序列，且满足

下面介绍详细设计过程，首先对变量p(k)和q(k)做如下分解：

因此，条件q(0)≤γ²[i-i]^ts[i-i]和与(34)等价。

为了对(29)中的不确定参数δ(k)进行估计，将(29)重写为如下等价式：

n(k)+h(k)f(k)e(k)+(h(k)f(k)e(k))t＜0(36)

其中

根据s-procedure引理，得到

很明显在(37)只存在不确定性参数δa(k)，δb(k)，δc(k)，为了估计它们，可以将上式重写为

其中

那么，由schurcomplement和s-procedure引理，(29)等效于(31)。类似地，可以得到(30)等效于(32)，至此完成了方差约束和h∞有限域非脆弱滤波器的设计。

(五)h∞和协方差有限域滤波器设计算法(nfd)概括如下

步骤1，对于给定的正实数γ＞0，正定矩阵s＞0,初始状态误差及方差约束矩阵ψ(0)。选取合适的初始值{q1(0)，q2(0)，q3(0)，p1(0)，p2(0)，p3(0)}满足初始条件(32)，令k＝0；

步骤2，时刻k，通过求解线性矩阵不等式组(31)-(33)得到矩阵及滤波器矩阵参数af(k)，bf(k)，cf(k)；

步骤3，令k＝k+1，调用更新表达式(35)得

步骤4，如果k＜n,那么执行步骤2,否则执行下一步；

步骤5，结束。

步骤三中验证方差约束和h∞有限域滤波器设计方法有效性的具体方式如下：

通过给出一个数值仿真实例，利用matlab/lmi工具箱对所设计的滤波器参数进行求解，并验证方差约束和h∞性能指标。

考虑如下离散系统：

零均值噪声ω(k)服从标准正态分布，非线性函数f(k，x(k))及外部扰动w(k)，v(k)如

下

随机变量r(k)，α(k)，β(k)，γ(k)的期望为

指数量化器h(·)的参数为χ1＝0.3，χ2＝0.6且不确定参数f(k)满足f^t(k)f(k)≤i。

滤波器增益变化中的已知矩阵参数ho(k)，eo(k)(o＝a，b，c)如下

hc(k)＝0.2，hc(k)＝[0.1+0.2exp(-k)0]

且不确定参数δo(k)满足

初始状态x(0)＝[0.26-0.4]^t，初始估计状态正实数γ＝0.5，正定矩阵s＝diag{8，8},ψ(0)＝[1.4036，-0.0144；-0.0144，1.4036]，

ψ(k)＝diag{0.2，0.2}（k＝1，…，n)，q1(0)＝diag{1，1},q2(0)＝diag{1，1},q3(0)＝0及

验证结果如图2-4所示，图2、3分别给出了状态变量x1(k)-x2(k)及它们的估计量图4给出了输出z(k)和它的估计量通过对仿真结果的计算得到h∞性能指标j＝-0.1236，并验证了对所有的k＝0，…，n成立。仿真结果说明了本发明所提出的滤波器设计方法的有效性。

综上所述，本发明给出了一类离散时变随机非线性系统的方差约束和h∞有限域滤波器设计方法，所设计的滤波器具有随机发生的滤波增益变化且受量化作用影响。随机非线性现象是由服从bernoulli分布规律的随机变量描述的在两种非线性扰动之间的二元切换；滤波增益的随机变化用来描述受网络带宽影响发生的滤波器参数的微小随机变化；量化器采用指数型，并通过一定的方法将量化不确定转化为扇形有界不确定以降低问题的复杂性。在非线性干扰和外部扰动为非零均值时，成功进行方差约束设计。通过解一组递推线性矩阵不等式，给出了使滤波误差系统同时满足方差约束和h∞性能指标的滤波器存在的充分条件。最后，通过一个仿真实例说明了所提出的滤波器设计方法的有效性。