于 鏑,白雪峰,任偉建,康朝海
(1.東北石油大學電氣信息工程學院,黑龍江大慶163318;2.哈爾濱醫(yī)科大學(大慶校區(qū))醫(yī)學信息學院,黑龍江大慶163319)
近年來,由于復雜網(wǎng)絡大量存在于自然界和人類社會生活中,如WWW、城市交通網(wǎng)、生物代謝網(wǎng)絡和大型電力網(wǎng)絡等,所以復雜網(wǎng)絡研究引起眾多領域學者的極大關注[1-5]。同步是一類十分重要的非線性現(xiàn)象,并且歷史悠久,而復雜網(wǎng)絡同步由于其廣泛的應用前景而倍受研究人員的青睞。目前,復雜網(wǎng)絡漸近同步取得了豐富的研究成果[6-14]。文獻[6]對線性耦合常微分方程所描述的系統(tǒng),考慮存在耦合與非耦合并存時,通過矩陣特征值分析得到局部和全局同步的充分條件。文獻[7]由若干相同結點線性耦合組成的復雜動力學系統(tǒng),通過同步狀態(tài)線性化,用主穩(wěn)定函數(shù)法判斷同步流形穩(wěn)定性實現(xiàn)復雜系統(tǒng)的局部和全局同步。文獻[8]利用鄰近結點信息與度值最大的結點相連耦合實現(xiàn)網(wǎng)絡同步優(yōu)化。文獻[9]考慮了無向網(wǎng)絡且耦合結點均為線性情況時的同步控制,并總結歸納出同步控制與網(wǎng)絡協(xié)調(diào)控制在本質上是一致的,進而為復雜網(wǎng)絡同步控制開辟了新的控制思路。牽引控制是復雜網(wǎng)絡同步控制的有效方法之一[10-15]。其中文獻[12]對一類線性耦合復雜網(wǎng)絡證明了通過足夠大的耦合度,僅牽制部分節(jié)點即可實現(xiàn)整個網(wǎng)絡的同步。更有趣的是,文獻[13]提出了只應用一個控制器的牽制控制方案實現(xiàn)網(wǎng)絡同步。文獻[14,15]則對一般耦合復雜網(wǎng)絡提出線性控制協(xié)議,給出網(wǎng)絡實現(xiàn)同步的充分條件,并引入耦合強度和控制增益的自適應控制機制,提出了更為靈活的控制方案。
以上成果大都是關于復雜網(wǎng)絡的漸近同步,即網(wǎng)絡在無限時間內(nèi)達到同步。在實際應用中,復雜網(wǎng)絡同步通常需要在有限時間內(nèi)實現(xiàn)。因此,對于復雜網(wǎng)絡有限時間同步的研究具有重要的理論意義和應用價值,但該方面的研究成果較少。受多智能體網(wǎng)絡有限時間協(xié)調(diào)控制[16-20]的啟發(fā),筆者針對復雜動態(tài)網(wǎng)絡有限時間同步問題,提出非光滑控制協(xié)議,并給出實現(xiàn)有限時間同步的充分條件。
考慮一類由N個節(jié)點構成的復雜網(wǎng)絡,第i個節(jié)點動態(tài)描述如下
復雜網(wǎng)絡同步解s(t)可以是一個平衡點、周期軌跡或混沌軌跡。復雜網(wǎng)絡(1)的有限時間同步問題描述如下。
定義1 若存在時刻t*>0且對于?i,j和t≥t*都有‖xi(t)-s(t)‖=0,則稱該復雜網(wǎng)絡實現(xiàn)了有限時間同步。
筆者的目標是設計合適的反饋控制律,使復雜網(wǎng)絡(1)在該控制下實現(xiàn)有限時間同步。在給出主要結果前,先給出如下假設和引理。對于網(wǎng)絡(1)各個節(jié)點動態(tài)滿足下面假設。
假設1 存在常值矩陣K,使
引理1[19]考慮如下微分方程
其中c>0,0<α<1,則函數(shù) x(t)將會在有限時間內(nèi)收斂到零點,且對于所有t≥t*,有x(t)=0。
引理 2[20]令xi∈ (i=1,2,…,N),0 < α <1。若xi≥0 且,則
為實現(xiàn)復雜網(wǎng)絡系統(tǒng)(1)各節(jié)點同步,提出非光滑連續(xù)控制協(xié)議,并基于有限時間穩(wěn)定性引理,給出復雜動態(tài)網(wǎng)絡(1)實現(xiàn)有限時間同步的充分條件。不失一般性,令第i=N個節(jié)點作為驅動節(jié)點,則受控復雜網(wǎng)絡描述為
其中ui為第i個節(jié)點的控制協(xié)議,在ui控制作用下使復雜動態(tài)網(wǎng)絡(1)在t≥t*時實現(xiàn)x1(t)=x2(t)=…=xN(t),即實現(xiàn)網(wǎng)絡的有限時間同步。
定義同步誤差為ei(t)=xi(t)-xN(t),i=1,…,N,則定義ε(t)=[e1(t),…,eN(t)]T為網(wǎng)絡誤差矢量。分別設計兩種不同的控制協(xié)議為
其中i=1,…,N-1,函數(shù)),z=[z1,…,zn]T,0 < α <1,且 sigα(z)=[sigα(z1),…,sigα(zn)]T。di為控制增益,則在控制協(xié)議(5)和協(xié)議(6)的作用下,網(wǎng)絡誤差動態(tài)分別為
其中i=1,…,N-1。
定理1 考慮由式(1)描述的復雜動態(tài)網(wǎng)絡,令假設1成立,節(jié)點之間無向連通,在控制協(xié)議(5)的作用下,若λmax[‖K‖IN+c(L-D)]≤0(其中λmax(·)為矩陣的最大特征值),則網(wǎng)絡(1)實現(xiàn)有限時間同步。其中D=diag{d1,…,dN-1,0},L為節(jié)點之間的外部耦合矩陣。
證明 考慮如下Lyapunov函數(shù)
第1個不等式由假設1得出,第2個不等式由引理2和矩陣性質‖K‖IN≥K,K∈n×n得到。則當λmax[‖K‖IN+c(L-D)]≤0時,由Γ正定和引理1可得,網(wǎng)絡誤差在有限時間
內(nèi)達到零,則網(wǎng)絡(1)可在有限時間內(nèi)達到同步。
定理2 考慮由(1)描述的復雜動態(tài)網(wǎng)絡,令假設1成立,節(jié)點間無向連通,在控制協(xié)議(6)的作用下,若λmax[‖K‖IN+c(L-D)]≤0,則網(wǎng)絡(1)同樣可實現(xiàn)有限時間同步,并且有限時間為
證明 考慮如下Lyapunov函數(shù)
前面與定理1證明過程類似,只有最后一步根據(jù)下面的推導得證。令λi為L的特征值,同時wi為相應于λi的特征向量,令{wi}為N中的一組標準正交基。令,可得,進一步有,則有。當時,由引理1可得,網(wǎng)絡誤差在有限時間
內(nèi)達到零,即網(wǎng)絡(1)可在有限時間內(nèi)達到有限時間同步。
Chen系統(tǒng)和Lü系統(tǒng)是典型的復雜網(wǎng)絡,并且其同步研究在保密通信、信號處理和生命科學等方面有著十分廣闊的應用前景,引起了廣大學者的極大關注。下面就針對這兩個復雜網(wǎng)絡分別進行有限時間同步仿真研究。
實例1 首先考慮由5個Chen系統(tǒng)所構成的復雜網(wǎng)絡(1)。其中非線性特性描述如下
其中 xi=[xi1,xi2,xi3]T,取第 5 個節(jié)點作為驅動點,令‖K‖ =30.934 2[14],D=diag{3.8,6,8,6,0},c=32,內(nèi)部和外部耦合矩陣分別取為
則在控制協(xié)議(5)和協(xié)議(6)的作用下,網(wǎng)絡同步誤差曲線分別如圖1和圖2所示。由仿真圖形可知,大約在t=0.6 s后網(wǎng)絡同步誤差趨于零,網(wǎng)絡(1)實現(xiàn)有限時間同步。
圖1 Chen網(wǎng)絡同步誤差曲線(控制協(xié)議(5)作用下)Fig.1 Synchronization error curves of Chen networks(under the action of control protocol(5))
圖2 Chen網(wǎng)絡同步誤差曲線(控制協(xié)議(6)作用下)Fig.2 Synchronization error curves of Chen networks(under the action of control protocol(6))
實例2 考慮由5個Lü系統(tǒng)所構成的復雜網(wǎng)絡(1)。其中非線性特性描述如下
取第5 個節(jié)點作為驅動點,且‖K‖ =16.278 7[21],D=diag{4,4,4,4,0},c=17,內(nèi)部和外部耦合矩陣同上。則在控制協(xié)議(5)和協(xié)議(6)的作用下,網(wǎng)絡同步誤差曲線分別如圖3和圖4所示。由仿真圖形可知,大約在t=0.2 s后網(wǎng)絡同步誤差趨于零,網(wǎng)絡(1)實現(xiàn)有限時間同步。
圖3 Lü網(wǎng)絡同步誤差曲線(控制協(xié)議(5)作用下)Fig.3 Synchronization error curves of Lü networks(under the action of control protocol(5))
圖4 Lü網(wǎng)絡同步誤差曲線(控制協(xié)議(6)作用下)Fig.4 Synchronization error curves of Lü networks(under the action of control protocol(6))
在上述復雜網(wǎng)絡有限時間同步仿真研究過程中,控制增益矩陣D的取值決定網(wǎng)絡有限時間收斂的快速性,但同時也影響收斂過程中同步誤差的大小。如果控制增益矩陣D中的元素選取較大數(shù)值,則能得到快速的有限時間收斂的同步誤差曲線,但此時同步誤差振蕩幅度較大。所以需要選取合適的控制增益得到既快速又振蕩適度的收斂過程。
筆者針對具有線性耦合的復雜動態(tài)網(wǎng)絡有限時間同步問題,提出非光滑連續(xù)控制協(xié)議,并基于有限時間穩(wěn)定性引理,給出了實現(xiàn)有限時間同步的充分條件。在復雜網(wǎng)絡同步過程中,網(wǎng)絡同步不僅受網(wǎng)絡拓撲和結點之間耦合影響,而且在控制協(xié)議中加入非線性函數(shù)使復雜網(wǎng)絡同步特性得到改善,通過仿真實例驗證了所提方案的可行性和有效性。下一步將研究不確定復雜網(wǎng)絡的有限時間同步問題。
[1]BELYKH I,DI BERNARDO M,KURTHS J,et al.Evolving Dynamical Networks[J].Physica D:Nonlinear Phenomena,2014,267(1):1-6.
[2]張志穎,馮秀琴,姚治海,等.運動光格玻色-愛因斯坦凝聚體的混沌同步[J].吉林大學學報:理學版,2012,50(4):779-784.ZHANG Zhiying,F(xiàn)ENG Xiuqin,YAO Zhihai,et al.Chaotic Synchronization of Bose-Einstein Condensed Body in Moving Optical Lattices[J].Journal of Jilin University:Science Edition,2012,50(4):779-784.
[3]趙明,周濤,陳關榮,等.復雜網(wǎng)絡上動力系統(tǒng)同步的研究進展Ⅱ——如何提高網(wǎng)絡的同步能力[J].物理學進展,2008,28(1):23-31.ZHAO Ming,ZHOU Tao,CHEN Guanrong,et al.A Review on Synchronization of Dynamical Systems on Complex NetworksⅡ:How to Enhance the Network Synchronizability[J].Progress in Physics,2008,28(1):23-31.
[4]姜玉秋.Lorenz混沌系統(tǒng)的線性同步控制[J].吉林大學學報:信息科學版,2011,29(6):571-575.JIANG Yuqiu.Linear Synchronization Control of Lorenz Chaos Systems[J].Journal of Jilin University:Information Science Edition,2011,29(6):571-575.
[5]梅蓉.時滯三角型混沌系統(tǒng)的滑模自適應同步控制方法[J].吉林大學學報:信息科學版,2013,31(1):18-24.MEI Rong.Sliding-Mode Adaptive Synchronization Control for Delayed Triangle Chaotic System [J].Journal of Jilin University:Information Science Edition,2013,31(1):18-24.
[6]LU Wenlian,CHEN Tianping.New Approach to Synchronization Analysis of Linearly Coupled Ordinary Differential Systems[J].Physica D:Nonlinear Phenomena,2006,213(2):214-230.
[7]樊炳航,孫鶴旭,張華,等.網(wǎng)絡型耦合復動力系統(tǒng)的同步動力學[J].復雜系統(tǒng)與復雜性科學,2012,9(1):72-77.FAN Binghang,SUN Hexu,ZHANG Hua,et al.Synchronization Dynamics in Networked Coupled Complex Dynamical Systems[J].Complex Systems and Complexity Science,2012,9(1):72-77.
[8]FENG Junliang,ZHAO Ming,LAI Choyheng.Synchronization Optimal Networks Obtained Using Local Structure Information[J].Physical A,Statistical Mechanics and Its Applications,2012,391(21):5280-5281.
[9]LI Zhongkui,DUAN Zhisheng,CHEN Guanrong,et al.Consensus of Multi-Agent Systems and Synchronization of Complex Networks:A Unified Viewpoint[J].IEEE Transactions on Circuits and Systems Ⅰ:Regular Papers,2010,57(1):213-224.
[10]YANG Huizhong,LI Sheng.Pinning Control of A Generalized Complex Dynamical Network Model[J].Journal of Control Theory and Applications,2009,7(1):1-7.
[11]ZHANG Qingzhen,LI Zhongkui.Pinning Control of Complex Lur'e Networks [J].Chinese Physics B,2009,18(6):2176-2181.
[12]QIN Jiahu,ZHENG Weixing,GAO Huijun.A Study of Synchronization of Complex Networks via Pinning Control[C]∥Proceeding of International Symposium on Circuits and Systems.Brazil:[s.n.],2011:1596-1598.
[13]GUO Wanli.Pinning Complex Networks by a Single Controller on One Element of a Node [J].Journal of Computational Information Systems,2012,8(19):7865-7872.
[14]YU Wenwu,CHEN Guanrong,Lü Jinhu.On Pinning Synchronization of Complex Dynamical Networks[J].Automatica,2009,45(2):429-432.
[15]WANG Weiping,LI Lixiang,PENG Haipeng.Adaptive Synchronization of Complex Dynamical Multi-Links Networks with Similar Nodes[J].Mathematical Problems in Engineering,2013(3):1-12.
[16]WANG Long,F(xiàn)ENG Xiao.Finite-Time Consensus Problems for Networks of Dynamic Agents[J].IEEE Transactions on Automatic Control,2010,55(4):950-955.
[17]SHANG Yilun.Finite-Time Consensus for Multi-Agent Systems with Fixed Topologies[J].International Journal of Systems Science,2012,43(3):499-506.
[18]YANG Shaofu,CAO Jinde,LU Jianquan.A New Protocol for Finite-Time Consensus of Detail-Balanced Multi-Agent Networks[J].Chaos:An Interdisciplinary Journal of Nonlinear Science,2012,22(4):043134.
[19]CORTéS J.Finite-Time Convergent Gradient Flows with Applications to Network Consensus[J].Automatica,2006,42(11):1993-2000.
[20]于鏑,伍清河,王垚.非線性網(wǎng)絡有限時間分布式跟蹤控制[J].控制與決策,2013,28(10):1491-1496.YU Di,WU Qinghe,WANG Yao.Finite Time Distributed Tracking Control of Nonlinear Networks[J].Control and Decision,2013,28(10):1491-1496.
[21]梁義,王興元.基于低階矩陣最大特征值的復雜網(wǎng)絡牽制混沌同步[J].物理學報,2012,61(3):532-539.LIANG Yi,WANG Xingyuan.Pinning Chaotic Synchronization in Complex Networks on Maximum Eigenvalue of Low Order Matrix[J].Acta Physica Sinica,2012,61(3):532-539.