数学模型-饮酒驾车模型

一、实验目的理解数据拟合根本内容，掌握Matlab软件求解数据拟合的根本方法二、实验用仪器设备、器材或软件环境Matlab软件三、求解问题据报载，2003年全国道路交通事故死亡人数为10.4372万，其中因饮酒驾车造成的占有相当的比例。针对这种严重的道路交通情况，国家质量监督检验检疫局2004年5月31日发布了新的《车辆驾驶人员血液、呼气酒精含量阈值与检验》国家标准，新标准规定，车辆驾驶人员血液中的酒精含量大于或等于20毫克／百毫升，小于80毫克／百毫升为饮酒驾车〔原标准是小于100毫克／百毫升〕，血液中的酒精含量大于或等于80毫克／百毫升为醉酒驾车〔原标准是大于或等于100毫克／百毫升〕。大李在中午12点喝了一瓶啤酒，下午6点检查时符合新的驾车标准，紧接着他在吃晚饭时又喝了一瓶啤酒，为了保险起见他呆到凌晨2点才驾车回家，又一次遭遇检查时却被定为饮酒驾车，这让他既懊恼又困惑，为什么喝同样多的酒，两次检查结果会不一样呢？请你参考下面给出的数据〔或自己收集资料〕建立饮酒后血液中酒精含量的数学模型，并讨论以下问题：1.对大李碰到的情况做出解释；2.在喝了3瓶啤酒或者半斤低度白酒后多长时间内驾车就会违反上述标准，在以下情况下答复：酒是在很短时间内喝的；酒是在较长一段时间〔比方2小时〕内喝的。参考数据1.人的体液占人的体重的65%至70%，其中血液只占体重的7%左右；而药物〔包括酒精〕在血液中的含量与在体液中的含量大体是一样的。2.体重约70kg的某人在短时间内喝下2瓶啤酒后，隔一定时间测量他的血液中酒精含量〔毫克／百毫升〕，得到数据如下：时间(小时)12345酒精含量306875828277686858515041时间(小时)678910111213141516酒精含量3835282518151210774数学模型1.模型假设1、假设每瓶啤酒内的酒精含量一定。2、假设食物不影响人体对酒精的吸收。3、假设人体血液体积一定。4、假设人体对酒精的吸收、消化、排泄功能正常。5、假设人体血液中的酒精量与喝入的酒精量成线性关系。6、假设酒精进入中心室以后直接排出体外。2.符号说明：进入吸收室中的酒精含量〔单位：mg〕。：吸收室中的酒精转移到中心室的速率系数。：中心室的酒精排出体外的速率系数。：长时间内喝酒，酒精进入中心室的速率。：酒精在时刻转移到中心室的速率。：时间〔单位：时〕。：时刻血液中的酒精浓度。：吸收室在时的酒精含量。：吸收室在时刻的酒精含量。：中心室在时刻的酒精含量。：长时间内喝酒中心室在时刻血液中的酒精含量。：人体中血液的体积。：中心室中血液的体积。模型一〔快速饮酒问题〕〔一〕、模型建立与求解快速饮酒与口服或肌肉注射药物的过程相似,这就相当于酒精进入中心室前先有一个将酒精吸收入血液的过程可以简化为有一个吸收室。如图一于是满足(1)而药物进入中心室的速率为(2)将方程〔1〕的解代入〔2〕式得(3)又由假设得出微分方程(4)与血液浓度、房室容积之间显然有关系式=〔5〕将〔5〕式代入〔4〕式得〔6〕用MATLAB解此微分方程，可得〔7〕即〔8〕因为都是常数故〔8〕式可以简化为都为待定参数，其中〔二〕、参数估计将血样分成两局部，对于比拟大时，可近似有两边取对数得利用MATLAB将数据拟合得出那么对于其余数据进行以下处理t=[0.250.50.7511.522.533.544.5567]t=Columns1through130.25000.50000.7500Column14>>c=[3068758282776868585150413835]c=3068758282776868585150413835>>polyfit(t,c4,1)ans=>>exp(4.8708)ans=同理得出从而得出快速喝入两瓶啤酒，血液中酒精含量与时间的关系式：经过小时后再喝入两瓶啤酒，可以进一步改良模型得血液中酒精含量与时间的关系式：模型二〔慢速饮酒问题〕〔一〕、模型建立与求解慢速饮酒与恒速静脉滴注过程相似，我们假设血液中酒精是比拟长时间〔两小时内〕恒速进入中心室的，其速率为，和初始条件为：如图三，可得微分方程〔9〕又因为原理图原理图求导得〔10〕把〔10〕代入〔9〕可得〔11〕即〔12〕可以利用MATLAB解微分方程得〔13〕又因为〔14〕所以〔15〕把〔15〕代入〔13〕得〔16〕两小时以后，中心室没有酒精注入即〔17〕把〔17〕代入〔9〕可得〔18〕所以就有即〔19〕可以利用MATLAB解微分方程得〔20〕为常数，其值，〔21〕把〔21〕代入〔20〕可得〔22〕所以慢速饮酒时血液中的酒精含量与时间的关系为〔23〕〔二〕、参数估计据我们调查后得到两瓶啤酒的酒精含量为：所以得〔24〕把〔24〕代入〔23〕可得血液中的酒精含量与时间的关系式五．实验步骤、编程及调试方法图一的程序：>>t=[8910111213141516];>>c=[282518151210774];>>c1=log(c);>>polyfit(t,c1,1)ans=>>exp(5.2477)ans=>>c=190.1285*exp(-0.2316*t)图二的程序：>>c2=190.1285*exp(-0.2316*t)c2=Columns1through13Column14实验数据处理及结果分析七．分析、检验和结论;1、对大李遇到的情况做如下解释大李喝过一瓶啤酒，经过6小时后，代入模型可得其血液中的酒精含量为毫克/百毫升，所以符合新的国家标准。随后又喝一瓶啤酒，经过7小时后，代入改良后的模型可得其血液中的酒精含量为毫克/百毫升，所以被定为饮酒驾车。2.在喝3瓶啤酒的情况下酒是在很短的时间内喝的程序：>>c3=c2-cc3=Columns1through13Column14根据所建立的模型一，画出其散点图〔图四〕如下根据上图可知，在11小时内驾车违反上述标准。〔2〕酒是在较长的时间内喝的程序：>>c4=log(c3)c4=Columns1through13Column14根据所建立的模型二，画出其散点图如下代入模型计算可得慢速喝酒以后在10小时内开车是违反规定的，也可以由上图得出。八、心得体会求解数据拟合问题，最主要的就是建模，模型的建立是解决问题的关键。本次建