matlab建立模型对传染病进行预期

假设在某一封闭空间内有5000人，潜伏者（可以传染疾病的）有1人，传染率为0.4%，预期在固定时间段内感染者、康复者等人数

N=5000;%人口数
E=1;%潜伏者
D=0;%死亡患者人数
I=1;%感染人数
S=N-I;%易感人数
R=0;%康复者人数
r=20;%感染者接触数量
B=0.4;%传染概率
a=0.17;%潜伏者转化为感染者的概率
r2=8;%潜伏者接触人数
B2=0.03;%潜伏者传染正常人的概率
y=0.05;%康复概率
k=0.0001;%日致死率
B3=0.1;%转阴率
T=1:200;
for it=1:length(T)-1
if it>=7
r=0.20;
r2=1;
y=0.1
end
if it<7
S(it+1)=S(it)-rBS(it)I(it)/N-r2B2S(it)E(it)/N;%易感人群迭代
E(it+1)=E(it)+rBS(it)I(it)/N-aE(it)+r2B2S(it)E(it)/N;%潜伏者迭代
I(it+1)=I(it)+aE(it)-(k+y)I(it);%感染人数迭代
R(it+1)=R(it)+yI(it);%康复人数迭代
D(it+1)=R(it)+kI(it);%死亡患者人数迭代
else
S(it+1)=S(it