数学模型饮酒驾车模型.docVIP

一、实验目的 理解数据拟合基本内容，掌握Matlab软件求解数据拟合的基本方法 二、实验用仪器设备、器材或软件环境 Matlab软件 三、求解问题 据报载，2003年全国道路交通事故死亡人数为10.4372万，其中因饮酒驾车国家质量监督检验检疫局2004年5月31日发布《车辆驾驶人员血液、呼气酒精含量阈值与检验》国家标准标准规定，车辆驾驶人员血液中的酒精含量大于或等于20毫克／百毫升，小于80毫克／百毫升为饮酒驾车于0毫克／百毫升，血液中的酒精含量大于或等于80毫克／百毫升为醉酒驾车大于或等于0毫克／百毫升。啤酒驾车啤酒驾车饮酒驾车血液中酒精含量3瓶啤酒或者半斤低度白酒驾车 参考数据 1. 人的体液占人的体重的65%至70%，其中血液只占体重的7%左右；而药物（包括酒精）在血液中的含量含量瓶啤酒血液中酒精含量毫克／百毫升酒精含量30 68 75 82 82 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 数学模型：进入吸收室中的酒精含量（单位：mg）。 ：吸收室中的酒精转移到中心室的速率系数。 ：中心室的酒精排出体外的速率系数。 ：长时间内喝酒，酒精进入中心室的速率。 ：酒精在时刻转移到中心室的速率。 ：时间（单位：时）。 ：时刻血液中的酒精浓度。 ：吸收室在时的酒精含量。 ：吸收室在时刻的酒精含量。 ：中心室在时刻的酒精含量。 ：长时间内喝酒中心室在时刻血液中的酒精含量。 ：人体中血液的体积。 ：中心室中血液的体积。 模型一（快速饮酒问题） （一）、模型建立与求解 快速饮酒与口服或肌肉注射药物的过程相似,这就相当于酒精进入中心室前先有一个将酒精吸收入血液的过程可以简化为有一个吸收室。如图一 于是满足 (1) 而药物进入中心室的速率为 (2) 将方程（1）的解代入（2）式得 (3) 又由假设得出微分方程 (4) 与血液浓度、房室容积之间显然有关系式 = （5） 将（5）式代入（4）式得 （6） 用MATLAB解此微分方程，可得 （7） 即 （8） 因为都是常数 故（8）式可以简化为 都为待定参数，其中 （二）、参数估计 将血样分成两部分，对于比较大时， 可近似有 两边取对数得 利用MATLAB将数据拟合 得出 则 对于其余数据进行以下处理 t=[0.25 0.5 0.75 1 1.5 2 2.5 3 3.5 4 4.5 5 6 7] t = Columns 1 through 13 0.2500 0.5000 0.7500 1.0000 1.5000 2.0000 2.5000 3.0000 3.5000 4.0000 4.5000 5.0000 6.0000 Column 14 7.0000 c=[30 68 75 82 82 77 68 68 58 51 50 41 38 35 ] c = 30 68 75 82 82 77 68 68 58 51 50 41 38 35 polyfit(t,c4,1) ans = -0.4775 4.8708 exp(4.8708) ans = 130.4252 同理得出 从而得出快速喝入两瓶啤酒， 血液中酒精含量与时间的关系式： 经过小时后再喝入两瓶啤酒，可以进一步改进模型得血液中酒精含量与时间的关系式： 模型二（慢速饮酒问题） （一）、模型建立与求解 慢速饮酒与恒速静脉滴注过程相似，我们假设血液中酒精是比较长时间（两小时内）恒速进入中心室的，其速率为，和初始条件为： 如图三，可得微分方程 （9） 又因为 求导得 （10） 把（10）代入（9）可得 （11） 即 （12） 可以利用MATLAB解微分方程得 （13） 又因为 （14） 所以 （15） 把（