具有区间时变时滞的线性系统稳定性分析
闫冬梅1, 钟辉1, 任丽莉2, 王若琳3, 李红梅1
1.吉林大学 大数据和网络管理中心,长春 130012
2.长春师范大学 网络中心,长春 130032
3.吉林交通职业技术学院 管理工程学院,长春 130012
通信作者:钟辉(1979-),男,副研究员,博士研究生.研究方向:通信与信息系统.E-mail:zhongh@jlu.edu.cn

作者简介:闫冬梅(1978-),女,高级工程师.研究方向:网络系统管理与控制.E-mail:yandongm@jlu.edu.cn

摘要

研究了给定区间时变时滞线性系统的稳定性问题。通过构造一个具有三重积分项的新的Lyaounov-Krasovskii泛函,利用一个最新提出的基于自由矩阵的积分不等式,得到了新的保守性更小的保证系统稳定的充分条件。最后,通过2个实例验证了所提出的方法比现有方法具有更小的保守性。

关键词: 自动控制技术; 时滞系统; 时滞依赖稳定性; 区间时滞; 线性矩阵不等式
中图分类号:TP273 文献标志码:A 文章编号:1671-5497(2018)05-1556-07
Stability analysis of linear systems with interval time-varying delay
YAN Dong-mei1, ZHONG Hui1, REN Li-li2, WANG Ruo-lin3, LI Hong-mei1
1.Management Center of Big Data and Network, Jilin University, Changchun 130012,China
2.Network Center, Changchun Normal University, Changchun 130032,China
3.College of Administrative Engineering,Jilin Communications Polytechnic, Changchun 130012,China
Abstract

This paper investigates the stability problem of a linear system with time-varying delay. The time-varying delay is assumed to belong to a given interval. A new Lyapunov-Krasovskii function is constructed which contains some tripe integral terms. Using some integral inequalities including the recent proposed free-matrix-based integral inequality, new sufficient stability conditions are obtained, which are presented in terms of Linear Matrix Inequality (LMI). Finally, two examples are given to demonstrate that the obtained results are less conservative than some existing method.

Keyword: automatic control technology; time-delay system; delay-dependent stability; interval delay; linear matrix inequality
0 引 言

时滞现象存在于许多实际系统中, 例如化工控制过程, 网络化控制系统[1, 2]。时滞是引起系统不稳定和使系统性能恶化的主要原因之一[3]。近年来, 由于其在理论和实际应用中的重要性, 具有区间时变时滞的线性系统的稳定性分析问题得到了越来越多的关注。一般来说, 研究时滞系统稳定性分析有两类方法。一类是频域法, 另一类是时域法。频域法通常只适用于定常时滞系统, 而在许多实际系统中时滞通常是时变的。例如, 在网络化控制系统中, 从传感器到控制器或者从控制器到执行器的时滞往往是时变的。因此, 时域法被广泛用于研究时变时滞系统的稳定性条件。通过建立一个Lyapunov-Krasovskii泛函或者一个Lyapunov-Razumikhin函数, 可以得到一些更容易验证的充分条件。

现有文献中提出了一些典型的基于Lyapunov泛函的稳定性分析方法。文献[4, 5]提出了广义系统模型变换的方法, 并且证明了广义系统模型变换方法比其他3种模型变换方法更有效。文献[6, 7, 8, 9, 10]提出了自由权矩阵的方法。这种方法比广义系统模型变换方法具有较小的保守性, 因为它不涉及任何模型变换和交叉项的缩放技术。文献[11, 12]提出了一个包含三重积分项的Lyapunov泛函。这个Lyapunov泛函恰当地利用了时滞的下界信息, 因而可以得到具有更小保守性的结果。文献[13]提出了Jesen不等式并将其应用于时滞系统的稳定性分析。与广义系统模型变换方法和自由权矩阵方法相比, Jesen不等式方法涉及的矩阵变量较少, 因此运算效率高。文献[14]提出了一个基于Wirtinger积分不等式的方法, 用该不等式可以得到比运用Jesen不等式保守性更小的稳定性条件。文献[15]提出了凸组合技术。文献[16]提出了倒凸组合方法用来估计Lyaounov泛函的导数, 并且得到了具有更小保守性的结果。除了上面提到的方法, 还有许多其他的方法, 例如文献[17]中的离散Lyapunov泛函方法, 文献[18, 19, 20]中的时滞分割方法等。这些方法的详情可参阅文献[21, 22, 23, 24, 25]。

