饮酒驾车模型及matlab实现要点详解.pptVIP

7.5 饮酒驾车 大李喝下啤酒后，酒精先从肠胃吸收进入血液和体液中，然后从血液和体液向体外排出。可以建立二室模型，将肠胃看成吸收室，将血液与体液看成中心室(见下图)。 酒精量是指纯酒精的质量，单位为毫克(mg)；酒精含量是指纯酒精的浓度，单位是毫克/百毫升(mg/100ml)； t:时刻(h)； x1(t):在时刻t吸收室(肠胃)内的酒精量(mg)； k1:酒精从吸收室进入中心室的速率系数； g0:在短时间内喝下1瓶啤酒后吸收室内的酒精量(mg)； y1(t):在时刻t中心室(血液和体液)的酒量(mg); K2:酒精从中心室向体外排出的速率系数; V:中心室的容积(100ml). 数学实验 微分方程模型及非线性拟合实验 问题 吸收室 x1(t) 中心室 c1(t),V k1 k2 吸收 排出 吸收和排出的过程都可以简化成一级反应来处理，加起来得到体液内酒精吸收和排出过程的数学模型。因为考虑到时短时间内喝酒，所以忽略喝酒的时间，可使初始条件得以简化。 根据上面问题要求，可归结为如下问题： 1、建立数学模型，并解释大李在中午12点喝1瓶啤酒后，在下午6点检查时体内血液中的酒精含量小于20mg/10ml，符合“驾车标准”。 2、建立数学模型，并解释大李在晚饭时再喝1瓶啤酒后，在凌晨2点检查时体内血液中的酒精含量不小于20mg/100ml,不符合“饮酒驾车”的标准。 （4）考虑到大李在下午6点接受检查，之后由于离开检查地点以及停车等待等原因耽误了一定时间，因此假定大李在晚8点吃晚饭 (即大李从第一次接受检查到第二次喝酒之间相隔了2个小时) 大李在短时间内喝下2瓶啤酒后，酒精先从吸收室（肠胃）进入中心室（血液与体液），然后从中心室向体外排出。忽略喝酒时间，并假设： （1）吸收室在初始时刻t=0时，酒精量立即为2g0，酒精从吸收室进入中心室的速率（吸收室在单位时间内酒精量的减少量）与吸收室的酒精量成正比，比例系数为k1. (2)中心室的容积V保持不变；在初始时刻t=0时，中心室酒精量为0;在任意时刻，酒精从中心室向体外排出的速率（中心室的单位时间内酒精量的减少量）与中心室的酒精量成正比，比例系数为k2. (3)在大李（体重为70kg）适度饮酒没有酒精中毒的前提下，假设k1和k2都是常数，与酒精量无关。 根据假设（1），吸收室的酒精量x1(t)满足微分方程初值问题 根据假设（2），中心室的酒精量y1(t)满足微分方程初值问题： 根据（7.5.1）和（7.5.2）得到微分方程组初值问题： 解上述微分方程组初值问题，其matlab程序如下： [x1,y1]=dsolve(Dx1=-k1*x1,Dy1=k1*x1-k2*y1,x1(0)=N*g(0),y1(0)=0) [y,how]=simple([x1,y1]) y = [N*g(0)*exp(-k1*t), k1*N*g(0)*(exp(-k2*t)-exp(-k1*t))/(k1-k2)] 程序运行结果： 程序求解结果整理为： 即:解为 返回 式（7.5.3）可以写成 当前任务就是，确定k,k1,k2 用MATLAB的函数nlinfit(非线性最小二乘拟合),根据赛题所给数据拟合式（7.5.4）的参数k1,k2和k.此问题的MATLAB程序如下： f=@(k,x)k(3).*(exp(-k(2).*x)-exp(-k(1).*x)); x=[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]; y=[30 68 75 82 82 77 68 68 58 51 50 41 38 35 28 25 18 15 12 10 7 7 4]; k0=[2,1,80];%参数的初值 k=nlinfit(x,y,f,k0) plot(x,y,r*,0:0.01:18,f(k,0:0.01:18),k) xlabel(‘时间(h)) ylabel(‘酒精含量) title(‘血液中酒精含量的拟合图) axis([0 18 0 90]) legend(‘原始数据’,‘拟合曲线) 参数k1，k2和k的拟合结果为: k =2.0079 0.1855 114.432524 参数的初值设定思路： fc=@(x)k(3).*(exp(-k(2).*x)-exp(-k(1).*x)); figure(2) plot(x,y-fc(x),ro,[0 18],[0 0],k) axis([0,18,-10,10]) 结果如图： 从原始数据拟合图、拟合误差图观察发现拟合效果比较好。在拟合误差图中只有一个误差在-10附近，其他误差值都在（-6,6）之内，且分布比较均匀，这说明引入的假设和建立的模型比