多層單向耦合星形網絡的特征值譜及同步能力分析?

孫娟李曉霞 張金浩申玉卓李艷雨

1)(省部共建電工裝備可靠性與智能化國家重點實驗室(河北工業大學),天津 300130)

2)(河北省電磁場與電器可靠性重點實驗室(河北工業大學),天津 300130)

多層單向耦合星形網絡的特征值譜及同步能力分析?

孫娟1)2)李曉霞1)2)?張金浩1)2)申玉卓1)2)李艷雨1)2)

1)(省部共建電工裝備可靠性與智能化國家重點實驗室(河北工業大學),天津 300130)

2)(河北省電磁場與電器可靠性重點實驗室(河北工業大學),天津 300130)

多層星形網絡,有向網絡,特征值譜,同步能力

1 引 言

近幾年,復雜網絡研究主要集中在網絡的拓撲結構與其動力學行為之間的關系[1?4]、擴散動力學[5?9]、同步[10?14]等方面.目前大多數復雜網絡的研究都集中在單個網絡[15?27],然而現實世界中各類系統之間并不是獨立的,而是相互聯系的.研究結果表明,單個網絡建模的方法一般不足以確定出現在實際復雜系統中的多層網絡間的真正的相互依存關系.比如,一些不同社交網絡的用戶關系[28]、運輸系統[29]以及大腦的學習型組織[30].盡管多層網絡的研究還處于起步階段,但已經涌現出一些有影響力的研究成果[31,32],這將是復雜網絡理論的重要研究方向.

研究復雜網絡的終極目標之一是理解網絡拓撲結構對物理過程的影響.網絡上的同步現象是一項重要的研究課題[33,34].同步現象大量出現在我們的生活中,如掌聲同步、路由器同步、螢火蟲同時閃爍、物理振子以及生物振子的同步等,在物理、化學、生物、工程技術、經濟以及社會科學等領域都觀察到了各種各樣的同步現象[35?39].復雜網絡作為一個載體展示了豐富多彩的網絡同步現象,不同領域的學者從不同的角度對復雜網絡同步進行了深入研究.

為了研究更接近于實際復雜網絡的同步能力,學者們提出了一種基于多層網絡的擴散動力學方程的超拉普拉斯矩陣結構,根據其特征值譜來分析網絡同步能力[40].根據主穩定性函數方法,復雜系統的非線性動力學行為與網絡的拉普拉斯矩陣的第二個特征值有關[41],在擴散過程中,收斂于最大熵狀態[42],而且耦合網絡的拉普拉斯矩陣的最大特征值對網絡的同步穩定性具有決定性影響[43,44].文獻[45,46]研究了典型的規則兩層星形網絡,提出并嚴格推導出兩層星形網絡的特征值譜,并分析了同步能力,但由于實際復雜網絡更接近于有向網絡,本文即在此基礎上研究了更具一般性的多層有向星形網絡的特征值譜及同步能力.首先,嚴格推導出三層單向耦合的星形網絡的特征值多項式,分別分析了在同步域無界和有界兩種情況下,網絡層數、每層網絡節點數、層內耦合強度、層間耦合強度(層間中心節點之間的耦合強度和層間葉子節點之間的耦合強度)與網絡同步能力的關系.其次,進一步嚴格推導出多層單向耦合的星形網絡的特征值多項式以及特征值譜,分析了在同步域無界和有界兩種情況下,網絡層數、每層網絡節點數、層內耦合強度、層間耦合強度(層間中心節點之間的耦合強度和層間葉子節點之間的耦合強度)與網絡同步能力的關系.最后,分析各個參數對網絡同步能力的影響并進行了仿真驗證.

2 模型與方法

2.1 模 型

M層網絡中第K層、第i個節點的動力學方程如下:

我們假設多層網絡為M層網絡,每層有N個節點.α層的層內連接強度矩陣W(α)∈RN×N,S(α)∈RN×N為層內點強度的對角矩陣,層的層內拉普拉斯矩陣為L(α)=S(α)?W(α), 假設多層網絡層間的所有節點的層間連接完全一致,我們定義WI∈RM×M為每兩層網絡的層間連接強度矩陣,SI∈RM×M為每兩層網絡的層間點強度對角矩陣,每兩層網絡的層間拉普拉斯矩陣為LI=SI?WI.

我們得到整個網絡的超拉普拉斯矩陣:L=LL+LI.我們把L分成兩部分,LL為層內的超拉普拉斯矩陣,LI為層間的超拉普拉斯矩陣.LL為各層內拉普拉斯矩陣的直和,

LI為層間拉普拉斯矩陣與M×N的單位矩陣的Kronecker積,

2.2 方 法

本文采用主穩定性函數法來判定網絡的同步能力的穩定性.一個由N個相同的動力學系統˙Xi=f(Xi)為節點所構成的連續時間耗散耦合動態網絡,其狀態方程如下:

