第1页共7页吸收室x1(t)中心室x2(t)k1k2编号:时间:2021年x月x日书山有路勤为径,学海无涯苦作舟页码:第1页共7页问题的提出:国家质量监督检验检疫局2004年5月31日发布了新的《车辆驾驶人员血液、呼气酒精含量阈值与检验》国家标准,新标准规定,车辆驾驶人员血液中的酒精含量大于或等于20毫克/百毫升,小于80毫克/百毫升为饮酒驾车(原标准是小于100毫克/百毫升),血液中的酒精含量大于或等于80毫克/百毫升为醉酒驾车(原标准是大于或等于100毫克/百毫升)。某人在中午12点喝了一瓶啤酒,下午6点检查时符合新的驾车标准,紧接着他在吃晚饭时又喝了一瓶啤酒,为了保险起见他呆到凌晨2点才驾车回家,又一次遭遇检查时却被认定为饮酒驾车,这让他既懊恼又困惑,为什么喝同样多的酒,两次检查结果会不一样呢?请你参考下面给出的数据建立饮酒后血液中酒精含量的数学模型,并讨论一下问题:(1)对此人碰到的情况做出解释。(2)假设酒是在很短的时间内喝的,在喝了3瓶啤酒或者半斤低度白酒后多长时间内驾车就会违反上述标准?(3)怎样估计血液中的酒精含量在什么时间最高。(4)根据你的模型论证:如果天天喝酒,是否还能开车?以下是某人喝了两瓶啤酒后血液酒精浓度(毫克/百毫升):t(h)0.250.50.7511.522.533.544.55浓度306875828477706858515041t(h)678910111213141516浓度3835282518151210774模型:常微分方程模型主要命令:polyfit,fzero问题分析:本题类似于药物在体内的分布与排除,故可使用二室模型。模型假设及符号说明:(1)机体分为吸收室和中心室,两室容积在整个过程中保持不变,忽略酒精在其他器官内的浓度。(2)酒精从一室向另一室转移速率及向体外排出的速率与该室的酒精浓度成正比。(3)只有中心室与体外有酒精转移。(4)酒在很短的时间内就全部进入胃内,即进入吸收室内。根据以上假设,可以给出以下的酒精转移流程图x1(t):吸收室(主要为胃部)内酒精含量第2页共7页第1页共7页编号:时间:2021年x月x日书山有路勤为径,学海无涯苦作舟页码:第2页共7页x2(t):体内的酒精含量V1:吸收室的体积V2:中心室的体积c1(t):吸收室的酒精浓度c2(t):中心室的酒精浓度k1:酒精由吸收室向中心室的转移速率k2:酒精由中心室向体外的转移速率d:两瓶啤酒中的酒精含量模型建立:又,则模型求解:由上述方程可解得:其中,A1、A2、k1、k2待定。系数的确定:表中所给数据为血液中酒精浓度,即c2(t)。酒精吸收速率大于酒精排出体内的速度,故当t充分大时,即取时的数据:第3页共7页第2页共7页编号:时间:2021年x月x日书山有路勤为径,学海无涯苦作舟页码:第3页共7页t22.533.544.5567c2777068585150413835lnc24.34384.24854.21954.06043.93183.91203.71363.63763.5553t8910111213141516c2282518151210774lnc23.33223.21892.89042.70812.48492.30261.94591.94591.3863做一阶线性拟合:则,,当t较小时,取时的数据,并加上初始值的数据点(0,0)t00.250.50.7511.5c203068758284b1.00000.70440.34510.24390.14480.0516lnb0-0.3504-1.0641-1.4110-1.9323-2.9643做一阶线性拟合:忽略在截距上的误差0.0426,综上可知:用Matlab做c2-t曲线如下所示:第4页共7页第3页共7页编号:时间:2021年x月x日书山有路勤为径,学海无涯苦作舟页码:第4页共7页02468101214160102030405060708090t/hc/(mg/100ml)酒精浓度随时间的变化原始数据点c2(t)拟合曲线又主要是对人体血液中的酒精浓度进行考察,故不去求解c1(t)中参数的具体数值。c2(t)中参数讨论:(1)对此人的情况的分析设中午12点时刻t=0,根据此人的情况分析可得:喝第一瓶酒时:喝第二瓶酒时:综上可得第5页共7页第4页共7页编号:时间:2021年x月x日书山有路勤为径,学海无涯苦作舟页码:第5页共7页下午6点,即t=6h时,,因而在下午6点检查时符合新的驾车标准。第二次喝酒时,,因而在凌晨两点检查时不符合驾车标准。(2)根据之前的分析,喝了3瓶啤酒或半斤低度白酒后,体内酒精浓度的变化规律为令,解得。所以在喝完酒约11.2368个小时之内驾车就会违反上述规则。(3),其中A取决于喝酒的多少。令,解得。由上式可见,血...