饮酒驾车

南昌航空大学数学与信息科学学院实验报告课程名称：数学模型实验名称：饮酒驾车实验类型：验证性□ 综合性 □ 设计性□实验室名称：数学实验室班级学号：学生姓名：任课教师（教师签名）：成绩：实验日期：2012年5月6日一、实验目的理解数据拟合基本内容，掌握Matlab软件求解数据拟合的基本方法。

二、实验用仪器设备、器材或软件环境Matlab软件，电脑三、问题据报载，2003年全国道路交通事故死亡人数为10.4372万，其中因饮酒驾车造成的占有相当的比例。

针对这种严重的道路交通情况，国家质量监督检验检疫局2004年5月31日发布了新的《车辆驾驶人员血液、呼气酒精含量阈值与检验》国家标准，新标准规定，车辆驾驶人员血液中的酒精含量大于或等于20毫克／百毫升，小于80毫克／百毫升为饮酒驾车（原标准是小于100毫克／百毫升），血液中的酒精含量大于或等于80毫克／百毫升为醉酒驾车（原标准是大于或等于100毫克／百毫升）。

大李在中午12点喝了一瓶啤酒，下午6点检查时符合新的驾车标准，紧接着他在吃晚饭时又喝了一瓶啤酒，为了保险起见他呆到凌晨2点才驾车回家，又一次遭遇检查时却被定为饮酒驾车，这让他既懊恼又困惑，为什么喝同样多的酒，两次检查结果会不一样呢？请你参考下面给出的数据（或自己收集资料）建立饮酒后血液中酒精含量的数学模型，并讨论以下问题：1. 对大李碰到的情况做出解释；2. 在喝了3瓶啤酒或者半斤低度白酒后多长时间内驾车就会违反上述标准，在以下情况下回答：1）酒是在很短时间内喝的；2）酒是在较长一段时间（比如2小时）内喝的。

参考数据1. 人的体液占人的体重的65%至70%，其中血液只占体重的7%左右；而药物（包括酒精）在血液中的含量与在体液中的含量大体是一样的。

2. 体重约70kg的某人在短时间内喝下2瓶啤酒后，隔一定时间测量他的血液中酒精含量（毫克／百毫升），得到数据如下：四、数学模型模型假设：1、假设每瓶啤酒内的酒精含量一定。

2、假设食物不影响人体对酒精的吸收。

3、假设人体血液体积一定。

4、假设人体对酒精的吸收、消化、排泄功能正常。

5、假设人体血液中的酒精量与喝入的酒精量成线性关系。

6、假设酒精进入中心室以后直接排出体外。

符号说明：D ：进入吸收室中的酒精含量（单位：mg ）。

0k ：吸收室中的酒精转移到中心室的速率系数。

k：中心室的酒精排出体外的速率系数。

2k ：长时间内喝酒，酒精进入中心室的速率。

()t f 0：酒精在t 时刻转移到中心室的速率。

t ：时间（单位：时）。

()t c ：t 时刻血液中的酒精浓度。

()00x ：吸收室在0=t 时的酒精含量。

()t x 0：吸收室在t 时刻的酒精含量。

()t x 1：中心室在t 时刻的酒精含量。

()t x ：长时间内喝酒中心室在t 时刻血液中的酒精含量。

V：人体中血液的体积。

1V ：中心室中血液的体积模型一（段时间饮酒问题）（一）、模型建立与求解快速饮酒与口服或肌肉注射药物的过程相似,这就相当于酒精进入中心室前先有一个将酒精吸收入血液的过程可以简化为有一个吸收室。

