Epidemic model: Improved SEIR model based on Matalb

clear; clc; % -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- % parameter Settings %-------------------------------------------------------------------------- N =10000;                                                               
E = 10;                                                                 
I = 0;                                                                                                                     
R = 0;
UE = 0;
Q = 0;
D = 0;
S = N-I-E-R-UE-Q-D; 

R0 = 3;                                                                    
a = 0.5;                                                                  
b1 = 0.05;
b2 = 0.95;
r1=0.90;
e=0.01;
o1=0.96;
r2=0.10;
o2=0.4;
y = 0.004;                                                                 

T = 1:100;
for idx = 1:length(T)- 1
    S(idx+1) = S(idx)-R0*a*S(idx)*E(idx)/N(1);
    E(idx+1) = E(idx)+R0*a*S(idx)*E(idx)/N(1)-b2*E(idx)-b1*E(idx);
    UE(idx+1) = UE(idx)+b1*E(idx)-r1*UE(idx)-r2*UE(idx);
    Q(idx+1) =  Q(idx)+b2*E(idx)+r1*UE(idx)-e*Q(idx)-o1*Q(idx);
    I(idx+1) = I(idx)+r2*UE(idx)+e*Q(idx)-o2*I(idx)-y*I(idx);
    R(idx+1) = R(idx)+o1*Q(idx)+o2*I(idx);
    D(idx+1) = D(idx)+y*I(idx);
end
Copy the code

Version: 2014 a