数学建模第二章微积分方法建模--2.10药物在体内的分布与排除

§10药物在体内的分布与排除药物进入机体后，在随血液输送到各器官和组织的过程中，不断地被吸收、分布、代谢，最终排出体外。药物在血液中的浓度mvg称血药浓度。血药浓度的大小直接影响到药物的疗效，浓度太低不能达到预期的治疗效果，浓度太高又可能导致中毒、副作用太强或造成浪费。因此，研究药物在体内吸收、分布和排除的动态过程，对于新药研制时剂量的确定、给药方案设计等药理学和临床医学的发展具有重要的指导意义和实用价值。为了研究目的，将一个机体划分成若干个房室，每个房室是机体的一部分，比如中心室和周边室。在一个房室内药物呈均匀分布，而在不同的房室之间按一定规律进行转移。模型假设1、药物进入机体后，全部进入中心室（血液较丰富的心、肺、肾等器官和组织），中心室的容积在给药过程中保持不变；2、药物从中心室排出体外，与排除的数量相比，药物的吸收可以忽略；3、药物排除的速率与中心室的血药浓度成正比。模型建立记tf0―给药速率tc―中心室血药浓度tx―中心室药量V―中心室容积k―排除速率系数中心室)(),(txtcVk排除给药)(0tf一、求解各种给药方式下血药浓度变化情况上述各量间有关系kxtfx0即tfkxx0又tVctx得方程Vtftkctc0（1）1、快速静脉注射设给药量D,则初始条件VDc0，00tf（1）的解为tkeVDtc（2）00cttc2、恒速静脉注射设持续时间为，注射速率为0k，则有当t0时，00ktf，初始条件00c，当t时，00tf，初始条件keVkkc10，（1）的解为teeVkkteVkktctkktk10100（3）t)(tc03、口服或肌肉注射在药物输入中心室之前先有一个将药物吸收入血液的过程，可以看作为有一个吸收室，药物由吸收室进入中心室，药物由吸收室进入中心室的转移速率系数记成1k，给药量D，吸收室药量tx0，则吸收室tx0中心室txktf010Dxxktx00010tkeDtx10于是tkeDktf110，初始条件00c，（1）的解为tktkeekkVDktc111，kk1（4）tct0二、各种给药方式下参数估计1、快速静脉注射下估计k在0t时刻快速注射剂量为D的药物后，在一系列时刻itni,,2,1从中心室取血样获得血药浓度itc。由（2）式反解出几个itk，取算术平均值就得到k的估计值：niitknk112、恒速静脉滴注下估计k方法与1、类似，it取在,0内，用（3）中第一式反解iitk。3、口服或肌肉注射下估计k和1k因为kk1，记kkVDkA11，于是当t充分大时（4）近似为tkAetc或ktAtclnln对于适当大的it和测得的相应itc，用最小二乘法估计出k和Aln，从而再由kkVDkA11就可以估计出DAVAVkk1模型分析当要求精度较高时,可采用二室甚至多室模型,这时的机理分析和参数估计都比一室模型难度更大。需要建立微分方程组来进行分析。例如:二室模型图示如下。tf0给药中心室111,Vtxtc周边室222,Vtxtc12k21k13k排除