《SI传染病模型》由会员分享,可在线阅读,更多相关《SI传染病模型(5页珍藏版)》请在金锄头文库上搜索。
1、可编辑SISI传染病模型传染病模型1. 模型的建立由题意知道:在此环境中仅存在健康者(即易感者)和已感者(即病人),且在 t 时刻人数分别为 S(t),L(t),不考虑人口的出生与死亡,此环境中的人口数量不变 N 即 K,于是在单位时间内每天每个病人感染的人数S(t)L(t) ,它是病人的增加率,所以有:dL= =*S*St*L*Lt L L0=L1=L1dt (1)在t时 刻 健 康 者 与 已 感 者 满 足 关 系 式 : S St+L+Lt= =(2)此模型满足 Logistic 模型,所以它的解为:L L(t t)=1/1+(1/L1)-1)*exp(-=1/1+(1/L1)-1)*
2、exp(-*t)*t)1.1.求平衡点求平衡点syms r S L K ysyms r S L K yy=r*L*(K-L);y=r*L*(K-L);solve(y)solve(y) ans = ans = 0 0SISSIS 传染病模型传染病模型-可编辑1. 模型假设 SIS 模型的假设条件 1.2 与 SI 模型相同,增加的条件为:每天被治愈的病人数占病人的总数为-可编辑m ,此称为日治愈率。病人治愈后仍然可以成为被感染的健康者,显然,平均传染期为 1/m 。2. 模型建立 此模型可以修整为 : (a 代表)dLt a*St*Ltm*LtLtSt KL0 L1dt求平衡点: (s, l ,
3、k 分别代表 S, L ,K)syms a t s l m k fsyms a t s l m k ff=a*l*(k-l)-m*l;f=a*l*(k-l)-m*l; solve(f) solve(f) ans = ans = -a*(-k+l) -a*(-k+l)1.大于时的图像a,a 10,b 0.8bdi/dt与 i的 变 化 关 系32.5di/dt(病人比例对时间的变化率)21.510.50-0.5-100.10.20.30.40.50.6I(病 人 比 例 )0.70.80.912.2.小于小于 1 1 时的图像时的图像a 0.2,b 0.8-可编辑di/dt与 i的 变 化 关
4、系0data1-0.2-0.4di/dt(病人比例对时间的比率)-0.6-0.8-1-1.2-1.4-1.6-1.8-200.20.40.60.811.2病 人 比 例 i1.41.61.82三三SIRSIR 模型模型模型假设:在 SIS 模型中我们增加:人群可分为健康者,病人,病疫免疫的移出者,且三种人群的数量分别为St,Lt,Rt;病人的日接触率和日治愈率分别为,m 所以传染期为1.模型建立mdLt a*St*Ltm*LtLtSt KL0 L1 (1)dtdSt a*St*LtS0 K L0 (2)dt求平衡点syms a t s l m ks,l=solve(a*l*(k-l)-m*l,
5、-(a*s*(k-s) s = a*k-a*l a*k-a*l l = 0k健康者与病人数量在总人数中的比例st,it对时间的变化关系图为:-可编辑健 康 者 与 病 人 各 占 比 例 随 时 间 的 变 化 关 系10.90.80.70.6病 人 数 量 占 总 人 数 的 比 例 i(t)健 康 者 占 总 人 口 的 比 例 s(t)比例0.50.40.30.20.100510152025时 间3035404550健康者与病人各自占总人数的比例间的相互关系:i-s的 图 形 ( 相 轨 线 )1健 康 者 与 病 人 各 自 占 总 人 数 的 比 例 间 的 相 互 关 系0.90.8健康者所占总人数的比例0.70.60.50.40.30.20.1000.050.10.150.20.25病 人 所 占 总 人 数 的 比 例0.30.350.4 . .-