在现有的文献中, 针对如何减少稳定性条件的保守性主要有两类方法。一个是构造适当的泛函, 另一个是找到新的方法来尽可能紧地估计Lyapunov泛函的导数。本文通过构造一个新的Lyapunov泛函, 并且利用最新提出的基于自由权矩阵的积分不等式方法来估计Lyapunov泛函的导数, 得到了保守性更小的稳定性条件。最后, 通过实例验证了本文提出方法的有效性。

1 理论推导

考虑如下具有区间时变时滞的线性系统:

x·(t)=Ax(t)+A1x(t-d(t)), t0x(θ)=φ(θ), θ[-h2, 0](1)

式中: x(t)Rn为状态向量; ARn×nA1Rn×n均为系统矩阵; 初始条件 φ(θ)为一个连续可微的函数; d(t)为可微函数且满足:

0< h1d(t)h2μ1d·(t)μ2(2)

式中: 0< h1h2μ1μ2为已知常数。

以下研究在时滞满足式(2)的情况下, 针对系统提出具有更小保守性的稳定性条件。

引理1[22] x是一个可微的向量函数。对于正定对称矩阵 RRn×n> 0, X1R3n×3n> 0, X3R3n×3n> 0, 以及任意适当维数的矩阵 X2R3n×3n, N1R3n×n, N2R3n×n, 使得:

X1X2N1* X3N2* * R0

那么下面的积分不等式成立:

-t-bt-ax·t(s)Rx·(s)dsφ(t)Ξφ(t)(3)

式中:

φ(t)=x(t-a)x(t-b)1b-at-bt-ax(s)ds,

Ξ=(b-a)(X1+13Z3)+N1Π1+ΠT1NT1+N2Π2+ΠT2NT2,

Π1=I-I0, Π2=-I-I2I

引理2[11, 13] 对于一个矩阵 ZRn×n> 0, 和参数 τ2> τ1> 0, 下面的两个积分不等式成立:

-τ12t-τ2t-τ1xT(s)Zx(s)ds-t-τ2t-τ1xT(s)dsZt-τ2t-τ1x(s)ds, -τ1222-τ2-τ1t+θt-τ1xT(s)Zx(s)dsdθ--τ2-τ1t+θt-τ1xT(s)dsdθZ-τ2-τ1t+θt-τ1x(s)dsdθ

式中: τ12=τ2-τ1

下面的定理1给出了保证系统(1)(2)渐近稳定的充分条件。

定理1 对于给定的标量 h1h2μ1μ2, 如果存在矩阵 PR7n×7n> 0; QiR2n×2n> 0; SiR2n×2n> 0; RiRn×n> 0; i=1, 2; ZjRn×n> 0; j=1, 2, 3, 4; 对称矩阵 X1, X3, G1, G3, L1, L3R3n×3n, 和任意矩阵 X2, G2, L2R3n×3n; N1, N2, M1, M2, W1, W2R3n×n, 使得下面的线性矩阵不等式对于 d(t)={h1, h2}d·(t)={μ1, μ2}成立, 则系统 (1)(2)渐近稳定。

Ψ(d(t), d·(t))< 0(4)X1X2N1* X3N2* * Z10(5)G1G2M1* G3M2* * Z20(6)L1L2W1* L3W2* * Z20(7)

式中:

Ψ(d(t), d·(t))=H1PH2+HT2PHT1+H3(S1- Q1+Q2)HT3+H4(1-d·(t))(S2-S1)HT4+ H5(-S2-Q2)HT5+H6Q1HT6+ATbΘAb+ e6ΛeT6+h1e1Z3eT1+h12e3Z4eT3+H7Ω1HT7+ H8Ω2HT8+H9Ω3HT9-h1e8Z3eT8- αe9Z4eT9-βe10Z4eT10-2(e1-e8)R1(e1- e8)T-2(e2-e10)R2(e2-e10)T-2(e3- e9)R2(e3-e9)T,

H1=e1e2e3e4αe9βe10h1e8,