主穩定性函數法是研究線性耦合同步穩定性問題的一種主要方法,并被廣泛采用.下面主要詳細介紹主穩定性函數法.我們將連續時間動態網絡(4)式關于同步流形S(t)進行線性穩定性分析.令δi為第i個節點狀態的變分,可以得到下面的變分方程:

這里Jf和JH分別是f(X)和H(X)關于同步流形S的Jacobi矩陣,通常要求為有界.令則(5)式可寫為矩陣形式:

下面根據Lyapunov指數法來判斷同步流形的穩定性.考慮到方程中只有λk和ηk與k有關,當外耦合矩陣W為非對稱時,其特征值可能是復數,故定義主穩定性方程為

其最大Lyapunov指數γmax是實變量α和β的函數,稱為動力學網絡的主穩定性函數.已知一個給定的耦合強度c,在(α,β)所在的復平面上對應一個點cλk,如果該點對應的最大Lyapunov指數γmax為正數時,則該特征態為不穩定態;反之,則該特征態為穩定態.如果對應所有的特征態λk(k=2,3,···,N), 其所有點對應的最大Lyapunov指數γmax均為負數時,那么認定在該耦合強度c下,整個網絡的同步流形是漸近穩定的.

根據主穩定函數方法,本文網絡模型的同步能力由L的非零最小特征值λ2和最大特征值λmax決定[47?50].當同步域無界時,網絡的同步能力由λ2決定,并與λ2成正比;當同步域有界時,網絡的同步能力由λmax/λ2決定,此時網絡的同步能力與r=λmax/λ2成反比.對于同步化區域為空集的情況,無論耦合強度c怎么改變都不能實現網絡的同步,所以無需討論該情況.

3 特征值譜及同步能力分析

3.1 三層單向耦合星形網絡的結構

本節首先嚴格推導三層單向耦合的星形網絡的特征值譜,其中單向耦合包含葉子節點向中心節點的單向耦合和中心節點向葉子節點的單向耦合,耦合模型如圖1(a)和圖1(b)所示,層間耦合方式如圖2所示.圖1(a)的每層星形網絡內,葉子節點向該層中心節點耦合強度為a;圖1(b)的每層星形網絡內,中心節點向該層葉子節點耦合強度為a;圖2中每兩層網絡間葉子節點與葉子節點一一對應耦合,中心節點與中心節點一一對應耦合,耦合方式為雙向耦合,葉子節點之間耦合強度為d,中心節點之間耦合強度為d0.

圖1 (網刊彩色)(a)葉子節點向中心節點的單向耦合;(b)中心節點向葉子節點的單向耦合;黑色虛線表示每兩層間葉子節點耦合,紅色虛線表示每兩層間中心節點耦合,層內耦合強度為a,層間葉子節點之間的耦合強度為d,層間中心節點之間的耦合強度為d0Fig.1.(color online)(a)The coupling model is that leaf nodes point to the hub node;(b)the coupling model is that the hub node point to leaf nodes.The black dotted line represents the coupling connection between the leaf nodes of every two layers of the network,the red dotted line represents the coupling connection between the hub nodes of every two layers of the network,the intralayer coupling strength a the inter-layer leaf nodes coupling strength d,the interlayer hub nodes coupling strength d0.

圖2 (網刊彩色)層間耦合方式,每兩層葉子節點之間一一對應雙向耦合,每兩層中心節點之間一一對應雙向耦合;黑色虛線表示層間葉子節點之間的耦合連接,紅色虛線表示層間中心節點之間的耦合連接Fig.2.(color online)The coupling model between the layers.The coupling between the leaf nodes of every two-layer of network is one-to-one bidirectional coupling,and the coupling between the hub nodes of every two-layer of network is one-to-one bidirectional coupling;the black dotted lines represent the interlayer leaf nodes,and the red dotted lines represent the coupling between the inter-layer hub nodes.

3.2 三層單向耦合星形網絡的特征值譜及同步能力分析

首先,這里我們假設三層網絡的每層網絡的拓撲結構是完全一致的星形網絡,每兩層網絡間葉子節點之間、中心節點之間分別都是一一對應耦合;每層網絡的節點數都為N,整個網絡層數為M,層內耦合強度為a,層間中心節點之間的耦合強度為d0,層間葉子節點之間的耦合強度為d.對于簡單有向圖G,其拉普拉斯矩陣L(G)=S(G)?W(G),S(G)為頂點出度的對角矩陣,W(G)=[aij]為其鄰接矩陣,aij為節點Vi到Vj的邊的權值[51].由此得到的三層葉子節點向中心節點的單向耦合星形網絡的超拉普拉斯矩陣為

三層中心節點向葉子節點的單向耦合星形網絡的超拉普拉斯矩陣為

三層中心節點向葉子節點的單向耦合星形網絡的超拉普拉斯矩陣的特征值多項式為故特征值為:

