幂指分段函数在饮酒数学模型中的叠加应用

1、幂指分段函数在饮酒数学模型中的叠加应用吴磊重庆大学B区城市建设与环境工程学院市政工程专业，重庆（400045）E-mail：摘 要：本文分别采用分段函数的数学思想建立了饮酒后血液中酒精含量y(t)随时间t的变化关系的动态数学模型，并对问题有了很好的解答。最后对驾驶人员和想喝一点酒的人提出了一些建议和忠告，希望交通事故不再发生，通过模型验证，得出模型具有很好的应用价值。 关键词：幂指分段函数；数学模型；叠加效果；延迟效应中图分类号：o 01 0291. 问题的提出据报载，2003年全国道路交通事故死亡人数为10.4372万，其中因饮酒驾车造成的占有相当的比例。针对这种严重的道路交通情况，国家质量

2、监督检验检疫局2004年5月31日发布了新的车辆驾驶人员血液、呼气酒精含量阈值与检验国家标准，新标准规定，车辆驾驶人员血液中的酒精含量大于或等于20毫克百毫升，小于80毫克百毫升为饮酒驾车（原标准是小于100毫克百毫升），血液中的酒精含量大于或等于80毫克百毫升为醉酒驾车（原标准是大于或等于100毫克百毫升）。大李在中午12点喝了一瓶啤酒，下午6点检查时符合新的驾车标准，紧接着他在吃晚饭时又喝了一瓶啤酒，为了保险起见他呆到凌晨2点才驾车回家，又一次遭遇检查时却被定为饮酒驾车，这让他既懊恼又困惑，为什么喝同样多的酒，两次检查结果会不一样呢？请你参考下面相关数据（或自己收集资料）建立饮酒后血液中酒

3、精含量的数学模型1，并讨论以下问题：1. 对大李碰到的情况做出解释；2. 在喝了3瓶啤酒或者半斤低度白酒后多长时间内驾车就会违反上述标准，在以下情况下回答：1）酒是在很短时间内喝的；2）酒是在较长一段时间（比如2小时）内喝的。3. 怎样估计血液中的酒精含量在什么时间最高。4. 根据你的模型论证：如果天天喝酒，是否还能开车？5. 根据做的模型并结合新的国家标准写一篇短文，给想喝一点酒的司机如何驾车提出忠告。2. 问题的分析2.1 参考数据1. 人的体液占人的体重的65%至70%，其中血液只占体重的7%左右；而药物（包括酒精）在血液中的含量与在体液中的含量大体是一样的。2. 体重约70kg的某人在

4、短时间内喝下2瓶啤酒后，隔一定时间测量他的血液中酒精含量（毫克百毫升），得到数据如下：表1 测定的喝酒后某人血液中的酒精含量表Tab1 after drinking Determination of alcohol Content of a person's blood in Table 时间(小时)酒精含量时间(小时)2.2 参考数据及问题的分析针对问题的提出，假设y(t)为t时刻血液中的酒精含量，由参考数据分析，我们设了一个幂指分段函数作为我们的目标函数作为我们的目标函数，这个分段函数由两部分组成，第一段为抛物线型幂函数，第二段为衰减的指数型，采取了数据拟合的思想方法，由参考数据，

5、我们假设了用y=at+bt(0t1.5)作为抛物线型而用y=ce2(t1.5)(t>1.5)(<0)作为衰减的指数型，我们在模型的推广中应用第二个衰减指数型，为了求参数，我们针对不同的情况采用了不同的实测计算数据，针对体重约为70kg的某人在短时间内喝下1瓶或3瓶啤酒的模型参数，求出三种不同情况下的血液中的酒精含量与时间t的函数关系式，然后根据血液中的酒精含量与人的体重成反比，得出了不同体重情况下酒精含量与时间t的关系式，即把70/m作为函数关系式的系数而得到。然后，根据我们建立的短时间内饮1瓶酒的模型，对大李的情况进行解释，并在短时间内喝完酒的模型下对问题2与问题3进行了求解。图

6、1 喝酒后人体血液中酒精含量随时间变化的拟合曲线图Fig1 After drinking human blood-alcohol content changing with time fitting curve对于估计血液中的酒精含量在什么时间最高的问题我们应用了微分法求极值的数学思想9。随后，我们对模型进行了推广，把3瓶酒分成n等份，相当于在任一小区间2/n小时内喝完3/n瓶，即每一个微小的时刻都在喝一个微小量，从而酒精摄入就成为连续的了，能更好的表示酒在较长时间内摄入的情况，但是由生活常识及叠加时计算的复杂性，我们只讨论等分为3次和10次的情况，通过计算机数学软件Mathematica4和

