1. 引言
滤波问题在航天、经济、军事等领域被广泛应用而一直受到人们的关注。Kalman和
滤波是目前比较常见的滤波方法。Kalman滤波适用于系统中的扰动是已知的高斯白噪声或谱密度的噪声 [1] 。Kalman滤波的缺点是如果对噪声模型或测量的噪声的了解不够准确,最终会造成滤波发散的情况。
滤波适用于测量的噪声是有界能量的任意信号,
滤波不必要求噪声是高斯白噪声。因此,
滤波在很多领域上得到了广泛应用。在
滤波器设计中被广泛使用的有Riccati方程方法和LMI方法。其中,Riccati方程方法主要用在带有不确定项是范数有界的系统上,解法比较简单,但是Riccati方程方法用的是迭代法,没有规范的解法,因此保守性比较大。但LMI方法实用性比较强,具有较小的保守性,因此得到了广泛使用 [2] 。
在实际工程中会遇到各种复杂的问题,这些问题不能仅用一个连续(离散)的系统来建模分析,而需要一些连续(离散)的子系统根据控制信号进行“切换”,并得到一种混合系统模型。因此,切换系统作为一种混杂动态系统,有着很重要的应用价值 [3] 。时滞是造成系统不稳定的重要原因之一,因此对时滞切换系统的同步和稳定性研究引起了许多学者的关注。文献 [4] 中研究了不确定时变时滞切换系统的控制与设计问题。文献 [5] 讨论了基于平均停留时间的不确定开关时滞系统的时滞相关稳定性分析。在 [6] 中,主要研究了时变时滞切换系统的时滞相关的有限时间和L2-增益分析问题。对于中立系统的稳定性研究是时滞系统稳定性研究的重点和难点,国内外很多学者对切换中立系统的稳定性研究进行了深入的研究。实际上,很多控制系统的模型可以转化为中立系统。针对切换中立系统的研究也取得了一定的成果 [7] [8] 。 [7] 给出了一类切换线性中立型系统的镇定问题。 [8] 讨论了具有混合时滞的不确定切换中立型系统的鲁棒有限时间
控制问题。
近几年
滤波研究引起了很多学者的兴趣并取得了丰富的研究成果 [8] [9] [10] [11] [12] 。在 [13] 中,利用Riccati方程方法研究了时变时滞离散系统的
滤波问题。虽然Riccati方程方法的解法比较简单,但是Riccati方程方法用的是迭代法,所以没有规范的解法,同时Riccati方程方法在给出滤波器参数的时候,需要提前给定待定参数,这些参数需要根据人为经验来给定待定参数,这些人为给定的参数对系统的稳定性分析带来比较大的保守性。在 [14] 中,利用LMI方法研究了一类切换线性中立系统的
滤波器设计问题。通过研究发现,不论是利用经典的最优控制研究系统稳定性,还是利用其设计控制器,都可能会导致控制器有脆弱现象出现,从而导致加入控制器后得到的闭环系统性能下降,甚至难以保持稳定。因此,非脆弱滤波问题引起了学者们的兴趣,取得了一定的研究成果 [15] [16] [17] [18] 。 [18] 研究了时变时滞不确定系统的非脆弱
滤波问题。
本文受到以上研究的启发,研究了不确定时变时滞切换中立系统的非脆弱
滤波器设计。对于所容许的不确定性和时变时滞,设计一个非脆弱
滤波器,使闭环系统渐近稳定并满足
性能指标。最后给出数值实验,证明了所提结果的有效性和可行性。本文提出的线性矩阵不等式具有决策变量少、参数多的特点,可以最大程度地降低结果的保守性。
本文的其余部分安排如下:第2节阐述了具有不确定时变时滞切换中立系统的非脆弱
滤波问题的表述和前言;第3节介绍了我们的主要结果以及非脆弱
滤波器设计;第4节给出了数值模拟;第5节为本文的结论。
2. 问题提出和前言
考虑如下形式的不确定时变时滞切换中立系统:
(2.1)
其中,
为切换系统的状态向量;
为输出向量;
为有限能量的外部干扰输入向量;
为带估计的信号向量;
和
为已知适维的常值矩阵;
为实值向量连续函数;
表示系统的切换规则,其中n表示为切换系统的n个切换子系统;
和
为系统的不确定项并满足如下形式:
(2.2)
其中,
、
、
和
是已知为矩阵,
是具有Lebesgue可测元素的未知矩阵且满足条件:
和
分别为时变时滞和中立时滞且为连续函数满足:
(2.3)
其中,
和
是已知常数;
是
上的一个连续实值初始函数。其中
。
现在我们考虑一个系统(2.1)的非脆弱滤波器,其形式如下:
(2.4)
其中,
是滤波器状态向量,
是滤波器输出,
是适维的滤波器参数矩阵,不确定项
表示不确定增益摄动且满足条件:
(2.5)
其中,
是适维常矩阵,
是未知矩阵且满足条件
。
令
则闭环系统可描述如下形式:
(2.6)
其中,
(2.7)
对于这种情况下的系统,本文要解决的问题是针对不确定时变时滞切换中立型系统(2.1)寻找满足给定
性能指标
的非脆弱
滤波器存在的充分条件。
引理2.1 (schur补引理)给定矩阵
,
和
是具有适当维数的对称矩阵,那么下列LMIs是等价的:
引理2.2 E、H和M是适当维数的实数矩阵,
,那么对于N中的所有
,
,我们有
。
当且仅当存在一个标量
,使得:
3. 主要结果
3.1. 稳定性分析
定义3.1.1 对任意给定常数
,不确定性变时滞切换系统(2.1)在存在加性摄动形式的非脆弱
输出反馈增益情况下,可以得到对应非脆弱
滤波器存在形式,且满足给定常数
。运用对应切换规则进行系统之间切换,进而系统(2.1)是可以镇定的。如果存在形如(2.4)式的非脆弱
滤波器以及切换规则
,对于给定性能指标
,使闭环系统响应满足如下条件:
1) 当外部扰动
时,构造对应切换规则
,使得系统(2.1)是渐近稳定的。
2) 当系统(2.1)在
时的初始状态为0时,对于所有非零的
,
有如下不等式成立:
定理3.1.1 对于任意给定的性能指标
,如果对于系统(2.1)存在对称正矩阵
、
、
和矩阵
、
以及正标量
和满足
的n个实数,使得以下LMI成立:
(3.1.1)
切换规则如下:
(3.1.2)
其中,
(3.1.3)
当切换系统(2.1)具有非脆弱
滤波器时,具有不确定时变时滞切换中立系统(2.1)是渐近稳定的,并满足给定的性能指标
。
证明:
对于系统(2.6),令
是在
上形成的与切换规则相对应的切换序列。考虑系统的内部稳定性,Lyapunov函数的形式如下:
令
:
沿系统(2.6)的时间求导数并结合条件(2.3),可以得到:
(3.1.4)
首先,考虑到系统的内部稳定性,设
,可以得到:
令
(3.1.5)
令
,有:
(3.1.6)
其中,
根据不等式(2.3),如果
,有
,那么系统(2.6)是渐近稳定的。
其次,当
时,在零初始条件下,
,则性能指标函数
如下:
如果
,则闭环系统(2.6)为鲁棒稳定的,并满足给定的
性能指标
。
要使
,只需:
(3.1.7)
如果
,那么系统(2.1)满足
性能指标
,并且系统(2.1)在非脆弱
滤波器(2.4)的作用下是渐近稳定的。
如果方程(3.1.6)成立,那么下列形式的不等式必须成立:
在没有扰动输入的情况下,
仍然可以满足要求。
注释1:切换规则(3.1.2)是通过单Lyapunov泛函和凸组合技术构造的,当时间t不断变化时,系统会被切换到
的最小值所对应的子系统,这样就被激活的子系统的轨迹始终可以保持下降的状态,可以保证系统(2.1)渐近稳定的。
注释2:在定理3.1.1中,在时滞上限和时滞导数上限同时受限的条件(2.3)下,利用单Lyapunov函数方法和凸组合技术构造了Lyapunov函数及其相应的切换规则(3.1.2)。对Lyapunov函数的导数进行缩放,消除不等式中的时变时滞项,然后引入
函数,得到满足
性能指标
的非线性矩阵不等式。我们利用引理2.1将非线性矩阵不等式转化为线性矩阵不等式(3.1.3),当
时,则系统(2.1)满足
性能指标
和系统(2.1)在非脆弱
滤波器(2.4)和切换规则(3.1.2)的作用下是渐近稳定的。
3.2. 非脆弱H∞滤波器设计
本小节主要讨论
滤波器设计问题。我们设计滤波器参数
和
,使闭环系统(2.6)渐近稳定,且满足非脆弱
性能指标
。
定理3.2.1 对于任何给定的对称正定矩阵:
矩阵
、
、
、
和标量
,
使得以下LMI成立:
(3.2.1)
其中,
式(2.4)给出的非脆弱
滤波器的参数矩阵如下:
(3.2.2)
则存在非脆弱滤波器(2.4),使得闭环系统(2.6)具有非脆弱
性能指标
,并且系统(2.1)是渐近稳定的。
证明:根据矩阵
、
、
的定义,利用引理2.2和公式(2.1)~(2.7),线性矩阵不等式(3.1.3)可以转化为:
其中,
是不等式(3.1.3)的左侧部分,
分别表示为:
(3.2.3)
其中,
(3.2.4)
(3.2.5)
令
根据方程(2.2)和(2.6),处理不确定性项
的结果分别为:
根据引理2.2,可得:
根据引理2.1,可以得到不等式(3.2.1)。
注释3:定理3.2.1中,利用引理2.2,对线性矩阵不等式(3.1.3)中的不确定项进行处理,再利用引理2.1,得到线性矩阵不等式(3.2.1)和非脆弱
滤波器参数
和
。
4. 数值实验
给定两个不确定时变时滞连续中立子系统,形成一个线性切换系统(1.1)。其中,
选取
。通过求解线性矩阵不等式(3.1.1),可以得到
。
将
,
带入到(3.2.2),得到:
切换规则为
在表1中,可以看到,当时滞导数的上限发生变化时,
性能指标
也会发生变化。可以验证,所提出的方法与时滞相关。
Figure 1. System state
and its estimation
图1. 系统状态
及其估计
Table 1. Comparison of performance index γ for change in upper bound of time delay derivative
表1. 时滞导数上限变化时的性能指标γ的比较
Figure 2. System state
and its estimation
图2. 系统状态
及其估计
Figure 3. System signal
to be estimated and its estimation value
图3. 系统待估计信号
,其估计值为
Figure 4. Filtering error
图4. 滤波误差
为便于数值模拟,我们假设初始条件为
,干扰信号为
。在上述滤波器增益下,图1表示
的轨迹,其估计值为
。图2表示
的轨迹,其估计值为
。图3表示待估计信号
的轨迹,其估计值为
。图4表示滤波误差
的轨迹图。从以上图中可以看出,所得到的滤波器能够平滑地跟踪真实的状态,并能随时间迅速趋近于零,保证了滤波误差系统在满足要求的情况下是渐近稳定的。
5. 结论
本文研究了不确定时变时滞切换中立系统的非脆弱
滤波器设计问题。首先,通过构造Lyapunov函数和利用线性矩阵不等式方法,得到了使闭环系统具有
性能指标和渐近稳定的充分条件,其中线性矩阵不等式与时滞导数的上界有关。在此基础上,对线性矩阵不等式中的不确定项进行处理,得到了非脆弱
滤波器存在的充分条件,并得到非脆弱滤波器参数。最后,通过数值实验验证了非脆弱
滤波器和切换中立系统形成的闭环系统渐近稳定并具有
性能指标
,这说明了理论结果的有效性和可行性。对于非脆弱
滤波器,本文中的定理是只对本系统有效,而对其他类型系统不适用。对于其它类型的系统的非脆弱滤波器或者其他增益形式的滤波器的设计还有待进一步研究。
参考文献