第5章微分方程模型2.ppt

5.6 人口预测和控制 年龄分布对于人口预测的重要性 只考虑自然出生与死亡，不计迁移 人口发展方程 人口发展方程 一阶偏微分方程 人口发展方程 ~已知函数（人口调查） ~生育率（控制人口手段） 0 t r 生育率的分解 ?~总和生育率 h~生育模式 0 * * 第五章 微分方程模型(续) 动态模型 描述对象特征随时间(空间)的演变过程 分析对象特征的变化规律 预报对象特征的未来性态 研究控制对象特征的手段 根据函数及其变化率之间的关系确定函数 微分方程建模 根据建模目的和问题分析作出简化假设 按照内在规律或用类比法建立微分方程 5.1 传染病模型 问题 描述传染病的传播过程 分析受感染人数的变化规律 预报传染病高潮到来的时刻 预防传染病蔓延的手段 按照传播过程的一般规律，用机理分析方法建立模型 已感染人数 (病人) i(t) 每个病人每天有效接触(足以使人致病)人数为? 模型1 假设 若有效接触的是病人，则不能使病人数增加 必须区分已感染者(病人)和未感染者(健康人) 建模 ? 模型2 区分已感染者(病人)和未感染者(健康人) 假设 1）总人数N不变，病人和健康 人的 比例分别为 2）每个病人每天有效接触人数为?, 且使接触的健康人致病 建模 ? ~ 日 接触率 SI 模型 模型2 1/2 tm i i0 1 0 t tm~传染病高潮到来时刻 ? (日接触率)? ? tm? Logistic 模型 病人可以治愈！ ? t=tm, di/dt 最大 模型3 传染病无免疫性——病人治愈成为健康人，健康人可再次被感染 增加假设 SIS 模型 3）病人每天治愈的比例为? ? ~日治愈率 建模 ? ~ 日接触率 1/? ~感染期 ? ~ 一个感染期内每个病人的有效接触人数，称为接触数。 模型3 i0 i0 接触数? =1 ~ 阈值 感染期内有效接触感染的健康者人数不超过病人数 1-1/? i0 模型2(SI模型)如何看作模型3(SIS模型)的特例 i di/dt 0 1 ? 1 0 t i ? 1 1-1/? i 0 t ? ?1 di/dt 0 模型4 传染病有免疫性——病人治愈后即移出感染系统，称移出者 SIR模型 假设 1）总人数N不变，病人、健康人和移出者的比例分别为 2）病人的日接触率? , 日治愈率?, 接触数 ? = ? / ? 建模 需建立 的两个方程 模型4 SIR模型 无法求出 的解析解 在相平面 上 研究解的性质 模型4 消去dt SIR模型 相轨线 的定义域 相轨线 1 1 s i 0 D 在D内作相轨线 的图形，进行分析 s i 1 0 1 D 模型4 SIR模型 相轨线 及其分析 传染病蔓延 传染病不蔓延 s(t)单调减?相轨线的方向 P1 s0 im P1: s01/? ? i(t)先升后降至0 P2: s01/? ? i(t)单调降至0 1/?~阈值 P3 P4 P2 S0 模型4 SIR模型 预防传染病蔓延的手段 ? (日接触率)? ? 卫生水平? ?(日治愈率)? ? 医疗水平? 传染病不蔓延的条件——s01/? ? 的估计 降低 s0 提高 r0 提高阈值 1/? 降低 ?(=?/?) ? ?, ? ? 群体免疫 模型4 SIR模型 被传染人数的估计 记被传染人数比例 xs0 i 0 P1 ? i0 ?0, s0 ?1 ? 小, s0 ? ?1 提高阈值1/??降低被传染人数比例 x s0 - 1/? = ? 5.4 药物在体内的分布与排除 药物进入机体形成血药浓度(单位体积血液的药物量) 血药浓度需保持在一定范围内——给药方案设计 药物在体内吸收、分布和排除过程 ——药物动力学 建立房室模型——药物动力学的基本步骤 房室——机体的一部分，药物在一个房室内均匀分布(血药浓度为常数)，在房室间按一定规律转移 本节讨论二室模型——中心室(心、肺、肾等)和周边室(四肢、肌肉等) 中心室 周边室 给药 排除 模型假设 中心室(1)和周边室(2),容积不变 药物在房室间转移速率及向体外排除速率，与该室血药浓度成正比 药物从体外进入中心室，在二室间相互转移,从中心室排出体外 模型建立 线性常系数非齐次方程 对应齐次方程通解 模型建立 几种常见的给药方式 1.快速静脉注射 t=0 瞬时注射剂量D0的药物进入中心室,血药浓度立即为D0/V1 给药速率 f0(t) 和初始条件 2.恒速静脉滴注 t T, c1(t)和 c2(t)按指数规律趋于零