0 引言

    近年來(lái),，復(fù)雜網(wǎng)絡(luò)引起了科學(xué)家的廣泛關(guān)注^[1],。復(fù)雜網(wǎng)絡(luò)最重要的統(tǒng)計(jì)特征是小世界效應(yīng)^[2]和無(wú)標(biāo)度特性^[3]，其中,，小世界網(wǎng)絡(luò)具有小的特征路徑長(zhǎng)度和大的簇系數(shù),，而無(wú)標(biāo)度特性指的是網(wǎng)絡(luò)的節(jié)點(diǎn)度(節(jié)點(diǎn)的度定義為與該節(jié)點(diǎn)相關(guān)聯(lián)的邊的條數(shù))分布服從冪律^[4]。

    復(fù)雜網(wǎng)絡(luò)上的同步現(xiàn)象是一種重要的研究課題^[5-7],。研究表明,，無(wú)標(biāo)度網(wǎng)絡(luò)和小世界網(wǎng)絡(luò)的同步能力好于規(guī)則網(wǎng)絡(luò)^[8]。對(duì)這一現(xiàn)象的普遍解釋是小世界網(wǎng)絡(luò)和無(wú)標(biāo)度網(wǎng)絡(luò)的平均最短距離小,，使得信息之間的交流更高效,。然而本文發(fā)現(xiàn)網(wǎng)絡(luò)的同步能力幾乎獨(dú)立于網(wǎng)絡(luò)的平均路徑長(zhǎng)度，即平均路徑相同的網(wǎng)絡(luò)同步能力卻存在較大的差異,。NISHIKAWA T,、MOTTER A E等的研究表明，對(duì)無(wú)標(biāo)度網(wǎng)絡(luò)而言,，網(wǎng)絡(luò)度分布的不均勻性抑制了網(wǎng)絡(luò)同步能力的增強(qiáng)^[9],；HONG H、KIM B J等的研究則表明,，對(duì)WS小世界網(wǎng)絡(luò)而言,，使節(jié)點(diǎn)度分布變得不均勻(即增大度分布的方差)增強(qiáng)了小世界網(wǎng)絡(luò)的同步能力^[10]，與文獻(xiàn)的結(jié)論貌似沖突,。本文對(duì)此給出了詳盡的解釋,。

    BARAHONA M和PECORA L M^[11]研究了線性耦合網(wǎng)絡(luò)的同步穩(wěn)定性問(wèn)題，給出了主穩(wěn)定函數(shù)判據(jù),。如果在網(wǎng)絡(luò)的每個(gè)節(jié)點(diǎn)上放置一個(gè)動(dòng)力學(xué)系統(tǒng)(該動(dòng)力學(xué)系統(tǒng)既可以是極限環(huán),，也可以是混沌的)，讓有邊相連的兩個(gè)節(jié)點(diǎn)的動(dòng)力學(xué)系統(tǒng)之間存在相互耦合作用,，就形成了一個(gè)動(dòng)力學(xué)網(wǎng)絡(luò),。設(shè)網(wǎng)絡(luò)有N個(gè)節(jié)點(diǎn),那么第i個(gè)節(jié)點(diǎn)所滿足的狀態(tài)方程是：

    WS小世界網(wǎng)絡(luò)的構(gòu)造算法為：(1)給定一個(gè)節(jié)點(diǎn)數(shù)為N的規(guī)則網(wǎng)絡(luò)，每個(gè)節(jié)點(diǎn)都與它最鄰近的K=2k個(gè)節(jié)點(diǎn)相連,，通常N≥K≥1,；(2)依次訪問(wèn)網(wǎng)絡(luò)中的每個(gè)節(jié)點(diǎn),，對(duì)每個(gè)節(jié)點(diǎn)的所有連線以概率p進(jìn)行重連，連接該節(jié)點(diǎn)的一端不變,，另一端隨機(jī)選擇節(jié)點(diǎn)連接,，但兩個(gè)節(jié)點(diǎn)之間不允許重復(fù)連接，也不允許節(jié)點(diǎn)與自身相連,。

    本文考慮了一種擇優(yōu)選擇的小世界網(wǎng)絡(luò)模型(Small-World Networks with Preferential Attachment,，PA小世界網(wǎng)絡(luò))。其構(gòu)造算法為：(1)步驟(1)與WS小世界網(wǎng)絡(luò)模型相同,；(2)在WS小世界網(wǎng)絡(luò)斷開重連時(shí)，不是隨機(jī)選擇一個(gè)節(jié)點(diǎn)連接,，而是按照正比與節(jié)點(diǎn)度的概率選擇節(jié)點(diǎn)連接,。

