均勻網(wǎng)絡(luò)上SIS傳染病模型幾種對(duì)逼近方法的比較

李 星，羅曉峰

(山西大學(xué) a.復(fù)雜系統(tǒng)研究所； b.數(shù)學(xué)科學(xué)學(xué)院；c.計(jì)算機(jī)與信息技術(shù)學(xué)院，太原 030006)

?

均勻網(wǎng)絡(luò)上SIS傳染病模型幾種對(duì)逼近方法的比較

李 星a,b，羅曉峰a,c

(山西大學(xué) a.復(fù)雜系統(tǒng)研究所； b.數(shù)學(xué)科學(xué)學(xué)院；c.計(jì)算機(jī)與信息技術(shù)學(xué)院，太原 030006)

網(wǎng)絡(luò)傳染病模型的精度是評(píng)價(jià)模型優(yōu)劣的標(biāo)準(zhǔn)，其中網(wǎng)絡(luò)對(duì)逼近傳染病模型的優(yōu)劣取決于所選逼近方法的精度。借助SIS傳染病模型比較均勻網(wǎng)絡(luò)上染病鄰居數(shù)服從泊松分布、多項(xiàng)分布以及基于平均場(chǎng)思想的3種逼近方法的精度，發(fā)現(xiàn)在泊松分布下模型的誤差最小，即泊松分布下的逼近方法精度最高。

均勻網(wǎng)絡(luò)；對(duì)逼近；傳染病；基本再生數(shù)；隨機(jī)模擬

傳染性疾病一直影響著人類(lèi)的生活，有的甚至威脅著人們的生命，如非典(SARS)、艾滋病(HIV)、埃博拉等，因此對(duì)傳染病傳播的研究顯得尤為重要[1]。由于在人群中進(jìn)行疾病實(shí)驗(yàn)具有不可行性，因此數(shù)學(xué)模型可作為強(qiáng)大的工具來(lái)模擬實(shí)際疾病在人群中的傳播機(jī)理，進(jìn)而預(yù)測(cè)疾病的發(fā)展趨勢(shì)，做出預(yù)防和控制策略。所以，傳染病數(shù)學(xué)模型的精度成了能否精準(zhǔn)預(yù)測(cè)疾病發(fā)展趨勢(shì)的關(guān)鍵[2-4]。另外，由于疾病在人群中傳播，人與人之間的接觸關(guān)系(人群結(jié)構(gòu))對(duì)疾病的傳播也起著至關(guān)重要的作用。隨著網(wǎng)絡(luò)科學(xué)的發(fā)展，研究者把人群映射到網(wǎng)絡(luò)上，通過(guò)網(wǎng)絡(luò)拓?fù)浣Y(jié)構(gòu)體現(xiàn)人群內(nèi)部的結(jié)構(gòu)，進(jìn)而建立基于網(wǎng)絡(luò)的傳染病模型來(lái)研究疾病在人群中的傳播[5-8]。在網(wǎng)絡(luò)中，節(jié)點(diǎn)代表人群中的個(gè)體，節(jié)點(diǎn)之間的連邊代表人與人之間的接觸關(guān)系。網(wǎng)絡(luò)對(duì)逼近模型是網(wǎng)絡(luò)傳染病模型的一種，它以節(jié)點(diǎn)和邊作為變量建立傳染病動(dòng)力學(xué)模型，進(jìn)而研究疾病的傳播機(jī)理[9-10]。這類(lèi)模型的精度取決于所選逼近方法的精度，因此本文主要借助均勻網(wǎng)絡(luò)(個(gè)體之間接觸關(guān)系的差異性不大)上的SIS傳染病模型比較了3種逼近方法(即染病個(gè)體數(shù)服從泊松分布、多項(xiàng)分布以及基于平均場(chǎng)理論)的精度。

1 均勻網(wǎng)絡(luò)模型

對(duì)于SIS對(duì)逼近傳染病模型,設(shè)個(gè)體總數(shù)為N且保持不變，即網(wǎng)絡(luò)節(jié)點(diǎn)總數(shù)不變，個(gè)體按狀態(tài)分為易感者S和染病者I，染病者恢復(fù)后變?yōu)橐赘姓撸粋€(gè)染病者和一個(gè)易感者每次接觸傳染的概率為τ，染病者的恢復(fù)率為γ。對(duì)于[S]、[I]、[SI]、[SS]、[II]，利用主方程可得

