仇智鵬，魯富榮，杜亞星，錢宇華+

1.山西大學(xué) 大數(shù)據(jù)科學(xué)與產(chǎn)業(yè)研究院，太原 030006

2.山西大學(xué) 計算智能與中文信息處理教育部重點實驗室，太原 030006

3.山西大學(xué) 計算機(jī)與信息技術(shù)學(xué)院，太原 030006

1　引言

隨著對復(fù)雜網(wǎng)絡(luò)的深入研究，復(fù)雜網(wǎng)絡(luò)的結(jié)構(gòu)、演化以及動力學(xué)問題已經(jīng)獲得許多重要的突破[1-4]。對于復(fù)雜網(wǎng)絡(luò)動力學(xué)性質(zhì)的研究，其最終目標(biāo)是如何控制一個網(wǎng)絡(luò)以達(dá)到期待的目標(biāo)狀態(tài)。復(fù)雜網(wǎng)絡(luò)的可控性問題是近年來提出的重要研究熱點之一。復(fù)雜網(wǎng)絡(luò)控制方面主要的工作包括牽制控制[5-6]和結(jié)構(gòu)控制[3,7]等。近期的工作主要包括復(fù)雜網(wǎng)絡(luò)的可觀測性[8]、目標(biāo)可控性[9]、控制能量[10]、精確可控性[11]及多層網(wǎng)絡(luò)的結(jié)構(gòu)和動力學(xué)性質(zhì)[12]等方面的研究[13-19]。對于復(fù)雜網(wǎng)絡(luò)而言，控制的目的是如何通過恰當(dāng)?shù)妮斎朐谟邢迺r間內(nèi)使得網(wǎng)絡(luò)在自身動力學(xué)機(jī)制下演化到任意想要的狀態(tài)[20-21]。2011年，Liu等人[3]通過將復(fù)雜網(wǎng)絡(luò)的完全可控問題轉(zhuǎn)化為圖論中的最大匹配問題，找到了控制網(wǎng)絡(luò)的最少驅(qū)動節(jié)點數(shù)。2013年，Yuan等人[11]基于可控性的PBH秩判據(jù)提出了對于復(fù)雜網(wǎng)絡(luò)的精確控制，從而簡化了對于無向復(fù)雜網(wǎng)絡(luò)驅(qū)動節(jié)點數(shù)的計算。然而，對于真實網(wǎng)絡(luò)而言，對整個網(wǎng)絡(luò)進(jìn)行完全控制有時是費時費力的，甚至是不可行的，例如微生物網(wǎng)絡(luò)和一些社會網(wǎng)絡(luò)，它們規(guī)模龐大且結(jié)構(gòu)復(fù)雜。出于這個考慮，Gao等人[9]基于結(jié)構(gòu)控制理論，利用貪婪算法提出了對網(wǎng)絡(luò)的目標(biāo)控制。該算法首先選擇一部分節(jié)點作為需要控制的目標(biāo)節(jié)點，再通過迭代使用匹配算法最終找到了驅(qū)動該部分目標(biāo)節(jié)點達(dá)到理想狀態(tài)的最小驅(qū)動節(jié)點及其數(shù)目。

度相關(guān)性[22]是刻畫網(wǎng)絡(luò)屬性的一個重要指標(biāo)，反映節(jié)點與其相鄰節(jié)點之間度的聯(lián)系。本文考慮的是有向網(wǎng)絡(luò)的度相關(guān)性，具體為入度-入度相關(guān)性、入度-出度相關(guān)性、出度-入度相關(guān)性、出度-出度相關(guān)性。如果網(wǎng)絡(luò)中兩個節(jié)點的邊連接情況與兩個節(jié)點的度值無關(guān)，稱網(wǎng)絡(luò)不具有度相關(guān)性。如果度值較大的節(jié)點傾向于連接度值較大的節(jié)點，那么就稱這個網(wǎng)絡(luò)是同配的。反之，如果度值較大的節(jié)點傾向于連接度值相對較低的節(jié)點，那么就稱這個網(wǎng)絡(luò)是異配的。

