《全国大学生数学建模竞赛2004优秀论文设计C、D题》由会员分享，可在线阅读，更多相关《全国大学生数学建模竞赛2004优秀论文设计C、D题（69页珍藏版）》请在金锄头文库上搜索。
1、wordC题之一全国一等奖酒精在人体内的分布与排除优化模型某某工学院，袁孟强，王哲，X莉指导教师：数模辅导组摘要：酒精进入机体后，随血液运输到各个器官和组织，不断的被吸收，分布，代谢，最终排除体外。为了研究酒精在体内吸收，分布和排除的动态过程，以与这些过程与人体反响的定量关系，本文建立了一个酒精在人体内的分布与排除优化模型，在药物动力学的一室模型的根底上，进展优化，改良，分别建立了酒精在人体内分布的房室模型和房室模型，以与酒精在人体内的静态排除模型和动态排除模型，导出模型的体液酒精浓度的状态函数，用常数交叉拟合方法，采用VB编写程序，得到两个重要系数和。根据此模型，计算的体液酒精浓度理论值与实
2、验值十分相符，并很好地解释了给出的所有问题，得到一些有价值的结论。关键词：房室模型，排除模型，体液酒精浓度，动态和静态的转换酒精在人体内的分布与排除优化模型一、 问题的重述国家质量监视检验检疫局2004年5月31日发布了新的车辆驾驶人员血液、呼气酒精含量阈值与检验国家标准，新标准规定，车辆驾驶人员血液中的酒精含量大于或等于20毫克百毫升，小于80毫克百毫升为饮酒驾车原标准是小于100毫克百毫升，血液中的酒精含量大于或等于80毫克百毫升为醉酒驾车原标准是大于或等于100毫克百毫升。大李在中午12点喝了一瓶啤酒，下午6点检查时符合新的驾车标准，紧接着他在吃晚饭时又喝了一瓶啤酒，为了保险起见他呆到凌
3、晨2点才驾车回家，又一次遭遇检查时却被定为饮酒驾车，这让他既懊恼又困惑，为什么喝同样多的酒，两次检查结果会不一样呢？参考下面给出的数据或自己收集资料建立饮酒后血液中酒精含量的数学模型，并讨论以下问题：1. 对大李碰到的情况做出解释；2. 在喝了3瓶啤酒或者半斤低度白酒后多长时间内驾车就会违反上述标准，在以下情况下回答：1） 酒是在很短时间内喝的；2） 酒是在较长一段时间比如2小时内喝的。3. 怎样估计血液中的酒精含量在什么时间最高。4. 根据你的模型论证：如果天天喝酒，是否还能开车？ 5. 根据你做的模型并结合新的国家标准写一篇短文，给想喝一点酒的司机如何驾车提出忠告。参考数据1. 人的体液占
4、人的体重的65%至70%，其中血液只占体重的7%左右；而药物包括酒精在血液中的含量与在体液中的含量大体是一样的。2. 体重约70kg的某人在短时间内喝下2瓶啤酒后，隔一定时间测量他的血液中酒精含量毫克百毫升，得到数据如下(表1)：时间(小时)12345酒精含量306875828277686858515041时间(小时)678910111213141516酒精含量3835282518151210774二、 模型假设1、 酒精的转移速率，与向体外排除的速率， 与该室的血酒浓度成正比。2、酒精的转移速率，与向体外的排除速率，与时间有关，与空间人体的各个局部无关。3、中心室与体外有酒精交换，与酒精从体
5、外进入中心室，最后又从中心室排出体外。与转移和排除的数量相比，酒精的吸收可以忽略。三、模型建立与求解房室模型在短时间内喝下酒精量为在短时间内喝下酒精量为，酒精进入胃，人体吸收酒精，然后排除出体外。吸收酒精的过程相当于酒精进入体液中心室的过程，全过程可以简化为如下图：胃室中心室排除体外建模过程：短时间内进入胃的酒精；为胃室吸收室进入中心室的转移速率系数由人体机能确定的常数；是t时刻胃室的酒精；其微分方程为： 是t时刻进入中心室的酒精，其微分方程为： 2 酒精进入中心室的速率为： 3将方程1的解代入3得： 4房室模型在较长一段时间内喝酒假设在较长的一段时间内喝下的酒是匀速进入胃室，如此简化如如下图
6、： 常数胃室排除中心室建模过程：为酒精进入胃的速率： ,为喝酒时间。为酒精从中心室排除体内的速率为酒精进入中心室的速率为胃室进入中心室的转移速率由人体机能确定的常数为是酒从中心室向外排除的速率系数。是t时刻胃室的酒精，微分方程为： 5 6 是t时刻进入中心室的酒精 将方程5的解代入6得： 7 8静态排除模型与房室模型I配套的静态酒精排除模型中心室的血酒浓度；V人体体液量和人体血液量；酒精进入中心室的速率： 中心室的酒精量；微分方程为： 9酒精从中心室向体外排除的速率系数由人体机能确定的常数由方程9得： 10对应的通解为：微分方程的解为：令 得：.根据参考数据表1，：短时间内进入胃的酒精,人体体
7、液量V和一批实验数据，()。用交叉常数拟合原理在VB环境中编写程序，利用该程序算出两个重要系数和，给出模型的状态函数.假如初始值设为，如此动态排除模型与房室模型配套建立动态酒精排除模型.中心室的血酒浓度；V人体体液量和人体血液量；酒精进入中心室的速率：中心室的酒精量；微分方程为： 11 12酒精从中心室向体外排除的速率系数由人体机能确定的常数.由方程11得： 13对应的通解为：.将与代入得： 14四、 酒过程的描述1、在短时间内喝下酒精量为用静态排除模型描述：.2、在较长一段时间内喝酒用动态排除模型描述喝酒的过程，用静态排除模型描述酒后的过程.即先用状态函数描述喝酒的过程，然后用状态函数描述酒
8、后的过程.3、天天喝下酒精量为用动态排除模型描述第一天喝酒喝一小时的过程，用静态排除模型描述酒后23小时内的的过程，用动态排除模型描述第二天喝酒喝一小时的过程，用静态排除模型描述第二天酒后23小时内的的过程，再用动态排除模型描述第三天喝酒喝一小时的过程，用静态排除模型描述第三天酒后23小时内的的过程五、参数的选择一瓶啤酒的酒精量：.人体体液：速率系数：两小时内慢慢喝下两瓶啤酒的输入率：六、拟合效果.中心室中酒精含量毫克百毫升浓度状态函数拟合效果图效果图显示拟合程度极高，说明参数的选择与客观情况相符合.七、问题分析问题1的分析：大李在中午喝了一瓶啤酒下午6点检查时，利用静态排除模型：由浓度状态函
9、数在时的浓度：,可知大李此时符合新的驾车标准.紧接着他在吃晚饭时又喝了一瓶啤酒，为了保险起见他呆到凌晨2点才驾车回家.检测时距晚饭喝的那瓶啤酒已过了八小时，其胃中酒精与体液中的酒精含量分别为：由浓度状态函数可知，检查时被定为饮酒驾车.之所以被定为饮酒驾车，关键是此时方程中的初始条件而不是第一次喝酒的.问题2的分析：在很短时间内喝了3瓶啤酒,问多长时间内驾车就会违反新的驾车标准. 根据静态排除模型：3瓶啤酒的酒精量：；浓度状态函数为 .代入数据得：. 15假如违反新的驾车标准,如此：,根据人体体液酒精浓度曲线图可知，喝酒后约在t=0.07小时42秒与t=11.24小时之内驾车会违反新标准。将t=
10、0.07，t=11.24分别代入方程15检验得出：通过验证，证明观测值根本接近实际值。问题3的分析：在较长一段时间比如2小时内喝3瓶啤酒，多长时间内驾车就会违反新标准.假设匀速喝酒，如此此过程分两阶段：1在喝酒过程中，多少时间后驾车会违反新标准。2喝酒之后，多长时间内驾车会违反新标准。根据动态模型： 16 1 由人体体液酒精浓度曲线观测出：当t大于0.62小时约37.2分钟时，体液酒精浓度大于毫克，把酒精数据，代入方程15检验得：通过数据验证，证明观测值根本接近实际值。2由阶段1的数据算出，当t=2小时即停止喝酒时人体酒精为： 此时，根据胃里剩余的酒精方程：得出： 以2小时作为零时刻，设状态初
11、始值为,如此取，此时静态排除模型为：根据此模型，可得曲线：可观测出：当进展检测得：毫克/百毫升通过数据验证，证明观测值根本接近实际值。结论：假如连续2小时均匀喝下3瓶啤酒，如此在0.62与13.67小时之内驾车会违反新标准。问题4的分析：估算出血液中酒精含量何时达到峰值。 最高值的估计要分两种情况讨论：1、 据静态排除模型一次喝完的情况 酒精进入机体之前，胃里的酒精浓度为0，当喝酒后，胃里的酒精浓度逐渐增大，因为存在着酒精进入中心室的过程，根据人体机理可知，。当胃里的酒精浓度等于中心室的酒精浓度时，此时血液中酒精含量达到最高值。由静态排除模型得：假如要使的值达到最高，如此要求即：代入数据得：即
12、：喝酒1.36小时后，血液中酒精含量最高。（1） 根据动态排除模型酒是在很长一段时间内喝完 喝酒时，由于酒量在不断的增加，根据人体机理可知，胃中的酒精浓度始终是大于中心室的酒精浓度，即，易知血液中的酒精浓度是一个递增的函数。那么在这阶段中，血液的酒精含量的最高值为停止喝酒时血液中的酒精含量。运用动态排除模型：可计算出最值。喝完酒后，胃中的酒精浓度仍然大于中心室的酒精浓度，此时，.随着时间的推移，胃中的酒精浓度在不断的减少，而血液的酒精浓度不断增加，这就使得当时，血液的酒精浓度达到最高值即，假如代入具体数据就可算出最大值。问题5的分析： 根据我们建立的模型论证：如果天天喝酒，是否还能开车？天天喝一瓶酒的数学模型见图一个人天天喝酒，如果每天只喝一瓶啤酒第二天即12个小时后体内的血酒浓度远远小于20毫克/百毫升所以仍能开车；天天喝两瓶酒的数学模型见图：一个人天天喝酒，如果每天只喝两瓶啤酒第二天即12个小时后体内的血酒浓度远远小于20毫克/百毫升所以仍能开车；天天喝三瓶酒的数学模型见图：如果每天喝三瓶啤酒，12小时后