1 網(wǎng)絡(luò)幾何特征量

1.1 簇系數(shù)

    簇系數(shù)是表征網(wǎng)絡(luò)節(jié)點(diǎn)集群程度的物理量。節(jié)點(diǎn)i的簇系數(shù)定義為所有相鄰節(jié)點(diǎn)之間連邊的數(shù)目與可能的最大連邊數(shù)目的比例,。假設(shè)節(jié)點(diǎn)i的度為k_i,，它的所有鄰居節(jié)點(diǎn)之間的邊數(shù)為E_i，則節(jié)點(diǎn)i的簇系數(shù)定義為C_i=2E_i/(k_i(k_i-1)),。對(duì)網(wǎng)絡(luò)中所有節(jié)點(diǎn)的簇系數(shù)求平均值就是整個(gè)網(wǎng)絡(luò)的簇系數(shù),。

1.2 特征路徑長(zhǎng)度

    在網(wǎng)絡(luò)中，兩點(diǎn)間的距離定義為連接兩點(diǎn)的最短路徑所包含的邊的數(shù)目,，對(duì)所有節(jié)點(diǎn)對(duì)的距離求平均就得到了網(wǎng)絡(luò)的特征路徑長(zhǎng)度,。

1.3 節(jié)點(diǎn)度的方差

    節(jié)點(diǎn)度的方差是衡量節(jié)點(diǎn)度分布均勻性的一種度量。度分布越均勻,，方差越?。环粗?，方差越大,。其定義為：

1.4 節(jié)點(diǎn)的介數(shù)

    節(jié)點(diǎn)k的介數(shù)定義為：

式中，g_ij是節(jié)點(diǎn)i和節(jié)點(diǎn)j之間最短路徑的條數(shù),，g_ikj表示節(jié)點(diǎn)i和節(jié)點(diǎn)j之間最短路徑中經(jīng)過(guò)節(jié)點(diǎn)k的條數(shù),。對(duì)網(wǎng)絡(luò)中的所有節(jié)點(diǎn)對(duì)(i，j)求和即可以得到節(jié)點(diǎn)k的介數(shù),。由此可以看出,，一個(gè)節(jié)點(diǎn)的介數(shù)反映了該節(jié)點(diǎn)在其他節(jié)點(diǎn)通信中的重要程度。