由于復(fù)雜網(wǎng)絡(luò)節(jié)點數(shù)目多，連邊復(fù)雜度高，不能像低維小系統(tǒng)一樣通過遍歷的方法來尋找控制輸入位置。因此，如何有效地確定所需要獨立控制的節(jié)點數(shù)目和控制輸入的位置，以滿足復(fù)雜網(wǎng)絡(luò)可控性的要求就成為網(wǎng)絡(luò)控制需要解決的首要問題。對于控制而言，除了探究其控制的機(jī)理，提出有效的控制手段之外，探究影響該控制方法效率的因素也是一個重要的方面。度相關(guān)性對于復(fù)雜網(wǎng)絡(luò)結(jié)構(gòu)可控具有重要的影響已經(jīng)被證實[23]。在真實網(wǎng)絡(luò)中，例如電力網(wǎng)絡(luò)、生物網(wǎng)絡(luò)，其度相關(guān)性呈現(xiàn)出一定的異配性，而社會網(wǎng)絡(luò)呈現(xiàn)出一定的同配性，網(wǎng)絡(luò)的同配性與異配性會影響這些網(wǎng)絡(luò)控制的效率。而目標(biāo)控制是采用貪婪算法，在結(jié)構(gòu)控制的基礎(chǔ)上實現(xiàn)對于網(wǎng)絡(luò)部分節(jié)點的控制，因此本文選取度相關(guān)性作為目標(biāo)控制的影響因素進(jìn)行研究。探究復(fù)雜網(wǎng)絡(luò)目標(biāo)控制的影響因素，能夠更好地認(rèn)識網(wǎng)絡(luò)的拓?fù)浣Y(jié)構(gòu)對于網(wǎng)絡(luò)目標(biāo)控制的影響，同時也能夠在此基礎(chǔ)上改進(jìn)網(wǎng)絡(luò)使得網(wǎng)絡(luò)更易于被目標(biāo)控制。

本文探究了各種度相關(guān)性對于目標(biāo)控制效率的影響。實驗結(jié)果表明在目標(biāo)控制的機(jī)制下，出度-入度相關(guān)性對網(wǎng)絡(luò)的目標(biāo)控制效率影響較大，而隨著出度-入度相關(guān)性的逐步增加，驅(qū)動節(jié)點比例呈逐步下降趨勢。

2　相關(guān)內(nèi)容

2.1　目標(biāo)控制

本文采用的目標(biāo)控制策略是由Gao等人[9]提出的，研究的是線性時不變系統(tǒng)的目標(biāo)控制：

其中，x∈RN,u∈RM,y∈RS，這3個變量分別代表了系統(tǒng)的狀態(tài)向量、輸入向量以及輸出向量。A∈RN×N，B∈RN×M，C∈RS×N，分別表示狀態(tài)矩陣、輸入矩陣以及輸出矩陣。A描述的是系統(tǒng)的線性連接關(guān)系；B決定了外部輸入信號與網(wǎng)絡(luò)節(jié)點之間的連接關(guān)系；u是隨時間變化的輸入信號；C是所需要控制的目標(biāo)節(jié)點。對于一個具有N個節(jié)點的網(wǎng)絡(luò)N={1,2,…,N}，所需要控制的目標(biāo)節(jié)點集合是C={c1,c2,…,cS}，其中，S=|C|=fN。表示輸出的矩陣C=[I(c1),I(c2),…,I(cS)]，其中I(i)表示一個N×N單位矩陣的第i行。如果給定一個時變的輸入向量u(t)=(u1(t),u2(t),…,uM(t))T，這樣一個輸入向量可以使得目標(biāo)節(jié)點變成任意期望的最終狀態(tài)，那么就說這個系統(tǒng)是目標(biāo)可控的。目標(biāo)控制可以被看作一種（A,B,C）特殊的輸出可控，而且在滿足如下條件的情況下：

如果滿足（2），對于任意的初始狀態(tài)x(0)，就可以計算出最優(yōu)化的輸入向量u(t)使得網(wǎng)絡(luò)可以在有限時間tˉ>0達(dá)到任意最終狀態(tài)yˉ，也就是說最終y(tˉ)≡yˉ。