7、Matlab2编程绘制出了其图形，并由相应的结论，对所提出的问题作出了更进一步的解释。3.模型的建立及求解3.1模型一 短时间内饮酒2at+bt(0t1.5) （拟合过程间检测他的血液中的酒精含量的分段函数关系式为y=(t1.5)(t>1.5)ce在Matlab中运行）对于指数衰减函数y=ce(t1.5)(t>1.5)(<0)可以转化为线性lny=lnc+(t1.5)， - 2 -令z=lny,lnc=b作线性最小二乘拟合，经Matlab运行结果求出参数值，故函数关系式为y=22.513e0.1798(t1.5)265.0509t+151.2597t(0t1.5) 所以总的函

8、数关系式为：y=0.1798(t1.5)(t>1.5)22.513e2at+bt(0t1.5)血液中酒精含量与时间t的分段函数关系式为：y= (t1.5)(t>1.5)ce仍然把指数衰减函数转化为线性关系lny=lnc+(t1.5)，z=lny,lnc=b，经Matlab运行，经计算：c=13.2447230.0618t+76.2581t(0t1.5)总的分段函数关系式为y= 0.1942(t1.5)(t>1.5)13.2447e290.1854t+219.77437t(0t1.5) 所以总的函数关系模型为y=0.1924(t1.5)(t>1.5)38.9586e总之，

9、通过对某人短时间内喝完两瓶，一瓶和三瓶啤酒的数据分析，通过线性最小二乘法的拟合思想，得出三种情况下的含量与时间t的函数关系。3.2 模型二 长时间内饮酒t0,0.67,m1(t)t0.67,1.33,m1(t)+m2(t)时间分配为：t1.33,1.5,m1(t)+m2(t)+m3(t)t1.5,2.17,n1(t)+m2(t)+m3(t)t2.17,2.81,n1(t)+n2(t)+n3(t)又由一瓶啤酒在很短的时间内喝完的拟合函数为：230.0618t+76.2581t(0t1.5)y= 0.1942(t1.5)(t>1.5)13.2447e2m1(t)=30.0618t+73.25

