驾车饮酒的数学模型-Read.DOCVIP

驾车饮酒的数学模型 [摘要] 本文利用药物动力学房室模型，结合微元分析法和微分方程理论，利用MATLAB软件进行曲线拟合，较准确方便的地求出了短时间饮酒、长时间饮酒和天天饮酒情况下，血液中酒精浓度变化关系式，以及达到峰值浓度和新安全标准的时间关系式，从而实现对实际情况的预测仿真，比较有力地解释了大李所遇到的情况,并为制定科学的检测标准提供依据。 [关键词] 饮酒驾车；数学模型；MATLAB; 药物动力学；房室模型; 曲线拟合 问题的重述 根据新颁布的《车辆驾驶人员血液、呼气酒精含量阈值与检验》的标准，解决2004年“高教社杯”全国数学建模竞赛中的C题[1] 2． 模型的假设 1．假设肠胃中酒精吸收率正比于肠胃中酒精含量； 2．忽略从饮酒酒精开始吸收的延迟时间； 3．忽略事物摄入对酒精吸收速率造成的影响； 4．由于呼吸和排泄对体液中酒精含量变化影响很小，因而忽略呼吸和排泄对酒精含 量的影响； 5．假设酒精在血液和体液中均匀分布； 6．忽略酒精在人体内含量对人体消化能力的影响 3． 问题的分析 通过比较，我们认为大李超标的主要是由两次饮酒前血液中酒精初始浓度不同，喝酒次数和喝酒时间长短不同造成的。当酒是在短时间内喝下的，，用房室模型中的一室模型建立线性微分方程处理；当酒是分多次等量等时间间隔在较长时间内喝下时，就形成一个多剂量给药的问题，可以用迭加的方法解决。酒精在血液中的含量必定存在一个峰值，我们使用药物动力学模型对峰值的大小和达到峰值的时间进行了计算。最后讨论的是每天定时等量喝一次酒的情况，主要讨论了饮酒的时间和饮酒量的变化对血液中酒精浓度的影响。 5．基本变量符号和模型的建立 酒精可以经自由扩散或主动运输作用由胃肠道进入血液，然后分布于人体体液，并逐渐被消耗。结合房室模型，酒精首先迅速进入吸收室（胃肠道），在该室中逐渐吸收，然后同时将吸收到的酒精逐渐转移到中心室（体液、血液），并在中心室中被分解消耗.可以用框图直观的表示酒精在体内吸收和消耗的过程： 吸收室 中心室 D K1 k c(t) 胃肠道酒精量 血液酒精含量 记吸收室中酒精浓度为c1(t), 中心室(体液) 中酒精浓度为c(t)。记吸收室的体积为V1，中心室体积为V。根据房室模型，药物（酒精）从吸收室向中心室传输速率以及中心室分解酒精的消耗速率都与该室中酒精含量成正比，比例系数分别为转移系数k1和排除系数k。 根据这些假设即可求解模型。 5．1 短时间喝酒血醇浓度的确定 设t=0时刻，一次喝的酒中乙醇含量为D，则c1(t)由方程确定 …………………………………………（1） 中心室的酒精浓度c（t）的变化率由两部分组成：与c(t)成正比的排出，与c1(t)成正比的吸收，比例系数分别为k与k1 。则可得c(t)满足方程： ………………………………………（2） 由（1）可解出………………………………………………… （3） 将（3）代入（2）中，可以解得， …… （4） 从上可以看出， 结合题目所给数据，用MATLAB中的类似lsqcurvefit命令拟合，可以得到： =114.4315 k=0.1855 k1=2.0080 c0=0 也就是说，短时间喝两瓶啤酒后，体液中酒精含量变化为： …………………………………（5） 试验数据和拟合曲线如下所示 ： 根据（4）（5）两式，可得短时间内饮2瓶、3瓶啤酒，血醇随时间变化关系： 相应的血醇变化曲线分别为： 5．2 长时间慢速饮酒血醇浓度的确定 假设较长时间T0内，喝下含酒精量D0的啤酒，简单起见，假设饮酒速度是恒定的，则在T时刻喝进就进酒精的微元为： 根据（4）式，在时刻T、饮进酒精微元在喝完酒t小时后，血液中，酒精浓度微元为： 于是在喝完酒t小时之后，血醇含量为： 要预测2小时内喝下3瓶啤酒时，血液中酒精含量达标所需等待时间，只要将 k=0.1855 k1=2.0080 T0=2 代入，即可利用MATLAB绘出曲线图，如下：