H2=Ab(1-d·(t))eT5eT6eT7eT3-(1-d·(t))eT2(1-d·(t))eT2-eT4eT1-eT3, H3=e3e6, H4=e2e5, H5=e4e7,

H6=e1ATb, H7=e1e3e8, H8=e3e2e9,

H9=e2e4e10, Θ=h1Z1+12h12R1, Λ=h12Z2+12h122R2,

Ab=AA100000000,

Ω1=h1(X1+13X3)+N1Π1+ΠT1NT1+N2Π2+ΠT2NT2,

Ω2=α(G1+13G3)+M1Π1+ΠT1MT1+M2Π2+ΠT2MT2,

Ω3=β(L1+13L3)+W1Π1+ΠT1WT1+W2Π2+ΠT2WT2,

α=d(t)-h1, β=h2-d(t),

eTi=00i-1I0010-i, (i=1, 2, , 10)

证明:选择如下Lyapunov泛函:

V(xt)=ζT(t)(t)+t-d(t)t-h1ηT(s)S1η(s)ds+t-h2t-d(t)ηT(s)S2η(s)ds+t-h1tηT(s)Q1η(s)ds+t-h2t-h1ηT(s)Q2η(s)ds+-h10t+θtx·T(s)Z1x·(s)dsdθ+-h2-h1t+θt-h1x·T(s)Z2x·(s)dsdθ+-h10t+θtxT(s)Z3x(s)dsdθ+-h2-h1t+θt-h1xT(s)Z4x(s)dsdθ+-h10θ0t+λtx·T(s)R1x·(s)dsdλdθ+-h2-h1θ-h1t+λt-h1x·T(s)R2x·(s)dsdλdθ(8)

沿系统(1)的轨迹对Lyapunov泛函(8)求导, 可得:

V·(xt)=2ζT(t)Pζ·(t)+ηT(t-h1)S1η(t-h1)- (1-d·(t))ηT(t-d(t))S1η(t-d(t))+ (1-d·(t))ηT(t-d(t))S2η(t-d(t))- ηT(t-h2)S2η(t-h2)+ ηT(t)Q1η(t)-ηT(t-h1)Q1η(t-h1)+ ηT(t-h1)Q2η(t-h1)-ηT(t-h2)Q2η(t- h2)+h1x·T(t)Z1x·(t)-t-h1tx·T(s)Z1x·(s)ds+ h12x·T(t-h1)Z2x·(t-h1)- t-h2t-h1x·T(s)Z2x·(s)ds+ h1xT(t)Z3x(t)-t-h1txT(s)Z3x(s)ds+ h12xT(t-h1)Z4x(t-h1)- t-h2t-h1xT(s)Z4x(s)ds+ h122x·T(t)R1x·(t)--h10t+θtx·T(s)R1x·(s)dsdθ+ h1222x·T(t-h1)R2x·(t-h1)- -h2-h1t+θt-h1x·T(s)R2x·(s)dsdθ(9)

式中:

h12=h2-h1,

ζ(t)=col  x(t), x(t-d(t)), x(t-h1), x(t-h2), t-d(t)t-h1x(s)ds, t-h2t-d(t)x(s)ds, t-h1tx(s)ds,

η(t)=colx(t), x·(t)

利用引理1, 可得:

-t-h1tx·T(s)Z1x·(s)dsξT(t)H7Ω1HT7ξ(t)(10)-t-h2t-h1x·T(s)Z2x·(s)ds=-t-h2t-d(t)x·T(s)Z2x·(s)ds-t-d(t)t-h1x·T(s)Z2x·(s)dsξT(t)(H8Ω2HT8+H9Ω3HT9)ξ(t)(11)

式中:

ξ(t)=col  x(t), x(t-d(t)), x(t-h1),  x(t-h2), x·(t-d(t)), x·(t-h1),  x·(t-h2), 1ht-h1tx(s)ds,  1d(t)-h1t-d(t)t-h1x(s)ds, 1h2-d(t)t-h2t-d(t)x(s)ds

利用引理2, 可得:

