- 在线时间
- 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年伯努利家族的丹尼尔.伯努利对天花传播规律的研究。
- G9 ~( {" `" r& g( C8 x% o$ N( j0 U" U( L7 Z: G
本题主要使用微分方程进行建模。4 B" u/ j3 E6 s$ |8 a
- \5 N- o, K) n(一)梳理题目9 [. B& V) k q6 l( B% e/ a
& ?1 Y6 p0 M; r6 j2 B
![]()
% F4 a% J- x$ \7 w" I
0 J1 S6 `8 a# h% E% s" p: `: k5 t( O/ q! _6 E# y, n6 Q, O
4 I6 O5 n/ \+ l; o
) y- O* s8 t/ s) Q3 _9 q9 d
(二)Highlights which makes this paper stands out! ?. w4 @# U% |4 r, q
(1)对早期模型拟合曲线的残差分析 E8 H, Y: @+ \, [* D9 l
拟合模型一定要用残差分析绘制残差图来分析拟合效果。比只是看图说话好。" S( c2 Q, t @& i# c( P
e( n1 u. H1 T0 h% y
![]()
& F* j; N5 \2 j% ], e% R, [7 Y& N* Y7 |+ w, @5 Z
e i是第i天的计算值和实际值的残差
6 R' y8 Q# K$ s' u* \' j/ ye∗i e_i^*e i∗是减去期望E(ei)=0 E(e_i)=0E(e i)=0,再除以残差的标准差得到的标准化残差- o- v/ o+ e3 ?9 X' z# k9 t
标准化残差服从标准正态分布
0 i' h) k" }& a; m. }美中不足的是!!!
/ N8 Q; K, b% g7 E- r% C0 ?: X# S没有解释为什么用这个式子作为残差的标准差的估计值。。。一般情况下,样本标准差的无偏估计应该是:1 Q5 E: T) h W$ ], f4 F9 h9 g$ r
0 n7 `- p9 ^( p
, N/ c9 o+ g2 L" P& ?8 `
8 D, z$ Q5 ~$ k9 f5 {6 k! `5 m% L
如读者朋友知道原因,请评论告知,非常感谢
" X% @! n) F8 N
) E2 Z$ u4 v% m9 A( i ' G8 p* A, _7 U y% _/ T
* ^" ~- C$ u4 M5 u' Z0 L8 y" Q
& p# r& c% s/ S4 M( E
论文绘制的残差图表明早期模型只有前期拟合效果较好,中后期都与实际情况偏离较大。. e. I, I: a2 X' S4 A4 i
5 G- J( U. n3 b( ]8 L1 L: e(2)模型假设和符号定义) q' {% ?" [- i) |6 ~! I3 g
这个假设写得简直太数学太专业了,为后面用微分方程建模埋下了十足的伏笔啊。
8 ~1 D+ Q& h$ S% H: J) d$ J$ A7 z) Z# w- ?7 v
![]()
" O/ e4 d! G- N# @0 B
: Y! L3 N! T* _& }! f0 f( T9 V9 d8 s" B
这6个关键变量的找出,是不容易的。
5 {, t3 z' T$ F& c$ |* f" p# |/ I2 N; O
![]()
. c3 C: s- y8 u+ }8 x+ X
. ^5 t* M" g7 C+ V5 @9 T$ ^. @(2)基于SIR模型建立新模型& J% m# j! j q
基于一个经典模型,成功率较高,又有更多可参考的资料。( e3 U3 ?1 P/ J5 G$ m/ G& C" \! Y
SIR简单地把一个城市的人口分为三类,三类的状态转移图精准地刻画了传染病的传播过程。
6 S5 b3 I) E" y: r/ I$ M
$ ~; Q: N( M8 ^ p+ O9 q9 ` , i; L9 h* f/ Y& H$ F# x- U9 P, M
|" Z$ J+ q3 L1 X9 W( g4 `! L
利用微分方程组建立数学模型,这也是对上图的数学描述:
, Q9 M% M9 i' F' B- J2 F6 }( m k* m3 l9 c8 }6 ]8 o
( C n; J5 m% j {
0 m+ a. E7 n, V. \
,因为S类(易感类,能被感染的人群)随疫情发展减少。& G; B& S# ?" g( k5 e6 Q
其它数学公式论文中很清晰4 a0 l% K" }, n" J! [3 A: q/ {' s
/ D2 r# }1 e: Y% ^; s/ w & @( |! Z% P T* ?3 X
1 |3 D: @0 p4 @: k/ o2 K2 ]; R; n6 ?, V& r6 e9 y: ?
(3)求解模型8 H+ n1 G0 C4 G5 c
求解可以说是很考验数学功底了。深入挖掘模型中方程的关系和隐含信息。1 ~0 }4 {0 H8 v
5 f+ v/ y D |9 M
$ j( b- u$ L; F7 |4 m/ R0 ]. |( l* |/ e) X& W& e+ r
4 }! \* ]+ x# \2 x$ Q
![]()
3 R. C% c4 C& \1 c: D) A6 X- D% A5 l v1 l8 u
然后根据实际数据就得到了σ 必须小于1的结论:+ w3 F. o2 _/ |
5 s/ t. H+ Y5 F6 s L
6 x/ F- Y$ B/ S0 e5 n" u! H4 N
6 _; Y0 l- U! Q+ o7 ?& v' v+ V(4)用导数为0划分疫情发展的四个阶段1 k4 G1 b0 b2 H# M* R {" }
+ k/ M) Q- ^1 Y0 n$ \/ b * ?3 \; s; [9 f3 o. W& P. n9 K
" b) ~8 S E* E# E: M2 y![]()
. g* g$ t4 p& J% b9 k2 M
6 @2 S+ o: ?& x z/ I" ~& e(5)根据实际设计三个关键函数
0 p1 e9 ^. _" O: x6 }4 B z这才是体现智商和拉开区分度的重要赛点!!!前面那些都是小亮点,这个是闪瞎眼睛的关键。
) v8 f9 y- S% d8 F论文也说了,疫情的发展要分阶段研究,各个参数在不同阶段的取值和变化规律(函数)是不同的,所以用分段函数来描述是符合实际的。! q8 r8 [/ @% v$ W# a/ b# z
/ r- a, F. C$ @; q, Q! b平均传染期函数:
; e* v& d6 j( I' d- a% \( \0 R; ~
9 {! F$ o4 K5 L! j9 X/ I) W 8 s0 S2 E6 u" R, C. S
# G9 d7 p" i/ O
就诊率函数:
: Y4 W& R( B$ @9 K" K9 n- `* `. p# r( W+ c% b4 B
![]()
( @& i3 O$ u, k5 f4 m6 W5 u7 h
5 Y% D' m8 X: h0 s. ~4 }/ s平均接触率函数:! H E. J/ c) e7 ]
# `3 h4 {; V. `# ^8 a$ A
![]()
* V* M' ^* o8 M z
- j/ L: I& `7 w/ d: V模型预测效果图:
! P1 T& k# g" W( r
7 ^% ~( U) w5 K0 w" L![]()
3 B7 G* Z1 ~$ F8 `8 P
& A7 S) @$ F/ X& v# |! E# ~————————————————* M! y; }3 T! `2 y* Z
版权声明:本文为CSDN博主「doubleslow;」的原创文章,遵循CC 4.0 BY-SA版权协议,转载请附上原文出处链接及本声明。
6 U4 t: Z& h' v1 E8 n7 z原文链接:https://blog.csdn.net/qq_36607894/article/details/92246947
4 `* L! ]" S' L. g2 t
( D7 X' V j0 a" {6 t
! n+ F. U* a$ p; E; P0 M/ F |
zan
|