%SEIR模型修正
%假設(shè)第5天開始采取隔離措施1月19日第一例確診；1月23號武漢開始封城，此時其他省市也相應(yīng)作出響應(yīng)隔離措施
N=140000000;%人口數(shù)
E=0;%潛伏者
D=0;%死亡患者人數(shù)
I=1;%感染人數(shù)
S=N-I;%易感人數(shù)
R=0;%康復(fù)者人數(shù)
r=20;%感染者接觸數(shù)量
B=0.6;%傳染概率
a=0.17;%潛伏者轉(zhuǎn)化為感染者的概率
r2=8;%潛伏者接觸人數(shù)
B2=0.03;%潛伏者傳染正常人的概率
y=0.05;%康復(fù)概率
k=0.0001;%日致死率
B3=0.1;%轉(zhuǎn)陰率
T=1:200;
for?idx=1:length（T）-1
????if?idx>=7
????????r=0.20;
????????r2=1;
????????y=0.1
????end
????if?idx<7
????????S（idx+1）=S（idx）-r*B*S（idx）*I（idx）/N-r2*B2*S（idx）*E（idx）/N;%易感人群迭代
????????E（idx+1）=E（idx）+r*B*S（idx）*I（idx）/N-a*E（idx）+r2*B2*S（idx）*E（idx）/N;%潛伏者迭代
????????I（idx+1）=I（idx）+a*E（idx）-（k+y）*I（idx）;%感染人數(shù)迭代
????????R（idx+1）=R（idx）+y*I（idx）;%康復(fù)人數(shù)迭代
????????D（idx+1）=R（idx）+k*I（idx）;%死亡患者人數(shù)迭代
????else
????????S（idx+1）=S（idx）-r*B*S（idx）*I（idx）/N-r2*B2*S（idx）*E（idx）/N+B3*E（idx-6）;%易感人群迭代
????????E（id