在實際網絡中,節點數N是非常大的,遠遠大于層內耦合強度a和層間節點之間的耦合強度d與d0.其同步能力的強弱根據同步域無界和有界兩種情形來考慮,即討論非零最小特征值λ2(值越大同步性越強)與r=λmax/λ2(值越小同步性越強).

三層葉子節點向中心節點單向耦合的星形網絡,當同步域無界時,網絡的同步能力由λ2=min{a,3d0}決定;當同步域有界時,網絡的同步能力由r=(a+3d)/min{a,3d0}決定,而無論是同步域無界還是有界,網絡的同步能力都與節點數無關.

三層中心節點向葉子節點單向耦合的星形網絡,當同步域無界時,網絡的同步能力由λ2=3d決定;當同步域有界時,網絡的同步能力由r=[(N?1)a+3d0]/3d決定.

3.3 多層單向耦合星形網絡的特征值譜及同步能力分析

本文進一步推導出多層單向耦合星形網絡的超拉普拉斯矩陣及特征值多項式,多層葉子節點向中心節點單向耦合的星形網絡的超拉普拉斯矩陣為

多層中心節點向葉子節點單向耦合的星形網絡的超拉普拉斯矩陣為

多層中心節點向葉子節點單向耦合的星形網絡的特征值多項式為:故特征值為

多層葉子節點向中心節點單向耦合的星形網絡的同步能力,按同步域為無界和有界兩種情形分別討論,當同步域無界時,網絡的同步能力由λ2=min{a,Md0}決定;當同步域有界時,網絡的同步能力由r=(a+Md)/min{a,Md0}決定,除了d,a,d0,多層網絡的層數對網絡的同步能力也會產生影響.

多層中心節點向葉子節點單向耦合的星形網絡,按同步域為無界和有界兩種情形分別討論,當同步域無界時,網絡的同步能力由λ2=Md決定;當同步域有界時,網絡的同步能力由r=[(N?1)a+Md0]/Md決定.除了N,d,a,d0,網絡的同步能力也受到多層網絡的層數的影響.

4 數值仿真及結果

4.1 三層葉子節點向中心節點單向耦合的星形網絡的同步能力仿真

同步域無界時,λ與a,3d0中的最小值有關;同步域有界時,同步能力由r=(a+3d)/min{a,3d0}來決定,網絡的同步能力與N無關.仿真情況歸納見表1.

1)同步域無界

同步域無界、a<3d0時,λ2與a有關;a>3d0時,λ2與d0有關.分別變化a和d0,所得到的λ2如圖3所示.λ2隨著a的增大而增大,隨著d0的增大而增大,同步能力增強,增大層間中心節點之間的耦合強度更有利于提高同步能力,這是因為三層網絡類似于三個聚類的社團結構,聚類塊間連接緊密,從而使同步能力變得更強.最后,從圖3我們可以得到:當同步域無界時,對于層內耦合強度和層間中心節點之間的耦合強度,同步能力依賴于較弱的一方;在兩種耦合強度同步改變時,層間中心節點之間的耦合強度對同步能力產生的影響更大一些.

表1 對于三層葉子節點向中心節點單向耦合的星形網絡,λ2,r= λmax/λ2隨d,a,d0的變化Table 1.When the coupling model is that leaf nodes point to the hub node for three-layer star networks,λ2 and r= λmax/λ2changes with d,a,d0.

圖3 (網刊彩色)對于三層葉子節點向中心節點單向耦合的星形網絡,λ2隨層內耦合強度a的變化(紅色實線),a<3d0;λ2隨層間中心節點之間的耦合強度d0的變化(藍色實線),a>3d0Fig.3.(color online)When the coupling model is that leaf nodes point to the hub node for three-layer star networks,λ2change with the intra-layer coupling strength a(red solid line),a<3d0;λ2change with inter-layer coupling strength between the hub nodes d0(blue solid line),a>3d0.

2)同步域有界

同步域有界、a<3d0時,同步能力由r=1+3d/a決定;a>3d0時,同步能力由r=(a+3d)/3d0決定.

圖4 (網刊彩色)對于三層葉子節點向中心節點單向耦合的星形網絡,r隨層內耦合強度a的變化(紅色實線);r隨層間葉子節點之間的耦合強度d的變化(藍色實線);其中a<3d0Fig.4.(color online)When the coupling model is that leaf nodes point to the hub node for three-layer star networks,r change with the intra-layer coupling strength a(red solid line)and the inter-layer coupling strength between the leaf nodes d(blue solid line),a<3d0.

2?a>3d0時,取d=1,d0=0.01,變化a,取a=1,d0=0.01,變化d;取d=1,a=1,變化d0.得到的r如圖5所示.r隨a的增大而增大,r隨d的增大而增大;r隨d0增大呈指數減小.層內耦合強度越小,層間葉子節點之間的耦合強度越小,層間中心節點之間的耦合強度越大,同步能力變得越強.最后,從圖5(a)我們可以得到:在上述假定條件下,當層內耦合強度和層間葉子節點之間的耦合強度值取1時,同步能力相等.從圖5(b)和圖4我們可以得到:當同步域有界時,網絡的同步能力依賴于層內耦合強度和層間中心節點之間的耦合強度中較弱的一方.