2 仿真結(jié)果

    仿真中,，選擇網(wǎng)絡(luò)節(jié)點(diǎn)數(shù)N=500(仿真表明,，節(jié)點(diǎn)數(shù)的多少并不影響結(jié)果的普遍性)，K=2k=12,，即網(wǎng)絡(luò)的平均度<k>=12,。仿真圖的數(shù)據(jù)為節(jié)點(diǎn)數(shù)N=500時(shí)程序運(yùn)行50次的平均結(jié)果,。

    首先研究了特征根比例R隨重連概率p的變化情況。如圖1所示,，對(duì)于PA小世界網(wǎng)絡(luò),，隨著的p增大，R值先是單調(diào)減小然后單調(diào)增大,，表明該網(wǎng)絡(luò)的同步能力隨著p先增強(qiáng)后減弱,；對(duì)于WS小世界網(wǎng)絡(luò)，R單調(diào)遞減,，這表明WS小世界網(wǎng)絡(luò)的同步能力隨著p的增大而增強(qiáng),。那么發(fā)生這種差異的原因是什么呢？為此,，本文還研究了網(wǎng)絡(luò)的其他幾何特征量,，如網(wǎng)絡(luò)的特征路徑長(zhǎng)度l和簇系數(shù)C(如圖2所示)、節(jié)點(diǎn)度的方差σ²(如圖3所示),、節(jié)點(diǎn)的最大介數(shù)B_max(如圖4所示)以及節(jié)點(diǎn)的最大度數(shù)k_max(如圖5所示)等隨著p的變化情況,。

    由圖2可以看出，對(duì)于PA小世界網(wǎng)絡(luò)和WS小世界網(wǎng)絡(luò),，特征路徑長(zhǎng)度l和簇系數(shù)C都隨著p的增加而單調(diào)減小,，且?guī)缀醣３忠恢拢@說(shuō)明圖1中R值的差異不是由特征路徑長(zhǎng)度l和簇系數(shù)C的不同引起的,；同時(shí),，也表明對(duì)于特征路徑長(zhǎng)度和簇系數(shù)都幾乎相同的網(wǎng)絡(luò)，其同步性能仍然可以展現(xiàn)出相當(dāng)?shù)牟町?，也說(shuō)明了兩者中的任何一個(gè)都不是造成這種差異的原因,。

    圖3是對(duì)PA小世界網(wǎng)絡(luò)和WS小世界網(wǎng)絡(luò)的節(jié)點(diǎn)度的方差仿真的結(jié)果。顯而易見,，PA小世界網(wǎng)絡(luò)和WS小世界網(wǎng)絡(luò)的節(jié)點(diǎn)度的方差都隨著p的增大而增大,，但在PA小世界網(wǎng)絡(luò)中，由于重連時(shí)另一端為擇優(yōu)選擇,，導(dǎo)致節(jié)點(diǎn)度的方差增大的速度明顯快于WS小世界網(wǎng)絡(luò),。HONG H^[10]等人認(rèn)為小世界網(wǎng)絡(luò)中節(jié)點(diǎn)度的方差的增大會(huì)導(dǎo)致其同步能力的增強(qiáng)。本文的仿真也表明,，在WS小世界網(wǎng)絡(luò)中,，隨著p的增大，R逐漸減小,，節(jié)點(diǎn)度方差σ²逐漸增大,，那么是否節(jié)點(diǎn)度的方差的增大是導(dǎo)致網(wǎng)絡(luò)同步能力增強(qiáng)的原因呢？答案是否定的,?？梢钥吹?，在PA小世界網(wǎng)絡(luò)中，隨著重連概率p的增大,，方差也逐漸增大,，且增加的速率明顯大于WS小世界網(wǎng)絡(luò)，但由于節(jié)點(diǎn)度的方差增加的太快,，使得這種度分布的不均勻性對(duì)網(wǎng)絡(luò)的同步能力起了抑制作用(如在圖1中當(dāng)p≥0.7時(shí)R值的上揚(yáng),，表明了PA小世界網(wǎng)絡(luò)同步能力的下降)。從圖1可以看出,，在整個(gè)p值的范圍內(nèi),，PA小世界網(wǎng)絡(luò)的同步能力明顯不及WS小世界網(wǎng)絡(luò)(因?yàn)镽_PA>R_WS)。

    另由圖3可知,，在整個(gè)p值的范圍內(nèi),，PA小世界網(wǎng)絡(luò)的節(jié)點(diǎn)度的方差始終大于WS小世界網(wǎng)絡(luò)。因?yàn)閮煞N網(wǎng)絡(luò)的網(wǎng)絡(luò)特征路徑長(zhǎng)度l和簇系數(shù)C幾乎相等(由圖2可知),，由此可以推斷,，節(jié)點(diǎn)度的不均勻性實(shí)際上對(duì)網(wǎng)絡(luò)的同步能力起抑制作用,。由此,，在WS小世界網(wǎng)絡(luò)中，增強(qiáng)網(wǎng)絡(luò)同步能力的主要因素不是節(jié)點(diǎn)度方差逐漸增大(節(jié)點(diǎn)度方差只是由于小世界網(wǎng)絡(luò)中邊的重連引起的一種很自然且必然的結(jié)果),，網(wǎng)絡(luò)的特征路徑長(zhǎng)度下降才是最主要的原因(當(dāng)p值較小時(shí),，盡管度分布的不均勻性抑制了網(wǎng)絡(luò)的同步能力，但由于網(wǎng)絡(luò)的特征路徑長(zhǎng)度大幅下降,，網(wǎng)絡(luò)的同步能力還是增強(qiáng)了,；但是，當(dāng)網(wǎng)絡(luò)的特征路徑長(zhǎng)度下降幾乎為常數(shù)時(shí),，如果節(jié)點(diǎn)度分布的方差大于一定的臨界值,，其抑制作用就比較明顯了，如圖1和圖3所示),。所以,，短的特征路徑長(zhǎng)度和均勻的節(jié)點(diǎn)度分布綜合作用才能導(dǎo)致網(wǎng)絡(luò)同步能力的增強(qiáng)。此外,，NISHIKAWA T^[9]和HONG H^[10]均認(rèn)為網(wǎng)絡(luò)節(jié)點(diǎn)的最大介數(shù)B_max是表征網(wǎng)絡(luò)同步能力的一種相對(duì)合理的因素,。降低網(wǎng)絡(luò)節(jié)點(diǎn)的最大介數(shù)B_max將導(dǎo)致網(wǎng)絡(luò)同步能力的增強(qiáng)。圖4的仿真結(jié)果也說(shuō)明了這個(gè)問(wèn)題,?？梢钥吹剑赑A小世界網(wǎng)絡(luò)中,，B_max隨p先是單調(diào)下降然后單調(diào)上升,，與R隨p的變化趨勢(shì)相同,。所以，網(wǎng)絡(luò)中節(jié)點(diǎn)的最大介數(shù)B_max的確可以作為表征網(wǎng)絡(luò)同步能力的一種度量,。

    由圖4可知,，在整個(gè)p值范圍內(nèi)，PA小世界網(wǎng)絡(luò)中節(jié)點(diǎn)的最大介數(shù)始終大于WS小世界網(wǎng)絡(luò)中節(jié)點(diǎn)的最大介數(shù),；而圖1中,，在整個(gè)p值范圍內(nèi)，PA小世界網(wǎng)絡(luò)的同步能力始終不及WS小世界網(wǎng)絡(luò),，這從另外一個(gè)側(cè)面反映了最大介數(shù)表征網(wǎng)絡(luò)同步能力的合理性,。在無(wú)標(biāo)度網(wǎng)絡(luò)中，網(wǎng)絡(luò)中最大節(jié)點(diǎn)度k_max和最大介數(shù)B_max是正相關(guān)的,，即k_max越大,，通常B_max也越大；在WS小世界網(wǎng)絡(luò)中,，變化呈相反的趨勢(shì),，隨著k_max的增大，B_max單調(diào)下降,，如圖4和圖5所示,，在PA小世界網(wǎng)絡(luò)中，k_max和B_max呈現(xiàn)一種相對(duì)微妙復(fù)雜的關(guān)系,。

