徐云程，胡 華，孫小軍

(1.寧夏大學數(shù)學統(tǒng)計學院，銀川 750021；2.寶雞文理學院數(shù)學與信息科學學院，陜西 寶雞 721013)

0 引言

其中

1 三層無標度關(guān)聯(lián)網(wǎng)絡(luò)協(xié)同傳播模型

1.1 三層關(guān)聯(lián)網(wǎng)絡(luò)結(jié)構(gòu)圖

首先，我們建立一個三層關(guān)聯(lián)網(wǎng)絡(luò)結(jié)構(gòu)圖，整個系統(tǒng)由3個子網(wǎng)A，B，C構(gòu)成，假設(shè)其分別含有N，M，H個節(jié)點。子網(wǎng)A到B、B到C之間有交叉連接，子網(wǎng)A到C之間沒有連接，且子網(wǎng)之間的連接是隨機的，記從A到B、B到A、B到C、C到B的層間網(wǎng)絡(luò)分別為AB，BA，BC，CB。具有3個子網(wǎng)的關(guān)聯(lián)網(wǎng)絡(luò)結(jié)構(gòu)圖如圖1所示。

圖1 三層關(guān)聯(lián)網(wǎng)絡(luò)

1.2 相關(guān)知識

在建立模型之前，對相關(guān)符號做以下說明(見表1)。

表1 三層關(guān)聯(lián)網(wǎng)絡(luò)符號說明

續(xù)表1

根據(jù)表1，可得3個子網(wǎng)A，B，C的易感節(jié)點總數(shù)、染病節(jié)點總數(shù)、節(jié)點總數(shù)分別為

聯(lián)合度分布分別為

邊界度分布分別為

平均度分別為

1.3 三層無標度關(guān)聯(lián)網(wǎng)絡(luò)協(xié)同傳播模型

由于無標度網(wǎng)絡(luò)給出的度分布信息是針對大部分節(jié)點的，并沒有準確描述出單個節(jié)點的情況，因此，基于單個節(jié)點去研究無標度網(wǎng)絡(luò)上流行病的動力學行為是很有現(xiàn)實意義的。本文用SIS倉室模型來描述3個子網(wǎng)上的動力學行為，在整個系統(tǒng)上，節(jié)點可以有兩種狀態(tài)：易感態(tài)(S)、感染態(tài)(I)。

接下來，對整個系統(tǒng)的動力學過程作如下假設(shè)：子網(wǎng)A中的一個易感節(jié)點通過一條層內(nèi)連邊被子網(wǎng)A中的一個染病節(jié)點感染的概率為λaa，通過一條層間連邊被子網(wǎng)B中的一個染病節(jié)點感染的概率為λba，對于非相關(guān)性無標度網(wǎng)絡(luò)，子網(wǎng)A中任意一條給定的邊與子網(wǎng)A中一個染病節(jié)點相連的概率為

(1)

子網(wǎng)A中任意一條給定的邊與子網(wǎng)B中一個染病節(jié)點相連的概率為

(2)

子網(wǎng)B中的一個易感節(jié)點通過一條層內(nèi)連邊被子網(wǎng)B中的一個染病節(jié)點感染的概率為λbb，通過一條層間連邊被子網(wǎng)A(C)中的一個染病節(jié)點感染的概率為λab(λcb)，子網(wǎng)B中任意一條給定的邊與子網(wǎng)A中一個染病節(jié)點相連的概率為

(3)

子網(wǎng)B中任意一條給定的邊與子網(wǎng)B中一個染病節(jié)點相連的概率為

(4)

子網(wǎng)B中任意一條給定的邊與子網(wǎng)C中一個染病節(jié)點相連的概率為

(5)

子網(wǎng)C中的一個易感節(jié)點通過一條層內(nèi)連邊被子網(wǎng)C中的一個染病節(jié)點感染的概率為λcc，通過一條層間連邊被子網(wǎng)B中的一個染病節(jié)點感染的概率為λbc，同樣可定義子網(wǎng)C中任意一條給定的邊與子網(wǎng)B中一個染病節(jié)點相連的概率為

(6)

子網(wǎng)C中任意一條給定的邊與子網(wǎng)C中一個染病節(jié)點相連的概率為

(7)

同時，子網(wǎng)A、B、C中的染病節(jié)點也會分別以恢復率μa，μb，μc變?yōu)橐赘泄?jié)點。

(8)

其中，e=0,1,…,naa;f=0,1,…,nba;g=0,1,…,nab;h=0,1,…,nbb;l=0,1,…,ncb;u=0,1,…,nbc;v=0,1,…,ncc。

2 閾值分析

2.1 全局傳播閾值

圖2 感染規(guī)模運動方向圖

(9)

