新冠肺炎SEIR模型(修正版)

初始模型 seir模型

clear;clc;
N=1400000000;
I=1;
R=0;
D=0; 
E=0;
S=N-I;
r=5;
r11=0;
a=0.11;
B=0.05;
B1=0.04;
count=0;
y=0.02;
r2=15;
k=0.025373;
T=20:1000;
for idx =1:length(T)-1
     if(idx>=25)
         r=1;
         r2=1;
         B=0.01;
         y=0.6;
         k=0.03;
         B1=0.01;
     end
     if(idx>=100)
         r=0;
         r2=0;
        y=0.9;
       k=0.05;
     end
   S(idx+1)=S(idx)-r*B*I(idx)*S(idx)/N(idx)-r2*B*E(idx)*S(idx)/N(idx);   
   E(idx+1)=E(idx)-a*E(idx)+r*B*I(idx)*S(idx)/N(idx)+r2*B*E(idx)*S(idx)/N(idx); 
   I(idx+1)=I(idx)+a*E(idx)-(y+k)*I(idx); 
   R(idx+1)=R(idx)+y*I(idx);  
   D(idx+1)=D(idx)+k*I(idx);
   N(idx+1)=N(idx)-D(idx);
end
plot(T,S,T,E,T,I,T,R,T,D);
hold on
grid on;
xlabel('日期 ');
ylabel(' 人数');
legend('易感者','潜伏者','传染者','康复者','死亡者');
title('模型')

相比与之前的模型,这个模型考虑到了随着疫情的发展,各项参数都会发生变化.现在只能说是再次对整个疫情的一次数学模拟,并没有考虑地区的自然死亡率,以及出生率.
最终结果是模拟康复人数为8.4万人.
参考文档1

参考文档2
(致敬所有的抗疫英雄)