-t-h1txT(s)Z3x(s)ds-h1ξT(t)e8Z3eT8ξ(t)(12)--h10t+θtx·T(s)R1x·(s)dsdθ-2ξT(t)(e1-e8)R1(eT1-eT8)ξ(t)(13)显见:-t-h2t-h1xT(s)Z4x(s)ds=-t-h2t-d(t)xT(s)Z4x(s)ds-t-d(t)t-h1xT(s)Z4x(s)ds--h2-h1t+θt-h1x·T(s)R2x·(s)dsdθ--h2-d(t)t+θt-d(t)x·T(s)R2x·(s)dsdθ--d(t)-h1t+θt-h1x·T(s)R2x·(s)dsdθ(14)

利用引理2, 则有:

--d(t)-h1xT(s)Z4x(s)ds-αξT(t)e9Z4eT9ξ(t)(15)--h2-d(t)xT(s)Z4x(s)ds-βξT(t)e10Z4eT10ξ(t)(16)--d(t)-h1t+θt-h1x·T(s)R2x·(s)dsdθ-2ξT(t)(e3-e9)R2(e3-e9)Tξ(t)(17)--h2-d(t)t+θt-d(t)x·T(s)R2x·(s)dsdθ-2ξT(t)(e2-e10)R2(e2-e10)Tξ(t)(18)

由式(9)~(18)可以得到:

V·(xt)ξT(t)Ψ(d(t), d·(t))ξ(t)(19)

由于 Ψ(d(t), d·(t))d·(t)是关于 d(t)d·(t)的线性函数, 可以得到如果

Ψ|d(t)=h1, d·(t)=μ1< 0, Ψ|d(t)=h1, d·(t)=μ2< 0, Ψ|d(t)=h2, d·(t)=μ1< 0, Ψ|d(t)=h2, d·(t)=μ2< 0,

那么 Ψ(d(t), d·(t))< 0, V·(xt)< 0因此, 系统是渐近稳定的。证毕。

需要注意的是, 本文构造了一个包含三重积分的新的Lyapunov泛函。如文献[11, 12]所指出, 这个Lyapunov泛函对于减小结果的保守性是相当有效的。另外, 本文应用文献中最新提出的基于自由矩阵的积分不等式估计Lyapunov泛函的导数, 故能更进一步地减少所得稳定性条件的保守性。正是由于利用了这些新技术, 本文得到的稳定性条件比现有文献中的结果具有更小的保守性。

当时滞变化率未知时, 基于定理1可以得到下面的结论。

定理2 对于给定的标量 h1h2, 如果存在矩阵 PR4n×4n> 0; QiR2n×2n> 0; RiRn×n; i=1, 2; ZjRn×n> 0; j=1, 2, 3, 4; 对称矩阵 X1X3G1G3L1L3R3n×3n, 和任意适当维数的矩阵 X2G2L2R3n×3n; N1N2M1M2W1W2R3n×n, 使得下面的线性矩阵不等式对 d(t)={h1, h2}成立, 则系统(1)(2)渐近稳定。

Ψ^(d(t))< 0(20)X1X2N1* X3N2* * Z10(21)G1G2M1* G3M2* * Z20(22)L1L2W1* L3W2* * Z20(23)

式中:

Ψ^(d(t))=H^1PH^2+H^T2PH^T1+H3(S1-Q1+ Q2)HT3-H5Q2HT5+H6Q1HT6+ATbΘAb+ e6ΛeT6+h1e1Z3eT1+h12e3Z4eT3+H7Ω1HT7+ H8Ω2HT8+H9Ω3HT9-h1e8Z3eT8-αe9Z4eT9- βe10Z4eT10-2(e1-e8)R1(e1-e8)T-2(e2- e10)R2(e2-e10)T-2(e3-e9)R2(e3-e9)T,

H1=[e1e3e4h1e8], H2=[AbeT6eT7eT1-eT3]T

证明:构造如下Lyapunov泛函:

V(xt)=ζ^T(t)Pζ^(t)+t-h1tηT(s)Q1η(s)ds+t-h2t-h1ηT(s)Q2η(s)ds+-h10t+θtx·T(s)Z1x·(s)dsdθ+-h2-h1t+θt-h1x·T(s)Z2x·(s)dsdθ+-h10t+θtxT(s)Z3x(s)dsdθ+-h2-h1t+θt-h1xT(s)Z4x(s)dsdθ+-h10θ0t+λtx·T(s)R1x·(s)dsdλdθ+-h2-h1θ-h1t+λt-h1x·T(s)R2x·(s)dsdλdθ(24)

