2020年初,一种新型冠状病毒袭击了中国武汉。
随着卫生设施的改善、医疗水平的提高以及人类文明的不断发展,诸如霍乱、天花等曾经肆虐全球的传染性疾病已经得到有效的控制,但是一些新的、不断变异着的传染病毒却悄悄向人类袭来。长期以来,建立数学模型以描述疾病的传播过程从而寻找抑制传播的方法一直是科学家们探索的方向。但是由于疾病传播的复杂性,仅能从一般的传播机理角度建立大致的传播模型,现介绍一种与实际情况较为接近的传染病传播模型——
SIR模型
大多数传染病如天花、流感、肝炎、麻疹等治愈后均有很强的免疫力,所以病愈的人既非健康者(易感染者),也非病人(已感染者),他们已经退出传染系统,这种情况比较复杂,下面将详细分析建模过程。
一、基本假设
1.在疾病传播期内所考察地区的总人数N不变,既不考虑生死,也不考虑迁移.人群分易感染者(Susceptible)、**已感染者(Infective)和病愈免疫的移出者(Removed)**三类,以下简称健康者,病人和移出者,t时刻人数分别记为S(t),I(t),R(t),t时刻这三类人在总人数中所占的比例分别记作s(t),i(t)和r(t)。
2.每个病人每天有效接触的平均人数是常数λ,称为日接触率.当病人与健康者有效接触时,使健康者受感染变为病人。
3.每天被治愈的病人占病人总数的比例为常数μ,称为日治愈率,病人被治愈后仍有可能被感染为病人,那么可认为1/μ是该疾病的平均传染期。
4.初始时刻,只有少数个体处于感染状态,其他都是易染状态。
5.假设疾病的时间尺度远小于个体生命周期,从而不考虑个体的出生和自然死亡。
6.完全混合(Fully mixed):每一个个体与其他个体接触的机会均等。
不难看出,考虑以上几条假设,可得:
-
每个病人每天可使λs(t)个健康者变为病人。
-
每天有μNi个病人被治愈
-
定义接触数(σ):
σ = λ / μ \sigma=\lambda / \mu σ=λ/μ
可见接触数σ是病人平均传染期内有效接触的人数
二、模型建立
由假设1可得,整个群体由健康者,病人,移出者构成
N ( t ) = S ( t ) + I ( t ) + R ( t ) ≡ N \mathrm{N}(t)=S(t)+I(t)+R(t) \equiv N N(t)=S(t)+I(t)+R(t)≡N
根据假设1:
s ( t ) = S ( t ) N i ( t ) = I ( t ) N r ( t ) = R ( t ) N s ( t ) + i ( t ) + r ( t ) ≡ 1 \begin{aligned} &\mathrm{s}(t)=\frac{S(t)}{N}\\ &\mathrm{i}(t)=\frac{I(t)}{N}\\ &r(t)=\frac{R(t)}{N}\\ &s(t)+i(t)+r(t) \equiv 1 \end{aligned} s(t)=NS(t)i(t)=NI(t)r(t)=NR(t)s(t)+i(t)+r(t)≡1
由此,我们可以得到关于S(t),I(t),R(t)的微分方程
d S d t = − λ I ( t ) S ( t ) N d I d t = λ I ( t ) S ( t ) N − μ I ( t ) d R d t = μ I ( t ) \begin{aligned} &\frac{d S}{d t}=-\lambda I(t) \frac{S(t)}{N}\\ &\frac{\mathrm{d} I}{d t}=\lambda I(t) \frac{S(t)}{N}-\mu I(t)\\ &\frac{d R}{d t}=\mu I(t) \end{aligned} dtdS=−λI(t)N