饮酒驾车问题

《饮酒驾车问题》由会员分享，可在线阅读，更多相关《饮酒驾车问题（18页珍藏版）》请在金锄头文库上搜索。

1、 1 饮酒驾车问题的数学模型 饮酒驾车问题的数学模型 摘摘 要要 随着现在酒后驾车导致的事故频发， 本文从分析饮酒后血液中酒精含量与时 间关系的角度，建立了一个关于饮酒驾车的数学模型。本文主要采用房室模型和 微分方程组的算法对其进行分析。本文首先针对吸收和代谢两方面的效应，对饮 酒后血液血液中的酒精含量进行预估。同时首先依靠房室模型，建立一个不考虑 时间因素，单纯关于吸收和代谢的一般模型，再根据文意需要分别建立短时饮酒 和长时饮酒两个模型。接下来用常数变易法对模型进行求解，用最小二乘法并借 助于 Matlab 软件， 对模型一所得的体液中酒精浓度关于时间的表达式进行拟合， 得到了模型中吸收能力

2、系数 1 k 和代谢能力系数 2 k 的值分别为1.915和0.1925， 由此我们解决问题中模型的解。 针对问题一：得到下午 6 点的体液中酒精的含量为 16.803 毫克/百毫升，小 于国家标准规定的 20 毫克/百毫升，而凌晨 2 点的为 21.170 毫克/百毫升，大于 国家标准规定的 20 毫克/百毫升， 如此通过定量的分析合理的解释了大李遇到的 情况。针对问题二，将其分为短时和长时饮酒两种情况，带入模型求解，得到其 分别在 11.3526 和 12.4003 小时后可以开车。针对问题三：我们首先假设只喝 3 瓶啤酒，2?T，我们可以用求 )(tC 的驻点的方法来求 )(tC 的最大

3、值点即酒精含 量最高的点，最终得到其浓度最大的时间点分别为 1.5 小时和 2.5 小时。针对问 题四：在本问中我们假设每天只喝一次，是短时间进酒，两次间隔 24 小时。因 此带入模型一最终得到当一瓶啤酒中酒精的量约为 24000 毫克，每天最多可 0.375 瓶。针对问题五，本文主要普及了酒后驾车的危险程度，如何自己判断是 否可以规避就酒驾风险，过度饮酒的危害等角度，对司机朋友提出我们的建议。 提醒广大司机朋友们适度饮酒，理性饮酒。 关键词：关键词：房室模型 微分方程组 t C 驻点法 吸收和代谢 2 1 问题重述问题重述 据报载，2003 年全国道路交通事故死亡人数为 10.4372 万，

4、其中因饮酒驾 车造成的占有相当的比例。 针对这种严重的道路交通情况， 国家质量监督检验检疫局 2004 年 5 月 31 号 发布了新的车辆驾驶人员血液、呼气酒精含量阈值与检验国家标准，新标准 规定，车辆驾驶人员血液中的酒精含量大于或等于 20 毫克百毫升，小于 80 毫 克百毫升为饮酒驾车（原标准是小于 100 毫克百毫升），血液中的酒精含量 大于或等于 80 毫克百毫升为醉酒驾车（原标准是大于或等于 100 毫克百毫 升）。 大李在中午 12 点喝了一瓶啤酒，下午 6 点检查时符合新的驾车标准，紧接 着他在吃晚饭时又喝了一瓶啤酒，为了保险起见他呆到凌晨 2 点才驾车回家，又 一次遭遇检查时

5、却被定为饮酒驾车， 这让他既懊恼又困惑， 为什么喝同样多的酒， 两次检查结果会不一样呢？ 我们将参考下面给出的数据并自己收集资料建立饮酒后血液中酒精含量的 数学模型，并讨论以下问题： 1. 对大李碰到的情况做出解释； 2. 在喝了 3 瓶啤酒或者半斤低度白酒后多长时间内驾车就会违反上述标准， 在以下情况下回答： 1）酒是在很短时间内喝的； 2）酒是在较长一段时间（比如 2 小时）内喝的。 3. 怎样估计血液中的酒精含量在什么时间最高。 4. 根据你的模型论证：如果天天喝酒，是否还能开车？ 5. 根据你做的模型并结合新的国家标准写一篇短文，给想喝一点酒的司机 如何驾车提出忠告。 2 问题分析问题

6、分析 当人们喝酒时,酒精(即乙醇)就会随着血液扩散到全身的器官和组织中,然 3 后不断地被吸收、 代谢和排出.在这里我们把乙醇在血液中的浓度(即每 100mL 血 液中乙醇含量(mg)称血醇浓度,它随时间和所处范围的变化而变化. 为了建立饮酒后血液中酒精含量的数学模型,我们需要讨论乙醇在体内吸收、 分布和排出的动态过程以及此过程中的定量关系.为了达到这一目的,我们借鉴 药物动力学中的房室模型.所谓房室,是指由具有相近的药物转运速率的器官、 组 织组合而成.同一房室内各部分的药物处于动态平衡.可以认为乙醇在体内各部 位间均有较高及相近的转运速率,可在体内迅速达到分布平衡,即血醇浓度具有 均衡性.

7、但是在不同的房室则可按一定的规律进行乙醇的转移,在这里我们简单 地认为机体是一个血液均匀分布的中心， 喝酒后乙醇的扩散与分布可看作在进入 中心室之前有一个将乙醇由肠胃吸收入血液的过程将之简化为一个吸收室这种 简化在医学界和药理学界已被接受和承认，我们称之为一室模型。 本题中因为乙醇在血液中的浓度与在体液中的浓度大致相等所以我们可以 将血液与体液合并为一个房室即中心室。 下面在吸收室和中心室之间建立一个血 醇浓度的微分方程来描述其动态特征。 3 模型假设模型假设 （1） 体液总体积保持不变 （2） 在较短时间内喝酒的情况下，酒精量是瞬间进入到胃里的。 （3） 体液的总体积不变。 （4） 酒精在血

8、液中的含量与在体液中的含量大至相同。 （5） 不管喝的是什么酒，只以涉入的酒精总量纳入计算。 （6） 假设整体过程中人没有摄入任何影响代谢的药类物质和作剧烈性运 动。 （7） 人的吸收速率和代谢速率是恒定的。 （8） 忽略不同人对酒精代谢能力的差异。 4 符号说明符号说明 4 变量 含义 单位 备注 1 k 吸收能力系数 1/小时 2 k 代谢能力系数 1/小时 0 C 开始饮酒时人体体液中的酒 精浓度 毫克/百毫升 Q 人喝下酒精的总量 毫克 一瓶啤酒中酒精的量 约为 24000 毫克 0 V 人体液所占的体积 百毫升 约为 450 百毫升 t 时间 小时 T 较长时间饮酒时的持续时间 小时

9、 酒精由肠胃进入人体体液的 速率 毫克/小时 )( 1 tf 酒精由口进入肠胃的速率 毫克/小时 )(tx t 时刻体液中酒精的量 毫克 )(ty t时刻肠胃中酒精的量 毫克 )(tC 在 t 时刻人体体液（或血 液）中的酒精浓度 毫克/百毫升 5 模型建立与求解模型建立与求解 5.1 模型的建立模型的建立 根据以上分析，我们可以建立一个房室模型，如图 1： 5 图 1 房室模型结构 又考虑到吸收总量为Q ，有： 在吸收室中我们考虑的是肠胃吸收酒精进入体液所引起的酒精量的增量， 其中 用 )(tf 代表t时刻酒精由肠胃进入人体体液的速率（单位：毫克/小时），则时间 段 ,ttt? 内体液中酒精

10、量的增量可表示为： ( ) tt t f t dt ? ? 。 再以Q表示人喝下酒 精的总量，以刚开始饮酒的时刻为记时的初始时刻，即t=0，由于喝下的酒精在 时间无穷长时总会被全部吸收，故有： ( )= tt t f t dt Q ? ? 再考虑中心室，根据有关文献可知设从t 时刻到tt?时刻中心室（血液） 中乙醇变化量为 2 ()( )( )( ) tttt tt x ttx tf t dtk x t dt ? ? ? ? 即血液中乙醇变化量=吸收量-代谢量。等式两边同时除以t? ，并令0t? ? 可得： 2 ( ) ( )( ) dx t f tk x t dt ? 此外还有初值条件 0

11、(0)CC? 5.1.1 一般模型的建立一般模型的建立 由此得到一般情况下的模型： 6 由题中可知，酒可以在很短时间内喝完，也可以在一段较长的时间内喝完， 这样将有两个不同的吸收速率， 即可以得到两个不同的 )(tf 表达式， 于是可以得 到两个具体模型。第一个模型是考虑短时饮酒效应，第二个模型考虑的是长时饮 酒效应。 5.1.2 具体模型一具体模型一 (快速饮酒快速饮酒) 针对具体模型一： 该模型针对短时的饮酒效应， 即在喝酒后的很短的时间内， 酒精就全部进入肠胃，以 )(ty 表示t时刻肠胃中酒精的量，即有 0 )0(Qy?，根据 文献资料判断，可知肠胃中酒精的减少速率应该和剩余量成正比，

12、并设其比例系 数为a，因此得到微分方程： 1 ( ) ( ) dy t k y t dt ? ? 在此我们假设酒精进入人体和肠胃中酒精的减少为一个过程， 即不考虑其在 转化过程中的酒精的损耗，从而得到以下方程： 1 ( )( )f tk y t? 再结合之前通过方式模型得到的一般模型， 我们在此得到短时饮酒效应的模 型如下所示： 7 2 0 0 0 1 1 ( )( ) ( ) ( ) (0) ( )( ) ( ) ( ) (0) dC tf t k C t dtV f t dtQ CC f tk y t dy t k y t dt yQ ? ? ? ? ? ? ? ? ? ? ? ? ? ?

13、 ? ? ? ? ? ? ? ? 5.1.3 具体模型二（慢速饮酒）具体模型二（慢速饮酒） 针对具体模型二：该模型针对长期饮酒效应，可将其近似认为在持续饮酒的 过程中酒精是匀速进入肠胃的，参照模型一可有0)0(?y，在此我们引入函数 )( 1 tf来表示酒精进入肠胃的速率（单位：毫克/小时），T表示饮酒时的持续总 时间，则酒精进入肠胃的速率与整个过程中喝入的酒精量有如下关系： ? ? ? ? ? ? Tt TtTQ tf ,0 ,/ )( 0 1 而肠胃里的酒精量的变化与机体对酒精的吸收和喝入机体的酒精量都有关 系，而机体对酒精的吸收速率可以与模型一一样用 1 ( )k y t? 来表示，则时

14、间段 ,ttt? 内有： 肠胃里的酒精量的变化量 = 喝入机体的酒精量 - 机体对酒精的吸收量 即得： 1 1 ( ), ()( ) ( ), tttt tt tt t Q dtk y t dttT T y tty t k y t dttT ? ? ? ? ? ? ? ? ? ? ? 两边同除t?后让0?t取极限得： 8 1 1 ( ),( ) ( ), Q k y ttTdy t T dt k y ttT ? ? ? ? ? ? ? 而与模型一一样有： 1 ( )( )f tk y t? 再结合一般模型我们可以得到模型二如下： 0 0 0 1 1 1 ( )( ) ( ) ( )( ) (0

15、) ( )( ) ( ),( ) ( ), (0)0 T T dC tf t bC t dtV f t dtf t dtQ CC f tk y t Q k y ttTdy t T dt k y ttT y ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? 5.2 模型求解模型求解 5.2.1 具体模型一的求解具体模型一的求解 模型一求解：根据具体模型一得： 1 1 ( ) ( ) ( )( ) dy t k y t dt f tk y t ? ? ? ? ? ? ? ? 将其整理并带入一般模型中求解得到 ( )C t 与t 的关系：

16、 1 0 2 0 ( ) ( ) k t aCdC t ek C t dtV ? ? 观察其方程满足非齐次一阶线性微分方程的一般形式，带入公式进行求解， 最终得到如下公式： 9 212 1 0 012 ( )() () k tk tk t k Q C teeC e V kk ? ? ? 5.2.2 具体模型二的求解具体模型二的求解 模型二求解：根据题设，我们取 2?T。 由 1 1 ( ),2( ) 2 ( ),2 Q k y ttdy t dt k y tt ? ? ? ? ? ? ? 及0)0( 0 ?x可得： 1 11 1 2 1 (1),2 2 ( ) (1),2 2 k t kk t Q et k y t Q eet k ? ? ? ? ? ? ? ? ? ? ? 又 1 ( )( )f tk y t?，可得： 1 11 2 (1),2 2 ( ) (1),2 2 k t kk t Q et f t Q eet