1、南 昌 航 空 大 学数 学 与 信 息 科 学 学 院实 验 报 告课程名称： 数学模型 实验名称： 醉酒驾驶的数学模型 实验类型： 验 证 性 综 合 性 设 计 性 实验室名称： 数学实验室 D208 班 级 学 号 ： 11071120 学 生 姓 名 ： 张 万 晴 任 课 教 师 ： 张 邻 成 绩 ： 实 验 日 期 ： 2013-11-13 至 2013-11-20 - 2 -目 录一、 问题的重述3二、 问题的假设4三、 符号说明5四、 模型的建立5五、 问题分析与模型建立7六、 结果分析与检验17七、 模型的评价与改进18- 3 -一 、 问 题 重 述 据报载，2003
2、年全国道路交通事故死亡人数为 10.4372万，其中因饮酒驾车造成的占有相当的比例。 针对这种严重的道路交通情况，国家质量监督检验检疫局 2004年 5月 31日发布了新的车辆驾驶人员血液、呼气酒精含量阈值与检验国家标准。 驾驶人员血液中的酒精含量大于或等于 20毫克百毫升，小于80毫克百毫升为饮酒驾车 血液中的酒精含量大于或等于 80毫克百毫升为醉酒驾车 大李在中午 12点喝了一瓶啤酒，下午 6点检查时符合新的驾车标准 在吃晚饭时他又喝了一瓶啤酒，为了保险起见他呆到凌晨 2点才驾车回家，又一次遭遇检查时却被定为饮酒驾车为什么喝同样多的酒，两次检查结果会不一样呢？ 问题 1 对大李碰到的情况做
3、出解释 问题 2 喝了 3瓶啤酒或者半斤低度白酒后多长时间内驾车就会违反上述标准，在下列情况下回答该问题： 1. 酒是在很短时间内喝的；2. 酒是在较长一段时间（比如 2小时）内喝的。参 考 数 据1. 人的体液占人的体重的 65%至 70%，其中血液只占体重的 7%左右；而药物（包括酒精）在血液中的含量与在体液中的含量大体是- 4 -一样的。2. 体重约 70kg 的某人在短时间内喝下 2 瓶啤酒后，隔一定时间测量他的血液中酒精含量（毫克百毫升） ，得到数据如下：时间(小时) 0.25 0.5 0.75 1 1.5 2 2.5 3 3.5 4 4.5 5酒精含量 30 68 75 82 82
4、 77 68 68 58 51 50 41时间(小时) 6 7 8 9 10 11 12 13 14 15 16酒精含量 38 35 28 25 18 15 12 10 7 7 4二 、 问 题 假 设1. 不考虑酒精进入体内随呼吸或汗液排出的量，及肠道细菌产生的酒精，只考虑饮入的酒全进入肠胃，再由肝脏等分解的过程。2. 假设体液中的酒精消耗(向外排出、分解或吸收)的速度，与体液中的酒精浓度(或含景)成正比。3. 不考虑酒精进入体内随呼吸或汗液排出的量，及肠道细菌产生的酒精，只考虑饮入的酒全进入肠胃，再由肝脏等分解的过程4. 设人体血液和体液中酒精浓度相等,酒精进入血液后瞬间混合均匀。5. 假
5、设整体过程中人没有摄入任何影响代谢的药类物质和做剧烈- 5 -三、 符号说明t时刻肠胃中的酒精含量：p(t):t时刻血液中的酒精浓度；X(t):t时刻血液中的酒精含量；r1(t)：饮入酒精的速率；r2(t):肠胃内酒精进入血液的速率；r3(t):肠胃内酒精进入（除体液外）其它地方的速率；r4(t)：血液中酒精被分解的速率。：所饮酒中含的酒精量；0pV：体液的体积；X(t)：四、 模型的建立一、符号说明及模型假设1.符号说明-酒精在人体中的消除速率常数 1k-t时刻血液中的酒精浓度tcF-酒精在人体中的吸收度V-人体的血液体积V 酒 -喝酒的体积t-饮用酒的时间- 6 -t时刻血液中的酒精量)(
6、tx-t时刻人体吸收的酒精量1M-人的体重-人的体液占人的体重的百分含量-酒中的酒精含量-饮酒持续时间-人体饮入酒精总量0x-人的血液占人体重的百分含量-酒精在人体中的吸收速率常数 11k2.基本假设1. 酒精按一级吸收过程进入体内；2. 正常情况下，酒精在各人体中的吸收和消除速率基本相同3. 酒精进入人体后，不考虑其他因素对酒精的分解作用；4. 每瓶啤酒的酒精含量、体积基本相同；5. 如果在很短时间内饮酒，认为是一次性饮入，中间的时间差不计；6. 酒精在血液中的含量与在体液中的含量大至相同；假设模型 一个人的血液中酒精含量取决于他的饮酒量、体内原有的酒精含量以及喝酒方式等。由科普知识知道，酒
7、精是经胃肠(主要是肝- 7 -脏)的吸收与分解进入体液的。根据假设可建立以下基本模型：)()()(4324321tyktrxttrdttytx从生物学可知，酒类进入人体后，胃及血液中酒精随注入酒精速率、浓度、时间等不同产生不同的代谢速率。下面根据饮酒速率及方式的不同建立三种实用模型。短时间内快速饮酒模型1在基本模型的基础上，由短时间快速饮酒的特点可得出初始值：， , ，将其代入基本模型可得模型0)(1tr0)(px)(y)()(423txktdtyt五、问题分析与模型建立1.快速饮酒模型同药物一样，酒精进入机体后，作用于机体而影响某些器官组织的功能；另一方面酒精在机体的影响下，可以发生一系列的
8、运动和体内过程：自用药部位被吸收进入血液循环；然后分布于各器官组织、组织间隙或细胞内；有部分酒精则在血浆、组织中与蛋白质结- 8 -合；或在各组织（主要是肝脏）发生化学反应而被代谢；最后，酒精可通过各种途径离开机体（排泄） ；即吸收、分布、代谢和排泄过程。它们可归纳为两大方面：一是酒精在体内位置的变化，即酒精的转运，如吸收、分布、排泄；二是酒精的化学结构的改变，即酒精的转化亦即狭义的代谢。由于转运和转化以致形成酒精在体内的量或浓度（血浆内、组织内）的变化，而且这一变化可随时间推移而发生动态变化。又因为酒精有促进血液循环的作用 2。而药物动力学模型中的一室模型 3是指给药后，药物一经进入血液循环
9、，即均匀分布至全身,故快速饮酒情况可通过建立一室模型求解。虽然酒精在体内的分布状况复杂，但酒精的吸收、分解等则都在系统内部进行，酒精进入人体后，经一段时间进入血液，进入血液后，当在血液中达最高浓度时，随后开始消除 3，把酒精在体内的代谢过程看为进与出的过程，这样便会使问题得到简化。用和 分别表示酒精输入速率和输出速率。由于单位时间内indtxout血液中酒精的改变即变化率 就等于输入与输出速率之差，所以其dtx动力学模型为：= - dtxinoutdx（1）又因为酒精在血液中的消除速率与当时血液内的药量成正比，所以=kt，代入（1）式得： = -kt outdx dtxin（2）- 9 -则由
10、（2）式可知X(t)的变化规律由饮酒速率而定。而酒精在人体内的代谢可简单的由图一表示：1k(图一)则 t时刻吸收室的药量为 X1(t)，又药物是按一级吸收过程进入体内的，对于吸收室有：=-k1X1 (3)dtx对于房室， = ，于是（2）式变为： （4）int1k kxdtx1（3） 、 （4）两式构成一阶线性方程组，当 t=0时， ，X(0)01)(F=0,解（2）式得：，将其代入（4）式得一阶线性非齐次方程：tkeFx101tkdt1解之得： tktekFxt110)(从而，人体内酒精含量为： tktekVxtC1)()10在这种情况下，酒精含量最大值出现的时间：使 时 t的值。0dtc一
11、般情况下，又因为酒精在血液中的含量与在体液中的含量大至相吸收室X1(t) X(t)VFX0 k- 10 -同。则有： (F为常数且MFx0F0F1)X0=V 酒 则人体内酒精含量与时间函数关系为:（一般情况）)()()11tkteKVktC酒1.1.模型的求解根据图一中酒精含量实测数据拟合,显然它无法化为线性最小二乘,我们直接作非线性最小二乘拟合 4。用MATLAB优化工具箱的Leastsq5计算,拟合参数 ),(01VFxk程序见:JM2004C1.m。拟合得： =2.0079mgml-1h-11,k=0.1855mgml-1h-1, =11.2423（毫克/毫升） ，又由x0，得到问题中隐
12、含的一瓶啤酒的酒精量约为：FVx243.027543.635毫克。1.2问题一的解答虽然大李喝等量的酒，并且相隔的时间也相同的情况下,两次检查的结果不一样是因为第一次喝下去的酒，在 6小时内并没有完全分解，还残留有相当一部分在血液中，并且这一部分在较长时间内不能完全分解。由图（二）喝一瓶啤酒的酒精含量随时间变化的- 11 -函数图像可知，6 小时后，第一次喝的酒的酒精含量约 19.5毫克/百毫升。也就是说此时血液中已有一定的酒精量，这样虽然第二次喝的是同样多的酒，由于第一次残留部分的存在，相当于涉入的酒精量已增大了，使其同样再过 6小时，酒精含将会大于 20毫克/百毫升。这样大李碰到的情况也就
13、很自然的解释了（图二）现通过实际计算证明：设 A为第一次喝酒在六个小时后，残留在血液中的酒精量。则第二次喝酒时，酒精含量 C（t）与时间 t的函数关系为：)()()110tktekVAFxt 代值用 MATLAB计算（程序见 JM2004C5.LOG） ，此时再过 6小时酒精含量为：27.4 毫克/百毫升。这就大于了第一次的值。- 12 -1.3问题二第一问的求解当一个人在很短时间内喝 3瓶啤酒时，相当于在瞬间使其吸收室的酒精浓度达最大值，运用模型并用计算机拟合得其函数图象如(图三)（程序见 JM20004C2.LOG） ，(图三)由图可得： t:0.0650.24小时内饮酒驾车；t:0.34
14、4.5 小时内醉酒驾车；t:4.512 小时内饮酒驾车。2.匀速饮酒二室模型）2.1模型的建立及求解由于在两小时内慢速喝酒有多种方式，为了便于计算我们考虑为匀速饮入等量酒精，即在时间 内匀速吸收酒精。若体内酒精含量不超过一级消除动力学范围，假设人的酒精含量未达到平衡状态，随着人体吸收次数增多，血液中酒精浓度逐渐升高，当在 时间内- 13 -饮完后，由于 时刻前一小段时间内饮入的酒精在短时刻内没有被吸收，故在 时刻后一段时间内，酒精浓度将继续升高，某一时刻将达到最大值。这一时刻后，酒精含量就会逐渐减小。考虑到一室模型在求慢速喝酒情况下的局限性，我们采用建立二室模型 6。对二室模型我们将建立两个关
15、于酒精浓度的动力系统模型来描述其动态特性。2.2假设：1、 机体分为中心室和周边室，两个室的容积在过程中保持不变。2、 药物从一室向另一室的转移速率，及向体外的排除速率，与该室的酒精浓度成正比。3、 只在中心室一体外有酒精交换，即酒精从体外进入中心室，最后又从中心室排出体外，与转移和排除的数量相比，酒精的吸收可以忽略。2.3模型建立二室模型的示意图如(图四)所 示：1kK12吸收室X1(t) X(t)VFX0 k- 14 -(图一)两个房室中酒精量 满足的微分方程。 的变化率由一室)(,21tx)(1tx向二室的转移 ，一室向体外排除 ，二室向一室的转移2k13k及酒精 组成；21xk)(0t
16、f)(2tx的变化率由一室向二室的转移 及二室向一室的转移 组成，12xk 21xk于是有：)(0211321 tfxkxkdt （1）2112tx与血液中酒精含量 、房室容积 显然有关系式)(i )(tci iV（2）,1),.( itxii将（2）式代入（1）式可得：（3）21122 102131 )()(ckVdtcVtf喝酒相当于在酒精进入中心室之前先有一个将酒精吸收入血液的过程，可以简化为有一个吸收室，如下图， 为吸收室的酒精，酒)(0tx精由吸收室进入中心室的转移速率系数为 ，于是 满足：1k)(0tx（4）001)(Dxkdt当 时， （3）可以化为：0)(,2110 cVctf
17、- 15 -（5）ttBeAtc)(1（ 图四） 酒精经 吸收室进入中心室是饮入的酒精量，而酒精进入中心室的速率为：0D)()(010txktf（6）将方程（5）式代入（6）式得：tkeDtf010)(（7）在这种情况下方程（3）的解的一般形式为：tktatEeBAetc01)(1 （8）此时，这种情况下，酒精含量最大值出现的时间为：使 时01dtc的时间。其中 132132k（9）01kDE3.2.1.3参数估计不妨设 ，于是当 t充分大时 可近似为： 对)(1tctAetc)(1吸收室 )(0tx中心室01xkf- 16 -于适当大的 和相应的 ，用最小二乘法估计出 的值。然后it)(1i
18、tc A和计算 （10）tAec)(1再利用（5）式得： 对于较小的 和（10）式算出的tBtcln)(l1 it，仍用最小二乘法可得到 。)(1itc 和由上可得参数值： A,24 模型的求解将题中所给参考数据代入，运用 MATLAB工具，对前参数估计进行代值计算得:程序见 JM2004C6-C8A=69.4908;a=0.3234;B=-48.0251;b= 2.77;将所求得的参数代入（8）当中可得，并用 MATABT画出匀速喝三瓶啤酒的酒精含量-时间曲线图如图（五）程序见 JM2004C9.M- 17 -（图五）由上函数图象可知：当 t为 24.5小时内为醉酒驾车；当 t为4.5-12
19、小时为饮酒驾车。六、结果分析与检验由于酒精不同于一般的药物，人饮用之后，酒精能迅速的进入人体，所以在快速饮酒情况下用一室模型来研究酒精的作用过程及代谢过程，与几室模型来研究酒精的作用过程及代谢过程并无太大的区别， （图七）酒精含量 C（t）的拟合曲线与实测数据的比较,可以看出,我们的结果是较为准确的，并且拟合的曲线是一条较为理想的曲线。- 18 -（图七）对于匀速饮酒情况，由于匀速饮酒要考虑饮酒过程中的代谢情况，与快速饮酒有一定差异，所以采用二室模型来建立匀速饮酒模型，达到减小误差的目的。七、模型的评价与改进1.模型的优点1.本模型从三种情况分别建立模型，模型稳定性高，适用性强。2.本模型计算
20、步骤清晰，引用了医药动力学的二室模型进行计算，可靠性较高- 19 -3.综合运用 MATLAB和 LINGO两个软件，准确求解，在运用MATLAB进行数据拟合时，得到了较理想化的曲线。在表示喝三瓶啤酒的人什么时候是饮酒驾车，什么时候是醉酒驾车时，运用 MATLAB准确的做出了函数据图像，使结果一目了然。2.模型的缺点1.模型为使计算简便，使所得的结果更理想化，忽略了一些次要的因素。如：酒进入身体后随着血液流动，人体对酒精的吸收率。对于这些问题，由于时间关系本模型还未能更好的研究，有待以后的改进和完善。2.本文的模型参数仅是依靠题中给出的一组数据拟合求解得出，可能有一点偏差。3模型的改进由于人体内部的复杂性，及各器官对酒精转化的多样性，用一室或二室都较为初级，三室或多室的情况模型更准确，但考虑起来会很复杂，又由于短期收集资料的局限，实行起来较为困难，可留着时间充裕时考虑。
