实验二饮酒驾车报告.docVIP

东南大学《数学实验》报告 学号: 姓名: 成绩: 实验题目（二）：习题 6.8“饮酒驾车的药物注射模型”求解 一、实验目的 运用药物注射模型，使用曲线拟合方法，解释饮酒驾车的一些实际问题。 实验题目 国家质量监督检验检疫局2004年5月31日发布了新的《车辆驾驶人员血液、呼气酒精含量阈值与检验》国家标准，新标准规定，车辆驾驶人员血液中的酒精含量大于或等于20毫克/百毫升，小于80毫克/百毫升为饮酒驾车（原标准是小于100毫克/百毫升），血液中的酒精含量大于或等于80毫克/百毫升为醉酒驾车（原标准是大于或等于100毫克/百毫升） 某人在中午12点喝了一瓶啤酒，下午6点检查时符合新的驾车标准，紧接着他在吃饭时又喝了一瓶啤酒，为了保险起见他待到凌晨2点才驾车回家，又一次遭遇检查时却被定为饮酒驾车，这让他既懊恼又困惑，为什么喝同样多的酒，两次检查的结果会不一样呢？ （1）某人中午12点喝了一瓶啤酒，下午6点检查合格，晚饭又喝一瓶，次日凌晨2点检查未通过，请对此情况做出解释。 （2）短时间内喝啤酒3瓶多长时间之后才能驾车？ （3）怎样估计血液中的酒精含量在什么时候最高？ （4）如果天天喝酒，是否还能开车？ 二、实验原理 酒精在进入肠道之前即可被食道吸收，而且胃部吸收酒精速率速度快，所以本题采用“快速静脉注射模型”进行解决。 酒精主要存在于血液之中，所以本题应计算吸收室的血药浓度：c1(t)=A1e-αt+B1e-βt 相关系数可以通过拟合法进行求解。 模型的假设与建立 假设1：酒精从胃部向体液的转移速率、向外排出速率分别与胃部和体液中的酒精浓度成正比。? 假设2：体液总体积不变。? 假设3：入胃的酒精全部扩散到体液中。? 假设4：酒精不会从体液反渗入胃。? 假设5：酒精只通过血液排出体外。? 假设6：在较短时间内喝酒时，酒精量瞬时入胃。? 假设7：在较长时间内喝酒时，酒精量匀速入胃。 假设8：酒精正常吸收和排出，呕吐等非正常排出情况暂不予考虑。 假设9：忽略不同人体对酒精的敏感度及分解能力上的差异。 假设10：一瓶啤酒600毫升，酒精浓度为5%。 吸收过程中： 代谢过程中酒精总量守恒： 代入得： 列方程组得： （1）在短时间内饮酒： 此时肠胃对酒精的吸收速率v(t)与肠胃中酒精量Y(t)成正比：v(t)=aY(t) 且v(t)=-dY(t)/dt，方程组可化为： （2）在较长一段时间内饮酒： 这种情况下酒精均匀进入肠胃，酒精的改变量等于喝入肠胃的酒精量减去呼吸酒精的量，即： 则方程组可化为： 三、实验内容与要求 代码实现结果分析 format short g % 题中提供的某人喝了两瓶啤酒后血液酒精浓度随时间变化表 t=[0.25;0.5;0.75;1;1.5;2;2.5;3;3.5;4;4.5;5;6;7;8;9;10;11;12;13;14;15;16]; c=[30;68;75;82;84;77;70;68;58;51;50;41;38;35;28;25;18;15;12;10;7;7;4]; % 根据此变化表拟合求解相关系数 ft=fittype(A1*exp(-a*x)+B1*exp(-b*x)); options=fitoptions(Method,NonlinearLeastSquares); options.StartPoint=[0 -1000 0 0]; cfit=fit(t,c,ft,options); plot(cfit,t,c,o); A1=cfit.A1 B1=cfit.B1 a=cfit.a b=cfit.b 由此解得：(拟合曲线见右图) A1 = 110.55 B1 = -151.46 a = 0.17949 b = 2.8243 %---1---% %问题：某人中午12点喝了一瓶啤酒，下午6点检查合格，晚饭又喝一瓶，次日凌晨2点检查未通过，请对此情况做出解释。 t11_1=6; c11_1=(A1*exp(-a*t11_1)+B1*exp(-b*t11_1))/2 %下午6点酒精浓度 除以2是因为此人只喝了一瓶，浓度减半 t11_2=13.2; c11_2=(A1*exp(-a*t11_2)+B1*exp(-b*t11_2))/2 %中午喝酒在次日凌晨2点残留酒精浓度 t12=7.2; c12=(A1*exp(-a*t12)+B1*exp(-b*t12))/2 %晚饭喝酒在次日凌晨2点残留酒精浓度 c12A=c12+c11_2 %次日凌晨2点总残留酒精浓度 由此解得： c11_1 =18.829 %下午6点酒精浓度20，故检测合格 c11_2 =5.1712 c12 =15.181 c12A =20.352 %次日凌晨2点总残留酒精