饮酒驾车的数学模型精讲.docVIP

饮酒驾车的数学模型 （CUMCM-2004C题） 一 、摘要 本题是关于一个饮酒驾车的数学模型。因为酒精在一个房呈均匀分布，从吸收室到中央室按照一定的规律进行吸收和排除。所以根据不同时刻的吸收与排除情况，为了研究酒精的吸收和排除的动态过程，我们对市场上酒的分析调查为参考资料。以传统的常微分方程理论来建立控制饮酒驾车模型方程与曲线拟合的模型,近似于房室模型来解决.通过matlab数学软件求解模型，得到相关结果。最后从模型方程跟实际对比分析中找出实际与理论的差异。 关 健 词 ：常微分方程 曲线拟合 房室模型 二、 问题的提出 在2003年全国道路交通事故死亡数字的10.4372万中饮酒所造成的事故占着相当大的比例。针对这一比例所造成的事故国家质量验检局与2004年5月31日发布的新的《车辆驾驶人员血液、呼气酒精含量阈值与检验》国家标准规定了驾驶人员血液中的酒精含量。新标准规定，车辆驾驶人员血液中的酒精含量大于或等于20毫克／百毫升，小于80毫克／百毫升为饮酒驾车（原标准是小于100毫克／百毫升），血液中的酒精含量大于或等于80毫克／百毫升为醉酒驾车（原标准是大于或等于100毫克／百毫升）。为了减少和预防事故发生,保证人民的生命财产的安全,我们建立模型对饮酒驾车进行分析,为政府提供一些相关资料的参考。 三、问题的分析与假设 (一) 问题分析 因为在1个小时以内酒精未达到机体最大消除力时，假设在吸收过程仍符合一级动力方式消除。因为按酒精的一般规律，酒精的清除符合零级动力学方式，所以我们可以假设在一开始喝酒的过程时，酒精的排除符合零级动力方学方式。 另一种情况就是酒在长时间内喝的，近似于口服药液。 根据表格数据我们可知，酒精在血液中的浓度随时间的变化而变化 （二） 问题假设 假设在酒精的吸收收速率及排除速率，与该室的酒精浓度成正比。 假设机体分为中心室和吸收室（如图1），且两个室的容积在过程中保持不变。 假设当酒精进入中心室时，吸收和排除的数量相比，吸收可以忽略。 假设酒精在体内作用没有性别差异,没有种族差异。 图1 一室模型 （三）参数说明 假设酒精在吸收室的浓度为 ，中央室的浓度c(t),吸收率为，消除速率为k; 假设人体中心室的容量为v和吸收室的容量为； 假设一个正常人的体重为70kg,每瓶啤酒的体积为600ml; 四：模型的建立 实验数据如表1： 时间(小时) 0.25 0.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 12 13 14 15 16 酒精含量38 35 28 25 18 15 12 10 7 7 4 表1 由上述数据分析作插值运算并绘图（见图2）：（程序见附件１） 图2 由此我们可做相应推断: 酒精在人体内的吸收和排除的曲线图中含有一个波峰,固可推出此时酒精在人体内的含量达到最大值，相应的人体内的酒精的排除速率也达到最大值。 浓度在持续快速增加了大约2个小时后,开始下降。 因为酒精的吸收类似于医学上药物动力学研究的药物代谢过程，因此我们按药物动力学建立相关模型。 酒精被人体吸收进入机体后随血液输送到全身，不断的被吸收， 分布代谢，最终排出体 。现在我们假设血液中的酒精浓度为c(t),建立一个一室模型，将整个机体看作一个整体。 （一）、酒精的吸收阶段 解得 （１） 对(1)式取对数得 记 ，， 由实验数据用最小二乘法拟合中的参数. 用matlab软件作最小二乘法拟合可得到 ， 所以得 （ 2 ）　 （二）、酒精的排出阶段 因为一开始喝酒时，就有排出过程，当酒精未达到机体最大消除能力时（主要是未超出催化生物转化的酶的饱和限时），都将按一级动力学方式消除当其量超过机体最大消除能力时，将只能按最大消除能力这一恒量进行消除，变为零级消除动力学方式，即出现消除动力学模型转换 解得 （ 3 ） 根据表1的实验数据,用最小二乘法拟合下面方程组中的参数： 得 所以有 （4） 图3 （三）、较长时间内喝酒 设t=0时刻酒精含量为d,则吸收室的酒精浓度为的微分方程为 , 其中 解得 (5) 中心室的酒精浓度与