為了解決目標(biāo)控制問題，采用“K-walk theory”。該理論的原理是：一個節(jié)點可以控制這樣的一組目標(biāo)節(jié)點，即從驅(qū)動節(jié)點到任意目標(biāo)節(jié)點的路徑長度是互不相同的。使用“K-walk theory”就可以確定這樣的一個受控節(jié)點的集合。對于需要超過一個控制輸入的網(wǎng)絡(luò)，采用貪婪算法（圖1），來找到最少的輸入節(jié)點。圖1（a）中對一個小型網(wǎng)絡(luò)施加目標(biāo)控制。網(wǎng)絡(luò)中一共7個節(jié)點，其中節(jié)點{x3,x4,x6,x7}為目標(biāo)節(jié)點（紅圈標(biāo)出）。在圖1（b）中使用貪婪算法來解決網(wǎng)絡(luò)的目標(biāo)控制問題。首先，將原網(wǎng)絡(luò)映射成一個二分圖的形式。目標(biāo)節(jié)點是{x3,x4,x6,x7}，通過紅線所示的第一次迭代之后，找到匹配節(jié)點{x2,x3,x5,x6}，將這些節(jié)點作為第二次迭代的目標(biāo)節(jié)點。通過3次迭代，最終得出{x1,x5}是{x3,x4,x6,x7}的驅(qū)動節(jié)點。

Fig.1　Progress of target control圖1　目標(biāo)控制過程

2.2　度相關(guān)性

度相關(guān)性可以度量節(jié)點之間的連接趨勢，而且它可以有多種表示形式。本文采用皮爾遜相關(guān)系數(shù)刻畫度相關(guān)性：

2.3　度相關(guān)性影響

本文所研究的是度相關(guān)性對于網(wǎng)絡(luò)目標(biāo)控制的影響。如圖2所示，從一個簡單的角度考慮一個網(wǎng)絡(luò)。在圖2（a）中，該矩陣為隨機(jī)重連之前網(wǎng)絡(luò)的連接矩陣，左半部分代表源節(jié)點（即有向邊的出節(jié)點），右半部分代表終節(jié)點（即有向邊的指入點）中N=7。圖2（b）展示的是隨機(jī)重連之后網(wǎng)絡(luò)的連接矩陣。

Fig.2　Connection matrix of networks before and after change圖2　網(wǎng)絡(luò)重連前后的連接矩陣

Fig.3　Effect of degree correlations on target control圖3　度相關(guān)性對目標(biāo)控制驅(qū)動節(jié)點的影響

在圖3（a）～（d）中，所要控制的是{x3,x4,x6,x7}4個節(jié)點。圖（a）是初始網(wǎng)絡(luò)的連邊圖，圖（b）通過目標(biāo)控制算法找出的網(wǎng)絡(luò)最終驅(qū)動節(jié)點{x1,x5}，圖（c）是隨機(jī)重連之后生成的網(wǎng)絡(luò)連邊圖，圖（d）為改變度相關(guān)性之后，通過目標(biāo)控制算法找出的網(wǎng)絡(luò)最終驅(qū)動節(jié)點變?yōu)閧x5}，從而使得驅(qū)動節(jié)點的數(shù)目降低。因此，網(wǎng)絡(luò)的目標(biāo)可控性會受到度相關(guān)性的影響。

3　度相關(guān)性影響的實驗分析

3.1　模擬退火算法

為了探索度相關(guān)性對于復(fù)雜網(wǎng)絡(luò)目標(biāo)控制的影響，本文采用模擬退火算法，使得在保持度分布不變的情況下，改變度相關(guān)性到任意期望的目標(biāo)值。

（1）設(shè)定了初始溫度T和每個溫度下的迭代次數(shù)iter，以及能量函數(shù)E(r)=|r-r*|；

（2）采用隨機(jī)擾動的方法，每次從網(wǎng)絡(luò)中任取兩條邊，交換兩條連邊的終節(jié)點，該操作的目的是使得網(wǎng)絡(luò)中每個節(jié)點的出度和入度保持不變；

