饮酒驾车的数学建模.docxVIP

饮酒驾车的数学模型

学院：数学学院

姓名：赵景阳

班级：15-数学四班

学号【摘要】

本文的目的在于，通过对人饮酒后体内酒精含量进行建模，然后根据所建模型，对相关问题进行分析和处理，并予以解决。本文主要根据假设合理条件，用常微分方程建立酒精在人体内的变化模型。以时间为变量，分类讨论酒精在人体内的变化。最后，根据国家酒驾标准，结合所建立的模型，给司机朋友发出忠告。

【关键词】房室系统、MATLAB、酒后驾车，常微分方程。

一、问题重述

小王，12点喝一瓶啤酒，18:00被检查合格，吃晚饭喝一瓶啤酒，夜里2点，开车回家。

讨论问题：

（1）如果小王凌晨2点驾车上路遇到酒驾检查，问他能否顺利通过?

（2）喝3瓶啤酒，隔多久开车会违反标准，并回答：

1）酒是在很短时间内喝的；

2）较长一段时间内喝的。（2小时内）

3）估计体内酒精含量达到MAX的确切时间。

4）根据你的模型论证：如果天天喝酒，是否还能开车？

5）提出忠告。

参考数据

1.国家标准：驾驶员血液的酒精含量≥20毫克／百毫升，80毫克／百毫升为饮酒驾车，≥80毫克／百毫升为醉酒驾车。

2.体液占人体重的65%至70%，

3.体重70kg人短时间内喝下2瓶啤酒后其体内酒精含量（毫克／百毫升），数据如下：

时间(小时) 0.250.5 0.75 1 1.5 2 2.5 3 3.5 4 4.5 5

酒精含量 30 68 75 82 82 77 68 68 58 51 50 41

时间(小时) 6 7 8 9 10 11 1213 14 15 16

酒精含量 38 35 28 25 18 15 12 10 7 7 4

二、模型假设

dy

因为f1t=（k

dytdt+k2y

解得：

y

=c2

（其中A2=k

浓度=酒精量/体液容积:

ct

（其中，c3=c

（三）模型的讨论：

当酒是在较短时间内喝时

此时有X0=D0=

因为：

计算得：A

代入上式：

c

=-

=A[e-

(其中A=-B3=

故ct≈Ae-k

现利用参考数据：

时间(小时) 0.250.5 0.75 1 1.5 2 2.5 3 3.5 4 4.5 5

酒精含量 30 68 75 82 82 77 68 68 58 51 50 41

时间(小时) 6 7 8 9 10 11 1213 14 15 16

酒精含量 38 35 28 25 18 15 12 10 7 7 4

通过Matlab进行曲线拟合可得：A=118.5459k

一瓶啤酒的酒精量:640ml，密度:810mg/ml酒精浓度:4.5%，两瓶啤酒的酒精量

体重为70kg,体液质量为体重的65%左右，体液密度为1.05mg/ml，所以可得体液的总体积为v

又A=k1D0

综上:此条件酒精含量与时间关系式：

ct=118.5459[

用matlab作图：

（图二：拟合曲线）

较长时间内时：

当t时，可得：

dxt

此时k0=

则：

y

(其中A2=k0k

由上可得：

ct

(其中c3

即：

c

=A

得出：

A3=

B

所以：cT=

=B

=Be

tT时，中心室酒精浓度与时间关系式：

ct=k

x

cT=k0

综上所述，当tT：

ct

五、问题求解

问题一：

设小王第一次喝酒在短时间内完成，我们可以得到关系式：

ct=k1

已求解：k

ct

当时，得ct=18.2778毫克/

图三

小王第二次喝酒满足在较长时间内喝酒的条件，关系式：

ct=A[e-k2

其中，A=59.27295,

t=14，ct=20.3618毫克/百毫升，

图四

问题二：

较短时间条件下，根据模型，关系式：

c

由于t十分小，x

则：

c

又c

ct

三瓶啤酒，D

v0=70×65%×

所以：

当时，可求得t=11.3。

故此条件下，经过11.3小时后开车,不是饮酒驾车。

（2）当较长时间内喝酒时，体液内洒精含量与时间关系式：

此时，，，

因为已知，，，百毫升，

，t=13.5小时，在较长时间（如二个小时）喝下三瓶啤酒后，只有在13.5小时后开车，就不会违反规定。

问题三：

短时间内喝酒时

根据模型可知：

当c`t=0时

当时，取得峰值

（2）、当在较长时间内喝酒时

时，

求导得：

c

c

由K1大于K2知，

中心室酒精浓度