如图一于是)(0t x 满足⎩⎨⎧=-='00000)0()(D x x k t x (1) 而药物进入中心室的速率为)()(000t x k t f = (2)将方程（1）的解代入（2）式得()tk eD k t f 0000-= (3)又由假设得出微分方程)()()(10'1t kx t f t x -= (4))(1t x 与血液浓度)(1t c 、房室容积1V 之间显然有关系式 )(1t x =)(1t c 1V （5）将（5）式代入（4）式得100'0)(V eD k kc t c tk -+-= （6）用MATLAB 解此微分方程，可得kttk k ec ek k V D k t c -+-+-=])([)(1)(01000 )(0k k ≠ （7）即kttk ec ek k V D k t c --+-=101000)()( （8）因为k V D k ,,,10都是常数故（8）式可以简化为t k kt Be Ae c 0---=B A ,都为待定参数，其中)(0101k k V D k B -=（二）、参数估计将血样分成两部分，对于t 比较大时， 可近似有ktAet c -=)(两边取对数得ktA t c -=ln )(ln利用MATLAB 将数据拟合 （计算程序见附录一）得出2316.0,2477.5==k A则t e t c 2316.01285.190)(-= 对于其余数据进行以下处理t k kt Be t c Ae t c0)()(~--=-= 同理得出4252.130=B 4775.00=k 从而得出快速喝入两瓶啤酒，血液中酒精含量与时间t 的关系式：tteet c 4775.02316.04252.1301285.190)(---=经过n 小时后再喝入两瓶啤酒，可以进一步改进模型得血液中酒精含量与时间t 的关系式：)(4775.0)(2316.04775.02316.04252.1301285.1904252.1301285.190)(n t n t tteeeet c -------+-=图二为模拟数据（“o ”）与所提供的数据（“+”）的散点对照图。

显而易见，两者基本吻合。

模型二（慢速饮酒问题）（一）、模型建立与求解慢速饮酒与恒速静脉滴注过程相似，我们假设血液中酒精是比较长时间（两小时内）恒速进入中心室的，其速率为2k ，)(0t f 和初始条件为：0)0(,)(20==c k t f如图三，可得微分方程)()()(0't kx t f t x -= （9）又因为),()(1t c V t x =求导得)()('1't c V t x = （10）把（10）代入（9）可得)()()(10't c kV t f t Vc -= （11）即)()()(10't kc V t f t c -=（12）可以利用MATLAB 解微分方程得ktec KV t f t c -+=110)()( （13）又因为0)0(=c （14）所以kV t f c 101)(-= （15）把（15）代入（13）得)1()()(10ktekV t f t c --=（16）两小时以后，中心室没有酒精注入即0)2(0=-t f （17） 把（17）代入（9）可得)2()2('--=-t kx t x （18）所以就有)2()2(1'1--=-t c kV t c V即)2()2('--=-t kc t c （19）k t c t c -=--)2()2('可以利用MATLAB 解微分方程得)2(2)2(--=-t k eC t c （20）2C 为常数，其值)2(2C C =，)1()()2(210kekV t f C --=（21）把（21）代入（20）可得)2(210)1()()2(----=-t k keekV t f t c （22）所以慢速饮酒时血液中的酒精含量与时间的关系为⎪⎪⎩⎪⎪⎨⎧≥-≤<-=----)2()1()()20()1()()()2(21010t e e kV t f t e kV t f t c t k k kt（23）（二）、参数估计据我们调查后得到两瓶啤酒的酒精含量为：)(384001000203.0640mg D =⨯⨯⨯=所以得48.2202316.04728138400)(0=⨯⨯⨯=Vkt f （24）把（24）代入（23）可得血液中的酒精含量与时间的关系式⎩⎨⎧≥-≤<-=----)2()1(48.220)20()1(48.220)()2(2316.04632.02316.0t e e t e t c t t （25）用MATLAB画出散点图如上五、实验方法、步骤、及编程掌握polyfit求解线性最小二乘问题的用法，及lsqcurvefit和lsqnonlin求解非线性最小二乘问题的用法。

掌握数据拟合运用matlab来编程。

程序代码：（1）x=[0.25,0.5,0.75,1,1.5,2.0,2.5,3,3.5,4,4.5,5,6,7,8,9,10,11,12,13,14,15,16];y=[30,68,75,82,82,77,68,68,58,51,50,41,38,35,28,25,18,15,12,10,7,7,4];plot(x,y,'*')plot(x,y,'*',x,y)（2）x=[0.25,0.5,0.75,1,1.5,2.0,2.5,3,3.5,4,4.5,5,6,7,8,9,10,11,12,13,14,15,16];y=[3.75,17,28.125,41,61.5,77,85,102,101.9,102,112.5,102.5,114,122.5,112,112.5,90 ,82.5,72,65,49,52.5,32];plot(x,y,'*')>> plot(x,y)>> x1=[0.25,0.5,0.75,1,1.5,2.0,2.5,3];y1=[3.75,17,28.125,41,61.5,77,85,102];>> a=polyfit(x1,y1,2)a =-6.4863 55.7627 -9.2847>> x2=0.25:0.5:3;>> y2=-6.4863*x2.^2+55.7627*x-9.2847;??? Error using ==> +Matrix dimensions must agree.>> x2=0.25:0.5:3;y2=-6.4863*x2.^2+55.7627*x2-9.2847;>> plot(x,y,'*',x2,y2)>> polyfit(x1,x2,3)??? Error using ==> polyfitX and Y vectors must be the same size.>> polyfit(x1,y1,3)ans =2.1562 -16.9706 69.6712 -13.6200>> y2=2.1562*x2.^3-16.9706*x2.^2+69.6712*x2-13.6200;>> plot(x,y,'*',x2,y2)>> x=[0.25,0.5,0.75,1,1.5,2.0,2.5,3,3.5,4,4.5,5,6,7,8,9,10,11,12,13,14,15,16]; y=[3.75,17,28.125,41,61.5,77,85,102,101.9,102,112.5,102.5,114,122.5,112,112.5,90 ,82.5,72,65,49,52.5,32];plot(x,y,'*')>> plot(x,y)>> x1=[0.25,0.5,0.75,1,1.5,2.0,2.5,3,3.5,4];y1=[3.75,17,28.125,41,61.5,77,85,102,101.9,102];>> a=polyfit(x1,y1,2)a =-7.7319 59.3789 -10.9764>> x2=0.25:0.5:4;>> y2=-7.7319*x2.^2+59.3789*x2.-10.9764;??? y2=-7.7319*x2.^2+59.3789*x2.-10.9764;|Error: "identifier" expected, "-" found.>> x2=0.25:0.5:4;y2=-7.7319*x2.^2+59.3789*x2-10.9764;>> plot(x,y,'*',x2,y2)>> polyfit(x1,y1,3)ans =-0.2926 -5.8683 56.1772 -9.7604>> y2=-0.2926*x2.^3-5.8683*x2.^2+56.1772*x2-9.7604;>> plot(x,y,'*',x2,y2)>> x3=[4.5,5,6,7];>> y3=[102,102.5,114,122.5];>> b=polyfit(x3,y3,2)b =0.9347 -1.9849 90.9799>> x3=4:0.5:7;>> y3=0.9347*x3.^2-1.9849*x3+90.9799;>> plot(x,y,'*',x2,y2,x3,y3)>> polyfit(x3,y3,3)ans =0.0000 0.9347 -1.9849 90.9799>> x4=4:0.5:7;>> y4=0.0000*x4.^3+0.9347*x4.^2-1.9849*x4+90.9799;>> plot(x,y,'*',x2,y2,x4,y4)>> x5=[7,8,9,10,11,12,13,14,15,16];>> y5=[122.5,112,112.5,90,82.5,72,65,49,52.5,16];>> c=polyfit(x5,y5,2)c =-0.3201 -3.4140 161.6318>> x6=7:0.5:16;>> y6=-0.3201*x6.^2-3.4140*x6+161.6318;>> plot(x,y,'*',x2,y2,x4,y4,x6,y6)>> polyfit(x6,y6,3)ans =0.0000 -0.3201 -3.4140 161.6318>> x6=7:0.5:16;>> y6=0.0000*x6.^3-0.3201*x6.^2-3.4140*x6+161.6318;>> plot(x,y,'*',x2,y2,x4,y4,x6,y6)六、计算结果1、对大李遇到的情况做如下解释大李喝过一瓶啤酒，经过6小时后，代入模型可得其血液中的酒精含量为20c毫=19<.97克/百毫升，所以符合新的国家标准。