（3）計算出E(r)的值，看是否小于閾值，如果小于閾值，則已經(jīng)達(dá)到所需要的r值范圍；

（4）繼續(xù)迭代，直到達(dá)到所需要的r值范圍，或者溫度T小于所設(shè)定的最小溫度t。

3.2　度相關(guān)性對E-R隨機(jī)網(wǎng)絡(luò)的實驗分析

圖4為在E-R隨機(jī)網(wǎng)絡(luò)中度相關(guān)性對目標(biāo)控制的影響(nD=ND/(N×p))。其中，nD為目標(biāo)控制驅(qū)動節(jié)點比例，ND為驅(qū)動節(jié)點數(shù)目，N為節(jié)點數(shù)，p為選擇的目標(biāo)節(jié)點比例。本文研究的是度相關(guān)性對網(wǎng)絡(luò)目標(biāo)控制的影響，因此每次隨機(jī)選取一半的節(jié)點(p=50%)為目標(biāo)節(jié)點。該E-R隨機(jī)網(wǎng)絡(luò)的平均度為=1(red),=2(green),=3(blue),=5(black)。圖中每一個數(shù)據(jù)點為5次實驗取平均值。圖4（a）考慮了入度-入度相關(guān)性(rin-in)和驅(qū)動節(jié)點比例(nD)的關(guān)系；圖4（b）考慮了入度-出度相關(guān)性(rin-out)和驅(qū)動節(jié)點比例(nD)的關(guān)系；圖4（c）考慮了出度-入度相關(guān)性(rout-in)和驅(qū)動節(jié)點比例(nD)的關(guān)系；圖4（d）考慮了出度-出度相關(guān)性(rout-out)和驅(qū)動節(jié)點比例(nD)的關(guān)系。

實驗表明對于目標(biāo)控制而言，在4種度相關(guān)性中，出度-入度相關(guān)性的影響最為明顯。

對于E-R隨機(jī)網(wǎng)絡(luò)，首先從入度-入度相關(guān)性角度來看，在網(wǎng)絡(luò)處于異配的情況(rin-in<0)，網(wǎng)絡(luò)中入度較大的節(jié)點連接網(wǎng)絡(luò)中入度較小的節(jié)點，而該入度較小的節(jié)點連接網(wǎng)絡(luò)中入度較大的節(jié)點。從形成匹配的角度來看，該結(jié)構(gòu)會在形成二階路徑的過程中，浪費大量的入邊，從而抑制了網(wǎng)絡(luò)中其他匹配的形成，加大了網(wǎng)絡(luò)目標(biāo)可控的難度。

Fig.4　Correlations betweennDand all kinds ofrin random network圖4　隨機(jī)網(wǎng)絡(luò)中驅(qū)動節(jié)點比例和各種相關(guān)性的關(guān)系

因此，網(wǎng)絡(luò)異配性較強(qiáng)(r<-0.3)時，網(wǎng)絡(luò)的驅(qū)動節(jié)點比例較大。通過模擬退火算法改變網(wǎng)絡(luò)的度相關(guān)性后，網(wǎng)絡(luò)在逐步變?yōu)橥|(zhì)網(wǎng)絡(luò)(rin-in>0)的過程中，相關(guān)性逐步上升，由二階路徑所帶來的對于匹配的抑制性降低，所以網(wǎng)絡(luò)驅(qū)動節(jié)點的比例逐步下降。而當(dāng)網(wǎng)絡(luò)的入度-入度相關(guān)性逐步增加時，網(wǎng)絡(luò)中入度較高的節(jié)點相互連接，使得網(wǎng)絡(luò)中入邊存在較大的浪費，所以在度相關(guān)性增高的后半階段，目標(biāo)控制的驅(qū)動節(jié)點比例逐步上升。因此，隨著網(wǎng)絡(luò)入度-入度相關(guān)性的增加，網(wǎng)絡(luò)的驅(qū)動節(jié)點比例首先會逐步降低，在網(wǎng)絡(luò)入度-入度相關(guān)性增高的后半階段，網(wǎng)絡(luò)驅(qū)動節(jié)點比例逐步上升。

