微分方程模型(2)

2020-06-16 23:13

1.模型Ⅰ

这是一个最简单的传染病模型。设时刻t的病人人数x(t)是连续、可微函数,并且每个病人每天有效接触(足以使人致病的接触)的平均人数是常数?。考察t到t??t这段时间内病人人数的增加,于是就有

x(t??t)?x(t)??

x(?t), tx(t??t)?xt()??x(t).

?t再设t?0时,有x0个病人。并对上式取?t?0时的极限,得微分方程 方程(15)的解为

?t x(t)?x0e. (16)

dxdt ??x,x(0?)x0. (15)

模型(16)表明,随着t的增加,病人人数x(t)无限增长,这显然是不符合实际的。

上述建模失败的原因是:

(1)在病人有效接触的人群中,有健康人也有病人,而其中只有健康人才可以被传染为病人。所以在下面改进的模型中必须区别这两种人。

(2)人群的总人数是有限的,不是无限的。并且随着病人人数的增加,健康人的人数在逐渐减少。因此,病人的人数不会无限地增加下去。

2.模型Ⅱ 模型假设:

(1)在疾病传播期内所考察地区的总人数N不变,既不考虑生死,也不考虑迁移。 (2)人群分为易感染者和已感染者两类,以下简称健康者和病人。并记时刻t这两类人在总人数中所占的比例分别为s(t)和i(t)。

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

根据上述假设,每个病人每天可使?s(t)个健康者变为病人。因为病人人数为Ni(t),所以每天共有?Ns(t)i(t)个健康者被感染。于是?Ns(t)i(t)就是病人人数Ni(t)的增加率,即有

Ndidt??Ns(t)it(. ) (17)

31

又因为s(t)?i(t)?1,再记初始时刻(t?0)病人的比例为i0,则

didt. ??i(t)[?1it(i)],?(i0) (18)0方程(18)是Logistic模型,它的解为 i(t)?1?(11i0?1)e??t . (19)

模型(19)和方程(18)的图形见图4-3和图4-4。

由(18)、(19)式及图4-3和图4-4知,第一,当i?这个时刻为

tm???1

12时,?di?达到最大值??,dt?dt?mdi1ln(?1. )i0这时病人增加得最快,可以认为是医院门诊量最大的一天,预示着传染病高潮的到来,?是医疗卫生部门关注的时刻。tm与?成反比,因为日接触率?表示该地区的卫生水平,

越小卫生水平越高。所以改善保健设施,提高卫生水平可以推迟传染病高潮的到来。第二,当t??时,i?1。即所有人终将被传染,全变为病人,这显然不符合实际情况。

其原因是模型中没有考虑到病人可以治愈,人群中的健康者只能变成病人,病人不会再变成健康者。

为了修正上述结果必须重新考虑模型的假设,在下面的两个模型中将讨论病人可以治愈的情况。

3.模型Ⅲ

32

有些传染病如伤风、痢疾等愈后免疫力很低,可以假定无免疫性。于是病人被治愈后变成健康者,健康者还可以被感染再变成病人。

模型假设:

模型Ⅲ的前三个假设与模型Ⅱ的假设相同。

(4)每天被治愈的病人数占病人总数的比例为常数?,称为日治愈率。病人治愈后成为仍可被感染的健康者,显然1?是这种传染病的平均传染期。

考虑到假设(4),模型Ⅱ中的(17)式应修正为 Ndidt (20) ??Ns(t)it(?)?Nit(. )由于s(t)?i(t)?1,所以(20)式化为 方程(21)的解为

1????(???)t?1[?(?e)]?,??,????i????0 i(t)?? (22)

1?(?t?)?1, ???.?i0?didt??i(t)[?1it(?)?]iti(),?i(00. ) (21)

定义

???? (23) 注意到?和1?的含义,可知?是整个传染期内每个病人有效接触的平均人数,称为接触数。

利用?,(21)式可改写为

didt???i(t)i[t(?)1?(1?) ] (24).

方程(24)和模型(22)的图形见图4-5和图4-6。

33

由图4-5和图4-6不难看出,接触数??1是一个阈值。当??1时,i(t)的增减性取决于i0的大小,但其极限值i(?)?1?1?随着?的增加而增加;当??1时病人比例

最终趋于零,这是由于传染期内健康者变成病人的人数不超过原来病人数i(t)越来越小,的缘故。

4.模型Ⅳ

大多数传染病如天花、流感、肝炎、麻疹等治愈后均有很强的免疫力,所以病愈的人既非健康者(易感染者),也非病人(已感染者),他们已经退出传染系统。这种情况比较复杂,下面将详细分析建模过程。

模型假设:

(1)在疾病传播期内所考察地区的总人数N不变,既不考虑生死,也不考虑迁移。人群分为健康者、病人和病愈免疫的移出者,三类人在总人数N中占的比例分别记作

s(t),i(t)和r(t)。

(2)病人的日接触率为常数?,日治愈率为常数?,传染期接触数为????。 由假设(1)知

s(t)?i(t)? (25) r(t?) 1由假设(2)知方程(20)仍成立。对于病愈免疫的移出者而言应有 Ndrdt??Ni(t). (26)

再记初始时刻的健康者和病人的比例分别是s0(s0?0)和i0(io?0),且不妨假设移出者的初始值r0?0,则由(20)、(25)和(26)式得

34

?di??s(t)it(?)?it(i), ?(t0),0??dt ? (27)

ds????s(t)it()s, (?0s0) . ??dt(27)式即为我们所要建立的数学模型。由于方程(27)无法求出s(t)和i(t)的解析解,因此只能采用数值计算,具体应用时,可使用数学软件来完成。

在模型Ⅳ中????是一个重要参数,由于方程(27)无解析解,因此?、?都很难估计。而当一次传染病结束后,可以获得s0和s?,这时可采用下式对?进行估计

??lns0?lns?s0?s?. (28)

当同样的传染病到来时,如果估计?、?没有多大变化,那么就可以用上面得到的?分析这次传染病的蔓延过程。

关于对模型Ⅳ进行相轨线分析的问题,请有兴趣的同学参阅姜启源、谢金星、叶俊主编的《数学模型》第三版139-144页。

35


微分方程模型(2).doc 将本文的Word文档下载到电脑 下载失败或者文档不完整,请联系客服人员解决!

下一篇:普通话课文60篇

相关阅读
本类排行
× 注册会员免费下载(下载后可以自由复制和排版)

马上注册会员

注:下载文档有可能“只有目录或者内容不全”等情况,请下载之前注意辨别,如果您已付费且无法下载或内容有问题,请联系我们协助你处理。
微信: QQ: