|
EDA365欢迎您登录!
您需要 登录 才可以下载或查看,没有帐号?注册
x
S=1400050000;%中国人口总数
0 ^6 s5 I5 K- q, `- HE=79852;%潜伏者5 R. x7 o( C2 F& Z3 U& |7 a
I=10462;%未接受隔离的感染者3 W4 p' D: _- V" q l! T
I1=10462+7795+5000;
: o8 Q. ^; \' g! w: ISq=25000;%尚在接受医学观察的人数( M& ? M2 I5 [ Y1 \
Eq=5000;%估计值,为正在被隔离的潜伏者
) O/ E, Z6 p8 S. b. t; S$ ?H=7795;%正在住院的患者,为感染者和被隔离的潜伏者之和
/ T$ K# s+ a! C5 y3 m$ [, L5 M4 `R=263;%康复者人数
+ H$ G8 M5 o w7 u ? ZR1=263;%官方公布出院人数
" I, c+ X9 v5 {- v; QD=218;%死亡人数
3 m: G% W! f) B" P( `%模型参数设定
9 Y5 y- X+ o% a1 `: Tdelta=0.9;%感染者的隔离速度8 }1 y% y" ]+ r9 {( G" j
gamma=1/14;%隔离易感者转化为易感者概率' ?' a" ^, i# @0 l. c* h7 ~: `" M
gammaH=0.02;%隔离感染者的恢复速率* c2 n5 J- Y4 o2 u
gammaI=0.001;%感染者恢复率, \ K6 H" g( l$ m$ S) ~% p! G
beta=1/1.4*10^(-9);%传染概率; ^& L! D `7 u# d- o
sigma=1/7;%潜伏者向感染者的转化速度,平均潜伏期为7天,为7天的倒数
0 `4 W7 R7 G+ m2 w3 |d=0.03;%死亡率
- @1 O, G8 x1 W%差分迭代方程
, A1 B, c" i) t" kT=1:100;
6 M C, _! g m9 zfor i=1:length(T)-1
9 T& q# h) D, F" O- x S(i+1)=S(i)-(I(i)+E(i))*beta*S(i)+gamma*Sq(i);%易感人数迭代
: c2 D# h8 i7 i: U/ N E(i+1)=E(i)+beta*S(i)*(I(i)+E(i))-sigma*E(i);%潜伏者人数迭代
* J4 H% R& {8 c1 M A+ s0 o I(i+1)=I(i)+sigma*E(i)-(delta+d+gammaI)*I(i);%感染者人数迭代: [2 n7 Q& u5 }; s8 } S. [
I1(i+1)=I1(i)+sigma*E(i);%感染总人数3 P) n* R: J, q. |
Sq(i+1)=Sq(i)+beta*S(i)*(I(i)+E(i))-gammaI*Sq(i);%隔离易感染着人数迭代
" ~* T; M/ e* Y9 {' L Eq(i+1)=Eq(i)+beta*S(i)*(I(i)+E(i))-sigma*Eq(i);%隔离潜伏者人数迭代- O1 Q% E3 N- l- i& A1 R* d/ c: E
H(i+1)=H(i)+gammaI*I(I)+sigma*Eq(i)-(d+gammaH)*H(i);%住院患者人数迭代
. L8 G* I9 V' H1 H* ~ R(i+1)=R(i)+gammaH*H(i)+gammaI*I(i);%康复人数迭代- V/ h4 S4 m) ^" [
R1(i+1)=R1(i)+gammaI*I(i)+gammaH*H(i);%官方公布治愈人数迭代
) \0 h0 U' ?# d0 ~5 u, ]0 C D(i+1)=D(i)+I(i)*d;%死亡人数迭代" L! d" }% {/ T
end
2 c c- o8 x {3 L- ]9 ~9 _没有显示哪里出现错误,只说了标题这句话,求大神告知, l9 G8 J9 q8 J& K* t* R
|
|