從入度-出度相關(guān)性的角度來看，在求網(wǎng)絡(luò)中目標(biāo)節(jié)點的驅(qū)動時，所考慮的是節(jié)點之間的匹配關(guān)系。以連邊兩端的節(jié)點考慮，對于源節(jié)點而言，考慮的是源節(jié)點的出度，而對于連邊另一端的節(jié)點，考慮的是其入度，以此衡量網(wǎng)絡(luò)中的匹配數(shù)目。這樣的一個匹配數(shù)目與網(wǎng)絡(luò)的入度-出度相關(guān)性幾乎無關(guān)，因此網(wǎng)絡(luò)的目標(biāo)可控性對于網(wǎng)絡(luò)的入度-出度相關(guān)性沒有體現(xiàn)出很強(qiáng)的依賴關(guān)系。

從出度-入度相關(guān)性的角度來看，因為有向網(wǎng)絡(luò)的控制問題可以轉(zhuǎn)化為匹配問題，而匹配本身從源節(jié)點的角度看就是出邊，從終節(jié)點的角度看就是入邊，所以網(wǎng)絡(luò)本身體現(xiàn)出對于該相關(guān)性很高的依賴性。當(dāng)網(wǎng)絡(luò)體現(xiàn)出很高的異配性時（網(wǎng)絡(luò)中節(jié)點的出度很高而入度很低時），因為該出度高的節(jié)點只能控制其中一個入度低的節(jié)點，所以想要控制其他入度低的節(jié)點就需要另外施加控制輸入，網(wǎng)絡(luò)的異配性對于網(wǎng)絡(luò)的目標(biāo)可控體現(xiàn)出較高的抑制性。因此，對于目標(biāo)控制而言，其驅(qū)動節(jié)點的比例會相對較高。反之，如果網(wǎng)絡(luò)有良好的同配性，則對于出度較低的節(jié)點，其入度也會較低，對于網(wǎng)絡(luò)中出邊或者入邊的利用率較高，在形成匹配的過程中，目標(biāo)節(jié)點易處于一條有向路徑中。因此，當(dāng)網(wǎng)絡(luò)同配性越來越高時，網(wǎng)絡(luò)的驅(qū)動節(jié)點比例會降低。

對于出度-出度相關(guān)性而言，網(wǎng)絡(luò)驅(qū)動節(jié)點的比例體現(xiàn)出和入度-入度相關(guān)性良好的吻合。因為如果將原網(wǎng)絡(luò)中所有的連邊反向，則網(wǎng)絡(luò)中的匹配性質(zhì)不會發(fā)生變化，即網(wǎng)絡(luò)的驅(qū)動節(jié)點比例不會發(fā)生變化，此時網(wǎng)絡(luò)的出度-出度相關(guān)性等價于原來的入度-入度相關(guān)性，所以網(wǎng)絡(luò)的驅(qū)動節(jié)點比例與入度-入度相關(guān)性圖所體現(xiàn)出來的規(guī)律一致。

3.3　度相關(guān)性對S-F無標(biāo)度網(wǎng)絡(luò)的實驗分析

如圖5所示，對于該無標(biāo)度網(wǎng)絡(luò)，其平均度=1,=2,=3,=5,γ=3,N=5 000。目標(biāo)節(jié)點的選擇也是每次隨機(jī)選擇總節(jié)點數(shù)目的一半(p=50%)。每一個數(shù)據(jù)點都是取了5次實驗的平均值。圖5表示了驅(qū)動節(jié)點比例和4種度相關(guān)性之間的關(guān)系。對于S-F無標(biāo)度網(wǎng)絡(luò)而言，度相關(guān)性對于其影響表現(xiàn)出和E-R隨機(jī)網(wǎng)絡(luò)類似的性質(zhì)。

Fig.5　Correlations betweennDand all kinds ofrin scale-free network圖5　無標(biāo)度網(wǎng)絡(luò)中驅(qū)動節(jié)點比例和各種相關(guān)性的關(guān)系