式中: ζ^(t)=col{x(t), x(t-h1), x(t-h2), t-h1tx(s)ds}, η(t)=col{x(t), x·(t)}证明过程与定理1类似, 故从略。证毕。

2 实例验证

本文将通过2个数值仿真示例来验证提出方法的有效性。

例1 考虑具有如下参数的时滞系统:

A=-200-0.9, A1=-10-1-1

对于变量 μ2=-μ1=μh1=1, 对于给定的时滞下界, 表1列举了使用新方法得到的最大时滞上界, 并与文献[11, 12, 14, 22]中的结果进行了对比。可见, 通过新方法得到的最大时滞上界明显大于文献[11, 12, 14, 22]中的结果, 从而验证了本文提出方法的有效性。

表1 不同μ 得到的最大时滞上界 Table 1 MUDBs for different μ

如果时滞变化率不在一个对称区间内, 对于 h1=1以及不同的时滞变化率区间的情况, 表2对比了本文得到的最大时滞上界以及文献[14, 22]的结果。在表2中, [a, b]是指 ad·(t)b表2可以看出, 使用新方法可以得到较大的最大时滞上界, 从而说明本文提出的方法具有较小的保守性。

表2 不同时滞变化率区间的最大时滞上界 Table 2 MUBDs for asymmetric delay-rate interval

例2 考虑具有如下参数的时滞系统:

A= 0 1-1-2, A1= 00-11

对于给定的 h1=0.5, μ2=0.3和不同的 μ1, 表3对比了利用本文提出方法得到的最大时滞上界以及文献[11, 12, 14]中的结果。从表3可以看出, 本文提出的方法与文献[11, 12, 14]相比具有更小的保守性。

表3 不同μ 1得到的最大时滞上界 Table 3 MUBDs for defferent μ 1

例3 考虑如下的直流电机系统 :

x·(t)=0-10-0.1x(t)+00.1u(t)

通过网络控制上面的直流电机, 系统的控制器设计为:

u(t)=-3.75-11.75x(t)

考虑到网络中存在的网络延时, 上述系统可以建模为一个时滞系统:

x·(t)=0-10-0.1x(t)+0 0.1T-3.75-11.5x(t-d(t))

假设网络延时有界, 假设网络延时下界 h1=0.1, 表4对比了利用本文提出方法得到的最大网络延时上界以及文献[11, 12, 14]中的结果。从表4可以看出, 本文提出的方法与文献[11, 12, 14]相比具有更小的保守性。

表4 不同h1得到的最大网络延时上界 Table 4 MUBNDs for different h1

假设系统的初始状态为 x(0)=0 1T, 网络延时满足 0.1< d(t)< 1.1, 如图1所示, 系统的状态响应曲线如图2所示。从图2可知, 系统渐近稳定, 从而证明了本文所提出方法的有效性。

图1 网络延时Fig.1 Network delay

图2 系统状态响应曲线Fig.2 System state response curve

3 结束语

研究了具有区间时变时滞的线性系统的稳定性问题。通过构造了一个新的增广Lyapunov泛函, 并利用最新提出的基于自由矩阵的积分不等式估计Lyapunov泛函的导数, 得到了一些保守性更小的稳定性条件。最后, 通过2个仿真示例验证了所提出方法有效。

The authors have declared that no competing interests exist.

