SEIR代码新冠肺炎模型.zip

建模比赛刚刚用完亲测有效可出图出数据有分析有数据可直接运行

%SEIR模型修正
%假设第25天开始采取隔离措施（12月31日第一例确诊；1月23号开始封城，此时其他省市也相应作出响应隔离措施）
N=7500000;%哈尔滨市人口数
E=40;%潜伏者
D=0;%死亡患者人数
I=1;%感染人数
S=N-I;%易感人数
R=0;%康复者人数
SQ=0;
EQ=0;
beta=0.034;%传染率
q=0.001;%隔离比例
c=28;%接触比例范围[1.6，4.8]
lamda=1/14;%隔离解除速率
sigma=0.001;%潜伏转化成感染速率
alpha=0.0306%死亡率
r=0.30;%治愈率
w=0.3;%潜伏转化成易感人群速率
T=21:100;
for idx=1:length（T）-1
    if idx<13
         S（idx+1）=S（idx）-c*beta*（I（idx）+E（idx））*S（idx）/N;
         E（idx+1）=E（idx）+c*beta*（I（idx）+E（idx））*S（idx）/N-sigma*E（idx）;
         I（idx+1）=I（idx）+sigma*E（idx）-（r+alpha）*I（idx）;
         R（idx+1）=R（idx）+r*I（idx）;
    end
    if idx>=13
        beta=0.003;%传染率
        q=0.2;%隔离比例
        c=1;%接触比例范围[1.6，4.8]
        lamda=1/14;%隔离解除速率
        sigma=0.001;%潜伏转化成感染速率
        alpha=0.0306%死亡率
        r=0.63;%治愈率
        w=0.3;%潜伏转化成易感人群速率
        EQ（1:13）=0;
        SQ（1:13）=0;
        S（idx+1）=S（idx）-c*beta*（I（idx）+E（idx））*S（idx）/N-q*c*（1-beta）*（I（idx）+E（idx））*S（idx）/N+w*EQ（idx）;
        E（idx+1）=E（idx）+c*beta*（I（idx）+E（idx））*S（idx）/N-sigma*E（idx）-q*EQ（idx）;
        I（idx+1）=I（idx）+sigma*E（idx）-（r+alpha）*I（idx）+sigma*EQ（idx）;
        SQ（idx+1）=SQ（idx）+q*c*（1-beta）*（I（idx）+E（idx））*S（idx）/N-lamda*SQ（idx）;
        EQ（idx+1）=EQ（idx）+q*c*beta*E（idx）-w*EQ（idx）;
        R（idx+1）=R（idx）+r*I（idx）;
    end
end
%B={‘12-31‘‘01-26‘‘02-09‘‘02-29‘‘03-20‘‘04-09‘‘04-29‘‘05-19‘‘06-08‘‘06-28‘‘07-17‘};

%subplot（211）
plot（TI）;

grid on;
hold on;
clc;
filename=‘statistic1.xls‘;
A=xlsread（filename3‘E2:E61‘）;
C=zeros（861）;
C（27:86）=A（1:60）;
Y=C（1:86）;
X=1:length（Y）;
plot（XY‘r*‘）;
%ha=gca;
%set（ha‘xlim‘[2690]）
grid on;
hold on;
% plot（[25 25][0 2000]）;
% plot（[39 39][0 2000]）;
% plot（[53 53][0 2000]）;
% set（gca‘XTickLabel‘B）
xlabel（‘日期‘）;
ylabel（‘人数‘）;
legend（‘拟合患者数量‘‘实际数量‘）;
title（‘哈尔滨市采取隔离措施的SEIR拟合模型‘）;

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     文件        2047  2020-07-01 21:29  F2_SEIR.m
     文件        6816  2020-06-24 18:33  F_SEIR.m
     文件       47104  2020-06-24 13:57  statistic1.xls