式中：[A]表示狀態(tài)為A的總節(jié)點(diǎn)數(shù)；[AB]表示狀態(tài)為A和狀態(tài)為B的節(jié)點(diǎn)構(gòu)成的所有二元組(對(duì))數(shù)量；[ABC]表示狀態(tài)為A、狀態(tài)為B和狀態(tài)為C的節(jié)點(diǎn)構(gòu)成的所有三元組數(shù)量。假設(shè)網(wǎng)絡(luò)的平均度為n1，顯然有如下保持量：

(1)

接下來(lái)通過(guò)不同的逼近方法封閉模型進(jìn)行研究。

1.1 泊松分布下的SIS對(duì)逼近模型(P-PW)

(2)

1.2 多項(xiàng)分布下的SIS對(duì)逼近模型(B-PW)

(3)

1.3 平均場(chǎng)SIS對(duì)逼近模型(MF-PW)

Kiss等在文獻(xiàn)[11]中應(yīng)用平均場(chǎng)的想法給出三元組的逼近公式：

可得到相應(yīng)的平均場(chǎng)SIS對(duì)逼近傳染病模型(MF-PW)：

(4)

接下來(lái)利用無(wú)病平衡點(diǎn)的穩(wěn)定性計(jì)算該模型的基本再生數(shù)。

令系統(tǒng)(4)的右邊3式為0，得到系統(tǒng)的無(wú)病平衡點(diǎn)Q0(N,0,n1N)。該系統(tǒng)在無(wú)病平衡點(diǎn)Q0處的Jacobian矩陣為

(5)

因此， Jacobian矩陣(5)的特征方程為|λE-J0|=0， 即

(6)

顯然，特征方程(6)有一個(gè)負(fù)實(shí)根-γ，當(dāng)且僅當(dāng)

的所有特征根有負(fù)實(shí)部時(shí)，無(wú)病平衡點(diǎn)Q0是局部漸近穩(wěn)定的。由于h(λ)的判別式為

因此，不妨設(shè)h(λ)的特征根分別為λ1、λ2。由圓盤(pán)定理可知，兩特征根必須有負(fù)實(shí)部滿足以下條件：

(7)

下面將針對(duì)各模型及其基本再生數(shù)進(jìn)行隨機(jī)和數(shù)值模擬對(duì)比，并進(jìn)行誤差分析，進(jìn)而比較3種逼近方法的精度。

2 模擬與誤差分析

這部分在ER隨機(jī)網(wǎng)絡(luò)(均勻網(wǎng)絡(luò))[12]上通過(guò)3種模型(式(2)(3)和(4))進(jìn)行數(shù)值與隨機(jī)模擬(100次求平均)比較3個(gè)模型的準(zhǔn)確性。誤差分析結(jié)果表明模型(2)誤差最小，即染病鄰居服從泊松分布的逼近方法更合理。

首先，構(gòu)造一個(gè)總節(jié)點(diǎn)數(shù)為N=5 000的度分布服從泊松分布pk=P(λ)的配置網(wǎng)絡(luò)，其中,λ=20。該網(wǎng)絡(luò)的總邊數(shù)為n1N=1.000 0×105，平均度n1=20.0。

圖1 在平均度為20、總節(jié)點(diǎn)數(shù)為5 000的ER隨機(jī)網(wǎng)絡(luò)下[SI]的時(shí)間序列

圖2 ER隨機(jī)網(wǎng)絡(luò)時(shí)誤差隨時(shí)間的變化情況

圖3 ER隨機(jī)網(wǎng)絡(luò)R0和隨τ，γ的變化情況

3 結(jié)束語(yǔ)

疾病在人群中的傳播不僅與疾病本身的傳播機(jī)理有關(guān)，且取決于人群個(gè)體的分布。網(wǎng)絡(luò)傳染病模型就是基于此，把人群的分布映射到網(wǎng)絡(luò)上，然后讓疾病在網(wǎng)絡(luò)上傳播。模型的高準(zhǔn)確性是建模者所追求的，既有利于數(shù)學(xué)分析又能應(yīng)用到實(shí)際疾病研究中。本文研究了均勻網(wǎng)絡(luò)考慮染病鄰居數(shù)服從泊松分布、多項(xiàng)分布以及基于平均場(chǎng)思想3種不同逼近方法的SIS模型，對(duì)模型求解，發(fā)現(xiàn)染病鄰居數(shù)服從泊松分布的逼近方法得到的疾病閾值略大于其他兩種。通過(guò)模擬與誤差分析，得到染病鄰居數(shù)服從泊松分布的逼近方法誤差最小。