圖5 (網刊彩色)對于三層葉子節點向中心節點單向耦合的星形網絡 (a)r隨層內耦合強度a的變化(紅色實線)以及r隨層間葉子節點之間的耦合強度d(藍色實線)的變化;(b)r隨層間中心節點之間的耦合強度d0的變化;其中a>3d0Fig.5.(color online)When the coupling model is that leaf nodes point to the hub node for three-layer star networks,(a)r change with the intra-layer coupling strength a(red solid line)and the inter-layer coupling strength between the leaf nodes d(blue solid line),respectively;(b)r change with the inter-layer coupling strength between the hub nodes d0;a>3d0.

4.2 三層中心節點向葉子節點單向耦合的星形網絡的同步能力的仿真

三層中心節點向葉子節點單向耦合的星形網絡,當同步域無界時,網絡的同步能力由λ2=3d決定;當同步域有界時,網絡的同步能力由r=[(N?1)a+3d0]/3d決定.仿真情況歸納于表2.

表2 對于三層中心節點向葉子節點耦合的星形網絡,λ2和r=λmax/λ2隨N,M,d,a,d0的變化Table 2.When the coupling model is that the hub node point to leaf nodes for three-layer star networks,λ2and r= λmax/λ2changes with N,M,d,a,d0.

1)同步域無界

λ2與d有關,變化d,如圖6所示,λ2隨d的增大而增大.層間葉子節點之間的耦合強度增大,同步能力越強.中心節點向葉子節點單向耦合的三層星形網絡的超拉普拉斯矩陣的非零最小特征值是一定的,不需要討論,它與層間葉子節點之間的耦合強度有關,而與層內耦合強度和層間中心節點之間的耦合強度無關.

圖6 (網刊彩色)對于三層中心節點向葉子節點耦合的星形網絡,λ2隨層間葉子節點之間的耦合強度d的變化Fig.6.(color online)When the coupling model is that the hub node point to leaf nodes for three-layer star networks,λ2change with the inter-layer coupling strength between the leaf nodes d.

2)同步域有界

網絡的同步能力由r=[(N?1)a+3d0]/3d決定.取N=200,d=d0=1,變化a;取N=200,a=d0=1,變化d;取N=200,d=a=1,變化d0;取a=d=d0=1,變化N;得到的r如圖7所示.r隨a的增大而增大;r隨d的增大呈指數減小;r隨d0的增大而增大;r隨N的增大而增大.層內耦合強度越小,層間葉子節點之間的耦合強度越大,層間中心節點之間的耦合強度越小,網節點數越少,網絡的同步能力變的越強.從圖7我們得到:在上述假定條件下,網絡的同步能力更依賴于層間葉子節點之間的耦合強度,而層間中心節點之間的耦合強度對網絡同步能力的影響微乎其微.

圖7 (網刊彩色)對于三層中心節點向葉子節點耦合的星形網絡,(a)r隨層內耦合強度a的變化(紅色實線),r隨層間葉子節點之間的耦合強度d的變化(藍色實線),r隨層間中心節點之間的耦合強度d0的變化(黑色實線);(b)r隨節點數N 的變化Fig.7.(color online)When the coupling model is that the hub node point to leaf nodes for three-layer star networks,(a)r change with the intra-layer coupling strength a(red solid line),the inter-layer coupling strength between the leaf nodes d(blue solid line),and the inter-layer coupling strength between the hub nodes d0(black solid line),respectively;(b)r change with the number of nodes N.

4.3 多層葉子節點向中心節點單向耦合的星形網絡的同步能力的仿真

多層葉子節點向中心節點單向耦合的星形網絡,當同步域無界時,網絡的同步能力由λ2=min{a,Md0}決定;當同步域有界時,網絡的同步能力由r=(a+Md)/min{a,Md0}決定;網絡的同步能力與N無關.仿真情況歸納于表3.

情況 3 v7,v8全染顏色1,v2,v6染1,v4只能染2,則可繼續用上述方法將窮點v3改染為1, 并用3來染v。

表3 對于多層葉子節點向中心節點單向耦合的星形網絡,λ2和r=λmax/λ2隨M,d,a,d0的變化Table 3.When the coupling model is that leaf nodes point to the hub node for multilayer star networks,λ2 and r= λmax/λ2changes with N,M,d,a,d0.

1)同步域無界

2)同步域有界

圖8 (網刊彩色)對于多層葉子節點向中心節點單向耦合的星形網絡,(a)λ2隨層間中心節點之間的耦合強度d0的變化;(b)λ2隨層數M的變化Fig.8.(color online)When the coupling model is that leaf nodes point to the hub node for multilayer star networks,(a)λ2change with the inter-layer coupling strength between the hub nodes d0;(b)λ2change with the number of layers M.

