实验目
运药物注射模型熟练MATLAB曲线拟合方法解释饮酒驾车实际问题
实验原理
酒精需进入肠道吸收胃吸收速率非常快题应采快速静脉注射模型 酒精存血液中例应计算吸收室血药浓度c1(t)A1eαt+B1eβt A1αB1β间关联提高精确度重新解微分方程题目应模型拟合计算
实验容
国家质量监督检查检疫局2004年5月31日发布新车辆驾驶员血液呼气酒精含量阀值检查国家新标准新标准规定车辆驾驶员血液中酒精含量等20毫克百毫升80毫克百毫升饮酒驾车(原标准100毫克百毫升)血液中酒精含量等80毫克百毫升醉驾车(原标准100毫克百毫升)
某中午12点喝瓶啤酒午6点检查时符合新驾车标准紧接着吃晚饭时候喝瓶啤酒保险起见呆凌晨2点驾车回家次遭遇检查时定饮酒驾车懊恼困惑什喝样酒两次检查结果会样呢?
(1) 某中午12点喝瓶啤酒午6点检查合格晚饭喝瓶次日凌晨2点检查未通请情况做出解释
(2) 短时间喝啤酒3瓶长时间驾车?
(3) 样估计血液中酒精含量什时候高?
(4) 果天天喝酒否开车?
解答:
建立常微分方程模型假设喝进酒精胃吸收转移速率胃里酒精含量成正血液代谢酒精速度浓度成正
图示:
X
胃
C
血液
代谢物
K1
K2
设胃里初始含量X0血液中初始含量C00
解
题目数C00时
MATLAB命令:
T[025 05 075 1 15 2 25 3 35 4 45 5 6 7 8 9 10 11 12 13 14 15 16]
C[30 68 75 82 84 77 70 68 58 51 50 41 38 35 28 25 18 15 12 10 7 7 4]
cftool开曲线拟合工具箱X data选择TY data选择C拟合方式选择Custom Equation拟合参数图
拟合:a2273 b01822 c1034
K12273K201822X01034发现拟合较
第题
喝瓶啤酒时X0517时
符合驾车标准
紧接着喝瓶时X0约517C0188367凌晨二点8时时
发现没2072时时略20属酒驾题目情况晚喝酒快速喝导致误差
第二题
短时间喝三瓶啤酒时X01551时
MATLAB命令:
T00124
C168616*exp(01878*T)168616*exp(1971*T)
plot(TC’r’)
hold on
plot([0 24][20 20]’g’)
发现C20相交1112间
T110112
C168616*exp(01878*T)168616*exp(1971*T)
输出:
C
1 7 列
213665 209690 205789 201960 198202 194515 190896
8 11 列
187344 183859 180438 177081
114时驾车会违反规定
题目三
观察图知高点12间
T1012
C168616*exp(01878*T)168616*exp(1971*T)
输出:
C
1 6 列
1162545 1178569 1187560 1190852 1189541 1184526
7 11 列
1176543 1166193 1153972 1140283 1125457
喝酒约13时血液中酒精含量高
第四题
简化问题假设天喝次次快速喝n瓶
求: 类推
考虑48时影响需数日符合认符合里取十天
考虑常酒量排呕吐等正常代谢酒精方式脏器超负荷工作带误差考虑等10瓶啤酒
MATLAB:
建立函数
function C24CalcC24(nC0i)
T00524
C(C0+562053*n)*exp(01878*T)562053*n*exp(1971*T)
C24(C0+562053*n)*exp(01878*24)562053*n*exp(1971*24)
plot((T+i*24)C'Color'[((1)^n+1)(2*n) ((1)^(n+1)+1)(2*n) ((1)^n+1)(2*n)])
hold on
end
输入命令:
hold off
clear
clc
C00
for n 110
for i010
tempCalcC24(nC0i)
C0temp
end
end
plot([0 264][20 20]’r’)
:
发现日喝十瓶啤酒情况18时血液中酒精浓度降低20mg100ml早八点喝十瓶啤酒凌晨两点驾车违反新交法
隔12时快速饮酒n瓶样考虑十天
MATLAB
建立函数:
function C12CalcC12(nC0i)
T00512
C(C0+562053*n)*exp(01878*T)562053*n*exp(1971*T)
C12(C0+562053*n)*exp(01878*12)562053*n*exp(1971*12)
plot((T+i*12)C'Color'[((1)^n+1)(2*n) ((1)^(n+1)+1)(2*n) ((1)^n+1)(2*n)])
hold on
end
输入命令:
hold off
clear
clc
C00
for n 15
for i020
tempCalcC12(nC0i)
C0temp
end
end
plot([0 264][20 20]'r')
:
发现喝完2瓶啤酒98时降低20mg100ml喝完三瓶仅喝完1194时降低20mg100ml仅十二时36分钟符合求忽略计
12时饮酒2瓶驾车违反新交法
综述12时2瓶啤酒24时喝10瓶啤酒违反新交法
模型评价改进
1解应模型综合运MATLAB软件准确求解运MATLAB进行数拟合时较理想化曲线
2模型引医药动力学二室模型进行计算性较高
3问题出发分析应该考虑种情况建立般数学模型进行实例验证证明建立数学模型较解决实际问题
模型缺点
1文模型参数仅组数拟合求解出定偏差
2模型计算简便结果更理想化忽略次素
3计算题目述相拟合较考虑身操作误差
— END —
文档香网(httpswwwxiangdangnet)户传
《香当网》用户分享的内容,不代表《香当网》观点或立场,请自行判断内容的真实性和可靠性!
该内容是文档的文本内容,更好的格式请下载文档