張曉磊，劉茂省

（中北大學 理學院，太原 030051）

如今，許多現(xiàn)實世界的系統(tǒng)可以用復雜動態(tài)網(wǎng)絡描述，如互聯(lián)網(wǎng)、神經(jīng)網(wǎng)絡，生物系統(tǒng)。研究并探索復雜網(wǎng)絡中出現(xiàn)的不同動態(tài)行為之間的相互作用是增強對其理解與運用的有效舉措［1］。在不同的動態(tài)行為中，網(wǎng)絡上個體應對傳染病引起的集體行為和傳染病傳播之間的關系已經(jīng)引起大多數(shù)人的關注。隨著傳染病的傳播，個體會收到來自政府、媒體或人與人之間等其他途徑得知的傳染病信息來適應性地改變各自的行為，即采取集體行為以避免被感染，如經(jīng)常用清水洗手，避免去擁擠地方等，以提高自我保護，這意味著在傳染病傳播的過程中，個體之間有關傳染病信息的傳播可以自發(fā)地誘發(fā)集體行為，這表明網(wǎng)絡上個體應對傳染病引起的集體行為和傳染病傳播，可以同時發(fā)生并自適應地相互作用。

在過去的幾年中，復雜動態(tài)網(wǎng)絡的同步問題受到廣泛關注，如在小世界［2］和無標度動態(tài)網(wǎng)絡［3］中，推導出了動態(tài)網(wǎng)絡完全同步的準則，其中指出網(wǎng)絡拓撲結構和節(jié)點動態(tài)之間的相互作用對于網(wǎng)絡實現(xiàn)完全同步十分重要。就傳染病傳播過程中的同步研究，李科贊等［4］集中在傳染病傳播如何影響個體的集體行為方面開展研究，提出了基于異質(zhì)平均場（HMF）理論的SIS（易感－感染－易感）和SIS（易感－感染－恢復）傳染病同步模型，得出傳染病的傳播會引發(fā)個體行為的變化，進而影響傳染病的演變。李科贊等［5］研究了在復雜網(wǎng)絡上個體行為的同步和傳染病傳播之間的相互作用，即不僅研究傳染病傳播對個體行為同步的影響，還研究個體行為同步對傳染病傳播的影響，更分析了傳染病傳播的控制問題。孫孟鋒等［6］在傳染病傳播過程中，引入了在個體之間傳染病信息傳播與傳染病傳播之間的耦合時滯，構建了3種不同（無時滯、耦合時滯和雙重時滯）的數(shù)學模型，并分別研究了在復雜網(wǎng)絡上個體行為的同步與傳染病傳播之間的相關性，分析了傳染病同步數(shù)學模型的局部和全局穩(wěn)定性。李俊民等［7］研究了具有未知周期時變耦合和隨機噪聲擾動的復雜動態(tài)網(wǎng)絡的同步問題，得到了復雜動態(tài)網(wǎng)絡同步的充分條件。

然而，目前對傳染病傳播過程中個體行為同步的研究大都是不受噪聲影響的復雜網(wǎng)絡模型，為了更好地評估復雜網(wǎng)絡中傳染病傳播與個體行為之間的相互作用，研究了在受隨機噪聲擾動影響的復雜網(wǎng)絡上，個體應對傳染病引起的集體行為與傳染病傳播之間存在的相互作用的問題。在受噪聲影響復雜網(wǎng)絡模型的基礎上，應用實際網(wǎng)絡中觀察到的自適應機制，通過向網(wǎng)絡上節(jié)點添加合適的控制器，構造傳染病同步的數(shù)學模型來研究受控噪聲影響復雜網(wǎng)絡的同步問題。

1 模型的建立

基于文獻［5］的復雜網(wǎng)絡模型（1），建立具有N個節(jié)點組成的受噪聲影響復雜網(wǎng)絡（4），

其中xi（t）＝［xi1（t），xi2（t），…，xin（t）］T∈Rn，代表第i個節(jié)點在t時刻的狀態(tài)變量，t∈［0，＋∞），f∶Rn×R→Rn是一個連續(xù)非線性向量值函數(shù)，描述節(jié)點的局部動態(tài)；ci（t）＞0表示耦合強度；Γ＝diag（γ1，γ2，…，γn）∈Rn×n代表內(nèi)部耦合矩陣，是一個正定對角矩陣；A＝（aij）N×N是網(wǎng)絡的鄰接矩陣，若節(jié)點i與節(jié)點j（i≠j）之間有連接，則aij＝aji＝1，否則aij＝aji＝0?；卩徑泳仃嘇，網(wǎng)絡的Laplacian矩陣L＝（lij）N×N可表示為