参考文献
[1] Sun J, Liu G P, Chen J, et al. Networked predictive control for Hammerstein systems[J]. Asian Journal of Control, 2011, 13: 265-272. [本文引用:1]
[2] Sun J, Chen J. Networked predictive control for systems with unknown or partially known delay[J]. IET Control Theory Appl, 2014, 8(18): 2282-2288. [本文引用:1]
[3] 李学军, 陈虹, 于树友. 基于时滞系统的无偏H滤波[J]. 吉林大学学报: 工学版, 2009, 39(2): 473-479.
Li Xue-jun, Chen Hong, Yu Shu-you. Unbiased Hfiltering for time delay system[J]. Journal of Jilin University(Engineering and Technology Edition), 2009, 39(2): 473-479. [本文引用:1]
[4] Fridman E, Shaked U. Delay-dependent stability and H control: constant and time-varying delays[J]. International Journal of Control, 2003, 76: 48-60. [本文引用:1]
[5] Fridman E, Shaked U. An improved stabilization method for linear systems with time-delay[J]. IEEE Transactions on Automatic Control, 2002, 47: 1931-1937. [本文引用:1]
[6] He Y, Wang Q G, Lin C, et al. Augmented Lyapunov functional and delay-dependent stability criteria for neutral systems[J]. International Journal of Robust and Nonlinear Control, 2005, 15: 923-933. [本文引用:1]
[7] He Y, Wang Q G, Xie L, et al. Delay-range-dependent stability for systems with time-varying delay[J]. Automatica, 2007, 43: 371-376. [本文引用:1]
[8] He Y, Wang Q G, Xie L, et al. Further improvement of free-weighting matrices technique for systems with time-varying delay[J]. IEEE Transactions on Automatic Control, 2007, 52: 293-299. [本文引用:1]
[9] He Y, Wu M, She J H, et al. Delay-dependent robust stability criteria for uncertain neutral systems with mixed delays[J]. Systems & Control Letters, 2004, 51: 57-65. [本文引用:1]
[10] He Y, Wu M, She J H, et al. New delay-dependent stability criteria for robust stability of time-varying delay systems[J]. Automatica, 2004, 40: 1435-1439. [本文引用:1]
[11] Sun J, Liu G P, Chen J, et al. Improved delay-range-dependent stability criteria for linear systems with time-varying delays[J]. Automatica, 2010, 46: 466-470. [本文引用:1]
[12] Sun J, Han Q L, Chen J, et al. Less conservative stability criteria for linear systems with interval time-varying delays[J]. International Journal of Robust and Nonlinear Control, 2015, 25: 475-485. [本文引用:1]
[13] Gu K. Integral inequality in the stability problem of time-delay systems[C]∥Proceedings of the 39th IEEE Conference on Decision and Control, Sydney, Australia, 2000, 3: 2805-2810. [本文引用:2]
[14] Seuret A, Gouaisbaut F. Wirtinger-based integral inequality: application to time-delay systems[J]. Automatica, 2013, 49: 2860-2866. [本文引用:1]
[15] Park P, Ko J W. Stability and robust stability for systems with a time-varying delay[J]. Automatica, 2007, 43: 1855-1858. [本文引用:1]
[16] Park P, Ko J W, Jeong C. Reciprocally convex approach to stability of systems with time-varying delays[J]. Automatica, 2011, 47: 235-238. [本文引用:1]
[17] Gu K, Han Q L, Luo A C J, et al. Discretized lyapunov functional for systems with distributed delay and piecewise constant coefficients[J]. International Journal of Control, 2001, 74: 737-744. [本文引用:1]
[18] Gouaisbaut F, Peaucelle D. Delay-dependent robust stability of time delay systems[C]∥Proceedings of IFAC Symposium on Robust Control Design, Toulouse, France, 2006, 5: 453-458. [本文引用:1]
[19] Han Q L. New stability criteria for linear systems with interval time-varying delay[J]. Automatica, 2008, 44: 2680-2685. [本文引用:1]
[20] Hui Jun-jun, Zhang He-xin, Kong Xiang-yu, et al. On improved delay-dependent robust stability criteria for uncertain systems with interval time-varying delay[J]. International Journal of Automation and Computing, 2015, 12(1): 102-108. [本文引用:1]
[21] Fridman E. Tutorial on Lyapunov-based methods for time-delay systems[J]. European Journal of Control, 2014, 20(6): 271-283. [本文引用:1]
[22] Zeng H B, He Y, Wu M, et al. Free-matrix-based integral inequality for stability analysis of systems with time-varying delay[J]. IEEE Transactions on Automatic Control, 2015, 60(10): 2768-2772. [本文引用:2]
[23] Sun J, Chen J. A survey on Lyapunov-based methods for stability of linear time-delay systems[J]. Frontiers of Computer Science, 2017, 11(4): 555-567. [本文引用:1]
[24] Liu K, Fridman E, Johansson K H, et al. Generalized Jensen inequalities with application to stability analysis of systems with distributed delays over infinite time-horizons[J]. Automatica, 2016, 69: 222-231. [本文引用:1]
[25] Kim J H. Further improvement of Jensen inequality and application to stability of time-delayed systems[J]. Automatica, 2016, 64: 121-125. [本文引用:1]