传染病模型试验

简介: 传染病模型试验

一.SIS模型


假设:


(1) t时刻人群分为易感者(占总人数比例的s(t))和已感染者(占总人数比例的i(t))。


(2) 每个病人每天有效接触的平均人数是常数称为日接触率,当健康者与病人接触时,健康者受感染成为病人。


(3) 病人每天被治愈的占病人总数的比例为,称为日治愈率,显然为这种传染病的平均传染期。 则建立微分方程模型为:

image.png



二.分别作图:


当sigma>1时


Step1:先定义函数

function y=pr1(i,lambda,sigma)
y=-lambda.*i.*(i-(1-1./sigma))

step2:作图

lambda=0.3;sigma=2;
i=0:0.01:1;
y=pr1(i,lambda,sigma)
plot(i,y)

image.png


sigma<1时


Step1:先定义函数

function y=pr1(i,lambda,sigma)
y=-lambda.*i.*(i-(1-1./sigma))

step2:作图

lambda=0.3;sigma=0.5;
i=0:0.01:1;
y=pr1(i,lambda,sigma)
plot(i,y)

image.png


sigma=1时


Step1:先定义函数


function y=pr1(i,lambda,sigma)
y=-lambda.*i.*(i-(1-1./sigma))

step2:作图

lambda=0.3;sigma=1;
i=0:0.01:1;
y=pr1(i,lambda,sigma)
plot(i,y)

image.png


sigma>1时


Step1:先定义函数

function di=crb(t,i,lambda,sigma)
di=-lambda*i*(i-(1-1/sigma))


step2:求解并作图

clc
clear
lambda=0.01;
sigma=2;
[t,i1]=ode45(@crb,[0,1000],0.9,[],lambda,sigma);
[t,i2]=ode45(@crb,[0,1000],0.2,[],lambda,sigma);
plot(t,i1,t,i2,t,1/2)
legend('\sigma>1')

image.png

sigma=1时


Step1:先定义函数

function di=crb(t,i,lambda,sigma)
di=-lambda*i*(i-(1-1/sigma))


step2:求解并作图

clc
clear
lambda=0.01;
sigma=1
[t,i1]=ode45(@crb,[0,1000],0.9,[],lambda,sigma);
plot(t,i1)
legend('\sigma=1')

image.pngsigma<1时


Step1:先定义函数

function di=crb(t,i,lambda,sigma)
di=-lambda*i*(i-(1-1/sigma))


step2:求解并作图

clc
clear
lambda=0.01;
sigma=0.5;
[t,i1]=ode45(@crb,[0,1000],0.9,[],lambda,sigma);
plot(t,i1)
legend('\sigma<1')

image.png


目录
打赏
0
0
0
0
181
分享
相关文章
五项举措!为了湖北1300万慢病病人
新冠疫情之外,湖北还有近1300万慢病人群,一旦停药将引发各种断药反应,严重的会危机生命!为此,阿里健康紧急联动各方,推出五项举措,希望与社会各界一道,保护这1300万人的安全。
五项举措!为了湖北1300万慢病病人
膝盖中了一箭之康复篇-一周年纪念
右膝盖前十字韧带在2014年12月份的一次常规踢球中不小心弄断。2015年5月份才确诊。随后2个月后的7月1日做了十字韧带重建术。到现在刚好是术后一年了。 总体来说恢复的情况我还比较满意。
1068 0
森林病虫防治系统 (结束)
明天开始答辩。 今天老师给我们整理了一些知识点。 第一章  java Web 工作原理 重点: 什么是HTTP?(了解) HTTP是一种超文本传送协议,是一套计算机在网络通信的一种规则。 特记:HTML超文本标记语言,以.htm或.html为扩展名, HTTP请求/响应包含哪些内容?(了解) HTTP请求包含了:请求行,请求头,空行与消息体。
976 0