资源描述：
《数学建模章绍辉版第四章作业.doc》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。
1、第四章作业第二题：针对严重的交通情况，国家质量监督检验检疫局发布的国家标准，车辆驾驶人员血液中的酒精含量大于或等于20mg/100ml,小于80mg/100ml为饮酒驾车，血液中的酒精含量大于或等于80mg/100ml的为醉酒驾车。下面分别考虑大李在很短时间内和较长时间内（如2个小时）喝了三瓶啤酒，多长时间内驾车就会违反新的国家标准。1、问题假设大李在短时间内喝下三瓶啤酒后，酒精先从吸收室（肠胃）吸收进中心室（血液和体液），然后从中心室向体外排除，忽略喝酒的时间，根据生理学知识，假设（1）吸收室在初始时刻t=0时，酒精量立即为;在任意时刻，酒精从吸收室吸收进
2、中心室的速率（吸收室在单位时间内酒精含量的减少量）与吸收室的酒精含量成正比，比例系数为；（2）中心室的容积V保持不变；在初始时刻t=0时，中心室的酒精含量为0；在任意时刻，酒精从中心室向体外排除的速率（中心室在单位时间内酒精含量的减少量）与中心室的酒精含量成正比，比例系数为；（3）在大李适度饮酒没有酒精中毒的前提下，假设和都是常量，与饮酒量无关。2、符号说明酒精量是指纯酒精的质量，单位是毫克；酒精含量是指纯酒精的浓度，单位是毫克/百毫升；时刻（小时）；在时刻吸收室（肠胃）内的酒精量（毫克）；两瓶酒的酒精量（毫克）；在时刻吸收室（血液和体液）的酒精含量（毫克/
3、百毫升）；在时刻中心室（血液和体液）的酒精含量（毫克/百毫升）；中心室的容积（百毫升）；酒精从吸收室吸收进中心室的速率系数（假设其为常数2.0079）；酒精从中心室向体外排除的速率系数（假设其为常数0.1855）；在短时间喝下三瓶酒的假设下是指短时间喝下的三瓶酒的酒精总量除以中心室体积，即；而在较长时间内（2小时内）喝下三瓶酒的假设下就特指.1、模型建立和求解（1）酒是在很短时间内喝的：记喝酒时刻为（小时），设，可用来计算血液中的酒精含量，此时为假设中所示的常数，而.下面用MATLAB程序画图展示血液中酒精含量随时间变化并且利用fzero函数和fminbnd
4、函数来得到饮酒驾车醉酒驾车对应的时间段，以及血液中酒精含量最高的时刻。MATLAB程序如下：k1=2.0079;k2=0.1855;k3=155.79;c=@(t)(k1.*k3)./(k1-k2).*(exp(-k2.*t)-exp(-k1.*t));f=@(t)c(t)-20;g=@(t)c(t)-80;h=@(t)-c(t);t1(1)=fzero(f,1);t1(2)=fzero(f,12),t2(1)=fzero(g,1);t2(2)=fzero(g,12)[t3,c3]=fminbnd(h,0,24)fplot(c,[0,20],'k')hold
5、onplot([0,20],[20,20],'k',[0,20],[80,80],'k')holdoffxlabel('时刻t（小时），从开始喝酒算起')ylabel('血液中的酒精含量（mg/100ml）')title('短时间喝下三瓶酒时，血液中酒精含量随时间的变化过程')gtext('(0.06891,20)')gtext('(11.589,20)')gtext('(0.38052,80)')gtext('(4.1125,80)')gtext('(1.307,122.25)')运行结果如下：t1=0.0689111.589t2=0.380524.112
6、5t3=1.307c3=-122.25所绘图形如下：结果分析：所以，当时，，属饮酒驾车。当时，属醉酒驾驶；当时，血液中的酒精含量最高为122.25毫克/百毫升。（1）酒是在2小时内喝的：可假设三瓶啤酒是在2小时内匀速喝的.同样记喝酒时刻为（小时），设，则吸收室的酒精量满足分段的初值问题解得于是中心室内的酒精含量满足分段的初值问题解得其中，，，，，，因为，以及，所以，，，，，，下面用MATLAB程序画图展示血液中酒精含量随时间变化并且利用fzero函数和fminbnd函数来得到饮酒驾车醉酒驾车对应的时间段，以及血液中酒精含量最高的时刻。MATLAB程序如下：k
7、1=2.0079;k2=0.1855;k3=155.79;k4=42.743;k5=462.66;k6=419.92;k9=2328.3;k10=207.82;c1=@(t)(k4.*exp(-k1.*t)-k5.*exp(-k2.*t)+k6).*(t>=0&t<=2)+...(k10.*exp(-k2.*t)-k9.*exp(-k1.*t)).*(t>2);f1=@(t)c1(t)-20;g1=@(t)c1(t)-80;h1=@(t)-c1(t);t1(1)=fzero(f1,1);t1(2)=fzero(f1,12),t2(1)=fzero(g1,1)
8、;t2(2)=fzero(g1,12),[t3,c3
