第5章 微分方程模型要点.ppt

模型假设 中心室(1)和周边室(2),容积不变. 药物在房室间转移速率及向体外排除速率与该室血药浓度成正比. 药物从体外进入中心室，在二室间相互转移, 从中心室排出体外. 模型建立 中心室 周边室 给药 排除 c1(t), x1(t) V1 c2(t), x2(t) V2 转移 线性常系数非齐次方程 对应齐次方程通解 模型建立 几种常见的给药方式 1.快速静脉注射 t=0 瞬时注射剂量D0的药物进入中心室,血药浓度立即为D0/V1 给药速率 f0(t) 和初始条件 2.恒速静脉滴注 t T, c1(t)和 c2(t)按指数规律趋于零 0?t ?T 药物以速率k0进入中心室 3.口服或肌肉注射 相当于药物( 剂量D0)先进入吸收室，吸收后进入中心室. 吸收室药量x0(t) 吸收室 中心室 D0 参数估计 各种给药方式下的 c1(t), c2(t) 取决于参数k12, k21, k13, V1,V2 t=0快速静脉注射D0 ,在ti(i=1,2,?,n)测得c1(ti) 由较大的 用最小二乘法确定A,? 由较小的 用最小二乘法确定B,? 参数估计 进入中心室的药物全部排除 建立房室模型,研究体内血药浓度变化过程,确定转移速率、排除速率等参数,为制订给药方案提供依据. 机理分析确定模型形式，测试分析估计模型参数. 药物在体内的分布与排除 房室模型： 一室模型 二室模型 多室模型 非线性(一室)模型 c1较小时近似于线性 ~ 一级排除过程 如 c1较大时近似于常数 ~ 零级排除过程 背景 年份 1625 1830 1930 1960 1974 1987 1999 人口(亿) 5 10 20 30 40 50 60 世界人口增长概况 中国人口增长概况 年份 1908 1933 1953 1964 1982 1990 1995 2000 人口(亿) 3.0 4.7 6.0 7.2 10.3 11.3 12.0 13.0 研究人口变化规律 控制人口过快增长 5.6 人口的预测和控制 做出较准确的预报 建立人口数学模型 指数增长模型——马尔萨斯1798年提出 常用的计算公式 x(t) ~时刻t的人口 基本假设 : 人口(相对)增长率 r 是常数 今年人口 x0, 年增长率 r k年后人口 随着时间增加，人口按指数规律无限增长. 与常用公式的一致 rt e x t x 0 ) ( = ? 指数增长模型的应用及局限性 与19世纪以前欧洲一些地区人口统计数据吻合. 适用于19世纪后迁往加拿大的欧洲移民后代. 可用于短期人口增长预测. 不符合19世纪后多数地区人口增长规律. 不能预测较长期的人口增长过程. 19世纪后人口数据 人口增长率r不是常数(逐渐下降) 阻滞增长模型——逻辑斯蒂(Logistic)模型 人口增长到一定数量后，增长率下降的原因： 资源、环境等因素对人口增长的阻滞作用, 且阻滞作用随人口数量增加而变大 假设 r~固有增长率(x很小时) xm~人口容量（资源、环境能容纳的最大数量） r是x的减函数 dx/dt x O xm xm/2 t x O x增加先快后慢 xm x0 xm/2 阻滞增长模型(Logistic模型) 指数增长模型 Logistic 模型的应用 经济领域中的增长规律(耐用消费品的售量). 种群数量模型 (鱼塘中的鱼群, 森林中的树木). S形曲线 参数估计 用指数增长模型或阻滞增长模型作人口 预报，必须先估计模型参数 r 或 r, xm . 模型的参数估计、检验和预报 指数增长模型 阻滞增长模型 由统计数据用线性最小二乘法作参数估计 例：美国人口数据(百万) t 1860 1870 1880 … 1960 1970 1980 1990 2000 x 31.4 38.6 50.2 … 179.3 204.0 226.5 251.4 281.4 r =0.2022/10年，x0 =6.0450 模型的参数估计、检验和预报 指数增长模型 阻滞增长模型 r=0.2557/10年，xm =392.0886 年 实际人口 计算人口 (指数增长模型) 计算人口 (阻滞增长模型) 1790 3.9 6.0 3.9 1800 5.3 7.4 5.0 … … … … 1960 179.3 1