《生物药剂学与药物动力学》课程教学资源（文献资料）血管外给药的非线性房室模型解的逼近_胡晓虎

应用数学和力学，第35卷第9期Applied Mathematics and Mechanics2014年9月15日出版Vol.35,No.9,Sep.15,2014文章编号：1000-0887（2014)09-1033-13应用数学和力学编委会，ISSN1000-0887血管外给药的非线性房室模型解的逼近唐三一胡晓虎，（陕西师范大学数学与信息科学学院，西安710062）（我刊编委唐三一来稿）摘要：药物动力学模型的解析求解公式在新药设计特别是药物动力学参数确定等方面具有非常重要的意义.近年来，由于非线性米氏消除速率过程确定的药物动力学模型解析求解公式的获得，使得大多数单房室模型的解析解基本确定.但是，由于刻画血管外给药的非线性米氏药物动力学模型是一个非自治系统，进而不可能寻求其解析求解公式.该文的目的是讨论一次性血管外给药和周期血管外给药下非线性药物动力学模型解的逼近问题.采用微分方程和脉冲微分方程的比较定理并借助LambertW函数的定义以及相关性质给出模型的不同上下界，估计模型解的逼近程度，并通过数值模拟进行验证关键词：房室模型：米氏消除速率；血管外给药：LambertW函数：药物动力学中图分类号：0241.8：0242文献标志码：Adoi:10.3879/j.issn.1000-0887.2014.09.009引言药物的主动转运，肝脏的代谢，胆汁的分泌以及肾小管的分泌都需要酶或载体参与.这些酶或载体参与的过程不但特异性强，而且受酶或载体本身数量的限制，也就是说其转运或消除速度有一定的限制.上述特性决定了一级消除、零级消除速率过程不能准确刻画药物在体内的消除过程，需要采用基于酶催化反应而提出的非线性米氏消除速率过程-来刻画，即有下面的方程：-VC(t)dc(t)D(1), C(0) = C = ,=- + C()dt其中C（t）表示血药浓度，Vmx表示最大反应速率，K是Michaelis-Menten常数，D是一次静脉注射给药的药物剂量，V是中心室表观分布容积.来氏消除速率过程在药物动力学中应用非常广泛，尤其是在酶动力学中-Lundquist和Wolthers在1958年将米氏消除速率过程应用在恒定静脉注射和周期静脉注射中，提出了模型（1）并得到了该模型如下积分形式的解析解：1[DD- C(t) +Kmln((2)=t-loVlv(Vc())收稿日期：2014-03-18；修订日期：2014-06-16基金项目：国家自然科学基金（11171199)作者简介：胡晓虎（1989一），男，湖北宜城人，硕士：唐三一（1970一），男，教授，博士生导师（通讯作者，E-mail：sytang@snnu.edu.cn）1033?1994-2014 China Academic Journal Electronic Publishing House.All rights reserved.http://www.cnki.net

书 文章编号:1000-0887(2014)09-1033-13  应用数学和力学编委会，ISSN 1000-0887 血管外给药的非线性房室模型解的逼近* 胡晓虎， 唐三一 (陕西师范大学 数学与信息科学学院，西安 710062) (我刊编委唐三一来稿) 摘要: 药物动力学模型的解析求解公式在新药设计特别是药物动力学参数确定等方面具有非常 重要的意义． 近年来，由于非线性米氏消除速率过程确定的药物动力学模型解析求解公式的获得， 使得大多数单房室模型的解析解基本确定． 但是，由于刻画血管外给药的非线性米氏药物动力学 模型是一个非自治系统，进而不可能寻求其解析求解公式． 该文的目的是讨论一次性血管外给药 和周期血管外给药下非线性药物动力学模型解的逼近问题． 采用微分方程和脉冲微分方程的比较 定理并借助 Lambert W 函数的定义以及相关性质给出模型的不同上下界，估计模型解的逼近程度， 并通过数值模拟进行验证． 关 键 词: 房室模型; 米氏消除速率; 血管外给药; Lambert W 函数; 药物动力学 中图分类号: O241． 8; O242 文献标志码: A doi: 10． 3879 /j． issn． 1000-0887． 2014． 09． 009 引 言 药物的主动转运，肝脏的代谢，胆汁的分泌以及肾小管的分泌都需要酶或载体参与． 这些 酶或载体参与的过程不但特异性强，而且受酶或载体本身数量的限制，也就是说其转运或消除 速度有一定的限制． 上述特性决定了一级消除、零级消除速率过程不能准确刻画药物在体内的 消除过程，需要采用基于酶催化反应而提出的非线性米氏消除速率过程［1-4］来刻画，即有下面 的方程: dC(t) dt = － VmaxC(t) KM + C(t) ，C(0) = C0 = D V ， (1) 其中 C(t) 表示血药浓度，Vmax 表示最大反应速率，KM 是 Michaelis-Menten 常数，D 是一次静脉 注射给药的药物剂量，V 是中心室表观分布容积． 米氏消除速率过程在药物动力学中应用非常 广泛，尤其是在酶动力学中［5-6］． Lundquist 和 Wolthers［7］在 1958 年将米氏消除速率过程应用在恒定静脉注射和周期静脉 注射中，提出了模型(1)并得到了该模型如下积分形式的解析解: 1 Vmax D V － C(t) + KM ln D VC(t [ ( ) ] ) = t － t0 ． (2) 3301 应用数学和力学，第 35 卷 第 9 期 2014 年 9 月 15 日出版 Applied Mathematics and Mechanics Vol． 35，No． 9，Sep． 15，2014 * 收稿日期: 2014-03-18; 修订日期: 2014-06-16 基金项目: 国家自然科学基金(11171199) 作者简介: 胡晓虎(1989—)，男，湖北宜城人，硕士; 唐三一(1970—)，男，教授，博士生导师(通讯作者． E-mail: sytang@ snnu． edu． cn)．

1034胡晓虎唐三一利用上述积分形式的解，作者结合试验数据进行了药物动力学模型的数据拟合和相应药物动力学参数估计的研究。后来，Beal,Godfrey和Fitch[i系统地研究了积分方程（2）关于药物浓度C(t）的可解性。进一步，Goudar等，Schnell和Mendozai2分别给出了模型具有一次静脉推注时如下的解析求解公式： exp(Co - V (t -t0))C(t) = KmW(e(3)KMKM其中，W（）表示LambertW函数，下同.有关LambertW函数定义和有关性质在第1节给出，Tang和Xiaoi首先给出具有周期静脉推注给药方式下模型（1）的解析求解公式，同时建立了具有治疗窗口的药物动力学模型并给出了模型的解析求解公式.特别地，文献首次对如下恒定给药下的非线性药物动力学模型的解析求解公式进行了系统研究：VmaxC(t)dc(t)D(4)ddo =r-k+Co). C(c) = Co = D,其中r代表常数输入率。Tang和Xiao同样借助LambertW函数的定义和性质给出了其相应的解析求解公式：C(t) =TKubK (-1, - exp(D(_rKm + bCo +b (t - t0)b>0,KuVaasKyVux-Km+/(Km+C)*+2Vm（t-t),b = 0,rKM-bKumwfo, -kthcexp(-k+Cet(-)b<0,KVmaxKuVaux其中b=r-Vmx近些年，针对上述解析求解公式的数值计算和近似计算，也得到深入研究[44S对于血管外给药（比如口服给药），由于药物从用药部位进入血液循环要经过吸收过程，因此，通过口服剂量为D的药物首先进入吸收室，在室内逐渐被吸收，同时将吸收到的药物逐渐向中心室输送，根据这一过程可以建立相应的血管外给药的单房室模型.如果记C，为吸收室药物的浓度，V，为吸收室表观分布容积，则药物浓度在吸收室的变化规律为dC, (t)D= - K,C,(t), C,(0) = C =dtV.其中，K，为药物在吸收室的消除速率常数，求解上述方程得到一次性血管外给药吸收室的浓度为C, (t) = Cge-ky.如果假设从吸收室吸收的药物以生物利用度F的比率输送到中心室，因此中心室的血药浓度C（t）的变化规律可以描述为?1994-2014 China Academic Journal Electronic Publishing House.All rights reserved.http://www.cnki.net

利用上述积分形式的解，作者结合试验数据进行了药物动力学模型的数据拟合和相应药物动 力学参数估计的研究． 后来，Beal［8-9］，Godfrey 和 Fitch［10］系统地研究了积分方程(2)关于药物浓度 C(t) 的可解 性． 进一步，Goudar 等［11］，Schnell 和 Mendoza［12］分别给出了模型具有一次静脉推注时如下的 解析求解公式: C(t) = KM W C0 KM exp C0 － Vmax(t － t0 ) K ( ( ) ) M ， (3) 其中，W() 表示 Lambert W 函数，下同． 有关 Lambert W 函数定义和有关性质在第 1 节给出． Tang 和 Xiao［13］首先给出具有周期静脉推注给药方式下模型(1)的解析求解公式，同时建 立了具有治疗窗口的药物动力学模型并给出了模型的解析求解公式． 特别地，文献［7］首次对 如下恒定给药下的非线性药物动力学模型的解析求解公式进行了系统研究: dC(t) dt = r － VmaxC(t) KM + C(t) ，C(t + 0 ) = C0 = D V ， (4) 其中 r 代表常数输入率． Tang 和 Xiao［13］同样借助 Lambert W 函数的定义和性质给出了其相应 的解析求解公式: C(t) = － rKM b + KM VmaxW － 1，－ rKM + bC0 KM Vmax exp － rKM + bC0 + b2 (t － t0 ) KM V ( ( ) ) max ， b ＞ 0， － KM + (KM + C0 )2 + 2Vmax(t － t 槡 0 )， b = 0， － rKM b + KM VmaxW 0，－ rKM + bC0 KM Vmax exp － rKM + bC0 + b2 (t － t0 ) KM V ( ( ) ) max ， b ＜ 0              ， 其中 b = r － Vmax ． 近些年，针对上述解析求解公式的数值计算和近似计算，也得到深入研究［14-15］． 对于血管外给药(比如口服给药)，由于药物从用药部位进入血液循环要经过吸收过程． 因此，通过口服剂量为 D 的药物首先进入吸收室，在室内逐渐被吸收，同时将吸收到的药物逐 渐向中心室输送，根据这一过程可以建立相应的血管外给药的单房室模型． 如果记 Ca 为吸收 室药物的浓度，Va 为吸收室表观分布容积，则药物浓度在吸收室的变化规律为 dCa(t) dt = － KaCa(t)，Ca(0) = Ca 0 = D Va ， 其中，Ka 为药物在吸收室的消除速率常数． 求解上述方程得到一次性血管外给药吸收室的浓 度为 Ca(t) = Ca 0 e －Kat ． 如果假设从吸收室吸收的药物以生物利用度 F 的比率输送到中心室，因此中心室的血药 浓度 C(t) 的变化规律可以描述为 4301 胡 晓 虎 唐 三 一

1035血管外给药的非线性房室模型解的逼近V.C(t)dc(t)= FK,C, (t)KM+C()=dtV.C(t)(5)FCgK,e-ka --K + C(t)明显地，上述方程是一个非自治的一阶微分方程，由于变量C和t是不可分离的，因此无法求出该方程的解析解，由于实际问题和药物动力学研究的需要，有必要知道上述血管外给药下药物在体内的浓度如何变化？因此采用微分方程和脉冲微分方程的比较定理并借助LambertW函数的定义以及相关性质给出模型（5）的不同上下界，估计模型解的逼近程度，并通过数值模拟分析解的可行性，1LambertW函数和特殊情形在介绍具体的方法之前，首先给出Lambert W函数的定义和一些基本性质[定义1LambertW函数是函数f(z）=ze的反函数，且满足：W(z)exp(W(z)) = z,为了方便使用，本文将LambertW函数记做W.注意函数f(z）=ze在z>-1时有正导数f（z）=z+1)e.在区间[-1,αo）上定义f(z）=ze的反函数为W(0,z)=W(z)：易得当z>0时有W（0,z）>0，且当z>-1时W(0，z）为增函数.类似地，可以在区间（-，-1）上定义ze的反函数为W（-1，z).通过计算容易得到其导数满足如下关系式：W(2)W (z) =z(1 + W(z))为了研究模型的解的性质，首先考虑以下特殊情形情形1米氏常数远远大于血药浓度时的情形，即K》C.假设K>C，则方程（5）可以简化成Vc(t).dc()= FCK,e-k .(6)dtKu上述方程是一个简单的一阶线性非齐次方程，求解得到血药浓度C（t）满足：, VmaxFCK,KM(e-Vmat/KMK≠e-kt),K,K- VKMC(t):Va[FCK,te-*,K, =KM由此可见，当K》C时方程（5）刻画了血管外给药时体内药物浓度与时间的关系，即药物浓度在时刻KMIn(Vmxtmax = - 7KKM-V.KKM时达到峰值，血药浓度在时间区间（0，t）内由小变大，然后从t开始逐渐变小，其数值结果如图1所示，情形2米氏常数远远小于血药浓度时，即有KC假设K≤C，则方程（5）可以简化成dc(t)= - V + FCgK,e-Ka.dt?1994-2014 China Academic Journal Electronic Publishing House.All rights reserved.http:/www.cnki.net

dC(t) dt = FKaCa(t) － VmaxC(t) KM + C(t) = FCa 0Ka e －Kat － VmaxC(t) KM + C(t) ． (5) 明显地，上述方程是一个非自治的一阶微分方程，由于变量 C 和 t 是不可分离的，因此无法求 出该方程的解析解． 由于实际问题和药物动力学研究的需要，有必要知道上述血管外给药下药 物在体内的浓度如何变化? 因此采用微分方程和脉冲微分方程的比较定理并借助 Lambert W 函数的定义以及相关性质给出模型(5)的不同上下界，估计模型解的逼近程度，并通过数值模 拟分析解的可行性． 1 Lambert W 函数和特殊情形 在介绍具体的方法之前，首先给出 Lambert W 函数的定义和一些基本性质［16］． 定义 1 Lambert W 函数是函数 f(z) = zez 的反函数，且满足: W(z)exp(W(z)) = z， 为了方便使用，本文将 Lambert W 函数记做 W． 注意函数 f(z) = zez 在 z ＞ － 1 时有正导数 f '(z) = (z + 1)ez ． 在区间［－ 1，∞ ) 上定义 f(z) = zez 的反函数为 W(0，z) = W(z)． 易得当 z ＞ 0 时有 W(0，z) ＞ 0，且当 z ＞ － 1 时 W(0，z) 为增函数． 类似地，可以在区间( － ∞ ，－ 1］上定义 zez 的反函数为 W( － 1，z)． 通过计算容易得到其导数满足如下关系式: W'(z) = W(z) z(1 + W(z))． 为了研究模型的解的性质，首先考虑以下特殊情形． 情形 1 米氏常数远远大于血药浓度时的情形，即 KM  C． 假设 KM  C，则方程(5)可以简化成 dC(t) dt = FCa 0Ka e －Kat － Vmax KM C(t)． (6) 上述方程是一个简单的一阶线性非齐次方程，求解得到血药浓度 C(t) 满足: C(t) = － FCa 0KaKM KaKM － Vmax e －Vmaxt/KM － e ( －Ka ) t ， Ka ≠ Vmax KM ， FCa 0Ka te －Kat ， Ka = Vmax KM { ． 由此可见，当 KM  C 时方程(5)刻画了血管外给药时体内药物浓度与时间的关系，即药物浓 度在时刻 tmax = － KM KaKM － Vmax ln Vmax KaK ( ) M 时达到峰值，血药浓度在时间区间 (0，tmax) 内由小变大，然后从 tmax 开始逐渐变小，其数值结果 如图 1 所示． 情形 2 米氏常数远远小于血药浓度时，即有 KM  C． 假设 KM  C，则方程(5)可以简化成 dC(t) dt = － Vmax + FCa 0Ka e －Kat ． 血管外给药的非线性房室模型解的逼近 5301

1036胡晓虎唐三求解上述微分方程得到血药浓度C(t）为C(t) = C, - Vt - FCfe-kt.容易看出，若Vm≥FCK，此时C（t）≤0，即在此时血药浓度始终是降低的8(ru/Bu)/)5410203040501/hCo=15mg/mL,Vmx=1 mg/(h-mL),Km =10mg/mL,K, =0.4mg/(hamL),F =0.8图1当KC时方程（5)刻画的血药浓度与时间的关系Fig. 1Plasncording to model (5) when Ky Ctratinnthiy2.01.5C-00.5C(0)C,(0)~0.51051520253040503545t/hCg=2mg/mL,Vmx=1 mg/(hmL),Km=1.2mg/mL,K,=0.1 mg/(hamL),F=1图2模型（5)血药浓度的一个粗略上下界Fig. 2A rough upper and lower bounds of model (5)情形3方程（5）刻画的血药浓度的一个粗略上下界情形1和情形2在两种非常特殊的情况下得到了模型（5）的解析求解公式.正如引言中所强调的那样，一般形式的模型（5）已经不可能得到其解析公式了.因此，笔者希望借助引言中讨论的两类可求解非线性米氏消除速率过程刻画的药物动力学模型，即模型（1）和（4），通过微分方程的比较定理得到模型（5）解的一个粗略下界和上界为此记C,（）和C,（）分别表示下述两个模型?1994-2014 China Academic Journal Electronic Publishing House.All rights reserved.http://www.cnki.net

求解上述微分方程得到血药浓度 C(t) 为 C(t) = C0 － Vmax t － FCa 0 e －Kat ． 容易看出，若 Vmax ≥ FCa 0Ka，此时 C'(t) ≤ 0，即在此时血药浓度始终是降低的． Ca 0 = 15 mg /mL，Vmax = 1 mg /(h·mL)，KM = 10 mg /mL，Ka = 0． 4 mg /(h·mL)，F = 0． 8 图 1 当 KM  C 时方程(5)刻画的血药浓度与时间的关系 Fig． 1 Plasma concentration vs． time according to model(5) when KM  C Ca 0 = 2 mg /mL，Vmax = 1 mg /(h·mL)，KM = 1． 2 mg /mL，Ka = 0． 1 mg /(h·mL)，F = 1 图 2 模型(5)血药浓度的一个粗略上下界 Fig． 2 A rough upper and lower bounds of model(5) 情形 3 方程(5)刻画的血药浓度的一个粗略上下界． 情形 1 和情形 2 在两种非常特殊的情况下得到了模型(5)的解析求解公式． 正如引言中 所强调的那样，一般形式的模型(5)已经不可能得到其解析公式了． 因此，笔者希望借助引言 中讨论的两类可求解非线性米氏消除速率过程刻画的药物动力学模型，即模型(1)和(4)，通 过微分方程的比较定理得到模型(5)解的一个粗略下界和上界． 为此记 C1 (t) 和 C2 (t) 分别表示下述两个模型 6301 胡 晓 虎 唐 三 一

1037血管外给药的非线性房室模型解的逼近dC,(t)VC,()dC, (t)Vm C, (t)=FCoK,-dtK+CdtKM+C,()的解，其解析求解公式引言中已经给出.由于VmC, (t)VC, (t)VC(t)V'则N，是正的.当g,=0时容易得到方程(9）在区间[n,n+1）上的解为C,(o) =-K + /(C,(n) + K)? +2u. (t -n).若g。+0时，模型（9）在区间[n，n+1）上的解可以转化为下述等价形式：(1 + Kug,)dc,() = g.d..u,+g,C,())将上述方程右端从n到t开始积分得+Kug. -"s dc,(t) = g, J.dt1+u.+g,C,(o))?1994-2014 China Academic Journal Electronic Publishing House.All rights reserved.http:/www.cnki.net

dC1 (t) dt = － VmaxC1 (t) KM + C1 (t) ，dC2 (t) dt = FCa 0Ka － VmaxC2 (t) KM + C2 (t) 的解，其解析求解公式引言中已经给出． 由于 － VmaxC1 (t) KM + C1 (t) ＜ FCa 0Ka e －Kat － VmaxC(t) KM + C(t) ≤ FCa 0Ka － VmaxC2 (t) KM + C2 (t) ． 根据微分方程的比较定理，模型(5)的解 C(t) 满足关系 C1 (t) ＜ C(t) ≤ C2 (t)，这样就得到 了模型(5)解的一个粗略下界和上界，如图 2 所示． 由于在寻找模型(5)解的一个上下界时放缩幅度过大，得到了如图 2 所示的精确解与上 下界相距甚远． 这样不能用这两个上下界确定或判别模型(5)解的一些行为特点． 为此，在下 面将采用更加技巧性的方法确定模型(5)的一个非常精确的上下界． 2 一次性血管外给药下模型(5)解的精确上下界 为了寻求模型(5)的一个非常精确的上下界，采用如下的分析办法:在每一个整数区间 ［n，n + 1］(n = 0，1，2，.) 上考虑模型(5)，寻找该模型在区间［n，n + 1］上的一个上下界， 然后推广到整个区间上，实现寻求模型(5)解的一个精确上下界的目标． 在区间［n，n + 1］上首先考虑如下微分方程: dC1 (t) dt = FCa 0Ka e －Ka ［t］ － VmaxC1 (t) KM + C1 (t) ， t ∈［n，n + 1］， (7) 其中［t］表示取整，即在区间［n，n + 1］上［t］ = n 或［t］ = n + 1． 为此先考虑［t］ = n 时的 情形，即在区间［n，n + 1］上先考虑如下的微分方程: dC1 (t) dt = FCa 0Ka e －Kan － VmaxC1 (t) KM + C1 (t) ， t ∈［n，n + 1)． (8) 将上述方程重新整理得 dC1 (t) dt = un + gnC1 (t) KM + C1 (t) ， t ∈［n，n + 1)， (9) 其中 un = FCa 0KaKM e －Kan ，gn = FCa 0Ka e －Kan － Vmax ． 下面根据 gn 的符号考虑几种不同情形． 首先考虑 gn = 0 的情形，此时关于 n 求解 gn = 0 得到一个根，记做 N1，且 N1 = 1 Ka ln FCa 0Ka V ( ) max ． (10) 因此，如果 FCa 0Ka ＞ Vmax， (11) 则 N1 是正的． 当 gn = 0 时容易得到方程(9) 在区间［n，n + 1) 上的解为 C1 (t) = － KM + (C1 (n) + KM )2 槡 + 2un ·(t － n)． 若 gn ≠ 0 时，模型(9) 在区间［n，n + 1) 上的解可以转化为下述等价形式: 1 + KM gn － un un + gnC1 (t ( ) ) dC1 (t) = gn dt． 将上述方程右端从 n 到 t 开始积分得 ∫ C1(t) C1(n) 1 + KM gn － un un + gnC1 (t [ ] ) dC1 (t) = gn ∫ t n dt， 血管外给药的非线性房室模型解的逼近 7301

1038胡晓虎唐三则C,(t) - C, (n) + L,In(u, + g,C, (t))= g.(t - n),n+g,C,(n))其中l。=Km-u,/g.将ung.分别代入l中，可得l。=-（VmKm）/gn由此可见l。的符号完全由g.的符号决定.上式进一步变形可得C, (0) + u,/g exp(C, (t) + u. /g.)InLC (n) + u,/gz exp(C (n) + u /g, +g. (t -n)(12)I.L.C, (n) + u./gn exp(C,(n) + u./g, +g, (t - n))2, (0)=令.In本文分g.>0和g0此时l0，从而C,(n) +u,/g.(C,(n) + u,/g.)z, (t) ≥I.已知对于函数f(x）=xe总有xe≥-e-l因此，当t≥n时，z,（t）≥-e"，故z（t）e[-e"，0)．此时W(ks,() -C(0 +n./ez 0.1.其中k=0，-1.由于W(k() +1 - G@+kN，时：模型(9)的解析解为C(0) =-+L.(-1,C(n) +u /gne)(C,(n) +u,/gn +g,(t-n)).(13)g..In当t→n+1时，"- +I,W(-1,z.),F(C,) = Cut = - "ga其中C, + u./g exFxn(C, + u./g, +g.).L通过上面的讨论可知W（-1,z，）0.+W(-1,z,) C,+u,/g.?1994-2014 China Academic Journal Electronic Publishing House.All rights reserved.http://www.cnki.net

则 C1 (t) － C1 (n) + ln ln un + gnC1 (t) un + gnC1 (n ( ) ) = gn ·(t － n)， 其中 ln = KM － un / gn ． 将 un，gn 分别代入 ln 中，可得 ln = － (VmaxKM ) / gn ． 由此可见 ln 的符号完 全由 gn 的符号决定． 上式进一步变形可得 C1 (t) + un / gn ln exp C1 (t) + un / gn l ( ) n = C1 (n) + un / gn ln exp C1 (n) + un / gn + gn ·(t － n) l ( ) n ． (12) 令 z1 (t) = C1 (n) + un / gn ln exp C1 (n) + un / gn + gn ·(t － n) l ( ) n ， 本文分 gn ＞ 0 和 gn ＜ 0 考虑下面两种情形: ? 假设条件(11)成立． 若 gn ＞ 0，此时 ln ＜ 0，即 n ＜ N1 ． 由于 dz1 (t) dt = gn［C1 (n) + un / gn］ l 2 (n) exp C1 (n) + un / gn + gn ·(t － n) l ( ) n ， 此时 z'1 (t) ＞ 0，从而 z1 (t) ≥ C1 (n) + un / gn ln exp C1 (n) + un / gn l ( ) n ． 已知对于函数 f(x) = xex 总有 xex ≥－ e －1 ． 因此，当 t ≥ n 时，z1 (t) ≥－ e －1 ，故 z1 (t) ∈［－ e －1 ， 0)． 此时 W(k，z1 (t)) = C1 (t) + un / gn ln ＜ 0， 其中 k = 0，－ 1． 由于 W(k，z1 (t)) + 1 = C1 (t) + KM ln ＜ 0， 即 W(k，z1 (t)) ＜ － 1． 因此，方程(12)的解可以由 Lambert W 函数的下支确定． 故当 n ＞ N1 时， 模型(9)的解析解为 C1 (t) = － un gn + lnW － 1， C1 (n) + un / gn ln exp C1 (n) + un / gn + gn ·(t － n) l ( ( ) ) n ． (13) 当 t → n + 1 时， F(Cn )  Cn +1 = － un gn + lnW( － 1，zn )， 其中 zn = Cn + un / gn ln exp Cn + un / gn + gn l ( ) n ． 通过上面的讨论可知 W( － 1，zn ) ＜ － 1，故有 F'(Cn ) = W( － 1，zn ) 1 + W( － 1，zn ) Cn + KM Cn + un / gn ＞ 0． 8301 胡 晓 虎 唐 三 一

1039血管外给药的非线性房室模型解的逼近则AC,=W(-1,z.)(Ca +Km)- (1+W(-1,z,))(C, +I.W(-1,z.) -C, - ",gn即Ca+ = I,W(-1,z.) _ "s(14)gn由于 C. + "./en exp(C. + u./ga) e [- e",0)Z.>IL和dw(-1,z)W(-1,z)2(1 + W(-1,2))Ig.gn即△C。>0.因此，根据0C。()假设条件（11)成立.若g.N，则有l>0.此时z（t）的符号依赖于C,（n）和C，=-u./g,的关系.为此分情况考虑(i)当C,（n）0,从而可以得到z,（t）E[-e-,0).由于方程（12)的右端小于0,从而左端也是小于0的，则（C,（t）+u,/g,）。0,W(k,z) (t)) +1 =In从而W(k,z(t)）>-1.因为W(k,z, (t))d,()=W(kz ()dW(k,z, (t))1+W(kz,() >0,dtz,()[+W(k，z())]dt根据C,（n） 0.由上述已知W(0,z）e（-1,0），由于C, +u,/g.(C. + u./g. +g.) >C, +u./gnexp(C, +u,/g.)e[-e-l,0]Z.I.In1.和?1994-2014 China Academic Journal Electronic Publishing House.All rights reserved.http://www.cnki.net

则 ΔCn = W( － 1，zn )(Cn + KM ) － (1 + W( － 1，zn )) Cn + un g ( ) n = lnW( － 1，zn ) － Cn － un gn ， 即 Cn +1 = lnW( － 1，zn ) － un gn ． (14) 由于 zn ＞ Cn + un / gn ln exp Cn + un / gn l ( ) n ∈［－ e －1 ，0) 和 dW( － 1，z) dz = W( － 1，z) z(1 + W( － 1，z)) ＜ 0， 从而 lnW( － 1，zn ) － Cn － un gn ＞ ln Cn + un / gn l ( ) n － Cn － un gn = 0， 即 ΔCn ＞ 0． 因此，根据 0 ＜ F'(Cn ) ＜ 1 和 ΔCn = Cn +1 － Cn，故始终有 Cn +1 ＞ Cn ． ? 假设条件(11)成立． 若 gn ＜ 0，即 n ＞ N1，则有 ln ＞ 0． 此时 z1 (t) 的符号依赖于 C1 (n) 和 C* n = － un / gn 的关系． 为此分情况考虑． ? 当 C1 (n) ＜ － un / gn，此时 z' 1 (t) ＞ 0，从而可以得到 z1 (t) ∈［－ e －1 ，0)． 由于方程(12) 的右端小于 0，从而左端也是小于 0 的，则(C1 (t) + un / gn ) /ln ＜ 0，即 W(k，z1 (t)) ＜ 0． 可知 W(k，z1 (t)) + 1 = C1 (t) + KM ln ＞ 0， 从而 W(k，z1 (t)) ＞ － 1． 因为 dW(k，z1 (t)) dt = W(k，z1 (t)) z1 (t)［1 + W(k，z1 (t))］ dz1 (t) dt = gn ln W(k，z1 (t)) 1 + W(k，z1 (t)) ＞ 0， 根据 C1 (n) ＜ － un / gn 和方程(9) 解的唯一性，可知式(12) 左边的 C1 (t) 是单调递增的而且趋 向于 C* ． 因此，有 C1 (t) + un / gn ln exp C1 (t) + un / gn l ( ) n 趋向于 0，并且方程(12)的解可以由 Lambert W 函数的实上支确定． 当 t → n + 1 时 F'(Cn ) = W(0，zn ) 1 + W(0，zn ) Cn + KM Cn + un / gn ＞ 0． 由上述已知 W(0，zn ) ∈ ( － 1，0)，由于 zn = Cn + un / gn ln exp Cn + un / gn + gn l ( ) n ＞ Cn + un / gn ln exp Cn + un / gn l ( ) n ∈［－ e －1 ，0) 和 血管外给药的非线性房室模型解的逼近 9301

1040胡晓虎唐三dw(0,z)W(0,z)>0,dzz(1 +W(0,z))从而有(C,+u,/g.)u=0,AC. = I.W(O,z.) - C, - "= > LCg.1.gn即△C。>0.因此，根据0C(i)当C（n）>-u,/g，此时z,（t）>0.由于方程（12）的右端是大于0的，左端同理也是大于0的，从而可知（C（t）+u/g,）/。>0，即W（k，z,（t））>0.因为dw(k,z, ())_ gnW(k,z,(o). 1+W(kz, (0) 0.由上述已知W（0，z）>0，由于C,+u./gh(C,+ua/g,+g.)exnI.L.C, +u,/g.ex(C. + u./g-)e (0, + 0)I.L和dw(0,z)W(0,z)>0,z(1 + W(0,z))dz从而有AC. = I,W(0,z.) - C, - "= > L.( (C, +u./ga)- C, -un= 0,g.B即AC>0.因此，根据0C故当nVm，ln（FCgK,/V）>K,时，同样地，可以得到C,（t）的解析解为?1994-2014 China Academic Journal Electronic Publishing House.All rights reserved.http://www.cnki.net

dW(0，z) dz = W(0，z) z(1 + W(0，z)) ＞ 0， 从而有 ΔCn = lnW(0，zn ) － Cn － un gn ＞ ln Cn + un / gn l ( ) n － Cn － un gn = 0， 即 ΔCn ＞ 0． 因此，根据 0 ＜ F'(Cn ) ＜ 1 和 ΔCn = Cn +1 － Cn，故始终有 Cn +1 ＞ Cn ． ? 当 C1 (n) ＞ － un / gn，此时 z1 (t) ＞ 0． 由于方程(12) 的右端是大于0 的，左端同理也是 大于 0 的，从而可知(C1 (t) + un / gn ) /ln ＞ 0，即 W(k，z1 (t)) ＞ 0． 因为 dW(k，z1 (t)) dt = gn ln W(k，z1 (t)) 1 + W(k，z1 (t)) ＜ 0， 则式(12)左边的 C1 (t) 是单调减少的． 因此方程(12)可以由 Lambert W 函数的实上支确定． 当 t → n + 1 时 F'(Cn ) = W(0，zn ) 1 + W(0，zn ) Cn + KM Cn + un / gn ＞ 0． 由上述已知 W(0，zn ) ＞ 0，由于 zn = Cn + un / gn ln exp Cn + un / gn + gn l ( ) n ＜ Cn + un / gn ln exp Cn + un / gn l ( ) n ∈ (0，+ ∞ ) 和 dW(0，z) dz = W(0，z) z(1 + W(0，z)) ＞ 0， 从而有 ΔCn = lnW(0，zn ) － Cn － un gn ＞ ln Cn + un / gn l ( ) n － Cn － un gn = 0， 即 ΔCn ＞ 0． 因此，根据 0 ＜ F'(Cn ) ＜ 1 和 ΔCn = Cn +1 － Cn，故始终有 Cn +1 ＞ Cn ． 故当 n ＜ N1 时，模型(9)的解析解为 C1 (t) = － un gn + lnW 0， C1 (n) + un / gn ln exp C1 (n) + un / gn + gn ·(t － n) l ( ( ) ) n ． 由于 t ≥［t］，根据比较定理，从而得到了 C(t) 的一个上界，即 C(t) ≤ C1 (t)． 为了得到 C(t) 的一个下界，考虑如下系统: dC2 (t) dt = FCa 0Ka e －Ka(n +1) － VmaxC2 (t) KM + C2 (t) ． (15) 可将上述方程写为 dC2 (t) dt = un +1 + gn +1C2 (t) KM + C2 (t) ， t ∈［n，n + 1)， (16) 其中 un +1 = FCa 0KaKM e －Ka ·(n +1) ，gn +1 = FCa 0Ka e －Ka ·(n +1) － Vmax ． 当 FCa 0Ka ＞ Vmax，ln(FCa 0Ka /Vmax) ＞ Ka 时，同样地，可以得到 C2 (t) 的解析解为 0401 胡 晓 虎 唐 三 一

1041血管外给药的非线性房室模型解的逼近Cz(n) +C,(n) ++gn+lua+1+l.gn+1(n+a4n> N,,C (t) =/(C,(n)+K)+2ua+ (t - n),n=N2,C,(n) +C,(n)(tgatl:Ue+l+l0Ia+1La+1gn+1n<N2,其中VmKMFCK-1, ln+1N2=TVga+1此时，根据微分方程的比较定理，得到模型（9）所刻画血药浓度C（t）的上下界为C,（t）<C(t）≤C（t).采用图3中给出的参数，得到相应的数值结果.由此可见，可以利用解析求解公式C,（t）和C,（t）来近似的刻画模型C(t）的解的动态行为2.01.81.61.4121000.80.60.40.2oL102040030501/hCg = 2 mg/mL, Vmx = 1 mg/(h-mL),Km = 1.2 mg/mL, K, = 0. 1 mg/(h-mL), F = 1图3一次性血管外给药精确解与其上下界之间的逼近关系Fig. 3Approximation between the exact solution and the upper and lower boundsin the case of single extravadministration3周期血管外给药下模型（5）解的精确上下界大多数药物需要重复多次给药才能使血药浓度达到并保持在有效范围内，从而产生预期的疗效，多次重复周期性地血管外给药对于合理用药和制剂设计都是十分重要的，当药物按定时间间隔T等量多次给予，每一次给药时体内药物浓度基线发生改变，药物在体内不断地得到积累，最终达到稳定水平749]，对于周期血管外给药，本文采用脉冲微分方程2021的思想来刻画，药物在吸收室和中心室的浓度变化规律可用如下的脉冲微分方程描述：?1994-2014 China Academic Journal Electronic Publishing House.All rights reserved.http://www.cnki.net

C2 (t) = － un +1 gn +1 + ln +1 ( W － 1， C2 (n) + un +1 gn +1 ln +1 ( exp C2 (n) + un +1 gn +1 + gn +1 ·(t － n) ln +1 ) ) ， n ＞ N2， (C2 (n) + KM )2 槡 + 2un +1 (t － n)， n = N2， － un +1 gn +1 + ln +1 ( W 0， C2 (n) + un +1 gn +1 ln +1 ( exp C2 (n) + un +1 gn +1 + gn +1 ·(t － n) ln +1 ) ) ， n ＜ N2              ， 其中 N2 = 1 Ka ln FCa 0Ka V ( ) max － 1，ln +1 = － VmaxKM gn +1 ． 此时，根据微分方程的比较定理，得到模型(9) 所刻画血药浓度 C(t) 的上下界为 C2 (t) ＜ C(t) ≤ C1 (t)． 采用图 3 中给出的参数，得到相应的数值结果． 由此可见，可以利用解析求解公 式 C2 (t) 和 C1 (t) 来近似的刻画模型 C(t) 的解的动态行为． Ca 0 = 2 mg /mL，Vmax = 1 mg /(h·mL)，KM = 1． 2 mg /mL，Ka = 0． 1 mg /(h·mL)，F = 1 图 3 一次性血管外给药精确解与其上下界之间的逼近关系 Fig． 3 Approximation between the exact solution and the upper and lower bounds in the case of single extravascular administration 3 周期血管外给药下模型(5)解的精确上下界 大多数药物需要重复多次给药才能使血药浓度达到并保持在有效范围内，从而产生预期 的疗效． 多次重复周期性地血管外给药对于合理用药和制剂设计都是十分重要的． 当药物按一 定时间间隔 T 等量多次给予，每一次给药时体内药物浓度基线发生改变，药物在体内不断地得 到积累，最终达到稳定水平［17-19］． 对于周期血管外给药，本文采用脉冲微分方程［20-21］的思想来 刻画． 药物在吸收室和中心室的浓度变化规律可用如下的脉冲微分方程描述: 血管外给药的非线性房室模型解的逼近 1401

1042胡晓虎唐三dc, (t)-K,C, (t),dtt+nT,V.C(t)dc(t)= FK,C, (t)(17)dtK+C()[C,(nT*) = C,(nT) +Co,t = nT,IC(nT*)=C(nT),其中n=1,2.".对VtE（（n-1)T，nT，方程可变形为V...C(t)dc(c)= FC (n - 1) T+)K,e-k.(-(a-)ndtKM+C()进一步化简Vm..C(t)dc(t)= 入(T)e-K -KM+C()dtA(T)= FC,((n-1)T*)K,eK.(a-1)T其中同上一节一样，可以构造两个相应的子系统，并记相应的解为C.（t）和C（).假设入（T）>Vmx且In（>（T)/V）>K在区间（(n-1)T,nT)上可以得到方程（17)的上下界为C,(t) M,,C,(t) =/(C,(m)+Km)2+2u.(t-m)m=M,,C, (m) + "mEC, (m) +gm.(t-Bua+l.W01.ImgmmM,,C,(t) =/(C,（m)+Km）2+2um+（t-m）m = M2,C, (m) + "m+lC,(m) ++gmt(tgm+1gm+um+l + ImW0,Im+gm+1(m+1m<M2,其中VmKum=KM入(T)e-Km,gm=入(T)e-Km-Vm,lm=-gm?1994-2014 China Academic Journal Electronic Publishing House.All rights reserved.http://www.cnki.net

dCa(t) dt = － KaCa(t)， dC(t) dt = FKaCa(t) － VmaxC(t) KM + C(t) { ， t ≠ nT， Ca(nT+ ) = Ca(nT) + Ca 0， C(nT+ ) = C(nT { )， t = nT          ， (17) 其中 n = 1，2 .． 对 t ∈ ((n － 1)T，nT］，方程可变形为 dC(t) dt = FCa((n － 1)T+ )Ka e －Ka(t －(n －1)T) － VmaxC(t) KM + C(t) ． 进一步化简 dC(t) dt = λ(T)e －Kat － VmaxC(t) KM + C(t) ， 其中 λ(T) = FCa((n － 1)T+ )Ka eKa(n －1)T ． 同上一节一样，可以构造两个相应的子系统，并记相应的解为 C1 (t) 和 C2 (t)． 假设 λ(T) ＞ Vmax 且 ln(λ(T) /Vmax) ＞ Ka，在区间((n － 1)T，nT］上可以得到方程(17)的上下界为 C2 (t) ＜ C(t) ≤ C1 (t)， 其中 C1 (t) = － um gm + lm ( W － 1， C1 (m) + um gm lm ( exp C1 (m) + um gm + gm·(t － m) lm ) ) ， m ＞ M1， (C1 (m) + KM )2 槡 + 2um (t － m)， m = M1， － um gm + lm ( W 0， C1 (m) + um gm lm ( exp C1 (m) + um gm + gm·(t － m) lm ) ) ， m ＜ M1              ， C2 (t) = － um +1 gm +1 + lm +1 ( W － 1， C2 (m) + um +1 gm +1 lm +1 ( exp C2 (m) + um +1 gm +1 + gm +1 ·(t － m) lm +1 ) ) ， m ＞ M2， (C2 (m) + KM )2 槡 + 2um +1 (t － m)， m = M2， － um +1 gm +1 + lm +1 ( W 0， C2 (m) + um +1 gm +1 lm +1 ( exp C2 (m) + um +1 gm +1 + gm +1 ·(t － m) lm +1 ) ) ， m ＜ M2              ， 其中 um = KM λ(T)e －Kam，gm = λ(T)e －Kam － Vmax，lm = － VmaxKM gm ， 2401 胡 晓 虎 唐 三 一