10、81t(0t1.5)其中 (0.1942(t1.5)n1(t)=13.2447e2m2(t)=30.0618(t0.67)+73.2581(t0.67)(0t1.5) (0.1942(t0.67)1.5)n2(t)=13.2447e2m3(t)=30.0618(t1.33)+73.2581(t1.33)(0t1.5) (0.1942(t1.33)1.5)n3(t)=13.2447e- 3 -在Mathematica软件中编制程序10进行叠加，得到叠加后的图形。230.0618t+76.2581t(0t1.5)，那么十分之三瓶的啤酒在短时函数变化关系：y=0.1942(t1.5)(t>1.

11、5)13.2447e间内喝完的人体血液中的酒精含量与时间t的函数关系为：32(0t1.5)30.0618t+76.2581t）10，这样模型的分析与计算就同三瓶啤酒在两y=313.2447e0.1942(t1.5)）(t>1.5)10个小时内分三次等量喝完的情况，针对这种情况，即把3瓶啤酒分10次喝完，并且每次喝完是在短时间内进行，经过在Mathematica软件中编程并运行，得出了人体中血液中酒精含量随时间的变化关系。经对比发现，叠加后的图形高峰点上升，但高峰点到来的时间推后。4.问题的解答4.1 问题一的分析与解答由于饮酒后酒精在人体血液内的积存，对人体血液中酒精浓度的峰值有个延迟效

12、应，所以虽然大李两次喝同样的酒，但是由于第一次喝的酒精仍在起作用，所以累积的结果是他第二次喝后超标，假设他是在晚上七点半时喝酒，经计算，第一次检查时，血液中的酒精含量为18.51（毫克/百毫升），小于20（毫克/百毫升），即第一次喝后符合标准，而第二次，还要考虑第一次喝后的影响，两次产生的共同影响为20.71（毫克/百毫升），大于20（毫克/百毫升），因此第二次时，被定为饮酒驾车。因此，他喝同样的酒，却产生不同的检查结果。图2 两次喝酒后叠加后的拟合曲线Fig2 After two drinkings superposed fitting curve- 4 -4.2 问题二的分析与解答（1）对

13、三瓶啤酒或半斤低度白酒是在很短时间内喝完的情况的分析与解答：由于我们分别建立了短时间内喝完一瓶、两瓶和三瓶的模型一，根据三瓶啤酒在短时间内喝完的模型及图像：290.185t+219.7743t(0t1.5)y= (0.1924(t1.5)(t>1.5)38.9586e因为从函数解析式与附图知，当t=1.5时，yt=1.5=130.83（毫克/百毫升）所以抛物线段所对应的时间区域内不能驾车，多长时间内驾车就会违反上述标准，只可能在指数衰减阶段求解，由38.9586e(0.1924(t1.5)<20可得t1>11.26小时，所以短时间内喝完3瓶，只能在t1>11.26小时以

14、后才能驾车。（2）对3瓶啤酒或半斤低度白酒是在较长时间（比如2小时）内喝的情况下的分析与解答。根据由3瓶啤酒在较长时间内（2小时内分3次）内喝的血液中酒精含量与时间t的函数关系式和3瓶啤酒在较长时间内（2小时内分10次）的血液中酒精含量与时间t的函数关系式推广到2小时内分n次喝完血液中酒精含量与时间t图像。然后由图算法得出结论。只考虑指数衰减阶段，当血液中酒精浓度不大于20毫克/百毫升，t2=12.4小时即司机在长时间内喝了3瓶酒，在酒后12.4小时内驾车就会违反上述标准。综合模型（1）与（2）我们得到以下结论：对于一定量的啤酒，通过对在很短时间内喝完t1<t2这是符合实际情况的，也就是

15、在饮酒较长时间的情况下，和在较长时间内喝完的对比，血液中酒精含量，不仅考虑其纵标的叠加，还要考虑其累积后的延迟效应，即达到新的驾车标准的规定时，较长时间饮酒后，取血液中酒精浓度与标准值浓度对应的时刻t都要延迟。4.3 问题三的分析与解答针对短时间内饮酒的情况，我们只考虑饮酒量为一瓶时，对其它饮酒量我们可以采取乘以相应的系数，短时间饮酒一瓶的函数解析式为：230.0618t+73.2581t(0t1.5)y= (0.1924(t1.5)(t>1.5)13.2447e可以用抛物线求极值方法也可以用数形结合的思想：法一：y=30.0618t2+73.2581t当t=1.218458，ymax=

16、43.63法二：数形结合地思想，从图中找出血液中酒精浓度达到最大值时对应的时间t.针对长时间内（不妨设2小时）内喝完啤酒时（3瓶啤酒），我们采用了图算法求峰值地方法，即抛物线顶点处，在两小时内将3瓶啤酒分3次等量喝完，用图像法求血液中酒精含量什么时间最高：（当y取最大值时，t=2.15小时）结论：与长时间内将3瓶啤酒喝完血液中酒精含量峰值比较，此时相对延迟。4.4 问题（4）的分析与解答根据假设，考虑了每天定时，定量喝一次酒，而且第一天这个时间喝的酒到第二天其影响可以忽略掉，还有每天定时定量喝的酒都是在短时间内喝完，每天定量喝x瓶，结合实际情况 - 5 -进一步考虑：对于喜爱喝酒的司机来说，最

17、理想的情况就是喝了酒后，随时都可以驾车，那么，只需饮酒后血液酒精浓度达到最大值时含量仍符合新的驾车标准即可，为此，列出了抛物线型函数。饮酒量1瓶：y=30.0618t+73.258t(0t1.5)饮酒量x瓶：y=x(30.0618t+73.258t)，既当t1.22时，取最大值，所以yt=1.22=x(44.63+89.26194056)<20所以x<0.4681瓶即x<0.47瓶所以根据假设条件，若司机每天定时定量喝一次酒，二十四个小时后酒的影响可以忽略，每天喝的酒都是在短时间内喝完，那么当饮酒量x<0.47瓶时，司机天天喝酒，且还能开车，从这一点上看，本模型很具有实

18、际意义，推广价值。 224.5 问题（5）的分析与解答酒精在血液中的含量随时间的变化是先升高而后减少，从而可以推算出当一名司机在饮酒后，经过多长时间后驾车不影响其正常驾驶，也就是当司机饮酒后血液中的酒精浓度经过t时间后，要小于20毫克/百毫升。对于司机饮酒驾车，针对饮酒时间，即是短时间内喝完还是在较长一段时间内喝完，建立了针对上述两种情况的幂指分段函数的数学模型。并对这模型进行分析，求解和作图。通过这个模型让广大司机朋友能对自己的饮酒情况有所了解，并以此判断自己是否能够驾车。以题中的大李为例，当他在中午12点喝酒后其身体内血液中酒精含量升高，经过6小时，血液中的酒精大部分已排出体外，即血液中酒

19、精含量低于20毫克/百毫升，符合新的驾车标准，但其体内血液中仍存在部分酒精，当其晚上再次喝酒时，血液中的酒精浓度又一次升高，使两次进入血液中的酒精相叠加，再经过8小时，血液中酒精浓度便高于20毫克/百毫升并超过国家标准。第二种情况，如果该司机想天天喝酒又能正常开车，首先是对于建立的短时间内饮酒的模型，采用了分段函数的概念，分为开口向下的递增抛物线和指数衰减曲线两部分，这就是说，司机在短时间内，饮酒以后，血液中的酒精含量是先增后减，所以对于想喝一点酒的司机来说，按照问题（4）的解答，以满足一次饮酒后血液中酒精浓度最大值符合新的国家标准，即天天喝酒，也能照样开车。对于长时间内饮酒模型首先也是假设该

20、司机每次喝酒的量都是一定的，则经计算可得出，只要该司机喝啤酒的总数不超过0.47瓶，他便可以天天喝酒，也不会影响开车。其次，对于建立的喝酒是在较长时间内完成的模型，根据图形分析，其血液中的酒精含量的最大值将会延迟，且血液中酒精含量将产生叠加并累积的效应，而使最大值发生变化，使最大值增大，最后，其叠加后的图形的最大值也将增大，这样以来，对于想喝一点酒的司机，如果天天喝酒，按问题（4）的解答，只有酌情减少每次饮酒量，这样才能保证血液中酒精含量符合新的国家标准。再次，对于想喝一点酒的司机，如果想喝了酒立即开车，那就要保证血液中酒精含量的最大值符合新的驾车标准，如果喝了酒以后不立即开车，而是休息或等待

21、一段时间后再驾车，那么这时就要保证这时血液的酒精含量符合新的驾车标准；总之，当一次饮酒量小于一定数值时，如问题（4）的解答，那么司机可以立即开车；当一次饮酒量大于这个定数时，司机不能立即驾车，只能在等待一段合适的时间后，才能驾车。利用相关的数学与人体生理知识，向各位驾驶人员提出以下通俗忠告和建议：- 6 - 如果是您想喝一瓶啤酒，并且是短时间喝, 6小时后方驾车上路； 如果是您想喝2瓶啤酒,并且是短时间喝, 10小时后方驾车上路； 如果是喝3瓶啤酒或半斤低度白酒，当天都不要开车； 如果是每餐喝一瓶啤酒，在14小时内都不要开车； 如果是一段时间内都在喝酒的话，当天也不要开车了。5. 模型的评价与

22、推广5.1模型的评价针对饮酒驾车的实际问题，建立了这个综合模型，模型里又有针对不同条件的小模型，我们不仅建立了体重约70kg的某人短时间内喝一瓶啤酒的模型，也建立了喝两瓶啤酒，三瓶啤酒的模型，并建立了此人在较长时间内喝完3瓶啤酒的模型，应用Matlab中线性最小二乘法拟合图像，得出了函数的参数并求出了函数；还运用了Mathematica软件3求解图像和Matlab2中优化工具箱中leastsp计算，拟合图像。用幂指分段函数的思想建模，能够确切的体现数学建模的思想，即把实际问题转化为数学语言，再把数学语言用数学式子表达出来，结合实际问题，用数学的方法解决了实际问题。所建立的分段函数叠加模型简洁、

23、明了，便于使用数学工具，降低了编程求解的难度，缩短了运行时间，提高了工作效率。这两个模型对于喝酒外出的司机很有指导意义，特别是想天天喝一点酒的司机，更具有借鉴意义，其实用价值很大。整个模型的结果最终由普通的数学关系表达式和几个图形表示，具有很好的实用性。在基于准确性的同时，对该模型的参数作了较为细腻的估计，在其表达式中有意识的对参数取了较多的有效数位，使其准确度更高，数据拟合很好，结果可信度更高。在酒量相同的条件下，较短时间内喝酒的人其自身酒精浓度达到最高峰的时间比较长时间内喝酒的要稍低一些，这与现实生活中出现的情况很吻合。5.2 模型的推广所列的模型是针对体重为70kg的人，那么对于质量为m

24、kg的人，模型推广为各项乘以系数70/m，这样就更具有普遍指导意义。题中只涉及到天天喝一次，且每次定时定量，如果天天喝n次，必须保证叠加后的峰值小于20，从而算出每次喝的酒量用以指导司机或饮酒人员。因为题所设为酒精含量，而在实际中通常是多少度的酒，二者存在一个转化关系，转化之后才更贴近于实际。模型是根据酒精在血液里的变化消耗减小规律得出了酒精在血液中的浓度变化类似药物，进入人体内起药里作用的浓度变化，所以该模型可进一步的推广到药物在人体中的浓度随时间的变化，对医学和临床的治疗，给药等提供了很好的参考价值。- 7 -参考文献1姜启源 数学模型M 北京：高等教育出版社 1993.2张志勇 精通MatlabM 北京：北京航空航天大学出版社 2003.3萧树铁 数学实验M 北京：高等教育出版社 1999.4嘉木工作室 Mathematica应用实例教程M 北京：机械工业出版社.5杨启帆 方道元，数学建模M 杭州：浙江大学出版社，1999.6姜启源 谢金星，数学建模M 北京：高等教育出版社，2003.7李 涛 贺勇军，应用数学篇M 北京：电子工业出版社，2000.8唐焕文、贺明峰，数学模型引论M 北京：高等教育出版社，2001年9月第1版.9徐全智、杨晋浩，数学模型M 北京：高等教育出版社，2003年7月第1版