药物中毒施救模型.docVIP

数学建模 （药物中毒施救模型） 学院：数学科学学院 专业：信息与计算科学 班级：12级信算本班 组别：第六小组 成员：刘慧杰 秦凯 刘家明 学号： 044 050 049  药物中毒施救模型 摘要： 药物都有最大服用剂量，服用过多会导致危险。本文以氨茶碱片为引例，通过建立药物中毒最小剂量模型求得服用氨茶碱片后血药浓度能达到的最大值来确定服用氨茶碱的剂量上限。 关键词：最小剂量 血药浓度 半衰期 正 文 1 问题的提出 氨茶碱片服用过多会使血药浓度（单位血液容积中的药量）过高，当血药浓度达到100μg/ml时，会出现严重中毒，当达到200μg/ml则可致命。血液系统对药物的吸收率和排除率可以由半衰期确定，从药品的说明书可知，药片的剂量为每片100mg，氨茶碱吸收的半衰期约5h，排除的半衰期为6h。现在确定对于孩子（血液总量为2000ml）及成人（血液总量为4000ml）服用氨茶碱能引起严重中毒和致命的最小剂量。 2 合理假设与变量说明 为了判断孩子（或成人）的血药浓度会不会达到危险的水平，需要寻求胃肠道和血液系统中的药量随时间变化的规律。记胃肠道中的药量为想x(t)，血液系统中的药量为y(t)，时间t以及孩子（或成人）误服药的时刻为起点(t=0)。我们可以做以下假设： 1：胃肠道中药物向血液系统的转移率与药量x(t)成正比，比例系数记作λ(λ0)，总剂量x mg的药量在t=0瞬间进入胃肠道,同时x mg也表示所使用药物的最小剂量。 2：血液系统中药物的排除率与药量y(t)成正比，比例系数记作μ(μ0)，t=0时血液系统中无药物。y mg表示血液系统中的药物的含量。 3：氨茶碱被吸收的半衰期约5h，排除的半衰期为6h。 4：孩子的血液总量为2000ml，成人的血液总量为4000ml。 3 模型建立 根据假设对胃肠道中药量x(t)和血液系统中y(t)先建立一个小孩出现药物严重中毒的最小剂量模型。（其余的情况类似） 由假设1，x(0)=x mg，随着药物从胃肠道向血液系统的转移，x(t)下降的速度与x(t)本身成正比（比例系数λ0），所以x(t)满足微分方程 (1) 由（1）式分离变量化简得： (2) 由题意可得：（假设3） (3) (4) 由（2）（3）（4）化简得： (5) 由假设2，y(0)=0，药物从胃肠道向血液系统的转移相当于血液系统对药物的吸收，y(t)由于吸收作用而增加的速度是λx，由于排除而减少的速度与y(t)本身成正比（比例系数μ0），所以y(t)满足微分方程 (6) 由（6）式分离变量化简得： (7) 由题意可得：（假设3） (8) (9) 由（5）（7）（8）（9）化简得： (10) 4 模型求解 将（5）（10）式代入（7）式中并化简得 (11) 要求的是小孩出现严重中毒的最小剂量，即求（11）式中的最大值，将（11）式求导得 (12) 令（12）式得0，得 (13) 已知能引起小孩严重中毒的时，血液系统中药物浓度达到100ug/ml，即小孩血液系统中药物含量为200mg。则由（11）式可得 (14) x=497.6640≈5片 根据假设4，同理可得：孩子的血液总量为2000ml,出现致命的血药浓度200μg/ml相当于血液中药量y达到400mg。成人的血液总量为4000ml,出现严重中毒的致命分别相当于血液中药量y达到400mg和800mg。 求得x分别为995,995,1990。 所以孩子服用氨茶碱能引起严重中毒和致命的最小剂量分别是498mg和995mg。 成人服用氨茶碱能引起严重中毒和致命的最小剂量分别是995mg和1990mg。 5：模型检验 为了检验模型的正确性，我们用matlab做出本模型中能引起小孩（以及成人）严重中毒和致命的最小药物剂量的药量变化图像。 Matlab程序代码： hold on t=[0:0.1:24]; plot(6*498*(exp(-0.1155*t)-exp(-0.1386*t))) plot(6*995*(exp(-0.1155*t)-exp(-0.1386*t))) plot(6*1990*(exp(-0.1155*t)-exp(-0.1386*t))) 可以基本确定我们的求解是正确的。 6 模