新冠肺炎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)-1if(idx>=25)r=1;r2=1;B=0.01;y=0.6;k=0.03;B1=0.01;endif(idx>=100)r=0;r2=0;y=0.9;k=0.05;endS(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('模型')