SIRS传染病模型求解及MATLAB实现

模型假设

1、总人数N不变。人群分为健康者、病人和移出者三类。t时刻三类人数量分别记为s(t),i(t)和r(t)。

2、病人的日接触率为\lambda,日治愈率为\mu

3、移出者康复后只有暂时免疫力,单位时间内将有\gamma的移出者丧失免疫而可能再次被感染。

模型构成

由假设1显然有 

s(t)+i(t)+r(t)=N          (1-1)

建立关于s(t),i(t)和r(t)的三个方程

\left\{\begin{matrix}i(t+\bigtriangleup t)-i(t)=\lambda s(t)\cdot i(t)\cdot \bigtriangleup (t)-\mu i(t)\cdot \bigtriangleup (t) \\ s(t+\bigtriangleup t)-s(t)=-\lambda s(t)\cdot i(t)\cdot \bigtriangleup (t)+\gamma r(t)\cdot \bigtriangleup (t) \\ r(t+\bigtriangleup t)-r(t)=\mu i(t)\cdot \bigtriangleup (t)-\gamma r(t)\cdot \bigtriangleup (t) \end{matrix}\right.          (1-2)

记初始时刻的健康人、病人和移出者人的比例分别是s_{0}(s_{0}>0)i_{0}(i_{0}>0)r_{0}(r_{0}>0),则SIRS模型的方程可以写作

\left\{\begin{matrix} \frac{\mathrm{d} s}{\mathrm{d} t}=-\lambda si+\gamma r \\\frac{\mathrm{d} i}{\mathrm{d} t}=\lambda si-\mu i \\\frac{\mathrm{d}r }{\mathrm{d} t}=\mu i-\gamma r \end{matrix}\right.           (1-3)

方程(1-3)中无法求出s(t),i(t)和r(t)的解析解,我们先作稳定性分析。

稳定性分析

由于s(t)+i(t)+r(t)=N是一个常数,所以令r=N-s-i,则式(1-3)可以降阶为

\left\{\begin{matrix} \frac{\mathrm{d} s}{\mathrm{d} t}=-\lambda si+\gamma N-\gamma s-\gamma i \\ \frac{\mathrm{d} i}{\mathrm{d} t}=\lambda si-\mu i \end{matrix}\right.          (1-4)

引理1 令

SIR模型是一种常见的传染病传播模型,其中S代表易感者,I代表感染者,R代表康复者。在SIR模型中,假设人群总数为N,易感者数量为S,感染者数量为I,康复者数量为R。传染病的传播是通过易感者和感染者的接触而发生的,传播速率由传染病的基本再生数R0决定。 SIR模型的平衡点是指在传染病传播过程中,传染病的传播速率等于康复速率时,感染者数量不再增加或减少,达到一个稳定的状态。这个状态称为平衡点。 SIR模型的平衡点分析可以通过以下步骤完成: 1. 计算基本再生数R0,R0是传染病传播速率的度量,它表示一个感染者在易感者中传播病毒的数量。如果R0大于1,则传染病会继续传播,如果R0小于1,则传染病会消失。在SIR模型中,R0的计算公式为:R0 = β / γ,其中β是传染率,γ是康复率。 2. 计算平衡点,当感染者数量不再增加或减少时,有:dI/dt = 0,即传染病的传播速率等于康复速率。根据SIR模型的微分方程,可以得到平衡点的计算公式为:I* = (γ / β) * (1 - (1 / R0)), S* = N - I* - R*,其中I*是感染者数量的平衡点,S*是易感者数量的平衡点,R*是康复者数量的平衡点。 3. 分析平衡点的稳定性,当传染病的传播速率小于康复速率时,感染者数量将减少,当传染病的传播速率大于康复速率时,感染者数量将增加。因此,当R0小于1时,平衡点是稳定的,当R0大于1时,平衡点是不稳定的。这是因为当感染者数量略微偏离平衡点时,将会产生一个正反馈机制,导致感染者数量迅速增加或减少,直到达到一个新的平衡点。 以上就是传染病SIR模型的平衡点分析的基本步骤和原理。
评论 29
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值