圖9 (網刊彩色)對于多層葉子節點向中心節點單向耦合的星形網絡,(a)r隨層內耦合強度a的變化(紅色實線),r隨層間葉子節點之間的耦合強度d的變化(藍色實線);(b)r隨層數M的變化Fig.9.(color online)When the coupling model is that leaf nodes point to the hub node for multilayer star networks,(a)r change with the intra-layer coupling strength a(red solid line)and the inter-layer coupling strength between the leaf nodes d(blue solid line),respectively;(b)r change with the number of layers M.

圖10 (網刊彩色)對于多層葉子節點向中心節點單向耦合的星形網絡,(a)r隨層內耦合強度a的變化(紅色實線),r隨層間葉子節點之間的耦合強度d的變化(藍色實線);(b)r隨層間中心節點之間的耦合強度d0的變化;其中a>Md0Fig.10.(color online)When the coupling model is that leaf nodes point to the hub node for multilayer star networks,(a)r changes with the intra-layer coupling strength a(red solid line)and the inter-layer coupling strength between the leaf nodes d(blue solid line),respectively;(b)r change with the inter-layer coupling strength between the hub nodes d0;a>Md0.

a>Md0時,網絡的同步能力還受到層數的影響.取a=1,d=1,d0=0.01,變化M,得到r(圖11).r隨著M的增大呈指數減小.層數越少,同步能力越強.我們從圖11可以得到:在上述假定條件下,與a

圖11 (網刊彩色)對于多層葉子節點向中心節點單向耦合的星形網絡,r隨層數M的變化,a>Md0Fig.11.(color online)When the coupling model is that leaf nodes point to the hub node for multilayer star networks,r change with the number of layers M,a>Md0.

4.4 多層中心節點向葉子節點單向耦合的星形網絡的同步能力的仿真

多層中心節點向葉子節點單向耦合的星形網絡,當同步域無界時,網絡的同步能力由λ2=Md決定;當同步域有界時,網絡的同步能力由r=[(N?1)a+Md0]/(Md)決定.仿真情況歸納于表4.

表4 對于多層中心節點向葉子節點耦合的星形網絡,λ2和r=λmax/λ2隨N,M,d,a,d0的變化Table 4.When the coupling model is that the hub node point to leaf nodes for multilayer star networks,λ2 and r= λmax/λ2changes with N,M,d,a,d0.

1)同步域無界

λ2與M,d有關,取d=1,改變M;取M=20,變化d.得到的λ2如圖12所示.λ2隨M的增大而增大,λ2隨d的增大而增大.層間葉子節點之間的耦合強度越大,層數越多,網絡的同步能力越強.

2)同步域有界

網絡的同步能力由r=[(N?1)a+Md0]/(Md)決定.取N=200,M=20,d=d0=1,變化a;取N=200,M=20,a=d0=1,變化d;取N=200,M=20,d=a=1,變化d0;取a=d=d0=1,N=200,變化M.得到的r如圖13所示.r隨著a的增大而增大,r隨著d的增大而減小,r隨著d的增大而增大,r隨著M的增大而減小.層內耦合強度越小,層間葉子節點之間的耦合強度越大,層間中心節點之間的耦合強度越小,層數越大,網絡的同步能力越強.從圖13我們得到:在上述假定的條件下,網絡的同步能力更依賴于層間葉子節點之間的耦合強度,層間中心節點之間的耦合強度對網絡同步能力的影響微乎其微.比較圖7(a)與圖13可以得到:層內耦合強度和層間耦合強度對網絡同步能力的影響可以由三層網絡推廣到多層網絡;對于多層網絡,網絡的同步能力也依賴于層數.

圖12 (網刊彩色)對于多層中心節點向葉子節點耦合的星形網絡,(a)λ2隨層數M的變化;(b)λ2隨層間葉子節點之間的耦合強度d的變化Fig.12.(color online)When the coupling model is that the hub node point to leaf nodes for multilayer star networks,(a)λ2change with the number of layers M;(b)λ2change with the inter-layer coupling strengthbetween the leaf nodes d.

圖13 (網刊彩色)對于多層中心節點向葉子節點耦合的星形網絡,(a)r隨層內耦合強度a的變化(紅色實線),r隨層間葉子節點之間的耦合強度d的變化(藍色實線),r隨層間中心節點之間的耦合強度d0的變化(黑色實線);(b)r隨層數M的變化Fig.13.(color online)When the coupling model is that the hub node point to leaf nodes for multilayer star networks,(a)r change with the intra-layer coupling strength a(red solid line)and the inter-layer coupling strength between the leaf nodes d(blue solid line)and the inter-layer coupling strength between the hub nodes d0(black solid line),respectively;(b)r change with the number of layers M.

