王治萍
(朔州師范高等??茖W(xué)校 數(shù)學(xué)與計算機系,山西 朔州 036000)
考慮靜態(tài)均勻網(wǎng)絡(luò)上的SIR對逼近傳染病模型,即網(wǎng)絡(luò)的規(guī)模N保持不變,網(wǎng)絡(luò)的結(jié)構(gòu)保持不變。其中節(jié)點的狀態(tài)有三種狀態(tài),易感狀態(tài)S和染病狀態(tài)I以及恢復(fù)狀態(tài)R; 并且[S]表示網(wǎng)絡(luò)中處于易感狀態(tài)的節(jié)點(即易感者)數(shù)量,[I]表示網(wǎng)絡(luò)中處于染病狀態(tài)節(jié)點(即染病者)的數(shù)量,[R]表示網(wǎng)絡(luò)中處于恢復(fù)狀態(tài)節(jié)點(即恢復(fù)者)的數(shù)量。[SI],[SS],[SR],[II],[IR]表示網(wǎng)絡(luò)中二元組(對)狀態(tài)為S-I,S-S,S-R,I-I,I-R的數(shù)量。[SSI],[ISI],[RSI]表示網(wǎng)絡(luò)中狀態(tài)為S-S-I,I-S-I,R-S-I三元組的數(shù)量。一個染病者和一個易感者的接觸傳染率為λ,染病者的恢復(fù)概率為γ0。應(yīng)用主方程可得關(guān)于[S],[I],[SI],[SS],[SR],[II],[IR]的微分方程組如下
(1)
[S]+[I]+[R]=N,
[SS]+2[SR]+2[IR]+[II]+2[SI]=n1N
(2)
下面應(yīng)用三種不同的逼近方法封閉模型進而研究這三種逼近方法的優(yōu)劣。
當(dāng)節(jié)點的染病狀態(tài)I鄰居數(shù)服從泊松(Poission)分布時,存在以下的三元組逼近公式
結(jié)合(2)式可以對系統(tǒng)(1)進行封閉和降維,然后得到染病鄰居數(shù)服從Poission分布時SIR對逼近傳染病模型(P-PW)如下
(3)
下面可得無病平衡點為P1(N,0,n1N,0,0,0),其中染病倉室為[I],[SI]運用下一代再生矩陣的方法可得
當(dāng)節(jié)點的染病狀態(tài)I鄰居數(shù)服從多項式分布時,存在以下的三元組逼近公式
同樣結(jié)合(2)式對系統(tǒng)(1)進行封閉和降維,然后得到染病鄰居數(shù)服從多項式分布是SIS對逼近傳染病模型(B-PW)。
(4)
可得無病平衡點為P1(N,0,n1N,0,0,0),其中染病倉室為[I],[SI]運用下一代再生矩陣的方法可得
Kiss基于平均場思想提出了三元組的另一種逼近公式
同樣結(jié)合(2)式可以對系統(tǒng)(1)進行封閉和降維,然后得到相應(yīng)SIR對逼近傳染病模型(MF-PW)。
可得無病平衡點為P2(N,0,n1N,0,0,0),其中染病倉室為[I],[SI]運用下一代再生矩陣的方法可得
由以上討論可知,節(jié)點的染病鄰居數(shù)服從泊松分布所得疾病的基本再生數(shù)大于染病鄰居數(shù),服從多項式分布和基于平均場逼近方法得到的基本再生數(shù);同時染病鄰居數(shù)服從多項分布和基于平均場逼近方法得到的基本再生數(shù)相同。