传染病模型[行稳教育]

1、微分方程建模微分方程建模 传染病模型传染病模型 1基本课堂 传染病模型传染病模型 目的目的 描述传染病的传播过程描述传染病的传播过程 分析受感染人数的变化规律分析受感染人数的变化规律 预报传染病高潮到来的时刻预报传染病高潮到来的时刻 预防传染病蔓延的手段预防传染病蔓延的手段 按照传播过程的一般规律，按照传播过程的一般规律， 用机理分析方法建立模型用机理分析方法建立模型 2基本课堂 一、微分方程建模一、微分方程建模 在研究实际问题时，常常会涉及到某些变 量的变化率或导数问题，这样所得到变量 之间的关系式就是微分方程模型。微分方 程模型反映的是变量之间的间接关系，因 此，要得到直接关系，就得求解微

2、分方程。 求解微分方程有三种方法： 1）求精确解；2）求数值解（近似解）；3）定性 理论方法。 3基本课堂 建立微分方程模型的方法 （1）根据规律列方程 利用数学、力学、物理、化学等学科中的定理或 经过实验检验的规律等来建立微分方程模型。 （2）微元分析法 利用已知的定理与规律寻找微元之间的关系式， 与第一种方法不同的是对微元而不是直接对函 数及其导数应用规律。 4基本课堂 （3）模拟近似法 在生物、经济等学科的实际问题中，许多现象 的规律性不很清楚，即使有所了解也是极其复 杂的，建模时在不同的假设下去模拟实际的现 象，建立能近似反映问题的微分方程，然后从 数学上求解或分析所建方程及其解的性质

3、，再 去同实际情况对比，检验此模型能否刻画、模 拟某些实际现象。 5基本课堂 二、问题重述二、问题重述 问题: 有一种传染病（如SARS、甲型H1N1）正 在流行。现在希望建立适当的数学模型，利用已 经掌握的一些数据资料对该传染病进行有效地研 究，以期对其传播蔓延进行必要的控制，减少人 民生命财产的损失。考虑如下的几个问题，建立 适当的数学模型，并进行一定的比较分析和评价 展望。 1、不考虑环境的限制，设单位时间内感染人数的 增长率是常数，建立模型求t时刻的感染人数。 2、假设环境条件下所允许的最大可感染人数 为 。单位时间内感染人数的增长率是感染人数 的线性函数，最大感染时的增长率为零。建立

4、模 型求t时刻的感染人数。 m x 6基本课堂 3、现有卫生防疫部门采集到的某地区一定 时间内一定间隔区间的感染人数数据（见 下表），利用该数据确定上述两个模型中 的相关参数，并将它们的预测值与实际数 据进行比较分析（计算仿真偏差）并对两 个模型进行适当的评价。（注：该问题中， 设最大可感染人数为2000人） 7基本课堂 4、假设总人口可分为传染病患者和易感染者，易 感染者因与患病者接触而得病，而患病者会因治愈 而减少且对该传染病具有很强的免疫功能，建立模 型分析t时刻患病者与易感染者的关系，并对传染 情况（如流行趋势，是否最终消灭）进行预测。 8基本课堂 三、问题分析三、问题分析 1、这是一

5、个涉及传染病传播情况的实际问题，其 中涉及传染病感染人数随时间的变化情况及一些 初始资料，可通过建立相应的微分方程模型加以 解决。 2、问题表述中已给出了各子问题的一些相应的假 设。 3、在实际中，感染人数是离散变量，不具有连续 可微性，不利于建立微分方程模型。但由于短时 间内改变的是少数人口，这种变化与整体人口相 比是微小的。因此，为了利用数学工具建立微分 方程模型，我们还需要一个基本假设：感染人数 是时间的连续可微函数。 9基本课堂 三、问题求解三、问题求解 3.1、问题、问题1的解答的解答模型一模型一 A、模型假设、模型假设 1)、感染人数是时间的连续可微函数； 2)、单位时间内感染人数

6、的增长率是常数，或单 位时间内感染人数的增长量与当时的感染人数成 正比。 10基本课堂 B、模型构成、模型构成 设t时刻的感染人数为 ，初始时刻( ) 的感染者人数为 ，感染者的增长率为r， 根据单位时间内感染人数的增长率是常数 的假设，t到 时间内感染人数的增量 为： 因此， 满足如下的微分方程： ( )x t 0 x tt ()( )( )x ttx trx tt ( )x t 0t 0 , (0) dx rx dt xx 11基本课堂 C、模型求解、模型求解 这是一个线性常系数微分方程，容易求得其解为： D、模型分析、模型分析 由上述解的形式，可以看出，感染人数将随着时 间的增长按指数规

7、律无限增长。特别地，当时间 趋向于无穷时，感染人数也将趋向于无穷大。这 显然是不符合现实的，说明该模型不可能用于传 染病的长期预报，同时也说明迫切需要对该模型 进行必要的修正。 00 ( )(1) rtt x tx exr 12基本课堂 E、改进方向、改进方向 单位时间内感染人数的增长率不是常数， 而是逐渐下降的。原因：感染人数增长到 一定数量后，环境条件、人口总数等因素 将对感染者数量的增长起阻滞作用，且阻 滞作用随感染者数量增加而变大。增长率 是感染人数的减函数：感染者越多，增长 率越低。 13基本课堂 3.2、问题、问题2的解答的解答模型二模型二 A、模型假设、模型假设 1)、感染人数是

8、时间的连续可微函数； 2)、感染人数受环境条件的限制，有一个最 大的可感染人数 。 3)、单位时间内感染人数的增长率和感染人 数有关，是其线性函数，最大感染时对应 增长率为零。 m x 14基本课堂 B、模型构成、模型构成 仍然设t时刻的感染人数为 ，初始时刻 ( )的感染者人数为 ，感染者人数为0 时，感染人数的增长率为 。根据单位时间 内感染人数的增长率和感染人数有关，是其 线性函数的假设，可得增长率关于感染者人 数的线性函数关系式： ( )x t 0 x0t 0 r 0 ( )r xrkx 15基本课堂 进一步，由最大感染时对应的增长率为零 可确定参数k的值为： 因此，在该模型的假设下，

9、感染人数 应满足如下的微分方程： 0 m r k x ( )x t 0 0 ( )(1) , (0) m dxx r x xrx dtx xx 16基本课堂 C、模型求解、模型求解 这是一个非线性微分方程，利用微分方程中的分 离变量法，求得其解为： 0 0 ( ) 11 m r t m x x t x e x 17基本课堂 D、模型分析、模型分析 a)、根据前述微分方程作 出dx/dtx的曲线图，见图 1-1，这是一条抛物线。由 该图可看出感染人数增长 率随感染人数的变化规律： 增长率随着感染人数的增 加而先增后减，在xm/2时 达到最大。这预示着传染 病高潮的到来，是医疗卫 生部门关注和需要密切注 意的时刻。因为感染人数 增长率在一定程度上代表 了医疗卫生水平，增长率 越小卫生水平越高。所以 改善保健设施、提高卫生 水平可以推迟传染病高潮 的到来。 18基本课堂 b)、根据模型求解得 到的结果作出xt曲线， 见图1-2，这是一条S 型曲线。由该图可看 出感染人数随时间的 变化规律：可以看出， 当时间趋于无穷时， x(t)趋于xm，且对一切 t， x(t)1/ x(t)先升后降至先升后降至0 P2: y01/ x(t)单调降至单调降至0 1/ 阈值阈值 P3 P4 P2 y0 00 0 1