Laplacian矩陣L＝（lij）N×N的對角元素滿足以下等式：

這里kai表示節(jié)點i的度。假設L是不可約矩陣，意味網(wǎng)絡是強連通，沒有孤立集群。根據(jù)文獻［8］可知0是L矩陣的重數(shù)為1的最小特征值，并且其他特征值都嚴格為正。

復雜網(wǎng)絡模型（1）可用Laplacian矩陣表示為

具有N個節(jié)點組成的受噪聲影響復雜網(wǎng)絡，如下所示：

對于（4）與（1）（3）相比，除多了最后一項gi（x1，x2，…，xN）dωi（t）以外，其他符號表示一致，其中gi∈C（Rn×… ×Rn，Rn×n）是噪聲強度函數(shù)矩陣，ωi（t）＝（ωi1，ωi2，…，ωin）T∈Rn是n維向量Wiener過程，當i≠j時，假設 ωi（t）和 ωj（t）是相互獨立的過程。

對于一個孤立節(jié)點的解s（t）滿足

由于傳染病的信息傳遞具有時滯或其他特殊原因?qū)е聜魅静”┌l(fā)時，人們不能在短時間內(nèi)對其做出相應的反應，此時，對人們的行為進行一些指導，包括在醫(yī)院獲得治療和隔離，或經(jīng)常洗手和休息等，以達到抑止傳染病爆發(fā)的目的。這種行為的指導都可在反饋控制器下實現(xiàn)，本文中，通過在每個節(jié)點上增加一個反饋控制器ui（t）來控制個體的行為，以便傳染病在網(wǎng)絡中傳播期間實現(xiàn)所期望的行為同步。受控噪聲影響復雜網(wǎng)絡描述如下：

如果對于任何給定的初始狀態(tài)向量xi（0），受控噪聲影響復雜動態(tài)網(wǎng)絡（4）的解滿足以下條件，則稱其在均方意義下實現(xiàn)漸近同步，

SIS模型是一種典型的傳染病傳播模型，在這種傳染病傳播模型中每個節(jié)點處于易感或感染這2種狀態(tài)中的一種。感染節(jié)點以單位速率恢復為易感節(jié)點，易感節(jié)點以感染率為λ被其鄰居感染，變?yōu)楦腥竟?jié)點?？紤]復雜網(wǎng)絡上標準SIS模型，其中ρk（t）表示在t時刻度為k的感染節(jié)點的密度，則度為k的節(jié)點的演化方程可描述為

這里θ（t）為從一節(jié)點出發(fā)與其他節(jié)點隨機連邊，與感染節(jié)點相連的概率表易感節(jié)點連接到感染節(jié)點，它被感染的概率。在實際中個體會根據(jù)傳染病信息適應性地改變各自的行為（減少個體間接觸的頻率，并更經(jīng)常地采取集體保護措施）以避免被感染，因此這樣會導致感染率變化，用 （t）來量化這種影響，感染率λ變?yōu)棣?（t）。對于個體根據(jù)傳染病信息適應性地改變各自的行為以避免被感染的情況，可以看作是個體受到同步信息（即當傳染病傳播開來，政府、媒體發(fā)出公示或個體互相交流信息：如果采取某種行為可以預防或減少被感染的風險）進而采取同步行為的影響，同時 （t）＝（1－α）E（t）＋α，α∈（0，1）也可看作是準入率［9］，作為同步的信息，可以被認為是個人意識（或風險認知），即如果一個人在公共場合做出同步行為來達到感染控制的目的，這時由于一個人的反應有一種集體化，就可能會使這種同步行為更容易被接受，進而有更多的人去效仿，最終有利于感染得到控制，在整個過程中體現(xiàn)了個體由于受到同步的信息影響進而引起行為的同步。所有個體實現(xiàn)同步時，t→∞，E（t）→0，準入率 （t）達到最小值α。參數(shù)α值越小，相對應感染率λ （t）則變小，說明對集體行為的感知程度越高，當α＝1表示傳染病傳播不會受到同步的影響。

在制定具體的傳染病同步模型之前，做出以下基本假設［4－5］：

