4.4 【实验】舱室模型.docx

【实验】舱室模型

4.4.1舱室模型建模方法

SIR模型是“舱室”传染病模型的基础，可以进一步考虑更复杂的人群划分和转移，如考虑出生率、死亡率、防疫措施的作用、潜伏期、抵抗能力、地域传播、传播途径（接触、空气、昆虫、水源等）等。比如，

SEIS模型：易感—暴露—感染—不免疫。

SEIR模型：易感—暴露—感染—移出（免疫）。

SIRS模型：易感—感染—短时免疫—易感。

回顾SIR模型，梳理舱室转移关系图与微分方程组模型之间的对应关系（稍作修改），如图4-13所示。

4-13修长版SIR模型示意图

●每个舱室（即图中的圆圈）代表一类人群（占比），对应一个未知函数，箭头及上方的标签表示舱室之间的转移关系和转移速度。

●一个舱室表示为一个微分方程，左端导数代表该类人群变化速度。

●微分方程右端式子：出的箭头就是“减”，入的箭头就是“加”，都是用转移速度乘以来自的舱室函数。

●整体流量（流入、流出）平衡，比如1式与2式中的两个同项，2式与3式中的两个同项。

注意，由于S(t)+I(t)+R(t)=1以及流量平衡关系，这样表示出来的微分方程组中有一个微分方程是多余的（可由其余微分方程表示），数值求解时需要任意去掉一个。

所以，可以根据具体传染病的实际情况，建立任何复杂的舱室模型，只要先绘制出舱室转移关系图，然后按上述规律表示为微分方程组，再用ode45()函数进行数值求解即可。

4.4.2SEIR模型

在SIR模型的基础上，将已感染但处于潜伏期的人群（称为暴露者E）也考虑进来，即构建SEIR模型。

SEIR模型考虑易感者、暴露者、感染者、移出者四类人群，适合有潜伏期、治愈后获得终身免疫的传染病，如带状疱疹。

先绘制舱室转移关系图，如图4-14所示。

4-14SEIR模型舱室转移图

易感类S先以λI(t)的速度转移到暴露类E，再以α的速度转移到感染类I，最后以γ的速度转移到移出类R。再按前文的表示规律写出微分方程组：

由于S(t)+E(t)+I(t)+R(t)=1和流量平衡，有一个微分方程是多余的，去掉最后一个，得到SEIR模型：

满足初值条件：S(0)=S0，E(0)=E0，I(0)=I0。

【例s4_5】：

%%SEIR模型

%y(1)=S(t),y(2)=E(t),y(3)=I(t)

f=@(t,y,lambda,alpha,gamma)[-lambda*y(3)*y(1);

lambda*y(3)*y(1)-alpha*y(2);

alpha*y(2)-gamma*y(3)];

lambda=0.6;

alpha=0.5;

gamma=0.3;

SEIRfun=@(t,y)f(t,y,lambda,alpha,gamma);

S0=0.98;

E0=0.01;

I0=0.01;

y0=[S0;E0;I0];

tspan=[0,50];

[t,y]=ode45(SEIRfun,tspan,y0);

R=1-y(:,1)-y(:,2)-y(:,3);

sol=ode45(SEIRfun,tspan,y0);

t1=1:0.2:2;

deval(sol,t1)%t必须在tspan范围内

plot(t,y(:,1),t,y(:,2),r-.,t,y(:,3),k.,t,R,g),gridon;

xlabel(t),ylabel(y(t))

legend(S(t),E(t),I(t),R(t))

4-15模拟SEIR模型的疾病传播规律

以上四种模型只是传染病传播的基本模型，还有很多没有考虑到的因素，如人口的出生与死亡、迁入和迁出等。还可以考虑更细致的因素，如人群流动速度、易感人群的年龄分布、不同人群对疾病的易感性、患病者的症状轻重、人口密度、医疗卫生程度、检验检疫手段、政府重视程度、隔离措施以及人群心理因素等。这些因素都对暴露数、发病率、治愈率、传染期长度有着直接或间接的影响。根据需要自行添加延伸即可，原理相通。