從入度-入度相關(guān)性和出度-出度相關(guān)性的角度來看，在目標(biāo)控制中，所體現(xiàn)的性質(zhì)與E-R隨機(jī)網(wǎng)絡(luò)相差不大。然而，因為網(wǎng)絡(luò)為S-F無標(biāo)度網(wǎng)絡(luò)，不考慮網(wǎng)絡(luò)中具體連邊的方向，網(wǎng)絡(luò)本身體現(xiàn)出較高的異配性，即使網(wǎng)絡(luò)的入度-入度相關(guān)性逐步增加，網(wǎng)絡(luò)本身也很難有較高的同配性，所以不容易構(gòu)成入度-入度相關(guān)性對于網(wǎng)絡(luò)形成匹配的抑制，從而網(wǎng)絡(luò)在入度-入度相關(guān)性由-1到+1的過程中都保持驅(qū)動節(jié)點比例逐步降低。而當(dāng)網(wǎng)絡(luò)的入度-入度相關(guān)性極高時，網(wǎng)絡(luò)中入度較高的節(jié)點相互連接，使得網(wǎng)絡(luò)中入邊存在較大的浪費，所以在度相關(guān)性接近+1的時候，目標(biāo)控制的驅(qū)動節(jié)點比例存在上升的趨勢。因此，隨著網(wǎng)絡(luò)入度-入度相關(guān)性的增加，網(wǎng)絡(luò)的驅(qū)動節(jié)點比例首先會逐步降低，在網(wǎng)絡(luò)入度-入度相關(guān)性極高時，網(wǎng)絡(luò)驅(qū)動節(jié)點比例存在逐步上升的趨勢。

從入度-出度和出度-入度相關(guān)性的角度來看，入度-出度相關(guān)性對于網(wǎng)絡(luò)目標(biāo)控制驅(qū)動節(jié)點比例幾乎沒有影響，而出度-入度相關(guān)性對S-F無標(biāo)度網(wǎng)絡(luò)影響較大，且隨著出度-入度相關(guān)性的增加，驅(qū)動節(jié)點比例逐步降低。

3.4　真實網(wǎng)絡(luò)

本節(jié)用真實網(wǎng)絡(luò)數(shù)據(jù)檢驗相關(guān)性對于網(wǎng)絡(luò)目標(biāo)控制效率的影響。所用的真實網(wǎng)絡(luò)的部分?jǐn)?shù)據(jù)如表1所示，表格內(nèi)從左到右的內(nèi)容依次為真實網(wǎng)絡(luò)的編號、類型、名字、節(jié)點數(shù)目以及連邊數(shù)目。

Table 1　Data of real networks表1　真實網(wǎng)絡(luò)的數(shù)據(jù)

一共考慮了7個真實網(wǎng)絡(luò)，依舊隨機(jī)選擇總節(jié)點的一半節(jié)點(p=50%)作為目標(biāo)節(jié)點，其中每個數(shù)據(jù)點重復(fù)5次實驗取平均值，所得出的關(guān)系如圖6所示。

對于選取的TRN-EC-Alon、TRN-EC-RDB64、TRNYeast-Alon而言，入度-入度相關(guān)性整體變化趨勢不明顯，但是在度相關(guān)性極高時驅(qū)動節(jié)點比例是增高的。而入度-出度相關(guān)性在該類型網(wǎng)絡(luò)上幾乎沒有影響。出度-入度相關(guān)性對于網(wǎng)絡(luò)驅(qū)動節(jié)點的比例影響較為明顯，隨著相關(guān)性的逐步增加，網(wǎng)絡(luò)的驅(qū)動節(jié)點比例逐步減少。由于真實網(wǎng)絡(luò)的限制，網(wǎng)絡(luò)在出度-出度相關(guān)性上無法有很大范圍的變化，然而在小范圍內(nèi)的變化是伴隨度相關(guān)性的增加先降低后增高。