1）當傳染病開始傳播時，動態(tài)網(wǎng)絡中的個體之間存在弱線性耦合；

3）當集體保護行為顯著增加時，個人之間的保護信息交流將因達成保護協(xié)議而變得飽和。因此，耦合強度變化同步間的比例關系始終保持有效。

在以上假設的基礎上，用模型（8）描述傳染病傳播的特征，用模型（6）描述傳染病傳播過程中個體保護行為的演化過程，然后可以構造如下SIS傳染病同步模型：

2 受噪聲影響的復雜網(wǎng)絡上傳染病同步的穩(wěn)定性

定義1 系統(tǒng)（9）的同步流形可定義為

設Ω＝｛（S1，ρ1，…，Sd，ρd）∈R2d＋0≤Sk≤1，Sk＋ρk＝1，k＝1，…，d｝，Ω為系統(tǒng)（8）是正不變集，Ω0在Ω的內(nèi)部。

從上式可以看出θ＝0是該方程的一個平凡解，下面導出該方程存在正解0＜θ＜1的條件，為了分析θ（t）的性質(zhì)，構造輔助函數(shù)F（θ）：

經(jīng)計算可以得到，

從R0的推導可知當R0＞1，則傳染病模型（8）存來，將考慮無病平衡點E0的全局穩(wěn)定性。

其中Nk（ρ）＝－kλ （t）ρk（t）θ（t），則傳染病模型（8）可寫為： ρ＝Bρ＋N（ρ），

令B＝－I＋珟B，I表示單位矩陣，

已知Ω為傳染病模（8）是正不變集，且Ω0在Ω的內(nèi)部，所以只需要考慮解在Ω0的全局漸近穩(wěn)定，對傳染病模（8）的兩邊同時乘以kp（k），并求和得：

構造Lyapunov函數(shù)

沿著傳染病模型（8）求導，得

為得到定理2，接下來給出以下數(shù)學預備知識。

假設1 在受噪聲影響復雜網(wǎng)（4）中，假設存在li＞0，滿足

假設2 存在非負常數(shù)rij，i，j∈1，2，…，N，使得

定義2 （It formula）對于n維隨機微分系統(tǒng)

V（x（t），t）∈C2，1（Rn×R＋；R＋），這里C2，1（Rn×R＋；R＋）表示所有非負函數(shù)的簇，它們關于x是2次連續(xù)可微的，關于t是一次可微的。算子LV（x（t），t）被定義為

引理1［11］考慮系統(tǒng)（14）設q1、q2、q3為正數(shù)，假設存在函數(shù)C2，1（Rn×R＋；R＋），使得

那么平衡點x＝0在均方意義下是全局隨機漸近穩(wěn)定的。此外，對于任何∞＞t＞t0≥0，只要積分的期望值存在，則有

考慮受噪聲影響復雜網(wǎng)絡（4）的同步，為了實現(xiàn)同步目標（7），每個節(jié)點上增加的自適應控制器ui（t），（ki（t）是自適應參數(shù)）如下所示：

將控制器（17）應用于（6），通過ei（t）＝xi（t）－s（t），則受控噪聲影響復雜網(wǎng)絡的誤差系統(tǒng)可表示為

以下定理描述了受控網(wǎng)絡（6）全局漸近同步的一個充分條件。

定理2 當R0＞1，傳染病模型（8）的地方病平衡點E ，它在Ω0內(nèi)是全局漸近穩(wěn)定的，在假設1、2下系統(tǒng)（4），通過控制器（17）和自適應律（18）可保證受控網(wǎng)絡（6）在均方意義下全局漸近同步，也就是傳染病同步模型（9）的動態(tài)行為網(wǎng)絡的同步流形是全局漸近穩(wěn)定的。

證明：因為Si（t）＋ρi（t）＝1，i＝1，2，…，d，則傳染病模型（8）可表示為

定義如下矩陣：

對于V1（t），沿系統(tǒng)（8）的解對t求導，相似的分析過程可參考文獻［13］，可以得到對于函數(shù)V1（t）＝

對于

從上兩式和假設1、2下，可得

從引理1和式（24），可知

根據(jù)引理1，可知誤差動態(tài)系統(tǒng)（19）在均方意義下是全局漸近穩(wěn)定的，網(wǎng)絡（4）在均方意義下實現(xiàn)漸近同步。

同步流形S也是全局漸近穩(wěn)定的。

3 數(shù)值模擬

