在线时间 791 小时 最后登录 2022-11-28 注册时间 2017-6-12 听众数 15 收听数 0 能力 120 分 体力 36316 点 威望 11 点 阅读权限 255 积分 13855 相册 0 日志 0 记录 1 帖子 616 主题 542 精华 12 分享 0 好友 225
TA的每日心情 开心 2020-11-14 17:15
签到天数: 74 天
[LV.6]常住居民II
群组 : 2019美赛冲刺课程
群组 : 站长地区赛培训
群组 : 2019考研数学 桃子老师
群组 : 2018教师培训(呼伦贝
群组 : 2019考研数学 站长系列
本题要建立传染病模型。上网一查,最经典的传染病模型是SIR模型,出自1760年伯努利家族的丹尼尔.伯努利对天花传播规律的研究。
3 ~. _! Q: a% e8 N. c
0 x q0 g/ |* L0 ~3 Q# M 本题主要使用微分方程进行建模。$ w2 B. Q% o( O) ]2 i: [
% I% \( a9 _8 g& d- K' @" c. i
(一)梳理题目
9 `5 C5 k; X. V5 A
' n0 d: T+ P7 _9 U1 u ; o8 f" J- a6 P* G1 I
4 G& U' p/ _4 k( X
% T) O7 `* | p; k . F8 r% C$ S' u) p. d
& O. O% t$ x8 M/ R; C
(二)Highlights which makes this paper stands out; l4 z" o1 h6 @
(1)对早期模型拟合曲线的残差分析2 W1 o9 G+ o, M4 G
拟合模型一定要用残差分析绘制残差图来分析拟合效果。比只是看图说话好。) L. n' x4 f" E/ M: n7 @
4 N4 c( J! Q, U) r. ^" E! D ! n& r% w1 ?* D
6 t, }* ^, D: H Y/ R f! w9 t e i是第i天的计算值和实际值的残差0 [0 N/ l3 [: V$ w0 ^
e∗i e_i^*e i∗是减去期望E(ei)=0 E(e_i)=0E(e i)=0,再除以残差的标准差得到的标准化残差
, ]7 U% K! f" m 标准化残差服从标准正态分布1 z$ Z2 @3 I& _; v! D6 ?0 P' J' W
美中不足的是!!!1 F$ N& W& E/ ?6 K8 J$ u
没有解释为什么用这个式子作为残差的标准差的估计值。。。一般情况下,样本标准差的无偏估计应该是:
7 Q$ O; ~+ K5 G& Y5 e
8 J+ n8 I/ w- r/ z6 W
z# D. l3 z9 }/ i
. J$ w" j* Y% `& ? M' w+ R- c! p
如读者朋友知道原因,请评论告知,非常感谢
' G9 I; d+ q$ F! A0 @
+ Q2 J0 W% W4 N) P+ M/ \% f ; u4 k* y }6 }$ t
! ^+ {: [, g; t& ]9 x
; l$ I/ K/ g8 M$ l+ U1 _ 论文绘制的残差图表明早期模型只有前期拟合效果较好,中后期都与实际情况偏离较大。9 a0 I. {% U9 ~$ q/ x; p
* j3 {+ D( N% o) q
(2)模型假设和符号定义; Y1 W" [1 y* h; \3 i
这个假设写得简直太数学太专业了,为后面用微分方程建模埋下了十足的伏笔啊。
8 I- W1 i: g9 L/ u/ H; T0 K' U- q 6 V! a- G& Z7 V5 ]
# n/ ~" R6 B0 D) z
3 }9 _5 E) ~- S- \+ u
( K7 K" B7 Q+ V 这6个关键变量的找出,是不容易的。
+ p( l5 ^* g, {: r& Q- Z5 v 2 Z9 {6 k* {2 A5 l! x/ y# w
0 V9 ?3 Y4 x5 m1 }6 E1 s1 U
' V9 W% ]4 H+ [0 r9 i6 N (2)基于SIR模型建立新模型) K- n0 P @" k3 ]
基于一个经典模型,成功率较高,又有更多可参考的资料。" j" w) Z# r" b+ [. Q7 o; p
SIR简单地把一个城市的人口分为三类,三类的状态转移图精准地刻画了传染病的传播过程。
. O1 P* a. b$ \2 d5 O* h, N1 B e 8 s; e, T+ b8 R: u: x; ~7 a
* K' [1 h1 Y* F. A9 D- f/ c # L& V; L# z- ]! o6 d7 o
利用微分方程组建立数学模型,这也是对上图的数学描述:+ L3 {+ |9 u+ ?
) E7 L8 L. s! M" T4 T, o
/ @+ K9 j/ ]: a& P. e& w
, w$ d4 }/ c# N3 Q% P ,因为S类(易感类,能被感染的人群)随疫情发展减少。
! j3 _( z$ T+ l( t) u 其它数学公式论文中很清晰
6 {& V6 H5 s8 A, y: v o* O 1 u6 V& G$ z( j. u* S' c
) A6 {4 t9 k' X( T8 M0 Y
. m! ^1 E" D& F* f t9 T
# [; Q" e! J5 _$ Y' V (3)求解模型1 I W4 w7 x0 X3 o: ]: F1 c
求解可以说是很考验数学功底了。深入挖掘模型中方程的关系和隐含信息。
9 q1 V m% [% X" K/ _ s# L/ M ! K, L Q" J9 [7 j8 ^" j1 R
+ N# Y2 {/ m6 V5 ]; o6 N
& `& e* \2 h! l6 [3 e) e; s
" P5 C" s( V4 d( R8 f
3 P5 T; P' g7 M. X' w5 y0 i( r2 L
3 q$ V0 D5 L- I0 R# ` 然后根据实际数据就得到了σ 必须小于1的结论:
( C/ s: x5 D, n) r; J6 H5 E 8 Z0 {! |' _8 h1 H/ Q
, I/ P9 l- w: s4 Z% a
! H" L. S, {- u; p) t
(4)用导数为0划分疫情发展的四个阶段# Q4 u* i% [7 ^) j/ M
8 C3 a& ^8 p; K
4 c3 A1 P( _/ c, h) U % o* [- v. b8 Z3 p4 b
# o; x6 ]% @- ?4 G j+ R8 W ! K! U, g& e; W. @7 M/ C0 V9 |6 _
(5)根据实际设计三个关键函数. P' i7 C' ]% _+ B: ~+ x+ i
这才是体现智商和拉开区分度的重要赛点!!!前面那些都是小亮点,这个是闪瞎眼睛的关键。
, |+ A: o% O9 Q8 J5 s/ m/ c 论文也说了,疫情的发展要分阶段研究,各个参数在不同阶段的取值和变化规律(函数)是不同的,所以用分段函数来描述是符合实际的。
}2 Y; H% h3 B' z b% U6 c & h' L9 c+ f; P& \2 w0 |# [
平均传染期函数:+ T, Y. E* s7 W4 o9 P
8 z4 y- U( _' N' J) L5 l
6 U3 [% G* b& g' u, _5 Y8 G
; D" Z/ ?5 i t+ W5 [- n$ \8 Q
就诊率函数:9 Y, R9 }6 M+ e+ }1 J
2 e& U( B$ p. N, T
6 ~ g* L1 \0 G1 S+ u 8 J3 m6 @! ^% @+ F( ~
平均接触率函数:
% R: h$ U: P0 ]) h' v4 [
( }" }; O( W; r: z9 i% ?5 |3 b7 @8 p 5 `9 k/ P( ]8 R7 \: N8 M
( x' T0 X1 R: @6 C! [5 T9 a. I
模型预测效果图:
: u* S. k5 P0 Q6 w
6 ^" O) R) u' c & z5 `6 x) V+ Y
, D, A0 ] y* } ————————————————
% {% ^! o6 w6 {* r) |2 Y3 S 版权声明:本文为CSDN博主「doubleslow;」的原创文章,遵循CC 4.0 BY-SA版权协议,转载请附上原文出处链接及本声明。, l% b. \; j3 O0 @
原文链接:https://blog.csdn.net/qq_36607894/article/details/92246947/ o6 k; N _8 @2 T3 z
3 b* @( C8 R( V0 a& b& e3 j
2 Q9 _6 @9 j2 k: n
zan