取a=d=d0=1,M=20,變化N,得到r(圖14).r隨著N的增大而增大.節點數越少,同步能力越強.比較圖7(b)和圖14可以得到:節點數對網絡同步能力的影響與層數無關;節點數對網絡同步能力的影響可以由三層網絡推廣到多層網絡.

圖14 (網刊彩色)對于多層中心節點向葉子節點耦合的星形網絡,r隨節點數N的變化Fig.14.(color online)When the coupling model is that the hub node point to leaf nodes for multilayer star networks,r change with the number of nodes N.

5 討論與結論

本文在理論上嚴格推導出三層、多層單向耦合星形網絡的超拉普拉斯矩陣,得到其特征值譜,并分析了同步能力.對于葉子節點向中心節點耦合的三層及多層星形網絡:當同步域無界時,隨著層內耦合強度和層間中心節點之間的耦合強度的增大會使同步能力增強,且該影響可以由三層網絡推廣到多層網絡;而對于多層網絡,層數的增大也會使同步能力增強.當同步域有界時,在層內耦合強度較弱的情況下,層內耦合強度的增大會使同步能力增強,層間葉子節點之間的耦合強度和層數的增大會使同步能力變弱.在層間中心節點之間的耦合強度較弱的情況下,層間中心節點之間的耦合強度和層數的增大會使同步能力增強,而層內耦合強度和層間葉子節點之間的耦合強度的增大反而會使同步能力變弱.對于中心節點向葉子節點耦合的三層及多層星形網絡:當同步域無界時,層間葉子節點之間的耦合強度和層數的增大會使同步能力增強;當同步域有界時,層間葉子節點之間的耦合強度和層數的增大會使同步能力增強,而層內耦合強度、節點數和層間中心節點之間的耦合強度的增大反而會使同步能力減弱.本文結論與文獻[45]的結論部分一致,單向耦合的多層星形網絡與雙向耦合的兩層星形網絡在結構連接方式上具有一定的聯系.相同之處在于:當同步域有界時,文獻[45]和本文的葉子節點向中心節點耦合的三層及多層星形網絡在層內耦合強度較弱的情況下,都得出層內耦合強度的增大會增強網絡的同步能力,層間葉子節點之間的耦合強度的增大會減弱網絡的同步能力.通過比較各參數對單向耦合的多層星形網絡和雙向耦合的兩層星形網絡同步能力的影響,也得出它們之間的一些不同之處:文獻[45]設定層間中心節點之間的耦合強度d0和層間葉子節點之間的耦合強度d相等,在同步域有界的情況下,文獻[45]討論的是層內耦合強度較弱的情況和層間葉子節點之間的耦合強度較弱的情況,而本文沒有做此設定.對于葉子節點向中心節點耦合的三層及多層星形網絡,我們討論的是層內耦合強度較弱的情況和層間中心節點之間的耦合強度較弱的情況;而對于中心節點向葉子節點耦合的三層及多層星形網絡的超拉普拉斯矩陣的最小非零特征值是確定的,不需要討論.相對于文獻[45]研究的兩層網絡,本文重點研究了多層網絡,并且得出了網絡的同步能力深受層數的影響.

由于多層網絡復雜龐大,本文對有向網絡進行了簡單一般化,得到當同步域無界時,網絡層數的增大使網絡的同步能力增強;在同步域有界的情況下,當層間中心節點之間的耦合強度較弱時,層數對兩種有向耦合方式的多層星形網絡的同步能力的影響相同,而當層內耦合強度較弱時,層數對兩種有向耦合方式的多層星形網絡的同步能力的影響不同;對于葉子節點向中心節點單向耦合的多層星形網絡,節點數的變化對網絡的同步能力沒有影響.由于簡單一般化的網絡與實際網絡的復雜無規則性還有一定的差異,目前,多層網絡的研究仍然還有很多問題沒有解決:譬如本文的多層單向耦合的星形網絡,層間在有向耦合的情況下,網絡的特征值譜以及同步能力問題;對于多層星形網絡,在層內耦合強度改變了網絡的同步能力的情況下,如何改變其他參數來使同步能力保持不變.特征值譜方法對研究規則網絡的同步能力問題有巨大的貢獻,而Newman-Watts小世界網絡介于最近鄰耦合規則網絡與全局耦合網絡之間,是否可以推廣到多層小世界網絡的研究也有待驗證.目前多層網絡的理論模型在應用層面還較為薄弱,主要由于理論模型與實際數據存在一定差距,未來隨著多層網絡理論越來越逼近現實系統,會產生較好的應用前景.

[1]Tang L,Lu J A,Wu X,Lü J H 2013Nonlinear Dyn.73 1081

[2]Wang X F 2002Int.J.Bifur.Chaos12 885

[3]Barrat A,Barthélemy M,Vespignani A 2004Phys.Rev.Lett.92 228701

[4]Wang W X,Wang B H,Hu B,Yan G,Ou Q 2005Phys.Rev.Lett.94 188702

[5]Ide Y,Izuhara H,Machida T 2016Physics457 331