[1] 靳禎，孫桂全，劉茂省.網(wǎng)絡(luò)傳染病動(dòng)力學(xué)建模與分析[M].北京：科學(xué)出版社，2014.

[2] ANDERSON R M,MAY R M,ANDERSON B.Infectious diseases of humans:dynamics and control[M].Oxford:Oxford university press,1992.

[3] HETHCOTE H W.The mathematics of infectious diseases[J].SIAM review,2000,42(4):599-653.

[4] KEELING M J,ROHANI P.Modeling infectious diseases in humans and animals[M].Princeton ：Princeton University Press,2008.

[5] KEELING M J,EAMES K T D.Networks and epidemic models[J].Journal of the Royal Society Interface,2005,2(4):295-307.

[6] MOTTER A E,MATíAS M A,KURTHS J,et al.Dynamics on complex networks and applications[J].Physica D:Nonlinear Phenomena,2006,224(1):vii-viii.

[7] FU X,SMALL M,WALKER D M,et al.Epidemic dynamics on scale-free networks with piecewise linear infectivity and immunization[J].Physical Review E,2008,77(3):036113.

[8] PASTOR-SATORRAS R,CASTELLANO C,VAN MIEGHEM P,et al.Epidemic processes in complex networks[J].Reviews of modern physics,2015,87(3):925.

[9] LUO X F,ZHANG X,SUN G Q,et al.Epidemical dynamics of SIS pair approximation models on regular and random networks[J].Physica A:Statistical Mechanics and its Applications,2014,410:144-153.

[10]MORRIS A J.Representing spatial interactions in simple ecological models[D].Warwick ：University of Warwick,1997.

[12]ERDDS P,WI A R.On random graphs I[J].Publ Math Debrecen,1959(6):290-297.

(責(zé)任編輯 陳 艷)

Comparison of Several Pair-Approximation Methods for SIS Epidemic Models in the Homogeneous Networks

LI Xinga,b, LUO Xiao-fenga,c

(a.Complex System Research Center; b.School of Mathematical Sciences; c.School of Computer & Information Technology, Shanxi University, Taiyuan 030006, China)

The accuracy of network-based epidemic models is a criterion of judging these models. Interestingly, judging the pair-approximation models in networks depends on the accuracy of approximation methods. Via SIS models, we compared the accuracy of three approximation methods in the homogeneous networks. The first is based on the number of infected neighborhoods of individuals following Poisson. The second is based the number of infected neighborhoods of individuals following multinomial distribution. The last one is based on the mean-field theory. Then we find the model under Poisson distribution with smallest errors, i.e. the accuracy of approximation methods under Poisson distribution is highest.

homogeneous network; pair-approximation; epidemic; basic reproduction number; stochastic simulation

2016-12-08 基金項(xiàng)目：國(guó)家自然科學(xué)基金青年科學(xué)基金資助項(xiàng)目(11501340)；山西省回國(guó)留學(xué)人員重點(diǎn)科研資助項(xiàng)目(2013-重點(diǎn)3)

李星(1990— )，女，碩士研究生，主要從事復(fù)雜網(wǎng)絡(luò)及其傳播，E-mail:lixing20469@gmail.com；羅曉峰(1987—)，男，博士研究生，主要從事網(wǎng)絡(luò)傳播動(dòng)力學(xué)研究，E-mail:luo_xf1988@163.com。

李星,羅曉峰.均勻網(wǎng)絡(luò)上SIS傳染病模型幾種對(duì)逼近方法的比較[J].重慶理工大學(xué)學(xué)報(bào)(自然科學(xué))，2017(5):169-173.

format：LI Xing, LUO Xiao-feng.Comparison of Several Pair-Approximation Methods for SIS Epidemic Models in the Homogeneous Networks[J].Journal of Chongqing University of Technology(Natural Science)，2017(5):169-173.

10.3969/j.issn.1674-8425(z).2017.05.028

O175

A

1674-8425(2017)05-0169-05