- 在线时间
- 791 小时
- 最后登录
- 2022-11-28
- 注册时间
- 2017-6-12
- 听众数
- 15
- 收听数
- 0
- 能力
- 120 分
- 体力
- 36484 点
- 威望
- 11 点
- 阅读权限
- 255
- 积分
- 13906
- 相册
- 0
- 日志
- 0
- 记录
- 1
- 帖子
- 616
- 主题
- 542
- 精华
- 12
- 分享
- 0
- 好友
- 225
TA的每日心情 | 开心 2020-11-14 17:15 |
|---|
签到天数: 74 天 [LV.6]常住居民II
 群组: 2019美赛冲刺课程 群组: 站长地区赛培训 群组: 2019考研数学 桃子老师 群组: 2018教师培训(呼伦贝 群组: 2019考研数学 站长系列 |
本题要建立传染病模型。上网一查,最经典的传染病模型是SIR模型,出自1760年伯努利家族的丹尼尔.伯努利对天花传播规律的研究。8 Y) G4 L B$ F. H2 `, Y1 _
; L# a! _/ ]) O/ y/ g本题主要使用微分方程进行建模。
" Y5 W+ F. H- ~+ C0 _/ l' F6 ^, _, Z
8 ]! _$ R7 c" {. p7 e(一)梳理题目
3 q _; n$ i6 ` E/ L0 |! S! S: E1 M8 P4 D( [5 U: H: k* f0 v
![]()
2 Z5 |: S% x0 u2 p( x) u3 h& g9 K% D8 I" @- Y$ v- L/ w
9 _' X8 Q; W8 i" Y
7 M% \* }; U7 p5 i9 o( M0 F
# q: i4 i! b# W7 F# I7 S) m Z3 H
(二)Highlights which makes this paper stands out9 Z8 C3 J4 Q) h. C5 v9 o
(1)对早期模型拟合曲线的残差分析
2 x. o/ P/ y/ k) r. V拟合模型一定要用残差分析绘制残差图来分析拟合效果。比只是看图说话好。7 z2 S7 d5 b k% M0 j# W$ A
' i4 J2 w$ O: L/ M1 g 0 G" @2 M0 N, N( e1 E8 R
% [1 @1 N+ W8 I8 o" n: q
e i是第i天的计算值和实际值的残差& [! ~) ^5 ~& B' }
e∗i e_i^*e i∗是减去期望E(ei)=0 E(e_i)=0E(e i)=0,再除以残差的标准差得到的标准化残差
- Q. k- ?/ v0 J* U; w, h; I& g标准化残差服从标准正态分布
& a! o" G; m) }5 `1 X美中不足的是!!!- Z9 p/ B2 e# W* ~2 X- A6 A
没有解释为什么用这个式子作为残差的标准差的估计值。。。一般情况下,样本标准差的无偏估计应该是:% M. H6 }2 ]* f9 M/ ~
( Z0 u+ H& x8 u2 m$ M/ u6 _, [+ l
& p: s5 ?4 i6 j$ H: p* i8 f+ Y# @: i. q( ]7 y3 c
; v. H$ ]! P( L+ A# n- S如读者朋友知道原因,请评论告知,非常感谢! f+ D" q/ q% H2 l& p: S6 M+ ^
& \* O* X4 [ w0 A![]()
" b5 S& d8 x- j+ x
& q" H$ m; I5 G$ @" [: b
5 `& ]7 p. L- M5 r0 a论文绘制的残差图表明早期模型只有前期拟合效果较好,中后期都与实际情况偏离较大。4 E+ y* S, I$ h8 h. p9 Y; U% F) H
' p3 K$ ^, ~& d$ |$ T(2)模型假设和符号定义( X& }2 C5 k0 a K" E7 q! i2 x0 I
这个假设写得简直太数学太专业了,为后面用微分方程建模埋下了十足的伏笔啊。
+ M% w$ c8 J3 ?( `) R$ Q
- M4 m( ^4 I! i" _: f / \- O* H7 C% c
. T# v6 U) ]8 ~- ?& l
- ^6 W/ j- _6 Q& V u" m这6个关键变量的找出,是不容易的。
$ k3 X ~# q/ Y) m, a9 l7 \5 N# ]7 g( p
![]()
6 N' b p; c: W4 P6 h. U7 S V+ B
8 \& b8 r- A4 ~(2)基于SIR模型建立新模型
0 t, W. w z! @ f基于一个经典模型,成功率较高,又有更多可参考的资料。; s$ f4 i$ W2 L
SIR简单地把一个城市的人口分为三类,三类的状态转移图精准地刻画了传染病的传播过程。" e& n8 |; k0 p+ |% v% w& ?
! s, L2 u. M3 Y( N, h4 P , l, t( U+ \% ]! M: }, m. ]
& A! H" |/ M; a7 _) Y8 V! i4 r利用微分方程组建立数学模型,这也是对上图的数学描述:
) K' H2 t" H- C7 G3 l: a3 ^8 T6 a/ o
8 k& k9 {- ]+ v+ T+ ^
9 I" {' s7 [* _1 B1 N,因为S类(易感类,能被感染的人群)随疫情发展减少。) x7 O- }/ b$ J" y$ g' T$ W; L9 \
其它数学公式论文中很清晰
5 j2 L8 Z+ C6 ]. J7 b" y; E) {6 F7 [, R3 Z
![]()
8 q# Q' C h/ Y% c! w y6 V9 _
5 ?1 v( D+ ^1 S- b) j+ T% y- K* }' \3 |# D: [( D: J1 s
(3)求解模型. |4 H0 }9 s+ e M- F6 O
求解可以说是很考验数学功底了。深入挖掘模型中方程的关系和隐含信息。
* _; H* j& u+ O; p6 e0 X# F1 z+ P. r9 M2 Z+ G
- z+ O' b3 _4 w% f ~
( O2 O+ T1 B4 G8 A; p
9 {% [4 U' J# X![]()
+ }7 ]1 _5 r- [5 C% Z3 n/ Z y# g) }; G @; M
然后根据实际数据就得到了σ 必须小于1的结论:$ d m# L2 Z X4 p6 R% c
/ Z1 p, L& b0 e4 g: b% G) w
6 v6 o5 h" M' H$ j6 E ^* F
2 k6 ~2 f/ m* n# m6 G
(4)用导数为0划分疫情发展的四个阶段" G3 r0 Q8 J# G m
, P4 S. B/ j3 C: k& y- B( J7 p q . y8 H8 T3 v3 s( a0 V
6 I+ m6 i) C1 [* v
![]()
$ _+ P8 P- @, t! W0 s/ y6 c8 y6 w9 [- V' d' f6 o7 T2 `
(5)根据实际设计三个关键函数; @ w1 `) v" Y( n2 o
这才是体现智商和拉开区分度的重要赛点!!!前面那些都是小亮点,这个是闪瞎眼睛的关键。
6 `+ k; Q$ X& V0 {/ e. K论文也说了,疫情的发展要分阶段研究,各个参数在不同阶段的取值和变化规律(函数)是不同的,所以用分段函数来描述是符合实际的。
0 Q0 B. f) s0 ~9 @: l& P+ N ^/ Q
平均传染期函数:8 T: _9 v0 J5 J# x" P0 r9 E, W6 Q
6 Q9 V# T" J' \2 s% U+ v# l![]()
4 K7 u$ s% ^4 r0 g3 q: m) b
+ ?* ]6 g0 x0 v/ X( A就诊率函数:6 Q/ r* J$ }. r4 E6 Q+ B3 f( |
7 ?. W7 f5 z, V- f' a7 W6 y![]()
N+ _" u9 Y2 j& {- O x, m9 Q, r8 P+ X0 n3 {
平均接触率函数:
% ]' O# M9 ~/ x, g( Y3 p/ o# h# I5 ~) H0 p
A7 p1 Q5 l+ r8 L- B3 I1 n% f3 B
2 L$ Z) y/ A, R% P5 P/ |
模型预测效果图:
) O' J6 w( x6 i' t* y! P5 h3 e8 x ?+ ~% t
+ M- @6 k1 r0 q8 b
6 l% N R% v$ }) W4 Q
————————————————/ v! X- H+ X, v% S* F0 Y
版权声明:本文为CSDN博主「doubleslow;」的原创文章,遵循CC 4.0 BY-SA版权协议,转载请附上原文出处链接及本声明。
3 [; [3 k: q/ W2 C' }, j; y原文链接:https://blog.csdn.net/qq_36607894/article/details/922469472 `. w# J$ D9 W3 G+ ^+ a9 g1 }
" x! L5 ^6 y9 j6 V- M b- g- N; d# U' P& Z6 C
|
zan
|