数学建模暑期训练题2.doc

饮酒驾车的数学模型
摘要
本模型目的是解决司机安全驾车与饮酒的问题，通过建立一个数学模型分析司机如何适量饮酒不会影响正常的安全驾驶。根据一定合理的假设，建立人体内酒精浓度随时间变化的微分方程模型，并通过拟合曲线对数据进行分析。在不同饮酒方式下进行分类讨论，得出体内酒精浓度随时间的变化函数。在讨论过程中，我们得到两个结论：在短时间喝酒形式下，达到最大值的时间为1.23小时，与喝酒量无关；在长时间喝酒形式下，喝酒结束时酒精含量最高。
关键词：微分方程、模型、房室系统。
一、问题重述
饮酒驾车问题主要是分析驾驶员在喝过一定量的酒后，酒精在体内被吸收后，血液中酒精含量上升，影响司机驾车，所以司机饮酒后需经过一段时间后才能安全驾车，国家标准新规定，车辆驾驶人员血液中的酒精含量大于或等于20毫克/百毫升，小于80毫克/百毫升为饮酒驾车，血液中酒精含量大于或等于80毫克/百毫升为醉酒驾车，司机大李在中午12点喝下一瓶啤酒，6小时后检查符合新标准，晚饭地其又喝了一瓶啤酒，他到凌晨2点驾车，被检查时定为饮酒驾车，为什么喝相同量的酒，两次结果不一样？讨论问题：
对大李碰到的情况做出合理解释；
在喝三瓶啤酒或半斤白酒后多长时间内驾车会违反标准，喝酒时间长短不同情况会怎样？
分析当司机喝酒后何时血液中的酒精含量最高；
如果该司机想天天喝酒还能否开车；
结合模型和国家新标准写一篇关天司机如何驾车的忠告。
二、模型假设
酒精从胃转移到体液的速率与胃中的酒精浓度成正比。
酒精从体液转移到体外的速率与体液中的酒精浓度成正比。
酒精从胃转移到体液的过程中没有损失。
测量设备完善，不考虑不同因素所造成的误差。
酒精在体液中均匀分布。
三、符号说明
k:酒精从体外进入胃的速率；
f(t):酒精从胃转移到体液的速率；
f(t):酒精从体液转移到体外的速率；
X(t):胃里的酒精含量；
Y(t)：体液中酒精含量；
V:体液的容积；
K:酒精从胃转移到体液的转移速率系数；
K:酒精从体液转移到体外的转移速率系数；
C(t):体液中的酒精浓度。
：短时间喝酒情况下进入胃中的初始酒精量。
T：较长时间喝酒所用的时间或达到浓度最大值所需时间。
四、模型的分析与建立
（一）、模型分析：
假设酒精先以速率进入胃中，然后以速率从胃进入体液，再以速率f(t)从体液中排到体外。根据假设可以建立如图一所示的带有吸收室的单房室系统，其中胃为吸收室，体液为中心室。

图一
（二）模型建立：
用x(t)与y(t)分别表示酒精在胃、体液中的酒精量，c(t)表示酒精在体液中的浓度。根据酒精从胃进入体液的速度f(t)与胃中的酒精量成正比，速率系数为K；酒精从血液中排出的速率f(t)与血液中的酒精量y(t)成正比，速率系数为K，可以建立方程如下：
（1）（2）（3）
将（1）式代入（3）式可得：（4）

通过移项，上式可以转化为；（5）
利用一阶线性常微分方程的常数变易法对（5）式求解，可以得到；
（6）
又因为，联合（6）式可得：（7）

又对中心室（即体液）可建立方程组如：（8）
将（2）式代入（8）式可得：
将上式转化为：
因为，将其代入上式可得到：
（9）
求解（11）式可得；（10）
（其中，，）
又酒精浓度为酒精量与体液容积之比，，即：（11）
（其中，，，）。
（三）模型的讨论：
当酒是在较短时间内喝时此时有，，。
因为有，，
所以经计算整理后可得：，，
将A，B，C代入式（11）可以得到：酒在较短时间内喝下去时，体液中的酒精浓度与时间的函数关系式如下所示：
（12）
（其中）当t比较大时，显然K1>>K2,因此可认为：

利用数表一：
时间(小时)0.250.50.7511.522.533.544.55酒精含量153437.5414138.534342925.52520.5时间(小时)678910111213141516酒精含量1917.51412.597.5653.53.52（表一）
通过Matlab进行曲线拟合可得：，
根据查阅资料可知：一瓶啤酒的酒精量一般为640ml，密度为810mg/ml酒精浓度为84.5%所以两瓶啤酒的酒精总量
由于体重为70kg,体重的65%左右，体液密度为1.05mg/ml，所以可得体液的总体积为毫克/百毫升。由可求得：。
可得短时间内喝下两瓶啤酒时血液中的酒精含量与时间的关系式如下；
（13）
用Matlab软件画出图形为：

（图二：拟合曲线）
本程序是对图二曲线的拟合Matlab程序：
>>t=[1.522.533.544.55678910111213141516];
y1=[82776868585150413835282518151210774];
y2=log(y1);
polyfit(t,y2,1);
ans=
-0.19404.77