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

最新推荐文章于 2024-10-14 23:34:35 发布

请叫我God磊

最新推荐文章于 2024-10-14 23:34:35 发布

阅读量8.2w

点赞数 62

文章标签：数学建模 SIRS传染病模型 MATLAB

本文链接：https://blog.csdn.net/qq_37694821/article/details/86578645

版权

模型假设

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）