药物动力学模型一般说来，一种药物要发挥其治疗疾病的作用，必须进入血液，随着血流到达作用部位。
药物从给药部位进入血液循环的过程称为药物的吸收，而借助于血液循环往体各脏器组织转运的过程称为药物的分布。
药物进入体以后，有的以厡型发挥作用，并以厡型经肾脏排出体外；有的则发生化学结构的改变--称为药物的代。
代产物可能具有药理活性，可能没有药理活性。
不论是厡型药物或其代产物，最终都是经过一定的途径(如肾脏、胆道、呼吸器官、唾液腺、汗腺等)离开机体，这一过程称为药物的排泄。
有时，把代和排泄统称为消除。
药物动力学(Pharmacokinetics)就是研究药物、毒物及其代物在体的吸收、分布、代及排除过程的定量规律的科学。
它是介于数学与药理学之间的一门新兴的边缘学科。
自从20世纪30年代Teorell为药物动力学奠定基础以来，由于药物分析技术的进步和电子计算机的使用，药物动力学在理论和应用两方面都获得迅速的发展。
至今，药物动力学仍在不断地向深度和广度发展。
药物动力学的研究方法一般有房室分析；矩分析；非线性药物动力学模型；生理药物动力学模型；药物药效学模型。
下面我们仅就房室分析作一简单介绍。
为了揭示药物在体吸收、分布、代及排泄过程的定量规律，通常从给药后的一系列时间(t)采取血样，测定血(常为血浆，有时为血清或全血)中的药物浓度( C )；然后对血药浓度——时间数据数据(C——t数据)进行分析。
一一室模型最简单的房室模型是一室模型。
采用一室模型，意味着可以近似地把机体看成一个动力学单元，它适用于给药后，药物瞬间分布到血液、其它体液及各器官、组织中，并达成动态平衡的情况。
下面的图(一)表示几种常见的给药途径下的一室模型，其中C代表在给药后时间t的血药浓度，V代表房室的容积，常称为药物的表观分布容积，K 代表药物的一级消除速率常数，故消除速率与体药量成正比，D代表所给剂量。
图(a)表示快速静脉注射一个剂量D，由于是快速，且药物直接从静脉输入，故吸收过程可略而不计；图(b)表示以恒定的速率K，静脉滴注一个剂量D；若滴注所需时间为丅，则K=D/丅。
图(c)表示口服或肌肉注射一个剂量D，由于存在吸收过程，故图中分别用F和K代表吸收分数和一级吸收速率常数。
1. 快速静脉注射在图(a)中所示一室模型的情况下，设在时间t，体药物量为x，则按一级消除的假设，体药量减少速率与当时的药量成正比，故有下列方程：dxKtdt(5.1)快速静脉注射恒速静脉滴注口服或肌肉注射0K FK图(一)初始条件为t=0，x=0，容易解得-Ktx=De……………………..(5.2)注意到房室的容积为V，故c=x/V；记t=0时血药浓度为C，因此0C=D/V，则有C=C Kte…………………….(5.3)这就是快速静脉注射(简称静注)一个剂量D时，符合一室模型的药物及其血药浓度随时间递减的方程。
对方程3两边取对数得ln lnC C Kt这表明在一室模型的情况下，将实测的C_t数据在以t为横轴，ln C 为纵轴的坐标系上作图，各个数据点应呈直线散布趋势。
据此，用图测法或最小二乘法拟合一条直线，其斜率为K，截距为0ln C，于是K和C便可求得。
当然，如果数据点的散布明显地不是呈直线趋势，则可断言不宜采用一室模型来解释该药物在快速静脉注射时的体动力学过程。
在实际应用中，表征药物消除快慢常用的参数是生物半衰期，记为1t /2，它是指药物浓度降至原定值的一半所需的时间。
在方程(3)中令t=1t /2，C =0C /2，可得1ln 20.6922t K K………………(5.4) 可见半衰期是常数，且与消除速率常数成反比。
例如，给一名志愿者一次静脉注射某药物100mg ，测得给药后一些时刻的血药浓度见下表，和在坐标系上作出各数据点，它们是呈直线散布趋势，故可采用一室模型。
一次静注100mg 所得数据5.42 5.32 4.80 4.10如用最小二乘法拟合如下的直线方程ln Ca bt……………..(5.5)利用实测的C 一t 数据计算直线斜率和截距的公式为：1112211111ln ln 11ln nn ni iiii i i nniii i nniii i t C t C n btt naC bt n……………..(5.6)其中n 为C 一t 数据点的个数。
将上表中的有关数据代入(6)式得 b=-0.02744 a=1.7386 于是，拟合数据点的直线方程为lnC=1.7386-0.02744与方程(4)对照，便得0C 和K 的估计值为15.689/,0.0274C g ml K h进而，可得该药物的生物半衰期1/2t 和表观分布容积V 为1/2t =0.69325.30.0274h 010017.65.689D Vl C2.恒速静脉滴注在图(b)所示一室模型的情况不，体药量x 随时间t 变化的微分方程如下：dxK Kxdt……………(5.7)在初始条件t=0，x=0之下，可得其解为1KtK xe K…………..(5.8)其中0,t T ，这里T 为滴注持续的时间。
利用x=VC ，由(8)式得1KtK Ce VK…………(5.9)这就是恒速静脉滴注期间，符合一室模型的药物浓度随时间递增的方程。
假如t=丅时，所给剂量D滴注完毕，则此后的血药浓度便按静注射时的规律下降(如图二)，不过此时初始浓度为01/KTK e VK，故滴注停止后的C一t方程(为区别起见，特记为''C t)如下：''1KTKtK eC eVK…………..(5.10)由此可见，我们可以从滴注停止后测得''C t数据，求得K和V的估计值(K和丅皆已知)假如滴注总是持续进行，则由(10)式可知，血药浓度将趋于一个极限，记作001lim KtsstK KC eVK VK…………..(5.11)这个血药浓度称为稳态浓度，又称坪水平。
记在时刻t的血药浓度达到坪水平的分数为ssf，则有1/20.69311ttKtssssCf e eC…………..(5.12)可见达到稳态的快慢取决于消除速率常数K 或半衰期，与滴注速率K 无关。
例如，当滴注持续时间等于5倍半衰期时，由(12)式算得0.969ssf ，此时血药浓度约为坪水平彻97％。
3.口服或肌肉注射在图（c ）所示一室模型的情况下，设在时刻t ，体药量为x ，吸收部位的药量为a x ，则可建立如下的微分方程组a a aa adxK x Kx dtdx K x dt……………(5.13) 在初始条件t=0，a x =FD ，x=0之下，可解得a K tKta a K FDxeeK K…………(5.14)从而血药浓度随时间变化的方程为a K tKta a K FDCeeV K K………… (5.15)令M=/a aK FD V K K，则上式可写为 a K tKtC M ee………… (5.16)在通常情况下，吸收比消除快的多，即aK K ，故对于足够大的t ，血药浓度实际上是时间的单项指数函数，为区别起见，记为Kt C Me *-=………… (5.17)或 ln ln C M Kt *=-………… (5.18)据此可得K 和M 的估计值，然后计算足够大的t 之前各个实测浓度与按 (5.17) 式推算的C 与C 值之差称为“剩余浓度”r C :a K trC C C Me………… (5.19)或ln ln r a C M K t =-………… (5.20)据此可得K 的估计值。
上述这种估计消除和吸收速率常数的方法称为剩余法。
(二)二室型二室模型是从动力学角度把机体设想为两部分，分别称为中央室和周边室。
中央室一般包括血液及血流丰富的组织(如心、肝、肾等)，周边室一般指血液供应少，药物不易进入的组织(如肌肉、皮肤、某些脂肪组织等)。
在快速静注的情况下常见的二室模型如图4-2 所示。
图中1V 代表中央室的容积，10k 代表药物从中央室消除的一级速率常数，12k 和21k 分别代表药物从中央室到周边室和反方向的一级转运速率常数，其余符号同前。
设在时刻t ，中央室和周边室中的药物量分别为1x 和2x ，则可写出下列微分方程组：()1212121012121212dx k x k k x dt dx k x k xdt⎧=-+⎪⎪⎨⎪=-⎪⎩………… (5.14) 在初始条件120,,0t x D x ===之下，可解得()()()21121122t tt t D k D k x e e Dk x e e αββααβαβαβαβ----⎧--=+⎪--⎪⎨⎪=-⎪-⎩…………(5.15)其中α和β由下列关系式决定：1012211021k k k k k αβαβ+=++⎧⎨=⎩………… (5.16) 通常规定α＞β。
由于11x V C =，故描述血药浓度随时间变化的方程为()()()()212111t tD k D k C e e V V αβαβαβαβ----=+--………… (5.17)令()()()()211211/,/A D k V B D k V ααββαβ=--=-- 则有at t C Ae Be β--=+………… (4.18)根据(4.18)式，利用实测C ——t 数据，用剩余法或电子计算机作曲线拟合，可得α、β、及A 、B 的值，而后按下列公式计算模型参数：1211021121021D V A B A B k A B k k k k k βααβαβ⎧=⎪+⎪+⎪=⎪+⎨⎪=⎪⎪⎪=+--⎩………… (4.19) 这组公式不难从(4.17)、(4.18)式及A 、B 的定义导出。
(三) 多次给药在临床药物治疗中绝大多数药物都需要多次给药，以使血药浓度在足够长的一段时间处于安全，有效的治疗围。
因此，认识多次给药下血药浓度的变化规律是拟订合理的给药方案的基础。
这里，我们只讨论一室模型多次重复静活的情况。
假定某药在快速静注下，符合一室模型的动力学规律，那末，每隔一段时间τ，静注一个剂量D 时，血药浓度C 随时间t 将如何变化呢？静注第一剂后，C —t 关系为0kt C C e -=其中0/,0C D V t τ=≤≤，显然，最高浓度为0C ，最低浓度为0k C e τ-，记为()10max C C =()10min k C C eτ-=不难理解，静注第二剂后，则有()()()()()2000max 220min max 11K K K K K C C C e C e C C e C e e τττττ-----=+=+==+ 静注n 剂后，就有()()100max 111nK n K K n K e C C e e C e ττττ-----⎛⎫-⎡⎤=+++= ⎪⎣⎦-⎝⎭…………… (5.21)()()0min max 11nK K K n n K e C C e C e e ττττ----⎛⎫-== ⎪-⎝⎭…………(5.22) 由此可知，重复静注n 剂后，血药浓度随时间的变化规律为011nK K n K e C C e e τττ---⎛⎫-= ⎪-⎝⎭0t τ≤≤………… (5.23) 假如n 充分大，使血药浓没达到稳态，那么，对(5.22)式取n →∞的极限，使得稳态浓度的变化规律为01K K C C e e ττ-∞-⎛⎫= ⎪-⎝⎭………… (5.24) 最高和最低稳态浓度分别为()()0max 11K K C D C e V e ττ∞--==--………… (5.25) ()()00min 111K K K K C C D C e e e V e ττττ-∞-⎛⎫=== ⎪---⎝⎭………(5.26) 在一个给药间隔时间，平均稳在浓度为00111T TKtss KC DC C dt e dte VKττττ-∞-⎛⎫===⎪-⎝⎭⎰⎰…………(5.27)图4-4表示每隔6小时重复静注一个剂量D产生的C——t曲线最后，我们举一个实例。