第四章作业。
第二题：针对严重的交通情况，国家质量监督检验检疫局发布的国家标准，车辆驾驶人员血液中的酒精含量大于或等于20mg/100ml,小于80mg/100ml为饮酒驾车，血液中的酒精含量大于或等于80mg/100ml的为醉酒驾车。
下面分别考虑大李在很短时间内和较长时间内（如2个小时）喝了三瓶啤酒，多长时间内驾车就会违反新的国家标准。
1、 问题假设。
大李在短时间内喝下三瓶啤酒后，酒精先从吸收室（肠胃）吸收进中心室（血液和体液），然后从中心室向体外排除，忽略喝酒的时间，根据生理学知识，假设。
1） 吸收室在初始时刻t=0时，酒精量立即为；在任意时刻，酒精从吸收室吸收进中心室的速率（吸收室在单位时间内酒精含量的减少量）与吸收室的酒精含量成正比，比例系数为；
2） 中心室的容积v保持不变；在初始时刻t=0时，中心室的酒精含量为0；在任意时刻，酒精从中心室向体外排除的速率（中心室在单位时间内酒精含量的减少量）与中心室的酒精含量成正比，比例系数为；
3） 在大李适度饮酒没有酒精中毒的前提下，假设和都是常量，与饮酒量无关。
2、 符号说明。
酒精量是指纯酒精的质量，单位是毫克；
酒精含量是指纯酒精的浓度，单位是毫克/百毫升；
时刻（小时）；
在时刻吸收室（肠胃）内的酒精量（毫克）；
两瓶酒的酒精量（毫克）；
在时刻吸收室（血液和体液）的酒精含量（毫克/百毫升）；
在时刻中心室（血液和体液）的酒精含量（毫克/百毫升）；
中心室的容积（百毫升）；
酒精从吸收室吸收进中心室的速率系数（假设其为常数2.0079）；
酒精从中心室向体外排除的速率系数（假设其为常数0.1855）；
在短时间喝下三瓶酒的假设下是指短时间喝下的三瓶酒的酒精总量除以中心室体积，即；而在较长时间内（2小时内）喝下三瓶酒的假设下就特指。
3、 模型建立和求解。
1） 酒是在很短时间内喝的：
记喝酒时刻为（小时），设，可用来计算血液中的酒精含量，此时为假设中所示的常数，而。
以及血液中酒精含量最高的时刻。
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 on
plot([0,20],[20,20],'k',[0,20],[80,80],'k')
hold off
xlabel('时刻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 =
t2 =t3 =
c3 =
所绘图形如下：
结果分析：所以，当时，，属饮酒驾车。当时，属醉酒驾驶；当时，血液中的酒精含量最高为122.25毫克/百毫升。
2） 酒是在2小时内喝的：
可假设三瓶啤酒是在2小时内匀速喝的。 同样记喝酒时刻为（小时），设，则吸收室的酒精量满足分段的初值问题。
解得。于是中心室内的酒精含量满足分段的初值问题。
解得。其中， ，
因为，以及，所以， ，
以及血液中酒精含量最高的时刻。
matlab程序如下：
k1=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);t2(2)=fzero(g1,12),t3,c3]=fminbnd(h1,0,20)
fplot(c1,[0,20],'k')
hold on
plot([0,20],[20,20],'k',[0,20],[80,80],'k')
hold off
xlabel('时刻t（小时），从开始喝酒算起')
ylabel('血液中的酒精含量（mg/100ml）')
title('短时间喝下三瓶酒时，血液中酒精含量随时间的变化过程')
gtext('(0.62321,20)')
gtext('(12.62,20)')
gtext('(1.6366,80)')
gtext('(5.1412,80)')
gtext('(2.6328,115.74)')
运行结果如下：t1 =
t2 =t3 =
c3 =
所绘图形如下：
结果分析：所以，当时，，属饮酒驾车。当时，属醉酒驾驶；当时，血液中的酒精含量最高，为115.74毫克/百毫升。
下面用图形比较两种不同假设下血液中酒精含量的变化过程。
matlab程序如下：
k1=2.0079;k2=0.1855;k3=155.79;
k4=42.743;k5=462.66;k6=419.92;k9=2328.3;k10=207.82;
c=@(t)(k1.*k3)./k1-k2).*exp(-k2.*t)-exp(-k1.*t));
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);
plot(0:0.01:
20,c(0:0.01:
20),'k',0:0.01:
20,c1(0:0.01:
20),'k',2,c1(2),'k')
xlabel('时刻t（小时），从开始喝酒算起')
ylabel('血液中的酒精含量（mg/100ml）')
title('短时间喝下三瓶酒时，血液中酒精含量随时间的变化过程')
legend('很短时间内喝三瓶啤酒','两小时内匀速喝下三瓶啤酒','函数的分段点')
所绘图形如下：
第四题：研究将鹿群放入草场后，草和鹿两个种群的相互作用，草的生长服从logistic规律，年固有增长率0.8，最大密度为3000个密度单位，在草最茂盛时，每只鹿每年吃掉1.
6个密度单位的草，若没有草，鹿群的年死亡率高达0.9,而在草最茂盛的时候草对鹿的死亡的补偿率为1.5。
1、建立差分方程组模型，比较将100只鹿放入密度为1000和密度为3000的两种草场的情况下，草和鹿两个种群的数量演变过程。
1）符号说明：
第k年草场的密度单位。
第k年草场上鹿的数量。
草场上草的年固有增长率。
由于捕食导致的草的密度单位减少的速度大小。
如果没有草，鹿群的年死亡率。
草对鹿群死亡的补偿率。
草的最大密度单位。
2）模型的建立与求解：
基于以上假设，由于草的生长服从logistic模型，建立差分方程组模型如下所示：
令，，与上述方程组联立得到平衡点为、、
以下用matlab实现差分方程组模型。
matlab程序如下：
n=20;r=0.8;a=1.6;b=1.5;d=0.9;
n=3000;t=1:n+1;
x1(1)=1000;y1(1)=100;
for k=1:n
x1(k+1)=x1(k)+r*x1(k)*(1-x1(k)/n)-a*x1(k)*y1(k)/n
y1(k+1)=y1(k)-d*y1(k)+b*x1(k)*y1(k)/n
endsubplot(2,2,1),plot(t,x1,'k^',t,y1,'kv'),axis([-1,21,0,3000]),xlabel('第k年'),ylabel('数量'),gtext('草密度'),gtext('鹿数量'),title('草和鹿随时间的演变');
subplot(2,2,2),plot(x1,y1,'ko'),axis([1000,3000,0,1000]),xlabel('草密度'),ylabel('鹿数量'),title('相平面图');
x2(1)=3000;y2(1)=100;
for k=1:n
x2(k+1)=x2(k)+r*x2(k)*(1-x2(k)/n)-a*x2(k)*y2(k)/n
y2(k+1)=y2(k)-d*y2(k)+b*x2(k)*y2(k)/n
endsubplot(2,2,3),plot(t,x2,'k^',t,y2,'kv'),axis([-1,21,0,3000]),xlabel('第k年'),ylabel('数量'),gtext('草密度'),gtext('鹿数量'),title('草和鹿随时间的演变');
subplot(2,2,4),plot(x2,y2,'ko'),axis([1000,3000,0,1000]),xlabel('草密度'),ylabel('鹿数量'),title('相平面图');
数学建模章绍辉版第四章作业
第四章作业。第二题 针对严重的交通情况，国家质量监督检验检疫局发布的国家标准，车辆驾驶人员血液中的酒精含量大于或等于20mg 100ml,小于80mg 100ml为饮酒驾车，血液中的酒精含量大于或等于80mg 100ml的为醉酒驾车。下面分别考虑大李在很短时间内和较长时间内 如2个小时 喝了三瓶啤酒...
章绍辉数学建模第四章
第四章。符号说明 酒精量是指纯酒精的质量，单位是毫克 酒精含量是指纯酒精的浓度，单位是毫克 百毫升 t 时刻 小时 x1 t 在时刻t吸收室 肠胃 内的酒精量 毫克 d0 在段时间内喝下2瓶啤酒后吸收室内的酒精量 毫克 c1 t 在时刻t吸收室 肠胃 的酒精含量 毫克 百毫升 c2 t 在时刻t中心...
第四章数学建模作业
严婉琳20110008015 邓南广 20110008005 11创新班 3,分析 设为t时刻的酵母菌数，且是连续函数并且可导或分段可导，设时刻酵母菌的增长率满足。是酵母菌固有增长率，称为酵母菌最大容量。得到阻滞增长模型 用分离变量法 1 式满足条件解为。由书本例题可知，散点图呈s型，故模型可简化为...