|
EDA365欢迎您登录!
您需要 登录 才可以下载或查看,没有帐号?注册
x
S=1400050000;%中国人口总数6 N7 b! {! `$ y. W( R4 I
E=79852;%潜伏者8 ~1 M+ E2 Z) l3 _/ [, W
I=10462;%未接受隔离的感染者
: J; t* L5 N g+ @7 JI1=10462+7795+5000;% N( t1 Y& g9 F) Y2 W% B
Sq=25000;%尚在接受医学观察的人数6 _0 p, ^0 p+ W! b5 Z
Eq=5000;%估计值,为正在被隔离的潜伏者4 D# u- [3 H- H3 `
H=7795;%正在住院的患者,为感染者和被隔离的潜伏者之和
$ c% _& t/ J% \8 v, cR=263;%康复者人数
! W: ?& E6 Y; m/ W, z% B+ h7 W7 eR1=263;%官方公布出院人数
7 `/ i" ?; N: F$ l% dD=218;%死亡人数4 b. G. m) h' U3 M
%模型参数设定' t, D3 {( t @4 l9 ]
delta=0.9;%感染者的隔离速度0 H& c; x, d, G& q, g2 \
gamma=1/14;%隔离易感者转化为易感者概率) F! o0 X# L/ n. S
gammaH=0.02;%隔离感染者的恢复速率
5 U% Q0 U! L$ ogammaI=0.001;%感染者恢复率
1 ?, N: m, h- E3 ?& L4 Vbeta=1/1.4*10^(-9);%传染概率
2 v( u1 { W. S+ Y0 \, ^sigma=1/7;%潜伏者向感染者的转化速度,平均潜伏期为7天,为7天的倒数9 X+ s0 y g L- U. ^7 q# z) M
d=0.03;%死亡率
" |% ^) P( s5 u- y% k%差分迭代方程: d" N D, E3 s: r% n
T=1:100;
- H8 p y# r# e( O- D! G0 [for i=1:length(T)-1$ W! \. S! F1 @
S(i+1)=S(i)-(I(i)+E(i))*beta*S(i)+gamma*Sq(i);%易感人数迭代
: U8 l0 T1 L9 K/ N E(i+1)=E(i)+beta*S(i)*(I(i)+E(i))-sigma*E(i);%潜伏者人数迭代3 W$ M: z# J, A$ J
I(i+1)=I(i)+sigma*E(i)-(delta+d+gammaI)*I(i);%感染者人数迭代* X1 \" k3 ?! `& W p3 m" z% c
I1(i+1)=I1(i)+sigma*E(i);%感染总人数
% L' C' G3 C& J: b! t3 k Sq(i+1)=Sq(i)+beta*S(i)*(I(i)+E(i))-gammaI*Sq(i);%隔离易感染着人数迭代, l# ?. E% n: @! u. M/ E
Eq(i+1)=Eq(i)+beta*S(i)*(I(i)+E(i))-sigma*Eq(i);%隔离潜伏者人数迭代0 O5 ^- n5 I4 M/ _: @2 V& k; @! m
H(i+1)=H(i)+gammaI*I(I)+sigma*Eq(i)-(d+gammaH)*H(i);%住院患者人数迭代
1 Z; D- q, g$ s2 v R(i+1)=R(i)+gammaH*H(i)+gammaI*I(i);%康复人数迭代
" I. ?% C9 ?2 q2 Y5 | R1(i+1)=R1(i)+gammaI*I(i)+gammaH*H(i);%官方公布治愈人数迭代$ |; _8 ^0 ?8 w! m9 L9 ~
D(i+1)=D(i)+I(i)*d;%死亡人数迭代
& k& Z/ Q1 M4 j) ~: bend- K# K( |$ a0 W- U' w8 c4 Q9 k
没有显示哪里出现错误,只说了标题这句话,求大神告知
; L" Y$ P6 ]0 r8 N* t |
|