微分方程模型.pptVIP

第五章 微分方程模型 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-McKendrick, P145图) SARS疫情的实证分析 与Kermack同样的方法 王铎,赵宵飞.SARS疫情的实证分析和预测[J].北京大学学报(医学版),2003,5(S):72-74. 一句话小结 不同的领域可以共享相同或类似的数学模型，但所关注的问题会有所不同； 不能求得解析解的方程仍可用相轨线办法分析解的性质。 进一步的问题 考虑出生和死亡因素的传染病模型 考虑潜伏期的传染病模型SEIR 考虑被动免疫的传染病模型MSIR 考虑随机接触率的传染病模型SSIR 参考/wiki/Epidemic_model 随机SIS: 参考A. Gray , D. Greenhalgh , L. Hu , X. Mao , and J. Pan, A Stochastic Differential Equation SIS Epidemic Model, SIAM J. Appl. Math. 71 , 2011, pp. 876-902 习题 P180ex1. 理论证明P144表下第1段“可以看出…” 。 在SIR 模型中考虑出生与死亡的因素。假设全体人群以相同出生率生育婴儿，且婴儿为易感人群。死亡率与出生率相等，从而人群总数不变。试建立数学模型描述疾病的流行特征，并分析传染病不蔓延的条件。 5.4 药物在体内的分布与排除（药物动力学之房室模型） 房室系统的概念 二房室模型的建立 模型求解 不同给药方式分析 参数估计技巧 进一步推广 药物动力学之房室系统 复习：常系数齐次线性方程组通解(n=2) 详解 因f0(t)=0, 为齐次方程，通解为 代入方程第一式（第二式也可） 比较两边e-?t, e-?t系数得到 由初始条件得 求解得 详解 方程非其次项是常数，所以设解非其次项也是常数，令方程的通解为 代入方程，比较两边e-?t, e-?t以及常数项系数得到6式，删去冗余两式得 由后两式解得 由前两式得 再利用初始条件c1(0)=0, c2(0)=0解得(可利用Matlab的符号工具箱) tT以后，静脉注射停止 房室模型建模小结 分析各房室的关联； 建立线性微分方程组模型； 写出微分方程组的通解； 用初始条件和代入方程求得特解； 用观测数据估计模型参数 参数估计可用分解技巧，简化计算，使结果更可靠。 进一步的问题 多房室系统模型 非线性房室模型 随机房室模型 房室模型在其他领域的应用 其他注射方式下的参数估计问题 (思考题：恒速静脉滴注情形的参数估计技巧？) 参考阅读 /models.php 周晓芳,陈小全,周鲁,生理房室模型药物动力学的研究进展 , 预防医学情报杂志 2002年06期 陈增敬, 关于血浆中放射性钙C^47α浓度的计算公式, 数理统计与应用概率，1995年 10卷 3期 习题 P180ex5。 P180ex6。 在北美的五大湖中，安大略湖处于伊利湖的下游，安大略湖不仅接受伊利湖来的水，还要接受非伊利湖流入的水，已知流入安大略的水有 5/6 是伊利湖流出的。试建模描述这两个湖的污染情况。进一步，假设除去控制不了由伊利湖自安大略湖的流动外， 流入伊利湖和安大略湖的所有污染都 暂时被停止了。试计算把安大略净化到 50%以及 5%所需要的时间。 5.6 人口预测和控制(偏微分方程模型) 人口模型介绍 ODE建模 PDE建模 人口预测 人口控制与计划生育 几个人口发展指数 参考文献 研究人口模型的意义 人口控制 人口系统工程 社会保障 寿险精算 种群生态学 人口模型概述 宏观模型: 总人口, 不考虑年龄， Malthus模型, Logistic模型（本节ODE） 微观模型：考虑年龄结构 1940‘s, Leslie差分方程模型(第6章) 1960s, Verhulst偏微分方程模型（本节PDE） 1970s, Pollard随机方程模型 (第9章) 中国人口增长情况 美国人口参数x0, r的估计 数据(课本P167，美国人口1790-2000) 线性化拟合(课本P167) 变换lnx=lnx0+rt，解线性方程组得lnx0和r。 短期拟合（1790-1900) r = 0.2743/10年，x17