SIR传染病模型（微分方程系列1）

一：基本参数

SIR模型是常见的一种描述传染病传播的数学模型，其基本假设是将人群分为以下三类：

S:(Susceptible):易感人群，指未得病者，但缺乏免疫能力，与感病者接触后容易受到感染。

I:(Infective):患病人群，指染上传染病的人，他可以传播给易感人群。

R:(Removed）：移除人群，被移出系统的人，因病愈（具有免疫力）或死亡的人。这部分人不再参与感染和被感染过程。

二：问题假设

1：没有天生携带抗体的人

2：总人数N不变且仅包含三种健康状态，健康人、病人和移出者的比例分别为s ( t ) , i ( t ) , r ( t )

3：病人的日接触率为 $\lambda$ ，日治愈率为 $\mu$ ，接触数 $\sigma = \lambda / \mu$

4：假设给定的数据真实可靠，不考虑人口的自然出生，自然死亡以及其他类型的人口流动

5: 假设严格隔离的新型冠状病毒感染病人都不再传染他人

三：模型建立

具体这个模型是咋写出来的，下面是我的推导过程（手稿）

实在无法理解的话，仅需要把SIR的固定模式牢记即可（上面的黄图)

这里我们设 $\lambda$ =1， $\mu$ =0.5，即每个患者可以传染一个易感染者，治愈率为0.5；

初值条件我们设定的初始易感染者S为0.99，初始患者I为0.01，那为啥没有R（治愈者）的初值条件呢，因为我们已经假设没有天生就携带抗体的人，所以R（0）=0，就不需要写，但在之后的过程中会出现携带抗体的人即治愈者，R=1-S-i

四代码实现

首先建立一微分方程函数

function   y = covid(t,x)    %因为我们的微分方程是X关于t的变化率的体现a=1;      %感染率miu的参数配置b=0.5;    %治愈率lamda的参数配置y=[a*x(1)*x(2)-b*x(1);-a*x(1)*x(2)];
end

主函数

ts=0:40;
x0=[0.01, 0.99];
[t,x]=ode45(@ill,ts,x0);
r=1-x(:,1)-x(:,2);
plot(t,x(:,1),t,x(:,2),ts,r),grid
legend('i(t)','s(t)','r(t)')

运行结果

SIR传染病模型（微分方程系列1）相关推荐

SIR传染模型Matlab代码,sir传染病模型 MATLAB代码运行不了,
问题描述: sir传染病模型 MATLAB代码运行不了, function y=ill(t,x) a=1;b=0.3; y=[a*x(1)*x(2)-b*x(1),-a*x(1)*x(2)]'; ts ...
最近疯传的SIR传染病模型是什么？
总第188篇/张俊红最近看到在网上传的一张SIR传染病模型的图,很多人应该对这个模型不是很了解,今天就讲一下这个模型.这一篇只讲学术,不讨论别的. SIR模型是传染病模型中最经典的一个,类似的还有S ...
SIR传染病模型介绍+R语言简单应用
此文用于整理回顾写论文时看的文献资料和学到的知识,也希望能带来一些参考. 什么是SIR模型? 1927年,Kermack 和 McKendrick 为了研究17世纪肆虐伦敦的黑死病和20世纪席卷孟买的 ...
Python学习：SIR传染病模型
SIR模型是传染病模型中最经典的一个,类似的还有SI和SIS两种.SIR是三个单词首字母的缩写,其中S是Susceptible的缩写,表示易感者:I是Infective的缩写,表示感染者,R是Remo ...
NetLogo学习笔记5 —— 物种与传染病模型
NetLogo学习笔记5 -- 物种与传染病模型 (模型经过一些修改,与标准SIR模型有些出入) 在上一篇文章,我们学习了随机选择.伪并发.ifelse和of语法.实现了用于解释隔离现象的谢林模型这 ...
数学建模清风第一次直播：传染病模型和微分方程拟合
学习网址:课件下载数学建模清风第一次直播:传染病模型和微分方程拟合数学建模清风第二次直播:模拟退火算法数学建模清风第三次直播:excel在数学建模中的应用博客笔记: 数学建模清风第一次直播:传 ...
数模学习第三天--微分方程（传染病模型）
微分方程定义包含关于未知变量.未知变量的导数以及自变量的方程第一类为可分离变量的微分方程. 在高数中我们学过相关内容,如: 建模准备 1.运用已知规律列方程利用数学.力学.物理.化学等学科中的 ...
微分方程（人口预测与传染病模型）
一.定义微分方程:含导数或微分的方程微分方程的阶数:所含导数或微分的最高阶数,如y'''+2y''-2x=0是三阶微分方程微分方程的解:使得微分方程成立的函数例如y'-2x=0的解可以为x²或 ...
matlab求解常微分方程组/传染病模型并绘制SIR曲线
看了很多关于传染病模型的matlab程序,大都是绘制出两条曲线(I.S)的,本文最大的不同是绘出SIR三条曲线. 先给出SIR微分方程组函数文件: run的程序:
传染病模型（1）——SI模型及matlab详解
前言常见的传染病模型按照具体的传染病的特点可分为 SI.SIS.SIR.SIRS.SEIR 模型.其中"S""E""I""R&q ...

SIR传染病模型（微分方程系列1）

SIR传染病模型（微分方程系列1）相关推荐

最新文章

热门文章