作者简介:胡萌萌(1990—),女,深圳大学硕士研究生.研究方向:动力系统在复杂网络中的应用.E-mail: 13728814178@163.com
中文责编:方 圆; 英文责编:淡 紫
深圳大学数学与统计学院,广东深圳518060
HU Mengmeng, FENG Jianwen, ZHAO Yi, and XU ChenCollege of Mathematics and Statistics, Shenzhen University, Shenzhen 518060, Guangdong Province, P.R.China
Lur'e network; quasi-synchronization; parameter mismatch; nonlinear coupling; time-varying delay; aperiodically intermittent control; pinning control
DOI: 10.3724/SP.J.1249.2018.01085
带有参数不匹配的异构网络是不能达到完全同步的,但可以实现网络中各节点的状态和目标状态间存在有界误差的拟同步.基于非周期间歇牵制控制策略,讨论一类非线性时变时滞耦合异构Lur'e网络的拟同步.通过构造简单的Lyapunov函数,经过严格理论分析,得到若干与参数相关的充分条件,保证了该异构网络实现拟同步.数值模拟验证了所得理论结果的正确性和有效性.
The heterogeneous networks with mismatched parameters cannot reach complete synchronization,but they can achieve quasi-synchronization,which means the state of each node in networks will keep a bounded error comparing with the target state.Quasi-synchronization is investigated for a class of nonlinear time-varying delay coupled heterogeneous Lur'e networks based on the aperiodically intermittent pinning control technique.By constructing a simple Lyapunov function,some sufficient conditions regarding the parameters are obtained to guarantee the heterogeneous networks to achieve quasi-synchronization through rigorous theoretical analysis.The validity and effectiveness of the theoretical results are verified by numerical simulation.
近二十年来,复杂网络同步控制由于其理论挑战性和应用广泛性受到了不同学科研究者的极大关注,为科学研究中的一个热点问题.研究讨论了各种不同类型网络的多种同步模式,如完全同步[1]、簇同步[2]、指数同步[3]及有限时间同步[4-5]等.
在现有关于复杂网络同步问题的研究中,大部分涉及的是恒同网络,即网络中的所有节点都有相同的动力学行为.然而,由非恒同节点构成的网络,在自然界和工程领域普遍存在.例如,因特网中的节点一般是不同的; 社交网络中的每个个体也存在差异,这就可能产生网络中的参数不匹配现象.一般来说,由参数不匹配的节点所构成的异构网络不能达到完全同步,拟同步则为异构网络的同步研究提供了全新思路.拟同步意味着各个节点的状态和目标节点状态之间存在一个有界误差.随着对网络同步问题的深入研究,异构网络的拟同步问题引起了越来越多研究关注,LIU等[6]通过反馈控制方法研究存在参数不匹配和随机扰动的异构复杂网络的拟同步; HE等[7]利用一个广义的Halanay不等式和矩阵测度得到了带有参数不匹配的时滞耦合系统的延迟拟同步判据; HE等[8]还基于分布式脉冲控制讨论了带有参数不匹配的异构复杂网络的拟同步,并对误差进行估计优化.WANG等[9]研究了由非恒同节点构成的非线性耦合网络的簇同步问题.以上研究忽略了节点间信息交换时的时滞,而时滞在现实中是不可避免的.另外,Lur'e系统是一类非常典型的非线性系统,如常见的蔡氏电路和Cohen-Grossberg神经网络均具有Lur'e系统的形式.近期出现了一些对Lur'e网络系统的同步研究,如SHI等[10]研究了带有数据采样反馈控制的混沌Lur'e系统的主从同步; TANG等[11]分析了由非连续Lur'e系统构成的非线性耦合复杂网络的有限时间簇同步.但对于有非线性时变时滞耦合的异构Lur'e网络的拟同步问题尚未涉及.
另一方面,控制方式在实现网络同步中起着重要作用,各种控制策略已成功用于网络同步研究中.按控制时间可分为连续状态反馈控制[12]和非连续状态反馈控制,后者包括脉冲控制[13]、样本采样控制[14]及事件激发控制[15]等.非连续牵制控制由于其低成本高效率的特点,在网络同步研究中受到青睐.特别是周期间歇牵制控制策略受到了关注,如LI等[16]基于周期间歇牵制控制策略实现了分数阶复杂动力网络的同步.然而,周期间歇牵制控制策略在现实中的应用并不合理,例如智能电网中风能发电呈现出非周期.因此,将非周期间歇牵制控制策略用于网络同步研究中具有现实意义.LIU等[17]成功地将非周期间歇牵制控制用于一类复杂网络的同步问题中; 2016年LIU等[18]进一步基于非周期间歇牵制控制方法研究了带有参数不匹配的复杂网络的拟同步问题.
受到以上相关研究的启发,本研究将考虑有非线性时变时滞耦合,且带有参数不匹配异构Lur'e网络的拟同步问题.基于非周期间歇牵制控制策略,通过构造合适的Lyapunov函数,经过严格理论分析得到实现该类网络拟同步的若干充分条件.并给出相应数值模拟论证所得理论结果有效.
首先,考虑一般Lur'e系统
s·(t)=As(t)+Bf(Cs(t))(1)
其中, s(t)=(s1(t),s2(t),…,sn(t))T∈Rn是t时刻的状态向量; A∈Rn×n、 B∈Rn×m 及C∈Rm×n是常数矩阵; f(Cs(t))=(f1(cT1s(t)),…,fm(cTms(t)))T∈Rm为系统的非线性向量值函数,满足 f(0)=0, 且 cTl表示矩阵 C的第l行, f(·)是非线性函数, l=1,…,m.
本文将在非周期间歇牵制控制下,研究有时变时滞的非线性耦合异构Lur'e网络的拟同步.不失一般性,只对网络中的第1个节点实施非周期间歇控制,这样受控制的网络模型可以描述为
{Zx·1(t)=A1x1(t)+B1f(Cx1(t))+c1∑Nj=1l1jΓxj(t)+c2∑Nj=1g1jΓh(xj(t-τ(t)))-c1dΓ(x1(t)-s(t))
x·i(t)=Aixi(t)+Bi f(Cxi(t))+c1∑Nj=1lijΓxj(t)+c2∑Nj=1gijΓh(xj(t-τ(t))),
i=2,…,N, t∈[tk,sk], k=0,1,2,…(2)
x·i(t)=Aixi(t)+Bi f(Cxi(t))+c1∑Nj=1lijΓxj(t)+c2∑Nj=1gijΓh(xj(t-τ(t))),
i=1,2,…,N, t∈(sk,tk+1), k=0,1,2,…
其中, xi(t)=(xi1(t),xi2(t),…,xin(t))T∈Rn是第i个节点的状态变量; 常数c1和c2>0表示耦合强度; Ai∈Rn×n, Bi∈Rn×m均是常数系数矩阵; 控制增益d>0; Γ=diag{ζ1,ζ2,…,ζ2}∈Rn×n是正定内部耦合矩阵; τ(t)是内部耦合时变时滞; L=(lij)∈RN×N和 G=(gij)∈RN×N分别为外部线性耦合矩阵和外部非线性耦合矩阵,反映网络拓扑结构,分别定义如下:若节点i和节点j之间有联系,则lij=lji>0, gij=gji>0, 否则lij=0, gij=0. 定义对角元素lii=-∑Nj=1, j≠i lij, gii=-∑Nj=1, j≠i gij.非线性向量值函数h(xi(t-τ(t))=(h1(xi(t-τ(t))),…,hn(xi(t-τ(t))))T∈Rn. 非周期间歇牵制控制策略是指将时间任意分为若干区间[tk,tk+1),k=0,1,2,…, 在每一区间内的部分时间施加控制, 即在第k个区间[tk,tk+1)内, [tk,sk]是第k个工作区间(其长度为sk-tk)和(sk, tk+1)是第k个休息区间(其长度为tk+1-sk).
系统(2)的初始条件为xi(t)=φi(t)∈C([-τ,0],Rn), i=1,2,…,N, 其中, C([-τ,0],Rn)表示从区间[-τ,0]映射到 Rn的连续函数集, =φi(κ)=2=supκ∈[-τ,0](∑nz=1φiz(κ)2)1/2, i=1,2,…, N为其范数的定义.假设存在一个正常数ε^-, 对任意时间t, 均有=s(t)=≤ε^-, 这表明 s(t)有界.
以下给出本文用到的定义、引理及假设.
定义1[19] 设Ω为相空间中包含系统(1)的混沌吸引子的相关区域,若存在常数T≥t0和 >0, 使得对任意初值 φi∈Ω, s(t0)∈Ω, 当t≥T时=ei(t)=≤(i=1,2,…,N)成立,则称系统(2)和系统(1)间实现了拟同步,其中,误差项 ei(t)=xi(t)-s(t).
假设1[17] 存在两个正数0<θ<ω,那么,对于k=0,1,2,…, 有
{infk(sk-tk)=θ>0,
supk(tk+1-tk)=ω<+∞.
注1 假设1表明每个控制区间的长度不会小于θ, 而休息区间的长度不会大于(ω-θ).
假设2 对于Lur'e网络中的非线性向量值函数满足f(·)以下扇形条件:
0≤(fl(a2)-fl(a1))/(a2-a1)≤δl, a1, a2∈R, a1≠a2, δl>0, l=1,2,…,m.
假设3 对于非线性向量值函数 h(·), 存在非负常数qij(i, j=1,2…,n), 对于任意 z1和 z2∈Rn, 有|hi(z1)-hi(z2)|≤∑nj=1qij|z1j-z2j|.
引理1[20] 若定义在[t0-τ,+∞)上的非负连续函数y(t)满足y ·(t)<-ay(t)+byτ(t)+ε, 其中, a>b>0, yτ(t)=sup-τ≤κ≤0(y(t+κ)), 则当t≥t0时有
y(t)≤yτ(t0)e<sup>-r(t-t0)+ε/r(3)
此处, r>0是a=berτ+r的唯一解.
引理2[20] 若定义在[t0-τ,+∞)上的非负连续函数y(t)满足y ·(t)<ay(t)+byτ(t)+ε, 其中, b>0, a+b>0, 那么当t≥t0时有
y(t)≤(yτ(t0)+ε/(a+b))e<sup>(a+b)(t-t0)-ε/(a+b)(4)
引理3[20] 若定义在[t0-τ,+∞)上的非负连续函数y(t), 满足对于k=0,1,2,…, 有
y(t)≤yτ(tk)e<sup>-r(t-tk)+α, tk≤t≤sk,
y(t)≤(yτ(sk)+β)e<sup>ξ(t-sk)-β, sk<t<tk+1.
若假设1成立,记 ρ=(ρ1-ρ2)>0,其中, ρ1=r(θ-τ), ρ2=ξ(ω-θ), 则有
y(t)≤yτ(0)eρe<sup>-(ρ/ω)t+v/(1-e-ρ)+α(5)
对于t≥0时成立,其中, v=(α+β)e<sup>ρ2-β.
为叙述方便,用ΔAi=Ai-A, ΔBi=Bi-B表示参数不匹配矩阵,此外还引入下述符号:
ηi(t)=f(Cei(t))=f(Cxi(t))-f(Cs(t)),
ηi(t)=(ηi1,ηi2,…,ηim)T,
其中,ηi1=f1(cTlxi(t))-fl(cTls(t)),
l=1,2,…,m,
η(t)=(ηT1(t),ηT2(t),…,ηTN(t))T,
ΔA=(ΔAT1,ΔAT2,…,ΔATN)T,
ΔB=(ΔBT1,ΔBT2,…,ΔBTN)T,
e(t)=(eT1(t),eT2(t),…,eTN(t))T,
M(t)=(mT1(t),mT2(t),…,mTN(t))T, 这里
mi(t)=ΔAis(t)+ΔBi f(Cs(t)), i=1,2,…,N,
h(ei(t-τ(t)))=h(xi(t-τ(t)))-
h(s(t-τ(t))),
h(ei(t-τ(t)))=(h1(ei(t-τ(t))),…,
hn(ei(t-τ(t))))T,
H(e(t-τ(t)))=(hT(e1(t-τ(t))),…,
hT(eN(t-τ(t))))T,
A ^=diag(A1,A2,…,AN),B ^=diag(B1,B2,…,BN),
D=diag(d,0,…,0),L ^=L-D,
Δ=diag(δ1, δ2,…,δm),
T=diag(τ1, τ2,…,τm).
则由系统(1)和(2)可得相应的误差系统为
{Ze·1(t)=A1e1(t)+B1f(Ce1(t))+c1∑Nj=1l1jΓej(t)+c2∑Nj=1g1jΓh(ej(t-τ(t)))-c1dΓe1(t)+m1(t)
e·i(t)=Aiei(t)+Bi f(Cei(t))+c1∑Nj=1lijΓej(t)+c2∑Nj=1gijΓh(ej(t-τ(t)))+mi(t),
i=2,3,…,N, t∈[tk,sk], k=0,1,2,…
e·i(t)=Aiei(t)+Bi f(Cei(t))+c1∑Nj=1lijΓej(t)+c2∑Nj=1gijΓh(ej(t-τ(t)))+mi(t),
i=1,2,…,N, t∈(sk,tk+1), k=0,1,2,…
则上述误差系统又可写为
{Ze·(t)=A ^e(t)+B ^η(t))+c2(GΓ)H(e(t-τ(t)))+c1(L ^Γ)e(t)+M(t),
t∈[tk,sk], k=0,1,2,…
e·(t)=A ^e(t)+B ^η(t))+c2(GΓ)H(e(t-τ(t)))+c1(LΓ)e(t)+M(t),
t∈(sk,tk+1), k=0,1,2,…(6)
定理1 若假设1至假设3成立,且存在一个正常数ψ使得max1≤i≤N(=ΔAi=+=ΔBi==ΔC=)≤ψ, 此外,还存在正数γ1, γ2, μ3>μ2>μ1, α1和矩阵 Q=(qij)n×n满足下列条件:
ⅰ. Π=
[A ^+A ^T-2α1(INΓ)B ^+(INCTΔT)
B ^T+(INTΔC)-2(INT)]≤0;
ⅱ.(γ2+μ1)I+2c1(L ^Γ)+c2γ1(GGTΓΓT)+2α1(INΓ)≤0;
ⅲ.(c2)/(γ1)λmax(QTQ)-μ2≤0;
ⅳ.(γ2-(μ3-μ1))I+2c1(LΓ)+c2γ1(GGTΓΓT)+2α1(INΓ)≤0;
ⅴ. ρ=ρ1-ρ2=r(θ-τ)-(μ3-μ1+μ2)(ω-θ)>0, r是方程r-μ1+μ2eγτ=0的唯一正解.
则存在一个正数T, 使得当t≥T时, =e(t)=≤ε2+(ε1)1/2, 这里ε2可以是任意小正数,即系统(1)和(2)达到拟同步,其中,ε1=v/(1-e-ρ)+α,v=(α+β)e<sup>(μ3-μ1+μ2)(ω-θ)-β, α=(Nε^-2ψ2)/(rγ2), β=(Nε^-2ψ2)/(γ2(μ3-μ1+μ2))
【证】定义以下Lyapunov函数
V(t)=∑Ni=1eTi(t)ei(t)=eT(t)e(t),
V(t)沿着误差系统(6)求导,
当t∈[tk,sk]时,
V ·(t)=2∑Ni=1eTi(t)e·i(t)=
2∑Ni=1eTi(t)[Aiei(t)+Biηi(t)+
c1∑Ni=1lijΓej(t)+c2∑Ni=1gijΓh(ej(t-τ(t)))+
mi(t)-c1dΓe1(t)]=
V1(t)+V2(t)+V3(t)+V4(t)+V5(t)(7)
根据假设2,得0≤(ηil)/(cTlei(t))≤δl, l=1,2,…,m, 上式等价于ηil(ηil-δlcTlei(t))≤0,(8)
由于τl≥0, l=1,2,…,m, 根据式(8)可得
∑ml=1τlηil(ηil-δlcTlei(t))≤0,
那么, ηTi(Tηi-TΔCei(t))≤0, i=1,2,…,N,
于是得
V1(t)=2∑Ni=1eTi(t)[Aiei(t)+Biηi(t)]=
2eT(t)A ^e(t)+2eT(t)B ^η(t)≤
eT(t)(A ^+A ^T)e(t)+2eT(t)B ^η(t)-
2ηT(t)(INT)η(t)+
2ηT(t)(INTΔC)e(t)-
2α1eT(t)(INΓ)e(t)+
2α1eT(t)(INΓ)e(t),
记y(t)=(eT(t),ηT(t))T, 注意到定理1中条件ⅰ, 有
V1(t)≤yT(t)Πy(t)+2α1eT(t)(INΓ)e(t)≤
2α1eT(t)(INΓ)e(t)(9)
V2(t)=2c1∑Ni=1eTi(t)∑Nj=1lijΓej(t)=
2c1eT(t)(LΓ)e(t)(10)
根据Cauchy不等式和假设3可得
V3(t)=2c2∑Ni=1eTi(t)∑Nj=1gijΓh(ej(t-τ(t)))=
2c2eT(t)(GΓ)H(e(t-τ(t)))≤
c2γ1eT(t)(GGTΓΓT)e(t)+
(c2)/(γ1)λmax(QTQ)eT(t-τ(t))e(t-τ(t))(11)
此外,由于max1≤i≤N(=ΔAi=+=ΔBi==ΔC=)≤ψ, 有
V4(t)=2∑Ni=1eTi(t)mi(t)≤
γ2∑Ni=1eTi(t)ei(t)+1/(γ2)∑Ni=1mTi(t)mi(t)≤
γ2eT(t)e(t)+N/(γ2)ε^-2ψ2(12)
最后
V5(t)=-2c1d∑Ni=1eTi(t)Γe1(t)=
-2c1eT(t)(DΓ)e(t)(13)
将式(9)至式(13)代入式(7)并结合定理1条件ⅱ和ⅲ,可得
V ·(t)≤eT(t)[2α1(INΓ)+c2γ1(GGTΓΓT)+
γ2I+2c1(L ^Γ)+μ1I]e(t)-
μ1eT(t)e1(t)+((c2)/(γ1)λmax(QTQ)-
μ2)eT(t-τ(t))e(t-τ(t))+
μ2eTi(t-τ(t))e(t-τ(t))+N/(γ2)ε^-2ψ2≤
-μ1eT(t)e(t)+μ2eT(t-τ)e(t-τ)+
N/(γ2)ε^-2ψ2=-μ1V(t)+μ2V(t-τ)+
N/(γ2)ε^-2ψ2
这样,根据引理1可得到
V(t)≤=Vτ(tk)=e-γ(t-tk)+(Nε^-2ψ2)/(rγ2).(14)
当t∈(sk,tk+1)时,用类似方法并注意到定理1中条件ⅰ、ⅲ及ⅳ有
V ·(t)=V1(t)+V2(t)+V3(t)+V4(t)≤
eT(t)[2α1(INΓ)+c2γ1(GGTΓΓT)+
γ2I-(μ3-μ1)I+2c1(LΓ)]e(t)+
(μ3-μ1)eT(t)e(t)+μ2eT(t-τ)e1(t-τ)+
((c2)/(γ1)λmax(QTQ)-μ2)eT(t-τ)e(t-τ)+
N/(γ2)ε^-2ψ2=(μ3-μ1)V(t)+μ2V(t-τ)+
(Nε^-2ψ2)/(γ2).
由引理2,可得
V(t)≤(=Vτ(sk)=+(Nε^-2ψ2)/(γ2(μ3-μ1+μ2)))×
e<sup>(μ3-μ1+μ2)(t-sk)-(Nε^-2ψ2)/(γ2(μ3-μ1+μ2))(15)
最后,由引理3,以及式(14)和式(15)可得
V(t)≤Vτ(0)eρe-(ρ/ω)t+ε1(16)
其中, α、 β、 ν、 ε1和r如定理1中所定义,
Vτ(0)=Nmax1≤i≤N sup-τ≤κ≤0(φi(κ)-s(0)))T ×
(φi(κ)-s(0)).
又因为V(t)≤=e(t)=2,
则由式(16)可得
=e(t)=≤(Vτ(0)eρe-(ρ/ω)t+ε1)1/2≤
(Vτ(0)eρ)1/2e-(ρ/(2ω))t+(ε1)1/2
故对任意小正数ε2, 存在正数T, 使得对任意t≥T, 有
=e(t)=≤ε2+(ε1)1/2,
由定义1可知,系统(1)和系统(2)取得拟同步,证毕.
当 Ai=A, Bi=B时,显然max1≤i≤N(=ΔAi=+=ΔBi==ΔC=)≡0, 即Lur'e网络(2)与系统(1)的参数匹配,相应地可得α=0、 β=0、 ν=0及ε1=0, 则得到以下推论.
推论1 若假设1至假设3成立,且存在正数γ1, μ3>μ2>μ1, a1和矩阵 Q=(qij)n×n, 满足条件
ⅰ. Π=[A+AT-2α1Γ B+CTΔT
BT+TΔC -2T]≤0,
ⅱ. μ1I+2c1(L ^Γ)+c2γ1(GGTΓΓT)+2α1(INΓ)≤0;
ⅲ.(c2)/(γ1)λmax(QTQ)-μ2≤0;
ⅳ. -(μ3-μ1)I+2c1(LΓ)+c2γ1(GGTΓΓT)+2α1(INΓ)≤0;
ⅴ. ρ=ρ1-ρ2=r(θ-τ)-(μ3-μ1+μ2)(ω-θ)>0, r是方程r-μ1+μ2eγτ=0的唯一正解.那么, =e(t)=≤(Vτ(0)eρe-(ρ/ω)t)1/2,当t→ +∞时, =e(t)=→0, 可见,系统(1)和系统(2)达到完全同步.
考虑由6个节点构成的异构Lur'e网络,每个节点的维数为3,模型由下面微分方程给出
s·(t)=As(t)+Bf(Cs(t)),
其中, s(t)=[s1, s2, s3]T;
A=[-18/7 9 0
1 -1 1
0 -14 0]; B=[14
0
0];
C=diag(1 0 0);
f(Cs(t))=[(|s1+1|-|s2-1|)/2 0 0]T.
内部耦合矩阵 Γ=diag(3, 3, 3);
外部耦合矩阵
G=[-14 4 0 10 0 0
4 -7 3 0 0 0
0 3 -3 0 0 0
10 0 0 -16 6 0
0 0 0 6 -11 5
0 0 0 0 5 -5], L=G,
选择控制增益d=6, τ=0.01. 考虑非周期间歇牵制控制,这里[tk,sk]定义为
[0,3]∪[3.010,6.105]∪[6.110,9.125]∪[9.195,12.195]∪[12.2,15.2]∪[15.250,18.305]∪[18.35,21.35]∪[21.405,24.425]∪[24.505,27.525]∪…
可见, θ=3; ω=3.1. 图1为节点状态图,图2为目标节点状态图.为满足定理中的条件,通过LMI计算,可得α1=4.29、 c1=10、 c2=0.35、 γ1=0.36、 γ2=0.32、 μ1=3、 μ2=1及μ3=58, 由r-μ1+μ2eγτ=0, 可得r=1.98. 取 Ai=(1+i/1 000)A, Bi=(1+i/1 000)B, i=1,2,…,6, 由图3可知,当t→+∞时, =e(t)=≤10, 通过计算max1≤i≤6=ΔAi=+=ΔBi==ΔC=≤0.108 8, 代入相关式子估算得到误差界为43.92,由此可知异构复杂网络(2)和系统(1)可以达到拟同步.
本研究基于非周期间歇牵制控制策略,讨论带有参数不匹配的非线性耦合时变时滞的异构Lur'e网络的拟同步问题.通过设计合理的Lyapunov函数,利用线性矩阵不等式以及严格的理论分析,实现了Lur'e网络的拟同步.最后,通过数值模拟验证所得结果的有效性.
深圳大学学报理工版
JOURNAL OF SHENZHEN UNIVERSITY SCIENCE AND ENGINEERING
(1984年创刊 双月刊)
主 管 深圳大学
主 办 深圳大学
编辑出版 深圳大学学报理工版编辑部
主 编 阮双琛
国内发行 深圳市邮电局
国外发行 中国国际图书贸易集团有限公司(北京399信箱)
地 址 北京东黄城根北街16号
邮 编 100717
电 话 0755-26732266
0755-26538306
Email journal@szu.edu.cn
标准刊号 ISSN 1000-2618
CN 44-1401/N