對于S420、S838網(wǎng)絡(luò)而言，從入度-入度的角度看，在度相關(guān)性為負(fù)時網(wǎng)絡(luò)的驅(qū)動節(jié)點比例浮動不大，而當(dāng)度相關(guān)性取正值的時候網(wǎng)絡(luò)的驅(qū)動節(jié)點比例隨度相關(guān)性的增加而逐步上升。入度-出度相關(guān)性對于網(wǎng)絡(luò)驅(qū)動節(jié)點比例的影響不大。而伴隨網(wǎng)絡(luò)出度-入度相關(guān)性的增加，網(wǎng)絡(luò)的驅(qū)動節(jié)點比例迅速下降，對于出度-出度相關(guān)性而言，網(wǎng)絡(luò)體現(xiàn)出和入度-入度相關(guān)性良好的匹配，在負(fù)值時網(wǎng)絡(luò)的驅(qū)動節(jié)點比例變化不大，當(dāng)網(wǎng)絡(luò)度相關(guān)性逐步變?yōu)?1的過程中，網(wǎng)絡(luò)驅(qū)動節(jié)點的比例逐步上升。

對于Dolphins、Prisoninmate網(wǎng)絡(luò)而言，伴隨著網(wǎng)絡(luò)入度-入度相關(guān)性的增加，網(wǎng)絡(luò)的驅(qū)動節(jié)點比例先是逐步下降，然后逐步上升。入度-出度相關(guān)性對其幾乎沒有影響，伴隨著出度-入度相關(guān)性的增加，網(wǎng)絡(luò)的驅(qū)動節(jié)點比例逐步降低。出度-出度相關(guān)性與入度-入度相關(guān)性的表現(xiàn)幾乎一致。

4　結(jié)論

通過對E-R隨機(jī)網(wǎng)絡(luò)、S-F無標(biāo)度網(wǎng)絡(luò)和真實網(wǎng)絡(luò)的實驗數(shù)據(jù)分析，本文最終能得出各種度相關(guān)性對網(wǎng)絡(luò)目標(biāo)控制效率影響的結(jié)論。

（1）對于絕大部分網(wǎng)絡(luò)而言，在入度-入度相關(guān)性由-1變?yōu)?1的過程中，網(wǎng)絡(luò)的驅(qū)動節(jié)點比例先是逐步降低，在度相關(guān)性為0的附近，網(wǎng)絡(luò)的驅(qū)動節(jié)點比例取得最低，并且在度相關(guān)性很小范圍內(nèi)幾乎保持不變。隨后，當(dāng)網(wǎng)絡(luò)的度相關(guān)性增大，網(wǎng)絡(luò)的驅(qū)動節(jié)點比例也隨之上升。

（2）入度-入度相關(guān)性，對于網(wǎng)絡(luò)目標(biāo)控制的驅(qū)動節(jié)點比例幾乎沒有影響。

（3）出度-入度相關(guān)性，伴隨著出度-入度相關(guān)性的增大，網(wǎng)絡(luò)的驅(qū)動節(jié)點比例逐步降低，且該相關(guān)性對于網(wǎng)絡(luò)驅(qū)動節(jié)點比例的影響最大。

（4）出度-出度相關(guān)性，該相關(guān)性對于網(wǎng)絡(luò)目標(biāo)控制驅(qū)動節(jié)點的比例影響與入度-入度相關(guān)性幾乎一致，都是隨著網(wǎng)絡(luò)度相關(guān)性的增大，驅(qū)動節(jié)點比例先降低后升高。

Fig.6　Correlations betweennDand all kinds ofrin real networks圖6　真實網(wǎng)絡(luò)中驅(qū)動節(jié)點比例和各種相關(guān)性的關(guān)系

5　展望

本文研究了在復(fù)雜網(wǎng)絡(luò)中，4種度相關(guān)性對于目標(biāo)控制的影響，其中網(wǎng)絡(luò)的出度-入度相關(guān)性對于網(wǎng)絡(luò)的目標(biāo)可控性具有較強(qiáng)的影響，并揭示了對于復(fù)雜網(wǎng)絡(luò)目標(biāo)控制的影響因素，及其帶來的一些問題。接下來的研究方向主要包括：如何通過調(diào)整網(wǎng)絡(luò)的結(jié)構(gòu)使得在目標(biāo)控制的情況下，網(wǎng)絡(luò)的驅(qū)動節(jié)點數(shù)目盡可能得少；探究各種度相關(guān)性影響的內(nèi)在機(jī)理；對于真實網(wǎng)絡(luò)，如何調(diào)整才能使其更易被控制等。

