%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