微分方程型—传染病.ppt

微分方程模型 May. 05, 2003 a disease that has rocked Asian markets, ruined the tourist trade of an entire region, nearly bankrupted airlines and spread panic through some of the worlds largest countries. 5.1 传染病模型 SIS的解析解 疫情实证分析(Kermack, P143图) SARS疫情的实证分析 与Kermack同样的方法 王铎,赵宵飞.SARS疫情的实证分析和预测[J].北京大学学报(医学版),2003,5(S):72-74. 一句话小结 不同的领域可以共享相同或类似的数学模型，但所关注的问题会有所不同； 不能求得解析解的方程仍可用相轨线办法分析解的性质。 进一步的问题 考虑出生和死亡因素的传染病模型 考虑潜伏期的传染病模型SEIR 考虑被动免疫的传染病模型MSIR 考虑随机接触率的传染病模型SSIR 参考/wiki/Epidemic_model 补充习题 理论证明P143第1~3行 。 在SIR 模型中考虑出生与死亡的因素。假设全体人群以相同出生率生育婴儿，且婴儿为易感人群。死亡率与出生率相等，从而人群总数不变。试建立数学模型描述疾病的流行特征，并分析传染病不蔓延的条件。 5.4 药物在体内的分布与排除（药物动力学之房室模型） 房室系统的概念 二房室模型的建立 模型求解 不同给药方式分析 参数估计技巧 进一步推广 药物动力学之房室系统 复习：常系数齐次线性方程组通解(n=2) 详解 方程非其次项是常数，所以设解非其次项也是常数，令方程的通解为 代入方程，比较两边e-?t, e-?t以及常数项系数得到6式，删去冗余两式得 由后两式解得 由前两式得 再利用初始条件c1(0)=0, c2(0)=0解得(可利用Matlab的符号工具箱) tT以后，静脉注射停止 房室模型建模小结 分析各房室的关联； 建立线性微分方程组模型； 写出微分方程组的通解； 用初始条件和代入方程求得特解； 用观测数据估计模型参数 参数估计可用分解技巧，简化计算，使结果更可靠。 进一步的问题 多房室系统模型 非线性房室模型 随机房室模型 房室模型在其他领域的应用 其他注射方式下的参数估计问题 (思考：恒速静脉滴注情形的参数估计技巧？) 参考阅读 /models.php 周晓芳,陈小全,周鲁,生理房室模型药物动力学的研究进展 , 预防医学情报杂志 2002年06期 陈增敬, 关于血浆中放射性钙C^47α浓度的计算公式, 数理统计与应用概率，1995年 10卷 3期 补充习题3 给出P155-156口服或肌肉注射情形下二房室模型的解。 提示： 补充习题4 在北美的五大湖中，安大略湖处于伊利湖的下游，安大略湖不仅接受伊利湖来的水，还要接受非伊利湖流入的水，已知流入安大略的水有 5/6 是伊利湖流出的。试建模描述这两个湖的污染情况。 假设除去控制不了由伊利湖自安大略湖的流动外， 流入伊利湖和安大略湖的所有污染都 暂时被停止了。试计算把安大略净化到 50%以及 5%所需要的时间。 5.6 人口预测和控制(偏微分方程模型) 人口模型介绍 PDE建模 人口预测 人口控制与计划生育 几个人口发展指数 参考文献 研究人口模型的意义 人口控制 人口系统工程 社会保障 寿险精算 种群生态学 人口模型概述 宏观模型: 总人口, 不考虑年龄， Malthus模型, Logistic模型(第一章) 微观模型：考虑年龄结构 1930s, Lotka 积分方程模型 1940s, Leslie差分方程模型(第七章) 1960s, Verhulst偏微分方程模型 1970s, Pollard随机方程模型 人口PDE建模和预测 补充习题5 验证P164(7)式为P163方程(5)的解 。 胡鞍钢 ：“一对夫妇一个孩儿”该结束了 “一对夫妇一个孩儿”该结束了 调整人口生育政策势在必行 少子化. 妇女总和生育率的过快下降，明显低于正常的人口生育更替水平。在1995年前后我国0至14岁少儿人口绝对数达到了最高峰，大约为3.34亿人，2008年的时候减少到2.52亿人。 老龄化。根据联合国人口署的预测，到2020年我国60岁以上人口将占到总人口的16.7%，2050年将进一步上升到31.1%，大大高于届时的世界平均水平(21.9%)。 劳动人口缺乏. 15至59岁劳动人口大约在2015至2020年之间也会达到最高峰，大约9.23亿人，而后开始持续下降。到2050年，中国的15至59岁劳动年龄人口数大约要比印度少2.44亿人。 “一对夫