[1]Newman M,BarabásiAL,Watts D J.The structure and dynamics of networks[M].Princeton:Princeton University Press,2011.

[2]Boccaletti S,Latora V,Moreno Y,et al.Complex networks:structure and dynamics[J].Physics Reports,2006,424(4):175-308.

[3]Liu Yangyu,Slotine J J E,Barabási A L,et al.Controllability of complex networks[J].Nature,2011,473(7346):167-173.

[4]Barabási A L,Albert R.Emergence of scaling in random networks[J].Science,1999,286(5439):509-512.

[5]Wang Xiaofan,Chen Guanrong.Pinning control of scale-free dynamical networks[J].Physica A:Statistical Mechanics and itsApplications,2002,310(3):521-531.

[6]Manaffam S,Seyedi A.Pinning control for complex networks of linearly coupled oscillators[C]//Proceedings of the 2013 American Control Conference,Washington,Jun 17-19,2013.Piscataway:IEEE,2013:6364-6369.

[7]Cowan N J,Chastain E,Vilhena D A,et al.Nodal dynamics,not degree distributions,determine the structural controllability of complex networks[J].PLOS ONE,2011,7(6):e38398.

[8]Liu Yangyu,Slotine J E,Barabási A L,et al.Observability of complex systems[J].Proceedings of the National Academy of Sciences of the United States of America,2013,110(7):2460-2465.

[9]Gao Jianxi,Liu Yangyu,D'Souza R M,et al.Target control of complex networks[J].Nature Communications,2014,5:5415.

[10]Yan Gang,Ren Jie,Lai Yingcheng,et al.Controlling complex networks:how much energy is needed?[J].Physical Review Letters,2012,108(21):218703.

[11]Yuan Zhengzhong,Zhao Chen,Di Zengru,et al.Exact controllability of complex networks[J].Nature Communications,2013,4:2447.

[12]Boccaletti S,Bianconi G,Criado R,et al.The structure and dynamics of multilayer networks[J].Physics Reports,2014,544(1):1-122.

[13]Gao Jianxi,Barzel B,Barabási A L.Universal resilience patterns in complex networks[J].Nature,2016,530(7615):307-312.

[14]Buldyrev S V,Parshani R,Paul G,et al.Catastrophic cascade of failures in interdependent networks[J].Nature,2010,464(7291):1025-1028.

[15]Liu Xueming,Stanley H E,Gao Jianxi,et al.Breakdown of interdependent directed networks[J].Proceedings of the National Academy of Sciences of the United States of America,2016,113(5):1138-1143.

[16]Zhang Xiyun,Boccaletti S,Guan Shuguang,et al.Explosive synchronization in adaptive and multilayer networks[J].Physical Review Letters,2014,114(3):038701.

[17]Menichetti G,Dall'Asta L,Bianconi G,et al.Network controllability is determined by the density of low in-degree and out-degree nodes[J].Physical Review Letters,2014,113(7):078701.

[18]D?rfler F,Chertkov M,Bullo F,et al.Synchronization in complex oscillator networks and smart grids[J].Proceedings of the National Academy of Sciences of the United States ofAmerica,2013,110(6):2005-2010.

[19]Zhang Y,Garas A,Schweitzer F,et al.Value of peripheral nodes in controlling multilayer scale-free networks[J].Physical Review E,2016,93(1):012309.

[20]Slotine J J E,Li Weiping.Applied nonlinear control[M].Upper Saddle River:Prentice Hall,1991.

[21]Haussmann U.Linear systems and optimal control[J].Siam Review,2012,31(4):696-698.

[22]Foster J G,Foster D V,Grassberger P,et al.Edge direction and the structure of networks[J].Proceedings of the National Academy of Sciences of the United States of America,2010,107(24):10815-10820.

[23]Pósfai M,Liu Yangyu,Slotine J J E,et al.Effect of correlations on network controllability[J].Scientific Reports,2013,3:1067.