QQ登录

只需要一步,快速开始

 注册地址  找回密码
查看: 3302|回复: 0
打印 上一主题 下一主题

Python小白的数学建模课-B2. 新冠疫情 SI模型

[复制链接]
字体大小: 正常 放大
杨利霞        

5273

主题

82

听众

17万

积分

  • TA的每日心情
    开心
    2021-8-11 17:59
  • 签到天数: 17 天

    [LV.4]偶尔看看III

    网络挑战赛参赛者

    网络挑战赛参赛者

    自我介绍
    本人女,毕业于内蒙古科技大学,担任文职专业,毕业专业英语。

    群组2018美赛大象算法课程

    群组2018美赛护航培训课程

    群组2019年 数学中国站长建

    群组2019年数据分析师课程

    群组2018年大象老师国赛优

    跳转到指定楼层
    1#
    发表于 2021-7-21 17:25 |只看该作者 |倒序浏览
    |招呼Ta 关注Ta

    # _, z' [1 y8 ~7 E0 TPython小白的数学建模课-B2. 新冠疫情 SI模型
    ! g) e- G/ W5 s0 @# t; U
    ) D( n7 W: Q# d5 w" ?% I
    2 `4 {$ e& `: I7 ]
    传染病的数学模型是数学建模中的典型问题,常见的传染病模型有 SI、SIR、SIRS、SEIR 模型。* J3 ]  ?* x5 {- \6 r
    9 h4 q1 _5 O+ }3 a; z" s- [
    ( x$ i$ @: G3 k* m' @, Y8 E* c
    SI 模型是最简单的传染病模型,适用于只有易感者和患病者两类人群。. e8 T: C' y4 ]& h; \* g( Z' z

    ; Y) g+ A8 E& Z7 u/ v2 h2 b
    / m* M( V) k* o7 ^
    我们就从 SI 模型开始吧,从模型、例程、运行结果到模型分析,全都在这个系列中。/ n+ I) ?# Q1 P/ F1 X3 {

    ; B1 I" _. \. p& M. w
    * N2 _3 x0 Y6 f' [+ x5 W
    『Python小白的数学建模课 @ Youcans』带你从数模小白成为国赛达人。4 q! H: `4 I, m! ~% w5 j6 m

    / }# O& k1 q9 D; V  [
    7 s: Z: X* i5 j+ B( O! M9 A5 n( N' w
    , w4 c% i4 d" d: u, r

    6 }* k6 V* o) w& r1. 前言
    . q  d6 v0 x  e5 @5 A7 B8 Y新冠疫情不仅严重影响到全球的政治和经济,深刻和全面地影响着社会和生活的方方面面,也已经成为数学建模竞赛的背景帝。8 ?8 D0 E$ Z0 J+ N3 O# h. Y3 V
    . B5 p* k+ d) s4 K9 q& S4 {
    ' ]' \* o0 O; Q* C+ R' p/ G# K. F
    传染病的数学模型是数学建模中的典型问题,标准名称是流行病的数学模型(Mathematical models of epidemic diseases)。建立传染病的数学模型来描述传染病的传播过程,研究传染病的传播速度、空间范围、传播途径、动力学机理等问题,以指导对传染病的有效地预防和控制,具有重要的现实意义。
    & e6 \1 Q/ z1 W$ w3 p  p8 f8 l2 _' d( h, }" x

    6 A# q) u( c5 s+ X2 j8 R: w# f, C+ P不同类型传染病的传播具有不同的特点,传染病的传播模型不是从医学角度分析传染病的传播过程,而是按照传播机理建立不同的数学模型。* h% V3 O+ L# u

    9 M4 V; O4 ~. j1 H4 N, n3 J" a

    / ?" `+ |. e; B( A' ?' j首先,把传染病流行范围内的人群分为 S、E、I、R 四类,具体含义如下:
    5 B6 g7 \: G2 H% N+ d( |& o7 j
    0 e' ^* ~( v. `5 I  l/ m
    7 _  j% ]" ]. r! t9 w
    S 类(Susceptible),易感者,指缺乏免疫能力的健康人,与感染者接触后容易受到感染;
    8 ?+ \" [, P8 o4 i/ H8 o4 j9 A% G6 B* l3 p

    4 }3 N4 B, R; L9 y- HE 类(Exposed),暴露者,指接触过感染者但暂无传染性的人,适用于存在潜伏期的传染病;9 @& H5 J$ {$ C- n( D
    1 C# s  A* |: a9 K
    . s3 T" O9 v- R1 G" k' y& r
    I 类(Infectious),患病者,指具有传染性的患病者,可以传播给 S 类成员将其变为 E 类或 I 类成员;
    3 p6 O  S; [; y) |2 }
    & s1 Y9 H( i) B+ E+ r7 |
    + X! `' N" a" S3 R9 s7 t1 w
    R 类(Recovered),康复者,指病愈后具有免疫力的人。如果免疫期有限,仍可以重新变为 S 类成员,进而被感染;如果是终身免疫,则不能再变为 S类、E类或 I 类成员。4 i& Q0 A/ o' c2 l  H9 y( ?
    # `( p: }6 _) X0 S/ ]
    ) R4 x9 x. ~, }
    常见的传染病模型按照传染病类型分为 SI、SIR、SIRS、SEIR 模型等,就是由以上四类人群根据不同传染病的特征进行组合而产生的不同模型。+ {' z+ z) P: W( u- k
    8 y$ w7 v0 F! W. ?5 O8 h

    8 z8 `: Z5 d& s9 {+ i9 n; a1 \% t% J1 _

    : ?$ c" `+ L" L; \0 mPython小白的数学建模课-A3.12个新冠疫情数模竞赛赛题及短评
    8 o! \0 T4 ~9 s, `2 E0 ~# a+ A' tPython小白的数学建模课-B2. 新冠疫情 SI模型
    7 X& J4 y! P# a9 aPython小白的数学建模课-B3. 新冠疫情 SIS模型, Z9 e3 ^8 w: e1 y" ?
    Python小白的数学建模课-B4. 新冠疫情 SIR模型
    - |! W0 z% x0 D1 zPython小白的数学建模课-B5. 新冠疫情 SEIR模型
    2 [. C0 K& ~1 |$ G8 YPython小白的数学建模课-B6. 新冠疫情 SEIR改进模型
    7 Z2 x5 @6 I7 t; V4 j/ DPython数模笔记-PuLP库
    " g) U5 C5 ~, w4 \  G/ ^: N) j9 x  x' c  ~2 T% c) C8 u

    5 P& {9 a& z. E7 L" X
      j% U* a. q; T* V

    * j4 s1 C: L  O( w) z; R4 R% W2. 疫情传播 SI 模型3 B; U7 |& ^' r0 b4 C
    2.1 SI 模型的适用范围
    " l0 q0 ~# i% WSI 模型适用于只有易感者和患病者两类人群,且无法治愈的疾病,例如 T型病、僵尸。
    0 Z" [/ e" z- f. W/ B4 u
    ; }+ G; Z7 |" w2 x
    ) f% R7 p4 {% e9 P" J6 F( ^
    * m, U& |. H- C8 P* A% {
    5 z6 I" l+ Z3 \. ]
    2.2 SI 模型的假设
    6 b+ m- L8 x9 b' Q3 K, P考察地区的总人数 N 不变,即不考虑生死或迁移;
    9 c! K: w! C6 N7 @) n人群分为易感者(S类)和患病者(I类)两类;6 q2 [+ J' u' C% ~' d6 i4 \7 u- b
    易感者(S类)与患病者(I类)有效接触即被感染,变为患病者,无潜伏期、无治愈情况、无免疫力;
    - B6 e; r8 @* b7 B- A每个患病者每天有效接触的易感者的平均人数(日接触数)是 λ \lambdaλ,称为日接触率;
    + E+ S( A. U# Y. a# ^2 L将第 t 天时 S类、I 类人群的占比记为 s ( t ) s(t)s(t)、i ( t ) i(t)i(t),数量为 S ( t ) S(t)S(t)、I ( t ) I(t)I(t);初始日期 t = 0 t=0t=0 时, S类、I 类人群占比的初值为 s 0 s_0s 4 e6 B, z) X2 H' s
    0
    - I' P1 b8 ^7 D+ P) a# v​       
    0 ?7 h, u8 W  c6 h7 }# A! S0 p 、i 0 i_0i " g" d3 ]; s9 E5 a- b: C
    0
    7 }" _* ~+ F2 ?8 I​       
    / b- x8 }: L0 H+ g9 a* Q1 _. C% W% z2 T- ]! D# U) J+ m5 ?
    2.3 SI 模型的微分方程
    ( }1 P" h( N$ C2 O- U
    / I) P- c7 X. U8 o* `1 I; pN d i d t = N λ s i N\frac{di}{dt} = N\lambda s i
    ! ?' ~. ?+ D  l( M  u! bN
    * x$ P5 _1 H0 b* Bdt3 @+ }3 Z3 v1 H3 T5 I. C" c$ H) {
    di5 _9 s/ w& k6 s" o; t  U
    ​       
    , v* t* o2 S# n0 H =Nλsi; {; x9 z) Z2 W" @
    * B: m" Y. t5 t% q+ X
    6 L, l% x' h% O. P, g# ~% H
    得:
    : V. d4 F( I4 ]3 Z& z" v  t) Jd i d t = λ i ( 1 − i ) ,   i ( 0 ) = i 0 \frac{di}{dt} = \lambda i (1-i),\ i(0) = i_0$ Z2 u: }  A! {2 y; L' j( A
    dt
    9 r6 q& b' J5 sdi
    3 m. P0 k$ e* h+ O​        " j3 @0 z# W3 I+ I0 v5 S, G  h: |. ]% J
    =λi(1−i), i(0)=i
    : R: _9 k) o& K  H" |& M0- \% ?" P- r* c7 W! d8 l$ t0 W6 |, v
    ​       
    9 \: u: `# p' S0 F$ Z/ D ; r5 b8 N& C: b- R7 {
    - u5 F3 U/ L) O- J$ r2 e) w$ [5 ]
    . F7 Q8 O4 K5 l, `; L
    这是 Logistic 模型,用分离变量法可以求出其解析解为:
    / m. x# {# o) r( T) V! A+ @i ( t ) = 1 1 + ( 1 / i 0 − 1 )   e − λ t I ( t ) = N   i ( t ) i(t)=\frac{1}{1+(1/i_0 - 1)\ e^{-\lambda t}}\\ I(t)= N\ i(t)
    - `; P" q. a9 O- ei(t)=
    " [& j, q$ j, N6 L* l# ~& ~1+(1/i
    8 f4 N! Q( ^+ i4 z. V8 `  [+ v0
    ' Y5 H# f) Y7 H2 G​       
    1 w" y( y0 I1 Q1 W, Q −1) e $ Q' Q) p8 H3 O: t
    −λt
    ! A; ^7 m) b$ G* }! c; e7 k! f9 u
    # }, r( Y& Q& V3 i1
      t: |. C: C' M* t3 _​       
    ) k; K# s+ N% ]. c+ ?+ q8 [/ J
    ' T* p) U* ]1 L8 G( mI(t)=N i(t)
    / r1 Q# e" m3 Y+ s5 ^
    - `/ C9 ^$ `! ^# }. x2 T

    ) i) F6 [3 A9 n( [, P! u1 w3 w+ z8 y$ X
    ( Q: g2 r1 O, L0 Q' ?- ]0 s
    ) }9 g6 w/ i# P5 m
    3. SI 模型的 Python 编程2 R( F9 l+ p( R" w
    3.1 SI 模型的解析解
    4 ~0 Y, @. g5 T; G( `上文已经得到 SI 模型的解析解,对此很容易通过 Python 编程实现,详见本文例程。! f+ U1 @9 A) o4 {
    6 s) d: ~- r8 ?
    2 \& R) C$ R! @) \
    虽然 SI 模型的解析解并不复杂,而且解的精度当然是最好的,但我们仍然不鼓励用解析解的方法。原因在于,一是对于小白求解析解的过程相对复杂困难,而且可能出错,二是对于更复杂的模型是没有解析解的,即便大神也只能用数值方法求解。既然如此,不如从一开始就学习、掌握数值求解方法,熟悉数值解法的编程实现。
      a1 f' W2 U% N' K/ b" L! x( e) k! [

    8 j- w- h: c/ M2 A$ t* F. E" _3.2 SI 模型的数值解
    , M  N- F/ @  U2 t7 F4 G( ]$ {SI 模型是常微分方程初值问题,可以使用 Scipy 工具包的 scipy.integrate.odeint() 函数求数值解,具体方法可以参考前文《Python小白的数学建模课-09 微分方程模型》。# @0 |2 x/ x) l0 g8 T8 m6 ]

    / _  `6 N' T2 Y$ C* Y3 V# _6 [& ~( X

    ; W  }/ S& L% F+ E0 g+ dscipy.integrate.odeint(func, y0, t, args=())- Q. D) i8 X0 ^
    1+ X$ @5 a$ l0 J! \& k0 O. h
    **scipy.integrate.odeint()**是求解微分方程的具体方法,通过数值积分来求解常微分方程组。) F, ?; B. j% L* e$ [

    0 a! i1 j! x, Q/ e/ ?# a, S5 k
    # U( q' h# q2 D) f2 D# u
    odeint() 的主要参数:+ _+ L/ z. K5 M; }* Z( J, F6 t

    + E2 X" n! \5 j: j# S3 i& Z
    . [$ [8 i1 U' _8 n. I
    func: callable(y, t, …)   导数函数 f ( y , t ) f(y,t)f(y,t) ,即 y 在 t 处的导数,以函数的形式表示3 L7 x: j5 v4 T3 K
    y0: array:  初始条件 y 0 y_0y
    5 k( t2 S" X5 d  Z0
    % \% ~' E( o0 F- S% I​        * \* `+ {& \0 O2 y* Z) S* ~- r2 z
    ,对于常微分方程组 y 0 y_0y " x5 Q2 G& t5 P8 h4 _7 O- Z
    0; D2 v' @  e) L2 k
    ​        9 D: h( W. _$ U: Q# S
      则为数组向量
    # F" j6 g! j9 s5 ht: array:  求解函数值对应的时间点的序列。序列的第一个元素是与初始条件 y 0 y_0y , Q3 I( C# N' K5 M4 r
    0# v' Y- ~7 Q  M0 ^8 Z+ l9 n7 S
    ​       
    + D6 _1 g* s; E" n* J& [3 n  对应的初始时间 t 0 t_0t 9 ^6 v! G, P8 ^- r  b0 i+ R+ `2 M) G
    0
    9 w+ F% P5 U. U. L, Z, Q7 q​       
    : V% C  N, P- E: P( ]) O ;时间序列必须是单调递增或单调递减的,允许重复值。
    . ^& d' h8 i0 X" V" b3 Zargs: 向导数函数 func 传递参数。当导数函数 f ( y , t , p 1 , p 2 , . . ) f(y,t,p1,p2,..)f(y,t,p1,p2,..) 包括可变参数 p1,p2… 时,通过 args =(p1,p2,…) 可以将参数p1,p2… 传递给导数函数 func。
    ' V$ K* F% c3 y( P* aodeint() 的返回值:0 ^2 k' U# F" A5 D# l! e0 [
    + U( B: b! o  b5 A$ L
    5 K) M2 c. J% Z% `* Z
    y: array   数组,形状为 (len(t),len(y0),给出时间序列 t 中每个时刻的 y 值。8 ^+ k0 K5 J$ Q1 h
    odeint() 的编程步骤:
    7 q9 e, h' c6 X. d
    9 z  n( i1 ^$ Y5 _5 D/ q: K$ w' d

    * x4 x4 w  K9 N7 B, X! N导入 scipy、numpy、matplotlib 包;7 p7 Q; S+ f' S+ \* A& e# p
    定义导数函数 f ( i , t ) = λ i ( 1 − i ) f(i,t)=\lambda i (1-i)f(i,t)=λi(1−i) ;: c: C  i: \6 ?# r) [% ?& v9 m
    定义初值 i 0 i_0i " N& z: c7 t' @# E0 Z
    0
    ' [( a! _* U0 a; D7 n​        + X/ P  U; W9 m8 _4 \/ D' W# I
      和 i ii 的定义区间 [ t 0 ,   t ] [t_0,\ t][t + H# u% p; U8 A% T( L
    0
    * L0 d) F4 p8 N​       
    # I2 f8 H9 y$ _, t) \$ J , t];( Z. q2 H4 r0 y: O0 e
    调用 odeint() 求 i ii 在定义区间 [ t 0 ,   t ] [t_0,\ t][t 4 F6 B- U; o- a1 M% o6 O
    0
    + ]* D1 q9 _  n! A3 x' P​       
    : |) h3 S* U0 Y , t] 的数值解。  v( M% _# B- C+ C4 e2 J
    0 K- @% \- L2 w$ a; e3 z2 A# F
    1 t1 j5 n& y" `4 w& w8 x8 f5 {
    3.3 Python例程:SI 模型的解析解与数值解
    % m# B: I+ w4 x/ R% O2 `: p# 1. SI 模型,常微分非常,解析解与数值解的比较5 `) m$ I3 x6 W7 ^& U; O
    from scipy.integrate import odeint  # 导入 scipy.integrate 模块0 }- i3 ~; [8 p  K& Y/ }% `4 z
    import numpy as np  # 导入 numpy包5 B) W. T9 E3 Y" R' f) m1 [
    import matplotlib.pyplot as plt  # 导入 matplotlib包
    8 h0 V! \% s$ r9 }0 b* Z
    / v2 f5 T: T& k0 W6 _
    9 b3 N+ Y! `  y8 N( i& a
    def dy_dt(y, t, lamda, mu):  # 定义导数函数 f(y,t)7 _) m6 c; P9 b' T
        dy_dt = lamda*y*(1-y)  # di/dt = lamda*i*(1-i)/ b9 M3 m3 L# n4 D! ?
        return dy_dt
    $ d: k1 a; c  _+ o" n- q9 Y4 U3 c# k
    5 ~# ]0 {2 T1 S% j6 F& Q

    ! m6 J4 k8 h" w& e( F# 设置模型参数
      p* |5 A1 h5 z( inumber = 1e7  # 总人数
    9 r4 q8 C- ]) S0 E( u6 e" b# clamda = 1.0  # 日接触率, 患病者每天有效接触的易感者的平均人数8 D- m, N7 ]+ `- o
    mu1 = 0.5  # 日治愈率, 每天被治愈的患病者人数占患病者总数的比例
    # I6 W) R) H- Y# f. xy0 = i0 = 1e-6  # 患病者比例的初值
    9 O; ]4 v) ?$ `! mtEnd = 50  # 预测日期长度) j( c! d5 B8 ^3 B2 `7 K0 v) P, A
    t = np.arange(0.0,tEnd,1)  # (start,stop,step)
    ' [: P; D+ i2 w1 y  u' y- E# s3 y/ Q' l/ \, _* C; d, X9 \( p6 |
    ; j: E1 ^+ d( I+ B6 d/ A
    yAnaly = 1/(1+(1/i0-1)*np.exp(-lamda*t))  # 微分方程的解析解
    ' Q' T5 K- i5 JyInteg = odeint(dy_dt, y0, t, args=(lamda,mu1))  # 求解微分方程初值问题6 I; ?- g  A) M& Q+ r: T; z, S' i
    yDeriv = lamda * yInteg *(1-yInteg)
    ; A5 m1 j, A, l* t7 ?" }- d
    2 h) t: a0 {+ |
    8 b% l$ n3 E" L& z
    # 绘图
    5 u! T# x' X9 ~7 B! q" lplt.plot(t, yAnaly, '-ob', label='analytic')
    ; Z' a/ a; j* [, hplt.plot(t, yInteg, ':.r', label='numerical')7 j  l1 \4 ^/ \6 Q1 V
    plt.plot(t, yDeriv, '-g', label='dy_dt')
    8 v% q! A. S, L$ mplt.title("Comparison between analytic and numerical solutions")' s: _8 q; Z" \4 g, D) I
    plt.legend(loc='right')- w0 Y0 {' N! {& Y: e
    plt.axis([0, 50, -0.1, 1.1])$ z- u. b% ?5 X, x6 S: b: c
    plt.show()
    - i5 H* p3 Y6 ?' r% A1
    5 }- @$ V# i. E9 c$ z  Y2
    + l4 e% G8 Q+ w4 g4 J& U3 }4 }3 |3
    3 z6 V0 x$ z6 d0 U6 p3 g5 U4
    - x2 u( Y9 [$ K% @: s4 m/ U5
    5 u& S. q7 n2 _. j# ]6) I$ X, m( O; p2 ^
    7
    / `- r+ _* ?: ~) s! t% ]" \88 M3 J7 f. D- h- e; c: q5 o9 q
    9# y6 w1 r8 S- {# @% N
    10
      F+ I/ _" c4 L" n11
    5 x0 x: q9 T& |1 W9 V12
    . C3 p; ]( t: o) U- Y2 ~13+ `" S3 ]+ r% F/ P, k4 v
    14
    " R) v, @1 \0 y0 Y' j15
    $ ~" _, n0 _7 q: `161 ~8 N* m; |0 F% D* p. b/ m& X
    17: ?! {% k4 T: k& ^. z! E
    18
    3 e* B# [# j% ?$ Q" |19) [* i; u9 o! _( W: `
    20
    5 U6 [! O$ B% z9 k5 C21
    8 I  _+ ]9 L! `22
    3 C  N$ o7 y$ u2 M) G( o" I8 f23
    3 }% Y  Q% E; z8 Y1 {24
    ! D3 R7 B% p6 A' c5 `: _) J25% e5 A, z) R6 k7 ~- p+ o
    26
    & M% j$ \  ?2 \8 g  ?278 z$ A3 |& @; R4 ]* p
    28
    5 `) G# |$ g) I5 j2 _7 w( H29
    8 z' ~7 s3 Q* ^, W( l1 G9 i
    # g- l, }" \1 s, l
      b. }8 N+ i; b$ W# F& c
    3.4 解析解与数值解的比较* H) R% M7 t- y  P4 g8 n
    3 R8 Y6 _5 w4 X* t+ L

    7 B  y5 d3 s' o3 _2 b, C  v" `本图为例程 2.3 的运行结果,图中对解析解(蓝色)与使用 odeint() 得到的数值解(红色)进行比较。在该例中,无法观察到解析解与数值解的差异,表明数值解的误差很小。- C0 x/ W9 |1 M7 h
    $ Z$ H8 A+ M2 s1 Q) i

    ( W; m" o" `6 p3 [) W" M图中 d i / d t di/dtdi/dt 具有最大值,最大值表示疫情增长的高潮,达到最大值后 d i / d t di/dtdi/dt 逐渐减小,但患病者比例很快增长到 100%,表明所有人都被感染成为患者。+ d0 i. \7 B( G! X
    , X0 V% S" o# [/ F* K: `
    ) P. R( U+ ?) m! f% _, G$ X
    这是特定参数的结果,还是模型的必然趋势,需要对参数的影响进行更详细的研究。+ l  m' U: c$ W) i7 o2 u4 G

    3 M! Z# ?3 |! l# M9 C1 r" h; {
    2 |: c  F3 _! M- w9 o1 c* @! S$ ?
    0 y& B8 a3 b' r- b9 `! [9 u2 Q
    : O) W& V; e9 N
    4. SI 模型参数的影响
    1 U4 w2 e6 p$ H9 e( t对于 SI 模型,只有日接触率 λ \lambdaλ 和患病者比例的初值 i 0 i_0i ) u$ l- u) B! i. f( K; c
    0+ q: y- f6 Z" |- G% t* F
    ​       
    $ }9 ]& ^, o! r& P0 D: c  会影响模型的结果,其它参数如总人数 N 并没有影响。
    $ X+ @! E9 R! K! L7 D
    " c8 p) y, v4 I6 ~. b, z7 G5 H

    $ M4 D  H8 [' p' }; P( T4.1 日接触率对 SI 模型的影响7 _" x" F/ J2 e6 L% G

    9 e# |5 i# ]3 e6 J
    - ?* s3 O8 k' a! M
    对不同日接触率的比较表明:" }* M6 G2 P. p: I  a  ?

    ' k, Q. s/ Y0 i7 |% Q% D8 E
      w2 H4 X. w$ o- r
    日接触率越大,疫情从发生到爆发的时间越短,爆发过程的增长速度也越快。
    ' L( Q# E5 X, j; u; |5 k" l不论日接触率多大,患病者的比例最终都会增长到 1,表明所有人都被感染成为患者。
    % q# _9 j6 Y; W6 y0 t不论日接触率多大,都具有缓慢发展、爆发、增长放缓 3 个阶段,进入爆发阶段后患病者的比例急剧增长,疫情就很难控制了。
    / U' K* ]4 Z% F& s# R; H" w" f  K- v4 t- I: Y6 Q

    # V. \* d4 D& u" o' o# M" w7 a4.2 患病者比例的初值对 SI 模型的影响/ H' m( n" w' o! v

    2 {3 J4 D* [4 I6 C  s4 `
    0 r/ e1 f+ _1 t! O% W1 K( @' W
    对患病者比例初值的比较表明,患病者初值的人数或比例只影响疫情爆发期到来的快慢,对疫情传播的过程和结果几乎没有影响。3 V2 p* H' b& c# }
    6 C6 i; B/ ~/ t
      V  y  \, W+ |2 Q& `% G# ?
    这与我们直观的经验不太一致,一个原因是 SI 模型本身存在不足,另一方面也说明如果对传染病不加控制,即使开始患病人数很少,经过一段时间的传播后也终将会引起爆发。, c* d/ [/ u/ D' @" X, n+ ?
    ; c: \. P$ M' {. A) a! z
    1 O4 _/ |# ]! N  x+ \! |) Z8 I
    4.3 SI 模型结果讨论6 m6 K& @; u( {
    在 i ( t ) = 0.5 ,   I ( t ) = N / 2 i(t)=0.5,\ I(t) = N/2i(t)=0.5, I(t)=N/2 时 $ di/dt$ 达到最大值,病人数目 I ( t ) I(t)I(t) 增加最快。由此可以预报传染病高潮的到来,即为医院的门诊量最大的一天,卫生部门要重点关注。
    # ~- r" W  d, G2 G2 b& \7 K6 ^9 B! Jt m t_mt ) |  X4 F# F8 P$ l
    m
    5 Y7 n3 j8 M8 d. `! a) ?5 f* G* u​        . i1 P% v6 _) g4 i) Y
      与 λ \lambdaλ 成反比。日接触率 λ \lambdaλ 反映卫生水平、防控手段,提高卫生水平、强化防控手段,降低病人的日接触率,可以推迟传染病高潮的到来。
    8 M* v; {7 g" _' `5 m0 ]" v6 Q当 t → ∞ t \to \inftyt→∞ 时 i → 1 i \to 1i→1 ,表明所有人最终都会被传染而变成病人。这完全不符合实际情况,表明该模型太不讲 politics 了,只能适用于美帝国家建模。, `' c& Z+ o+ l& p
    SI 模型非常明显而严重的缺陷,是该模型没有考虑患病者可以治愈,因此只能是健康人患病,而患病者不能恢复健康(甚至也不会死亡,而是不断传播疫情),所以终将全部被传染。
    : a9 o0 Y9 w2 e: M! l2 L/ H————————————————
    5 k4 J( s* ?1 Y1 b- V' c& n版权声明:本文为CSDN博主「youcans」的原创文章,遵循CC 4.0 BY-SA版权协议,转载请附上原文出处链接及本声明。
    & e7 j9 ]* G9 f5 k原文链接:https://blog.csdn.net/youcans/article/details/117740466/ d% h0 [0 g, ^
    1 c" _$ h0 x9 w, U

    * U' I; ]6 Z% |8 c
    zan
    转播转播0 分享淘帖0 分享分享0 收藏收藏0 支持支持0 反对反对0 微信微信
    您需要登录后才可以回帖 登录 | 注册地址

    qq
    收缩
    • 电话咨询

    • 04714969085
    fastpost

    关于我们| 联系我们| 诚征英才| 对外合作| 产品服务| QQ

    手机版|Archiver| |繁體中文 手机客户端  

    蒙公网安备 15010502000194号

    Powered by Discuz! X2.5   © 2001-2013 数学建模网-数学中国 ( 蒙ICP备14002410号-3 蒙BBS备-0002号 )     论坛法律顾问:王兆丰

    GMT+8, 2026-4-12 02:58 , Processed in 0.414476 second(s), 51 queries .

    回顶部