清除;关闭全部的;中电;
装载('数据.mat','表格恢复','桌面死亡','表格确认','时间')
indLocation=查找(包含(tableRecovered.CountryRegion,'意大利')=1个)
恢复的表格(indLocation,1:2)
indLocation=indLocation(1)
恢复=表2排列(表恢复(indLocation,5:end));
死亡人数=table2array(TableDeasures(indLocation,5:end));
Confirmed=table2array(tableConfirmed(indLocation,5:end));
最小值=40;
恢复(确认<=minNum)=[];
死亡人数(确认人数<=minNum)=[];
时间(确认<=minNum)=[];
确认(确认<=minNum)=[];
Npop=120e6;
猜测=[0.06,1.2,1/5,1/40,0.01,0.02,0.01,0.02];
E0=确认(1);
I0=确认(1);
Q0=确认(1);
R0=恢复(1);
D0=死亡人数(1);
[字母1、beta1、gamma1、delta1、Lambda1、Kappa1]=...
fit_SEIQRDP(确认恢复的死亡,恢复的,死亡,Npop,E0,I0,时间,猜测);
dt=0.1;
time1=日期时间(time(1)):dt:日期时间(2020,3,25,0,0,0);
N=纽梅尔(时间1);
t=[0:N-1].*dt;
[S,E,I,Q,R,D,P]=SEIQRDP(字母1,beta1,gamma1,delta1,Lambda1,Kappa1,Npop,E0,I0,Q0,R0,D0,t);
图形
符号学(时间1,Q,“r”,时间1,R,“b”,时间1,D,“k”);
持有在
符号学(时间,确认恢复的死亡,“滚”,时间,恢复,“波”,时间,死亡,'ko');
伊拉贝尔('案例数')
xlabel公司('时间(天)')
支腿={“已隔离(确认感染)”,'恢复','死亡'};
图例(leg{:},'位置',“南外”)
设置(gcf,'颜色',“w”)
网格在
轴紧的
设置(gca,'yscale',“林”)
0条评论
登录发表评论。