数模学习第三天--微分方程（药物分布与排除）

前言

人们常常采用一种叫“房室系统”的观点来考察问题.根据研究对象的特征或研究的不同精度要求，我们把研究对象看成一个整体（单房室系统）或将其剖分成若干个相互存在着某种联系的部分（多房室系统）。

此问题涉及到的方面。
1.药物进入机体形成血药浓度(单位体积血液的药物量)
2.血药浓度需保持在一定范围内——给药方案设计。
3.药物在体内吸收、分布和排除过程 ——药物动力学。
4.建立房室模型——药物动力学的基本步骤。
5.房室——机体的一部分，药物在一个房室内均匀分布，在房室间按一定规律转移

建立模型

体内药物在任一时刻都是均匀分布的，设t时刻体内药物的总量为x(t)；系统处于一种动态平衡中，即成立着关系式：

药物的分解与排泄（输出）速率通常被认为是与药物当前的浓度成正比的，即：
(dx/dt)出=kx;
且药物的输入规律与给药的方式有关
情况1 快速静脉注射
情况2 恒速静脉点滴
情况3 口服药或肌注

情况1 快速静脉注射

在快速静脉注射时，总量为D的药物在瞬间被注入体内。设机体的体积为V，则我们可以近似地将系统看成初始总量为D，浓度为D/V，只输出不输入的房室，即系统可看成近似地满足微分方程：
C（t）为药物浓度

情况2 恒速静脉点滴

药物似恒速点滴方式进入体内

情形3 口服药或肌注

口服药或肌肉注射时，药物的吸收方式与点滴时不同，药物虽然瞬间进入了体内，但它一般都集中于身体的某一部位（吸收室），靠其表面与肌体接触而逐步被吸收．设药物被吸收的速率与存量药物的数量成正比，记比例系数为k1，即若记t时刻吸收室残留药物量为y(t)，则y满足：

解得：

利用MATLAB解得三种结果

二室模型

上述是将机体看成一个均匀分布的同质单元，故被称单房室模型，但机体事实上并不是这样．药物进入血液，通过血液循环药物被带到身体的各个部位，又通过交换进入各个器官．因此，要建立更接近实际情况的数学模型就必须正视机体部位之间的差异及相互之间的关联关系，这就需要多房室系统模型。
k12表示由室I渗透到室II的变化率前的系数，而k21则表示由室II返回室I的变化率前的系数，它们刻划了两室间的内在联系，I中心室(心、肺、肾等)和II周边室(四肢、肌肉等)

模型假设
1.中心室(1)和周边室(2),容积不变
2.药物从体外进入中心室，在二室间相互转移,从中心室排出体外。
3.药物在房室间转移速率及向体外排除速率，与该室血药浓度成正比。


1.快速静脉注射

2.恒速静脉滴注

代入原始方程可解。
3.口服或肌肉注射
相当于药物( 剂量D0)先进入吸收室，吸收后进入中心室。