饮酒

饮酒驾车一、摘要本文建立了一个为确定司机酒后驾车时酒精含量是否超标提供参考的数学模型。首先，将酒精含量视为血药浓度，借用药物动力学的房室模型，将酒精在肠胃的吸收过程和在血液中的分解过程抽象为吸收室和中心室里所发生的作用，运用微分方程理论推导出了吸收速率和分解速率随时间变化的规律（,），并用回归分析方法结合题述经验数据具体导出一人在未喝过酒的情况下，饮入2瓶啤酒的血液中酒精含量与时间的关系模型()，以此为基础回答题述的五个问题。本文的最大特点是将房室模型灵活运用于司机各种不同的饮酒方式，使模型的应用范围和解释力都得到了加强。二、问题的重述据报载，2003年全国道路交通事故死亡人数为10.4372万，其中因饮酒驾车造成的占有相当的比例。针对这种严重的道路交通情况，国家质量监督检验检疫局2004年5月31日发布了新的《车辆驾驶人员血液、呼气酒精含量阈值与检验》国家标准，新标准规定，车辆驾驶人员血液中的酒精含量大于或等于20毫克／百毫升，小于80毫克／百毫升为饮酒驾车（原标准是小于100毫克／百毫升），血液中的酒精含量大于或等于80毫克／百毫升为醉酒驾车（原标准是大于或等于100毫克／百毫升）。大李在中午12点喝了一瓶啤酒，下午6点检查时符合新的驾车标准，紧接着他在吃晚饭时又喝了一瓶啤酒，为了保险起见他呆到凌晨2点才驾车回家，又一次遭遇检查时却被定为饮酒驾车，这让他既懊恼又困惑，为什么喝同样多的酒，两次检查结果会不一样呢？请你参考下面给出的数据（或自己收集资料）建立饮酒后血液中酒精含量的数学模型，并讨论以下问题：1.对大李碰到的情况做出解释；2.在喝了3瓶啤酒或者半斤低度白酒后多长时间内驾车就会违反上述标准，在以下情况下回答：酒是在很短时间内喝的；酒是在较长一段时间（比如2小时）内喝的。3.怎样估计血液中的酒精含量在什么时间最高。4.根据你的模型论证：如果天天喝酒，是否还能开车？5.根据你做的模型并结合新的国家标准写一篇短文，给想喝一点酒的司机如何驾车提出忠告。参考数据1.人的体液占人的体重的65%至70%，其中血液只占体重的7%左右；而药物（包括酒精）在血液中的含量与在体液中的含量大体是一样的。2.体重约70kg的某人在短时间内喝下2瓶啤酒后，隔一定时间测量他的血液中酒精含量（毫克／百毫升），得到数据如下：时间(小时)0.250.50.7511.522.533.544.55酒精含量306875828277686858515041时间(小时)678910111213141516酒精含量3835282518151210774三、问题分析.该部分包含于五、模型的建立与分析中.四、模型假设与符号约定假设人体密度是均匀的，而且酒精进入体液的时候马上均匀分布，且血液和体液的酒精浓度是一样的，都以c2表示。为了方便计算，保守计啤酒的酒精浓度为5%（g/ml）。机体分为吸收室和中心室，两个室的的容积(v1,v2),即肠胃容积和体液（血液是其中的一部分）的体积在过程中保持不变。两室的酒精总量分别为x1,x2。两室的酒精含量分别以浓度c1，c2表示。酒精从一室向另一室的转移速率（速率系数k1），及中心室分解酒精的速率（速率系数k2）,与该室的酒精浓度成正比。吸收室从外界补给酒精，并将其吸收和散布在体液中；中心室把体液中的酒精分解并将其产物排出中心室外的环境。两室某一时刻的酒精吸收速率或酒精的分解速率分别为,。g0为酒精的补给速率，是啤酒补给量对时间的导数。7、喝入2瓶啤酒的酒精总量为D0。在这些假设下的一种二室模型示意图如图所示。中心室（2室）x2(t),c2(t),v2吸收室（1室）x1(t),c1(t),v1K2K1G0五、模型的建立与分析5．1．下面建立一人在未喝过酒的情况下，短时间内饮入2瓶啤酒的酒精含量与时间的关系模型。由假设条件与上图，可以写出两个房室中酒精量所满足的微分方程。的变化率由1室向2室的转移,1室的酒精补给速率组成；的变化率由由1室向2室的转移及2室向体液外的转移组成。于是有-------------------------------------------（1）--------------------------------------------（2）-------------------------------------------------（3）由（1）、（2）、（3）式，得---------------------------（4）-----------------------------（5）（4）、（5）两式构成线性常系数非齐次微分方程，一人在未喝过酒的情况下，短时间内饮入2瓶啤酒，g0=0，即一下子喝完啤酒，（4）式可化为=------------------------------------（6）同时，有c1(0)=，c2(0)=0，由（5）、（6）联合，通过MATLAB可解得c1(t)=---------------------------------(7)c2(t)=---------------(8)令A=，则(8)可化为c2(t)=A—A-----------------------------（9）不妨设，第一、当t很大时，（9）可看成是c2(t)=A-------------------------------(10)两边取对数,使得(10)变为线性,得lnc2(t)=lnA-k2t-------------------------------(11)通过运用MATLAB进行线性回归分析,得lnA=4.9078A=135.34k2=0.20507第二、由（9）得A=A—c2(t)，在第一步的基础上代入k2，A，t可得出A—c2(t)的一系列具体的值，所以ln（A）=lnA-k1t=ln（A—c2(t)）就可以通过MATLAB做线性回归分析得出k1=1.4976，lnA=4.9425，A=140.1201由第一步和第二步的出的lnA有点误差，但从它们的线性残差分析可见，这两个结果都可用，把上面得出的结果代入（9）式，得c2(t)=---------------（12）这是一人在未喝过酒的情况下，短时间内饮入2瓶啤酒的酒精含量与时间的关系模型，（12）式作图如下：下面是题目给出的数据作出的散点图，对比一下上图，可见上图拟合的结果是比较精确的。5.2下面建立一人在未喝过酒的情况下，短时间内饮入n瓶啤酒的酒精含量与时间的关系模型.---------------------------（6）-----------------------------（5）c1(0)=，c2(0)=0c1(t)=-------------------------------------(13)c2(t)=---------------------------(14)其中,A=135.34,k2=0.20507,k1=1.49765.3.当一人喝酒前肠胃（吸收室）里没有酒精或酒精已被吸收完了，而体液里（中心室）的酒精仍有残余，这时他立即在短时间内喝完n瓶酒，吸收室酒精总量为nD0/2,此时体液酒精浓度（设为初始浓度）c2(0)=a,吸收室酒精初始浓度为c1(0)=nD0/2/v1，结合以下两式---------------------（5）=-----------------------（6）c2(0)=a,c1(0)=nD0/2/v1用MATLAB解之并化简得-------------------------------(13)---------------(15)其中,A=135.34,k2=0.20507,k1=1.49765.4当一人吸收室酒精含量为0,中心室亦为0的时候,先用2小时持续地以均匀速率给酒n瓶,在2小时后给酒停住.在这种情况下,c1(0)=0,c2(0)=0。率为g0=nD0/4.此时的微分方程组的形式为:-----------------------------（5）---------------------------（16）c1(0)=0,c2(0)=0用MATLAB解之并化简得到:-----------------------------------(17)---------------(18)其中,A=135.34,k2=0.20507,k1=1.4976六、问题的解答6．1.对大李碰到的情况做出解释大李在中午12点喝了一瓶啤酒，下午6点检查时符合新的驾车标准。因此符合5.1.的模型的条件。运用5.2.的模型，代入t=6，n=1c2(6)=19.76320,符合大李在中午12点喝了一瓶啤酒，下午6点检查时符合新的驾车标准。紧接着他在吃晚饭时又喝了一瓶啤酒，为了保险起见他呆到凌晨2点才驾车回家，又一次遭遇检查时却被定为饮酒驾车。根据5.1.的模型c1(t)=-------------------------------------(7)代入k1=1.4976，=25（根据经验，可假设大李开始喝酒那一刻的吸收室酒精含量为普通啤酒酒精浓度的一半）得c1(6)=0.003130忽略不记。因此他的情况是下午6点吃饭时吸收室酒精浓度为0。按正常情况，大李可能属于5.3模型，但计算5.3模型的c2（8）=16.95020，不符合超标的实际情况。所以，大李的情况必然是晚饭时酒非一饮而尽，且其体液（即中心室）中有酒精残余。但这种情况不包含在第五点所建立的四种情况的模型中。设他喝2小时酒，我们可以忽略着两个小时体液内的酒精消耗，设他从第8个小时开始一下子喝完，用5.3.得c2（6）=23.5955〉20所以超标6．2在喝了3瓶啤酒或者半斤低度白酒后多长时间内驾车就会违反上述标准，在以下情况下回答：（1）酒是在很短时间内喝的（2）酒是在较长一段时间（比如2小时）内喝的。对（1），运用5.2.的模型，代入n=3，c2(t)=---------------------------(14)其中,A=135.34,k2=0.20507,k1=1.4976解不等式得同理，解不等式得因此，若酒是在很短时间内喝的，以开始喝酒那一时刻计时间，内驾车就会违反上述“80mg”标准，内驾车就会违反上述“20mg”标准。对（2），运用5.4的模型，代入n=3，----------------(18)其中,A=135.34,k2=0.20507,k1=1.4976解不等式得同理，解不等式得因此，若酒是在较长一段时间（比如2小时）内喝的，以开始喝酒那一时刻计时间，内驾车就会违反上述“80mg”标准，内驾车就会违反上述“20mg”标准。6.3估计血液中的酒精含量在什么时间最高，需要对各种具体的饮酒方式的血液中酒精含量与时间的关系模型()求最植，下面列举一些简单的结果。对符合5。1模型的情况，时候酒精含量最高；对符合5。2模型的情况，时候酒精含量最高；对符合5。3模型的情况，何时酒精含量最高与喝酒前血液（中心室）中酒精残余含量有关。对符合5。4模型的情况，时候酒精含量最高。6.4随着生活的快节奏，我们假设酒都是一下子喝完的，利用5.2c2(t)=20代入t=12得出n=4即每天不得超过4瓶,否则不能开车。可见,只要掌握好了饮酒的方式和驾车的时机,天天喝酒还是能开车的.6.5短文:掌握好了饮酒的方式和驾车的时机-------给想喝一点酒的司机如何驾车提出忠告司机朋友们,掌握好了饮酒的方式和驾车的时机,既有利于公民对国家法规的遵守,又有利于公民人身安全的有效保障.在此,我们向各位司机朋友门提出几点忠告.若你想喝2瓶啤酒,并且是短时间喝,我们劝您10小时后方驾车上路.若你想喝1瓶啤酒,并且是短时间喝,我们劝您6小时后方驾车上路.为了您和他人的安全与幸福,请您掌握好了科学的饮酒的方式和适当的驾车的时机.七.参考文献[1]梁炼数学建模广州:华南理工大学出版社,2003.[2]姜启源谢金星叶俊数学模型(第三版)北京:高等教育出版社,2003[3]韩伯棠管理运筹学北京:高等教育出版社,2003[4]WILLAMF.LUCAS主编微分