数学建模社区-数学中国
标题:
传染病传播模型(含matlab代码)
[打印本页]
作者:
杨利霞
时间:
2020-4-18 16:00
标题:
传染病传播模型(含matlab代码)
传染病传播模型(含matlab代码)
! H0 W. I N- L4 G8 T' \
! a' [0 [. M1 H+ J" H# z
2020年初,一种新型冠状病毒袭击了中国武汉。
# C( n$ l1 }/ D1 E6 C* `. A0 i6 e* A
/ D. r$ }5 @. ?0 d
随着卫生设施的改善、医疗水平的提高以及人类文明的不断发展,诸如霍乱、天花等曾经肆虐全球的传染性疾病已经得到有效的控制,但是一些新的、不断变异着的传染病毒却悄悄向人类袭来。长期以来,建立数学模型以描述疾病的传播过程从而寻找抑制传播的方法一直是科学家们探索的方向。但是由于疾病传播的复杂性,仅能从一般的传播机理角度建立大致的传播模型,现介绍一种与实际情况较为接近的传染病传播模型——
# {' N8 D5 j8 P+ Y' R. {
+ S8 t. b4 |# S( Y. w; G
SIR模型
# @3 L% j) z" N, S9 ^' O0 U, k
4 v+ d* O, i' p$ ], Z6 Z
大多数传染病如天花、流感、肝炎、麻疹等治愈后均有很强的免疫力,所以病愈的人既非健康者(易感染者),也非病人(已感染者),他们已经退出传染系统,这种情况比较复杂,下面将详细分析建模过程。
. u9 k: ?8 u) q3 G
8 Z* P4 B6 c1 e5 l
一、基本假设
& m% u( Z) t! D, s. ` y: F
' w; R0 c" [, L7 X# W3 }/ Q8 I
1.在疾病传播期内所考察地区的总人数N不变,既不考虑生死,也不考虑迁移.人群分易感染者(Susceptible)、**已感染者(Infective)和病愈免疫的移出者(Removed)**三类,以下简称健康者,病人和移出者,t时刻人数分别记为S(t),I(t),R(t),t时刻这三类人在总人数中所占的比例分别记作s(t),i(t)和r(t)。
! {- C8 W& q7 l" ^4 D% b0 y
% D3 [ K1 T8 b* d! ], E& F
2.每个病人每天有效接触的平均人数是常数λ,称为日接触率.当病人与健康者有效接触时,使健康者受感染变为病人。
7 o4 e7 @/ t+ T
- d. |5 F( C @
3.每天被治愈的病人占病人总数的比例为常数μ,称为日治愈率,病人被治愈后仍有可能被感染为病人,那么可认为1/μ是该疾病的平均传染期。
$ Z4 R6 _7 w+ p* w
; @% F0 i# q- Y" }( b, }- G
4.初始时刻,只有少数个体处于感染状态,其他都是易染状态。
& ~1 r$ |- D% r
- d; q8 ?7 U- Y9 v! N
5.假设疾病的时间尺度远小于个体生命周期,从而不考虑个体的出生和自然死亡。
- z& P! L# m% c9 p' b3 O( h
5 t; m" z. Z P
6.完全混合(Fully mixed):每一个个体与其他个体接触的机会均等。
3 [ | G0 s+ |, W3 a" Y# i6 D
9 g# x& }' R. Y) U8 m$ s& O
不难看出,考虑以上几条假设,可得:
. X' C0 L) {, E6 G4 U4 }
3 m3 R9 L1 u7 R7 d- D, U2 Y
每个病人每天可使λs(t)个健康者变
. e* Z. x7 H- J6 ?! E
6 u0 r% b* z$ C, T) P2 B5 Z
访问原文:https://blog.csdn.net/wmhsjtu/article/details/104084414?ops_request_misc=%257B%2522request%255Fid%2522%253A%2522158719553519726869009124%2522%252C%2522scm%2522%253A%252220140713.130102334.pc%255Fall.%2522%257D&request_id=158719553519726869009124&biz_id=0&utm_source=distribute.pc_search_result.none-task-blog-2~all~first_rank_v2~rank_v25-14
; R5 p- ^- w' _# Z
————————————————
, c1 x5 b9 Y1 O% d! i/ O$ i
版权声明:本文为CSDN博主「wmhsjtu」的原创文章,遵循CC 4.0 BY-SA版权协议,转载请附上原文出处链接及本声明。
1 n+ y3 v0 [) v) E+ m1 k6 v/ U2 i
原文链接:https://blog.csdn.net/wmhsjtu/article/details/104084414
, J% Y# _6 P* F1 D* H
0 v9 V6 j8 t6 j- W$ B
, T- X3 M8 s3 X3 g7 g0 e, f
欢迎光临 数学建模社区-数学中国 (http://www.madio.net/)
Powered by Discuz! X2.5