方程(9)表示整個系統(tǒng)已經(jīng)進入了穩(wěn)定狀態(tài)，所以λaa,λbb,λcc將非常逼近λd，故可以用λd替換λaa,λbb,λcc，得:

(10)

(11)

基于式(11)將式(10)改寫為矩陣形式，有

(12)

即

(13)

令

2.2 孤立子網(wǎng)的傳播閾值

考慮各孤立子網(wǎng)上的情況，即僅考慮層內(nèi)傳播，不考慮層間傳播。孤立子網(wǎng)A,B,C的動力學模型分別為

(14)

在穩(wěn)定狀態(tài)，基于式(11)將上述模型寫為矩陣形式有

(15)

(16)

其中，Λmax(DA),Λmax(DB),Λmax(DC)分別為孤立子網(wǎng)A，B，C對應(yīng)的鄰接矩陣DA，DB，DC的最大特征值。

2.3 結(jié)論

為對比異質(zhì)網(wǎng)絡(luò)耦合的三層串狀關(guān)聯(lián)網(wǎng)絡(luò)的全局傳播閾值與各孤立子網(wǎng)傳播閾值的大小關(guān)系，這里將分兩種情形進行討論。

1)各孤立子網(wǎng)的恢復率均相等，即μa=μb=μc。

引理1(柯西交錯定理)[15]設(shè)H是一個n×n的對稱矩陣，而M是H的一個n-1主子陣，如果λ1≥λ2≥…≥λn和μ1≥μ2≥…≥μn-1分別是矩陣H和M的特征值，那么λ1≥μ1≥λ2≥…≥λn-1≥μn-1≥λn。

Λmax(P)≥Λmax(DA),Λmax(P)≥Λmax(DB),Λmax(P)≥Λmax(DC)，

2)各孤立子網(wǎng)的恢復率至少有兩個不相等，即μa,μb,μc中至少有兩個不相等。

證明：由式(12)有

λdDAρA+DBAρB-μaIρA=0

DABρA+λdDBρB+DCBρC-μbIρB=0

DBCρB+λdDCρC-μcIρC=0

(17)

基于擾動理論對各孤立子網(wǎng)的閾值進行研究，對傳播閾值及感染密度施加擾動(18)：

(18)

將式(18)代入式(17)，同時忽略二階項，有：

(19)

3 數(shù)值仿真

圖3和圖4分別展示了關(guān)聯(lián)網(wǎng)絡(luò)及對應(yīng)孤立網(wǎng)絡(luò)平均最終感染規(guī)模的三維切平面等高線圖，兩幅圖中最外側(cè)兩條曲線之間的區(qū)域代表疾病消亡，其余部分代表疾病變成了地方病將在人群中永遠流行。顯然圖3最外側(cè)兩條曲線之間區(qū)域的面積比圖4小，故關(guān)聯(lián)網(wǎng)絡(luò)全局閾值比孤立網(wǎng)絡(luò)閾值要小[16]，即協(xié)同傳播減小了傳播閾值，促進了疾病的傳播。

圖3 關(guān)聯(lián)網(wǎng)絡(luò)平均最終感染規(guī)模

圖4 孤立網(wǎng)絡(luò)平均最終感染規(guī)模

4 結(jié)論

本文基于異質(zhì)平均場理論構(gòu)建了三層無標度關(guān)聯(lián)網(wǎng)絡(luò)協(xié)同傳播模型，利用擾動理論對比分析了各子網(wǎng)的傳播閾值與全局閾值的相對大小。通過嚴格的數(shù)理方法分析得到：對三層無標度關(guān)聯(lián)網(wǎng)絡(luò)上的協(xié)同傳播而言，當各子網(wǎng)的恢復率均相等時，全局傳播閾值小于等于單個子網(wǎng)的傳播閾值；當各子網(wǎng)的恢復率至少有兩個不相等時，全局傳播閾值只能小于單個子網(wǎng)的傳播閾值。總之，協(xié)同傳播會使閾值減小，也即協(xié)同傳播促進了疾病的傳播。

文中所建模型不僅可以分析部分媒介傳染病的傳播，還可以用來分析艾滋病在男人、女人和兒童之間的傳播。所得結(jié)論有助于理解媒介傳染病的傳播機理、解釋某些流行病快速威脅人類的原因，且理論結(jié)果將為衛(wèi)生部門更好地控制和預防流行病提供科學依據(jù)。事實上，在現(xiàn)實生活中也存在部分能在人與家禽之間傳播的媒介傳染病，同時，媒體報道在疾病的預防和控制中也發(fā)揮著至關(guān)重要的作用，下一步將對受媒體報道影響的三層環(huán)狀無標度關(guān)聯(lián)網(wǎng)絡(luò)協(xié)同傳播模型進行研究。