為了驗證上述結果，對SIS傳染病同步模型（9）進行數(shù)值研究，嵌入模型（9）中的網(wǎng)絡被看成大小為N＝200的BA無標度（優(yōu)先連接）網(wǎng)絡，這個網(wǎng)絡是從m0＝4的初始網(wǎng)絡演變而來的，我們給每一個新節(jié)點添加m＝3新邊。不失一般性，假設模型（9）中的f為混沌洛倫茲振蕩，盡管從實際傳染病傳播過程的角度來看，這一假設很難證明是正確的，只是將其用于數(shù)值模擬。這種振蕩可以描述為

在這里a1＝10，a2＝28，a3＝8／3，?？醋魇菃挝痪仃嚕渌麉?shù)為 δ＝0．001，α＝0．5，di＝0．01。噪聲強度函數(shù)矩陣gi（x）＝diag（xi1－xi＋1，1，xi2－xi＋1，2，xi3－xi＋1，3），其中xN＋1，3＝x11，i＝1，…，N。xi的狀態(tài)初值是從服從均勻分布的［0，1］中隨機選取的，初始耦合強度為ci＝0．001，初始感染密度為 ρ1＝ρ2＝0．01，ρi＝0，i＝1，…，d，初始自適應參數(shù)為ki＝di＝0．01。

圖1中（a）、（b）表示模型（9）的 λ分別等于0．2，0．8時，同步誤差E（t），感染密度 ρ（t）和耦合強度

圖1 同步誤差E（t），感染密度 ρ（t），耦合強度ci（t）分別在 λ＝0．2，0．8時的變化

從圖1（a）可知：當λ＝0．2時，傳染病不會爆發(fā)，沒有形成地方病，此時E（t）不等于0，意味著同步?jīng)]有實現(xiàn)，個體沒有表現(xiàn)出集體行為，但是從圖上可以看出當傳播過程中， （t）滿足以下條件α≤ （t）＜1， （t）＝α時，個體之間存在同步信息，此時傳染病感染密度收斂為0的速度比 （t）＝1時個體之間不存在同步信息要快，說明個體之間同步信息的存在會加速傳染病的滅絕。從圖1（b）可知：當λ＝0．8時， （t）＝1傳染病感染密度曲線表明個體之間由于不存在同步信息傳染病會爆發(fā)，但是由于模型中存在自適應控制器導致傳染病生成的規(guī)模不是很明顯，但此時可以明顯看到個體之間存在同步信息，傳染病感染密度收斂為0，同時E（t）等于0，意味著同步實現(xiàn)，同步信息的存在會使個體表現(xiàn)出集體行為，說明個體的集體行為可以抑制流行病傳播行為，反之，這種流行病傳播行為可以加速個體的集體行為。

圖2是當λ＝0．8時，自適應控制器ui的變化曲線和自適應參數(shù)ki的變化曲線。

圖3 （a）、（b）表示模型（9）的 λ分別等于1．2，2時，同步誤差E（t），感染密度 ρ（t）和耦合強度ci（t）的變化。

從圖3（a）可知：當λ＝1．2時，由于模型中存在自適應控制器ui的原因?qū)е聜魅静∩傻囊?guī)模不明顯，從圖3（b）可知：當λ＝2時，盡管模型中存在自適應控制器ui，若個體之間不存在同步信息，傳染病會爆發(fā)，說明模型中的自適應控制器只在感染率在一定范圍內(nèi)才有效，但此時可以明顯看到若個體之間存在同步信息，傳染病感染密度收斂為0，同時E（t）等于0，意味著同步實現(xiàn)，同步信息的存在會使個體表現(xiàn)出集體行為，從而抑制傳染病的傳播。

4 結論

研究了受噪聲影響復雜網(wǎng)絡上由傳染病動力學引起的個體同步行為，構建了能刻畫這類現(xiàn)象的傳染病同步數(shù)學模型，研究了模型的感染率與同步穩(wěn)定性之間的關系，得到了同步模型的全局穩(wěn)定性條件。數(shù)值模擬結果表明：集體行為可以抑制傳染病的傳播行為，反之，這種傳染病的傳播行為可以加速集體行為，這一結論與實際復雜網(wǎng)絡上傳染病的傳播特征吻合。因此，在研究復雜網(wǎng)絡上的傳染病同步問題時，本研究可為更好地理解和控制這類復雜網(wǎng)絡系統(tǒng)提供基本框架。