药物代谢问题数学模型作业.docVIP

药物代谢问题 一 摘要 药物经口服而进入体内，因代谢而逐步消除药物。此建模根据t=0时含有X剂量的药物后血内药物量y（纳克）与时间t的关系，以及问题中的数据可推算出口服一定药物后的最高血药浓度Cmax。随代谢的进行，体内的血药浓度会逐渐降低，当血药浓度降低到一定值时，为保证药效需要再次服药，模型中给出第二次服药与第一次服药的时间间隔以及第二次服药的药剂量，为病人的再次服药提供了好的参考。模型中还给出了24小时之内的血药浓度曲线，可以直观地看出血药浓度的变化情况。 二问题重述 设表示时刻体内药量，药物经口服吸收而进入血内，因代谢而逐步消除药物(排泄). 已知在t = 0时口服含X(克)剂量的药物后血内药物剂量 (纳克) (1纳克=克)与时间t (小时)的关系为 ， 其中为未知的吸收速度常数，F为未知的吸收比例常数，K为未知的消除速度常数. 问题： 1. 问一体重60千克的人第一次服药X=X1=0.1克剂量后的最高血药浓度Cmax(纳克/毫升); 2. 为保证药效, 在血药浓度降低到437.15纳克/毫升时应再次口服药物, 其剂量应使最高浓度等于Cmax(纳克/毫升). 求第二次口服的时间与第一次口服的时间的间隔T2(小时)和剂量X2(克). 3. 画出符合2的二次服药情况下在24小时之内的血药浓度曲线(将所要求的三个量Cmax, T2，X2的数值的最后结果皆舍入到4位数字, 且要保证4位数字都是有效数字). 三参考数据与参数说明 现有一体重60千克的人在t =T1= 0时, 第一次口服某药(含剂量X=0.1(克))，经3次检测得到数据如下: t =3(小时)时血药浓度为763.9(纳克/毫升)(血药浓度 (纳克/毫升), V表示未知血液容积(毫升). t = 18(小时)时血药浓度为76.39纳克/毫升，t = 20(小时)时血药浓度为53.4(纳克/毫升). ：表示时刻体内药量 ：未知的吸收速度常数 F ：未知的吸收比例常数 K ：未知的消除速度常数. C(t): 表示时刻体内血药浓度 单位（纳克/毫升） V：表示身体内的血液的容积 X :表示服药的剂量 单位（克） 四问题分析 本题要解决人体服药后，血液中药物的含量是多少及下一次服药的时间间隔与服药剂量的问题。属于微分方程类型。因为药物是一次服入体内的，可根据药物动力学知识中的“二室模型”。题目已经给出了微分方程，所以由微分方程得到数学模型。最后可以利用MATLAB软件对所提供的数据进行拟合，得到血液中药物含量C(t)与时间的一个函数关系式。然后求出最高血药浓度Cmax第二次口服的时间与第一次口服的时间的间隔T2(小时)和剂量X2(克)，最后做出符合2的二次服药情况下在24小时之内的血药浓度曲线。 五建模假设 1设相同体重的人的药物代谢的情况相同. 2假设人体血液体积一定。 3假设人体对药物的吸收、消化、排泄功能正常 六数学建模与求解 （一）模型建立与求解 ---------------------------------------------（1） --------------------------------------------------（2） 用matlab求解微分方程（1）可得 y(t)=C2/exp(K*t) + (Ka*F*X)/(exp(Ka*t)*(K - Ka)) 由t=0,y(t)=0 可得 C2= (Ka*F*X)/ (Ka-K) -------------------------------------------------------------------（3） 由（2）（3）可得 C(t)= (Ka*F*X)（exp(-K*t)- exp(-Ka*t)）/ (Ka-K)/V --------------------------------(4) 令A= （4）式可以化为 C(t)=Aexp(-Kt)-Aexp(-Kat)--------------------------------------------------------------（5） 不妨设KKa, 第一、当t很大时，（5）可看成是 C(t)= Aexp(-kt) ---------------------------------(6) 两边取对数,使得(6)变为线性,得 ln C(t)=lnA-Kt----------------------------------------------(7) 通过运用MATLAB进行线性回归分析, 得 在命令窗口键入： x=[3,18,20]; x=[ones(3,1) x]; y=[6.638,4.336,3.978