[6]Dahan M,Levi S,Luccardini C,Rostaing P,Riveau B,Triller A 2003Science302 442

[7]Raquel A B,Borgeholthoefer J,Wang N,Moreno Y,Gonzálezbailón S 2013Entropy15 4553

[8]Boukobza E,Chuchem M,Cohen D,Vardi A 2009Phys.Rev.Lett.102 180403

[9]Weber S,Hütt M T,Porto M 2008Europhys.Lett.82 28003

[10]Luan Z 2008Phys.Rep.469 93

[11]Arenas A,Díaz-Guilera A,Kurths J,Moreno Y,Zhou C 2008Phys.Reports469 93

[12]Timme M,Wolf F,Geisel T 2004Phys.Rev.Lett.92 074101

[13]Motter A E,Zhou C,Kurths J 2005Phys.Rev.E71 016116

[14]Xue M,Yeung E,Rai A,Roy S,Wan Y,Warnick S 2011Complex Syst.21 297

[15]Cai S,Zhou P,Liu Z 2014Nonlinear Dyn.76 1677

[16]Chen Y,Yu W,Tan S,Zhu H 2016Automatica70 189

[17]Boccaletti S,Latora V,Moreno Y,Chavez M,Hwang D U 2006Phys.Rep.424 175

[18]Massah S,Hollebakken R,Labrecque M P,Kolybaba A M,Beischlag T V,Prefontaine G G 2004Phys.Rev.Lett.93 218701

[19]Song Q,Cao J,Liu F 2010Phys.Lett.A374 544

[20]He P,Jing C G,Fan T,Chen C Z 2014Complexity19 10

[21]Arenas A,Díazguilera A,Pérezvicente C J 2006Phys.Rev.Lett.96 114102

[22]Zhang Q,Zhao J 2012Nonlinear Dyn.67 2519

[23]Zhang Q,Luo J,Wan L 2013Nonlinear Dyn.71 353

[24]Ling L,Li C,Wang W,Sun Y,Wang Y,Sun A 2014Nonlinear Dyn.77 1

[25]Pacheco J M,Traulsen A,Nowak M A 2006Phys.Rev.Lett.97 258103

[26]Zhang J,Small M 2006Phys.Rev.Lett.96 238701

[27]Gómez-Garde?es J,Moreno Y,Arenas A 2007Phys.Rev.Lett.98 034101

[28]Murase Y,T?r?k J,Jo H H,Kaski K,Kertész J 2014Phys.Rev.E90 052810

[29]Cardillo A,Zanin M,Gómez-Garde?es J,Romance M,Amo A J G D,Boccaletti S 2013Eur.Phys.J.Special Topics215 23

[30]Bassett D S,Wymbs N F,Porter M A,Mucha P J,Carlson J M,Grafton S T 2011Proc.Nat.Acad.Sci.USA108 7641

[31]Li Y,Wu X,Lu J,Lu J 2015IEEE Trans.Circ.Syst.II Express Briefs63 206

[32]Kivel? M,Arenas A,Barthelemy M,Gleeson J P,Moreno Y,Porter M A 2014J.Complex Netw.2 203

[33]Lü J H 2008Adv.Mech.38 713(in Chinese)[呂金虎2008力學進展38 713]

[34]Lu J A 2010Complex Syst.Complex Sci.7 19(in Chinese)[陸君安2010復雜系統與復雜性科學7 19]

[35]Baptista M S,Garcia S P,Dana S K,Kurths J 2008Eur.Phys.J.Special Topics165 119

[36]Lee T H,Ju H P,Wu Z G,Lee S C,Dong H L 2012Nonlinear Dyn.70 559

[37]Qin J,Yu H J 2007Acta Phys.Sin.56 6828(in Chinese)[秦潔,于洪潔2007物理學報56 6828]

[38]Wang J,Zhang Y 2010Phys.Lett.A374 1464

[39]Zhao M,Wang B H,Jiang P Q,Zhou T 2005Prog.Phys.25 273(in Chinese)[趙明,汪秉宏,蔣品群,周濤 2005物理學進展25 273]

[40]Gómez S,Díazguilera A,Gómezgarde?es J,Pérezvicente C J,Moreno Y,Arenas A 2013Phys.Rev.Lett.110 028701

[41]Almendral J A,Díazguilera A 2007New J.Phys.9 187

[42]Granell C,Gómez S,Arenas A 2013Phys.Rev.Lett.111 128701

[43]Pecora L M,Carroll T L 1998Phys.Rev.Lett.80 2109

[44]Liang Y,Wang X Y 2012Acta Phys.Sin.61 038901(in Chinese)[梁義,王興元 2012物理學報 61 038901]

[45]Xu M M,Lu J A,Zhou J 2016Acta Phys.Sin.65 028902(in Chinese)[徐明明,陸君安,周進 2016物理學報 65 028902]

[46]Aguirre J,Sevillaescoboza R,Gutiérrez R,Papo D,Buldú J M 2014Phys.Rev.Lett.112 248701

