人体二室模型.docVIP

人体血药浓度模型 摘 要 本文讨论了药动学中的血药浓度随时间变化及相关药动学参数的问题。 针对问题一，在已知药动学二室模型的条件下，运用微分的数学思想，建立了单次给药时血药浓度关于时间的微分方程组。按照快速静脉注射、恒速静脉滴注、口服或肌肉注射3种给药方式给药速率的不同，带入不同的初值，分别得到3种给药方式下血药浓度与时间的微分方程组，通过线性微分方程组求解方法求出3种给药方式下血药浓度的解析表达式。结合具体实例并利用残数法，解出相关药动学参数，运用matlab软件画出药时曲线图。 针对问题二，采用问题一中单次给药方式下的血药浓度随时间变化的解析式，并且运用迭代的方法直接计算出多次给药方式下，血药浓度随时间变化的解析式，利用matlab软件画出药时曲线图。 对非线性模型进行初步讨论，提出血药浓度模型未来的发展趋势。 关键字 二室模型 血药浓度 药动学参数 残数法 一、问题重述 药动学通常用房室模拟人体，只要体内某些部位接受或消除药物的速率相似，即可归入一个房室。多数情况下二室模型能够准确地反映药物的体内过程特征。现解决以下问题： 问题一：根据已知的二室模型，讨论在快速静脉注射、恒速静脉注射、口服或肌肉注射3种给药方式下血药浓度的解析表达式，并画出浓度曲线图。 问题二：根据已知的二室模型，按照固定的时间间隔、每次给予固定的剂量的多次重复给药方式。为了维持药品的疗效和保证机体的安全，要求血药浓度控制在确定的合理范围内。讨论了在快速静脉注射、恒速静脉注射、口服或肌肉注射3种多次重复给药方式血药浓度的解析表达式，并画出浓度曲线图。 二、问题分析 对问题一，已知二室模型，根据快速静脉注射、恒速静脉滴注、口服或肌肉注射3种给药方式下药物的转移、消除与血药浓度成线性关系，建立两室血药浓度与时间的微分方程组。其中，快速静脉注射可简化为的瞬时将剂量的药物输入中心室，血药浓度立即上升为且在中心室分布均匀；恒速静脉滴注的速率为常数，经过时间后药物注射完成且中心室中药物分布均匀；口服或肌肉注射相当于在药物输入中心室之前先有一个将药物吸收的过程，可简化为存在一个吸收室，药物由外界进入吸收室，再由吸收室转运到中心室。利用线性微分方程组求解方法可以求解出3种给药方式下血药浓度的表达式，结合具体事例并利用残数法可计算出相关药动学参数。 对问题二，采用问题一中单次给药方式下的血药浓度随时间变化的解析式，并运用迭代的方法可以得出固定时间间隔及固定剂量的多次给药方式下，不同时段血药浓度随时间变化的解析式。 三、问题假设 1.机体分为中心室（1室）和周边室（2室），且两个室的容积（即血液体积或药物分布容积）在整个讨论过程中保持不变。 2.药物从一室向另一室的转运速率，以及向体外的排除速率，与该室的血药浓度成正比。 3.只有中心室与体外有药物交换，即药物从体外进入中心室，最后又从中心室排除体外，与转运和排除的数量相比，药物的吸收可以忽略。 4.静脉注射后中心室的药物浓度分布均匀。 四、符号表示 符号 变量说明 时间 第室的药量（=1,2，3分别表示中心室、周边室、吸收室） 第室的血药浓度（=1,2） 第室的容积（=1,2） 药物从中心室向周边室的转移速率系数 药物从周边室向中心室的转移速率系数 药物从中心室向体外排除的速率系数 药物从吸收室向中心室的转移速率系数 恒速静脉滴注的给药速率 给药速率 给药量 五、模型建立与求解 模型一： 1、三种不同给药方式的模型建立与相应求解过程 根据假设条件和已知二室模型（即图1），写出二室模型中，满足的微分方程组： （1） 与血药浓度，房室容积之间的关系式为 ； （2） 将（2）式代入（1）式解出： （3） 这是线性非齐次方程组，它的解由对应的齐次方程组的通解和非齐次方程组的特解构成。对应齐次方程组的通解为： （4） 其中为矩阵的特征根，且满足： （5） 根据3种不同的给药方式代入不同的初始值，求出特解，从而解出3种给药方式下血药浓度与时间的关系式，并利用matlab软件画图。下面将对三种不同的给药方式分别进行讨论： 1.1快速静脉注射 在时刻，将剂量的药物注入静脉，中心室血药浓度立即上升为，初始条件为： 时 （6） 则方程组（3）在（6）式的条件下的解为： （7） 其中， （8） 1.2恒速静脉滴注 开始注射后，以速率向静脉中滴注药物，初始条件为： 时，