3 結(jié)論

    本文研究了一種擇優(yōu)選擇的小世界網(wǎng)絡(luò)模型,，基于這種模型，詳盡地研究了網(wǎng)絡(luò)的同步能力與網(wǎng)絡(luò)中的各種幾何特征量(如網(wǎng)絡(luò)平均路徑,、簇系數(shù),、度分布以及介數(shù)等)之間的關(guān)系。研究結(jié)果表明,，網(wǎng)絡(luò)的特征路徑長(zhǎng)度和簇系數(shù)都不能單獨(dú)決定網(wǎng)絡(luò)的同步能力,，因?yàn)椋趦烧叨枷嗟鹊那闆r下網(wǎng)絡(luò)仍然可以有很多不同的配置,，因而網(wǎng)絡(luò)的同步能力仍然可以存在較大的差異,。節(jié)點(diǎn)度分布的不均勻性從本質(zhì)上抑制了網(wǎng)絡(luò)的同步能力。仿真結(jié)果表明,，網(wǎng)絡(luò)中節(jié)點(diǎn)的最大介數(shù)是表征網(wǎng)絡(luò)同步能力的一種相對(duì)統(tǒng)一的物理量,。

參考文獻(xiàn)

[1] NEWMAN M E J.The structure and function of complex networks[J].SIAM Review，2003,，45(2)：167-256.

[2] WANG X F,，CHEN G R.Synchronization in small-world dynamical networks[J].Bifurcation and Chaos，2002，12(1)：187-192.

[3] LI Z,，CHEN G R.Global synchronization and asymptotic stability of complex dynamical networks[J].IEEE Transactions Circuit and System,，2006，53(1)：28-33.

[4] WANG X F, CHEN G R.Synchronization in scale-free dynamical networks：robustness and fragility[J].IEEE Transactions on Circuits and Systems I,，2002,，49(8)：54-62.

[5] ZHAO M，ZHOU T,，WANG B H,，et al.Relations between average distance, heterogeneity and network synchronizaility[J].Physica A，2006,，371(2)：773-780.

[6] WU X,，WANG B H，ZHOU T,，et al.Synchronizability of highly clustered scale-free networks[J].Chinese Physics Letters,，2006，23(4)：1046.

[7] 魏民,，王莉,，張強(qiáng).基于復(fù)雜網(wǎng)絡(luò)理論的典型電力電子電路復(fù)雜性研究[J].電子技術(shù)應(yīng)用，2011,，37(9)：118-121.

[8] MOTTER A E,，ZHOU C，KURTHS J.Network synchronization,，diffusion,，and the paradox of heterogeneity[J].Physical Review E,，2005,，71(1)：016116.

[9] NISHIKAWA T，MOTTER A E,，LAI Y C,，et al.Heterogeneity in oscillator networks: are smaller worlds easier to synchronize?[J].Physical Review Letters，2003,，91(7)：014010(1-4).

[10] HONG H,，KIM B J，CHOI M Y,，et al.Factors that predict better synchronizability on complex networks[J].Physical Review E,，2004，69(6)：067105(1-4).

[11] BARAHONA M,，PECORA L M.Synchronization in small-world systems[J].Physical Review Letters,，2002，89(5)：054101(1-4).

作者信息：

呂  瑛，李穎華

(西北工業(yè)大學(xué)明德學(xué)院 信息工程學(xué)院,，陜西 西安710124)