[47]Dabrowski A 2012Nonlinear Dyn.69 1225

[48]Johnson G A,Mar D J,Carroll T L,Pecora L M 1998Phys.Rev.Lett.80 3956

[49]Sun J,Bollt E M,Nishikawa T 2008Europhys.Lett.85 60011

[50]Xu M,Zhou J,Lu J A,Wu X 2015Euro.Phys.J.B88 240

[51]Xu W G,Wang L G 2016Acta Math.Appl.Sin.39 801

Synchronizability and eigenvalues of multilayer star networks through unidirectionally coupling?

Sun Juan1)2)Li Xiao-Xia1)2)?Zhang Jin-Hao1)2)Shen Yu-Zhuo1)2)Li Yan-Yu1)2)

1)(State Key Laboratory of Reliability and Intelligence of Electrical Equipment,Hebei University of Technology,Tianjin 300130,China)

2)(Key Laboratory of Electromagnetic Field and Electrical Apparatus Reliability of Hebei Province,Hebei University of Technology,Tianjin 300130,China)

1 April 2017;revised manuscript

16 May 2017)

Previous studies on multilayer networks have found that properties of multilayer networks show great differences from those of the traditional complex networks.In this paper,we derive strictly the spectra of the Supra-Laplace matrix of three-layer star networks and multilayer star networks through unidirectionally coupling by using the master stability method to analyze the synchronizability of these two networks.Through mathematical analyses of the eigenvalues of the Supra-Laplace matrix,we explore how the node number,the intra-layer coupling strength the inter-layer coupling strength,and the layer number in fluence the synchronizability of multilayer star networks through unidirectionally coupling in two different ways.In particular,we focus on the layer number and the inter-layer coupling strength between the hub nodes,and then we conclude that the synchronizability of networks is greatly affected by the layer number.We f i nd that when the synchronous region is unbounded,the synchronizability of the two different coupling multilayer star networks is related to not only the intra-layer coupling strength or the inter-layer coupling strength between the leaf nodes of the entire network,but also the layer number.If the synchronous region of two different coupling multilayer star networks is bounded,and the intra-layer coupling strength is weak,the synchronizability of the two different coupling multilayer star networks is different with the changing of the intra-layer coupling strength and the inter-layer coupling strength between the leaf nodes and the layer number.If the synchronous region of two different coupling multilayer star networks is bounded,and the inter-layer coupling strength between the hub nodes is weak,the two different coupling multilayer star networks are consistent with the changing of the intra-layer coupling strength and the layer number while different from the inter-layer coupling strength between the leaf nodes and the inter-layer coupling strength between the hub nodes.We find that the node number has no effect on the synchronizability of multilayer star networks through coupling from the hub node to the leaf node.The synchronizability of the network is directly proportional to the layer number,while inversely proportional to the inter-layer coupling strength between the hub nodes.Finally,the effects of the coupling strength,the layer number and the node number on the synchronizability of the two different coupling star networks can be extended from three-layer network to multilayer networks.

*Project supported by the Natural Science Foundation of Hebei Province,China(Grant No.E2011202051).

?Corresponding author.E-mail:lixiaoxia@hebut.edu.cn

multilayer star networks,unidirectional network,eigenvalue spectrum,synchronization

PACS:89.75.Fb,89.75.Hc,89.75.–k,89.70.aDOI:10.7498/aps.66.188901

(2017年4月1日收到;2017年5月16日收到修改稿)

隨著復雜網絡同步的進一步發展,對復雜網絡的研究重點由單層網絡轉向更加接近實際網絡的多層有向網絡.本文分別嚴格推導出三層、多層的單向耦合星形網絡的特征值譜,并分析了耦合強度、節點數、層數對網絡同步能力的影響,重點分析了層數和層間中心節點之間的耦合強度對多層單向耦合星形網絡同步能力的影響,得出了層數對多層網絡同步能力的影響至關重要.當同步域無界時,網絡的同步能力與耦合強度、層數有關,同步能力隨其增大而增強;當同步域有界時,對于葉子節點向中心節點耦合的多層星形網絡,當層內耦合強度較弱時,層內耦合強度的增大會使同步能力增強,而層間葉子節點之間的耦合強度、層數的增大反而會使同步能力減弱;當層間中心節點之間的耦合強度較弱時,層間中心節點之間的耦合強度、層數的增大會使同步能力增強,層內耦合強度、層間葉子節點之間的耦合強度的增大反而會使同步能力減弱.對于中心節點向葉子節點耦合的多層星形網絡,層間葉子節點之間的耦合強度、層數的增大會使同步能力增強,層內耦合強度、節點數、層間中心節點之間的耦合強度的增大反而會使同步能力減弱.

10.7498/aps.66.188901

?河北省自然科學基金(批準號:E2011202051)資助的課題.

?通信作者.E-mail:lixiaoxia@hebut.edu.cn