MATLAB求解SEIR模型.zip

预测意大利疫情数据，MATLAB求解SEIR模型。建立的理论模型对意大利疫情在2月14日至3月10日的实际数据温度都较好（图？），上述模型的数据采用的2月14的数据进行预测的，可以很很明显的看出，后续的预测数据基本一致。3月9号，意大利官方宣布全国人员管控。根据管控的需要，启动了必要的防空隔离措施，如限制大范围聚集情况，对流动人口采取预防工作，并升级隔离措施，这实际上是降低接触率c。我们对3月9号的接触率进行一定的降低。因此对3月9号意大利的接触率进行一定程度的降低，分别取1，0.9，0.7，0.5；并作出了模拟情况。

clc;clear;close all% 计算意大利
hubei_data=xlsread（‘意大利数据.xlsx‘）;
hubei_data=hubei_data（13:end:）; %过滤2月14号以前数据
%% 参数
Q=1;    % 潜伏者相对于感染者传播能力的比值
A=1/14;  %隔离时间14天
G=1/7 ; % 潜伏期向感染者的转化速率潜伏期为7d
c=3.2;  %接触率
p=1;   %有效接触系数
Si=0.13;  %感染者隔离速率
Sq=0.13;  % 隔离
YI=0.007;YH=0.014; B=2.05*10^（-9）;
q=10^（-6）;a=0.0399;
h=0.01 ;  %求解的步长
T=60 ;  %预测50天的
t=0:h:T;
S=60430000 ;E=3;I=3*1.5;SQ=9;EQ=2;H=I+EQ;R=42;
DS=zeros（1length（t））;     DS1=zeros（1length（t））;     % 参数保存地方
DE=zeros（1length（t））;   DE1=zeros（1length（t））;
DI=zeros（1length（t））;   DI1=zeros（1length（t））;
DSQ=zeros（1length（t））;  DSQ1=zeros（1length（t））;
DEQ=zeros（1length（t））;   DEQ1=zeros（1length（t））;
DH=zeros（1length（t））;   DH1=zeros（1length（t））;
DR=zeros（1length（t））;   DR1=zeros（1length（t））;
DS（1）=S;%参数初始化
DE（1）=E;DI（1）=I;DSQ（1）=SQ;DEQ（1）=EQ;DH（1）=H;DR（1）=R;
for i=2:length（t）
    DS1（i-1）=-（p*c*B+p*c*q*（1-B））*（ DS（i-1） ）*（ DI（i-1）+1*DE（i-1）  ）+A*DSQ（i-1）;
    DS（i）=DS（i-1）+DS1（i-1）*h;
    DE1（i-1）=p*c*B*（1-q）*DS（i-1）*（ DI（i-1）+1*DE（i-1） ）-G*DE（i-1）;
    DE（i）=DE（i-1）+DE1（i-1）*h;
    DI1（i-1）=G*DE（i-1）-（Si+a+YI）*DI（i-1）;
    DI（i）=DI（i-1）+DI1（i-1）*h;
    DSQ1（i-1）=p*c*q*（1-B）*DS（i-1）*（ DI（i-1）+1*DE（i-1） ）-A*DSQ（i-1）;
    DSQ（i）=DSQ（i-1）+DSQ1（i-1）*h;
    DEQ1（i-1）=p*c*B*q*DS（i-1）*（ DI（i-1）+1*DE（i-1） ）-G*DE（i-1）;
    DEQ（i）= DEQ（i-1）+ DEQ1（i-1）*h;
    DH1（i-1）=Si*DI（i-1）+Sq*DEQ（i-1）-（a+YH）*DH（i-1）;
    DH（i）=DH（i-1）+DH1（i-1）*h;
    DR1（i-1）=（YI）*DI（i-1）+YH*DH（i-1）;
    DR（i）=DR（i-1）+DR1（i-1）*h;
end

figure（1）;plot（tDI）; hold on;plot（t-DR）; hold on ; 
plot（hubei_data（:1）-12hubei_data（:3）‘*‘）; hold on;
plot（hubei_data（:1）-12hubei_data（:4）‘*‘）; hold on;
xlabel（‘时间-初始时间是2月14日‘）;ylabel（‘人数（人）‘）;
legend（‘当前实际感染人数‘‘当前实际治愈人数‘‘理论预测-当前感染‘‘理论预测-当前治愈‘）;
title（‘意大利新冠肺炎感染趋势‘）;grid on

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     目录           0  2020-03-12 16:50  MATLAB求解SEIR模型\
     文件        1965  2020-03-11 17:18  MATLAB求解SEIR模型\Untitled4.m
     文件        2671  2020-03-11 17:06  MATLAB求解SEIR模型\Untitled5.m