医药数学模型.ppt

,单击此处编辑母版标题样式,单击此处编辑母版文本样式,第二级,第三级,第四级,第五级,*,欢迎大家！,医药数学模型,1.,药物动力学的房室模型,2.,药物服用模型,3.,传染病的传播模型,中 国 药 科 大 学,言方荣,等 编制,I know that it will happen,Because I believe in the certainty of chance.,药物动力学的房室模型,一、问题的背景与提出,在医学中有关药物的作用过程一般认为，从给药到产生药效，须经历三个主要过程,:,药剂学过程、药物动力学过程和药效动力学过程。,(,如下图所示,),通过药剂学过程，药物转化为可吸收的状态。接着，经历药物动力学过程，药物被吸收进入体循环，并在体内分布、代谢和排泄，使血液中有一定的药物浓度。当药物依度达到一定水平时，药物就可能产生效应。,图,1,药物作用的过程图,药剂学过程,剂 剂型的崩解,量 活性物质的,溶出,药物动力学过程,药物的吸收、,供吸收的药物 分布、代谢,及排泄,药剂动力学过程,药物与受体,血药浓度 相互作用 效应,药物动力学,药物动力学：研究药物、毒物及其代谢物在体内的吸收、分布、代谢和排泄的动态过程及这些过程与药理反应间的定量规律的学科分支。,对于新药研制、剂量确定、给药方案设计等药理学和临床医学的发展都具有重要的指导意义和实用价值。,房室模型,(Compartment model),房室模型是药物动力学研究上述动态过程的基本步骤之一。,药物在人体的分布过程中，可近似地把人体看成由有限个部分组成，每个部分称为一个房室。,房室具有以下特点,:,(1),每个房室有固定的容量，并且每一时刻的药物浓度都是均匀分布的。,(2),各房室间及各房室与外部环境间均可按照一定的规律进行药物交换。,房室说明,将一个机体分为几个房室，要看不同药物的吸收、分布、排泄过程的具体情况以及研究对象所要求的精度而定。,为了讨论方便，这里以,二,房室模型,为例。,即将机体分为,血液较丰富的中心室,(,包括心、肺、肾等器官,),血液较贫乏的周边室,(,如肌肉组织等,),。,二、模型假设,(1),机体分为中心室和周边室，两室容积,(,即血液体积或药物分布容积,),在过程中保持不变。,(2),药物从一室向另一室的转移速率及向体外的排泄速率与该室的血药浓度成正比。,(3),只有中心室与体外有药物交换，即药物从体外进入中心室，最后又从中心室排出体外。,(4),相对于转移和排泄的数量来说，忽略掉药物的吸收数量。,三、建模与分析,对于二房室系统来说，中心室用,l,标记，周边室用,2,标记，周围环境用,0,标记,x,i,(t),：第,i,室的药量,C,i,(t),：第,i,室的血药浓度,V,i,：,第,i,室的容积,(i=1,，,2),，,k,12,和,k,21,：,两室之间的药物转移速率系数，,k,10,：,1,室向体外排泄的速率系数，,f,0,(t):,体外给药速率,图,2,二房室系统模型示意图,血液较丰富 血液较贫乏,f,0,(t),中心室,k,12,周边室,C,1,(t),x,1,(t)C,2,(t),x,2,(t),V,1,k,21,V,2,排泄,k,10,x,i,(t)=C,i,(t)V,i,(i=1,2),药物二房室系统数学模型,根据物质平衡法则，从,t,到,t,十,t,时刻，第,i,房室的药量增加量,x,i,=x,i,(t+t),x,i,(t),应等于其余各室和环境流入,i,房室的药量之和再减去从第,i,房室流向环境和其余各室的药量之和。因此，我们有,x,1,(t+t),x,1,(t)=tk,21,x,2,(t),f,0,(t),k,12,x,1,(t),k,10,x,1,(t),x,2,(t+t),x,2,(t)=tk,12,x,1,(t),k,21,x,2,(t),两边除以,t,，再令,t0,可得,dx,1,/dt=,k,12,x,1,k,10,x,1,k,21,x,2,f,0,dx,2,/dt=k,12,x,1,k,21,x,2,由于,x,i,(t)=C,i,(t)V,i,i=1,2,整理可得到,dC,1,/dt=,(k,12,k,10,)C,1,+(V,2,/V,1,)k,21,C,2,+f,0,/V,1,dC,2,/dt=(V,1,/V,2,)k,12,C,1,k,21,C,2,(1),这是线性常系数的非齐次微分方程组，,它对应的齐次方程的通解,C,1,(t)=A,1,e,-at,B,1,e,-bt,，,C,2,(t)=A,2,e,-at,B,2,e,-bt,其中,a,，,b,满足,a+b=k,12,k,10,+k,21,，,ab=k,21,k,10,不同给药方式的讨论,下面，针对几种常见的给药方式和初始条件，具体给出方程,(1),的解。,l.,快速静脉注射模型,这种注射可理解为在初始时刻快速将剂量,D,的药物注入中央室，于是初始条件为,:,f,0,(t)=0,C,1,(0)=D/V,1,C,2,(0)=0,，,公式,(1),的解为,C,1,(t)=D(a,k,21,)/V,1,(a,b)e,-at,D(b,k,21,)/V,1,(b,a)e,-bt,，,C,2,(t)=Dk,12,/V,2,(b,a)(e,-at,e,-bt,),2.,恒速静脉注射,(,点滴,),模型,某些药物因溶解度低、刺激性过大，不宜快速注射。则采用恒速静脉注射,(,点滴,),。,初始条件为,:f,0,(t)=k,0,C,1,(0)=0,C,2,(0)=0,，,那么,(1),的解为,C,1,(t)=A,1,e,-at,B,1,e,-bt,+k,0,/k,10,V,1,，,C,2,(t)=A,2,e,-at,B,2,e,-bt,+k,12,k,0,/k,21,k,10,V,2,其中,A,1,=,k,0,(k,10,b)/V,1,k,10,(a,b),B,1,=,k,0,(a,k,10,)/V,1,k,10,(a,b),A,2,=V,1,(k,12,+k,10,a)/k,21,V,2,A,1,B,2,=V,1,(k,12,+k,10,b)/k,21,V,2,B,1,一旦在,t=T,时停止滴注,那么,C,1,(t),、,C,2,(t),在,t T,后将按指数规律衰减并趋于零。,3.,口服或肌肉注射模型,口服或肌肉注射与静脉注射的重要区别之一，就表现在药物的吸收过程，即在给药部位和药物进入中心室之间有一个将药物吸收入血液的过程。这就相当于有一个吸收室，如图,3,所示。,图,3,口服或肌肉注射时药物的吸收模式,吸收室 中心室,x,0,(t)f,0,=k,01,x,0,x,1,(t),药物进入中心室的速率,记,x,0,(t),为吸收室的药量，,k,01,为药物由吸收室进入中心室的转移速率系数，,于是有,:,dx,0,/dt=,k,01,x,0,x,0,(0)=D,其中,D,是给药量。,此时药物进入中心室的速率,f,0,(t)=Dk,01,exp(,k,01,t),那么,(1),的解为如下形式,:,C,1,(t)=Ae,-at,Be,-bt,+Eexp(-k,01,t),C,2,(t)=Fe,-at,Ge,-bt,+Hexp(-k,01,t),借助待定系数法和初始条件,C,1,(0)=C,2,(0)=0,，,可得系数,A,、,B,、,E,、,F,、,G,、,H,的值。,四、参数估计,从上可知，中心室的血药浓度,C,1,(t),与转移速率系数,k,12,、,k,21,，排除速率系数,k,10,、房室容积,V,1,、,V,2,以及输入参数,D,、,k,0,等因素有关,而房室模型恰恰是通过对,C,1,(t),的测量来确定一些对药理学及临床医学最为重要的参数，如转移速率系数,k,ij,，特别是以中心室向体外排除的速率系数,k,10,，显然这是微分方程的反问题，或称为系统辩识的问题。,下面以,快速静脉注射给药方式,为例来介绍估计诸参数的方法。,参数估计过程,从中心室采取血样，并获得血药浓度,C,1,(t,1,),C,1,(t,2,),C,1,(t,3,),，,，,C,1,(t,n,),我们将参数估计过程分成两步,:,(1),先计算,a,、,b,、,A,、,B,，其中,:,A=D(a,k,21,)/V,1,(a,b),，,B=D(b,k,21,)/V,1,(b,a),(2),再确定,k,12,、,k,21,、,k,10,。,1.,计算,a,、,b,、,A,、,B,因为,C,1,(t)=Ae,-at,Be,-bt,，不妨设,a0,为常数,k,值取决于药品的种类。,其解为,y(t)=y,0,e,-kt,，,t,0,，,T),当,t=T,时，由于经过时间间隔,T,，患者第二次服药，剂量仍为,y,0,，所以,t=T,时,y(T)=y,0,+y,0,e,-kT,则当,t,T,，,2T),时，体内药品浓度,:,y(t)=(y,0,+y,0,e,-kT,)e,-kt,t,T,，,2T),当,t=2T,时，患者第三次服药仍为,y,0,，所以,y(2T)=y,0,十,(y,0,+e,-kT,)e,-kT,=y,0,(1+e,-kT,+e,-2kT,),则,y(t)=y,0,(1+e,-kT,+e,-2kT,)e,-kt,t,2T,，,3T),以此类推，则当,t=nT,时，体内药品浓度,y(nT)=y,0,(1+e,-kT,+e,-2kT,+e,-nkT,),上式右边为一等比数列之和，求和得,y(nT)=y,0,(1-e,-(n+1)kT,)/(1-e,-kT,),则,y(t)=y,0,(1-e,-(n+1)kT,)/(1-e,-kT,),t,nT,，,(n+1)T,当,n,时,lim,n,y(nT)=y,0,/(1-e,-kT,),药物剂量的确定,如果治疗患者病情所需药物剂量水平接近,y,c,，我们近似地有,y,c,=y,0,/(1-e,-kT,),。,如果间隔时间,T,为确定量，那么剂量,y,0,可由,y,0,=(1-e,-kT,)y,c,所确定，,体内药品浓度的分布，可由图,1,说明，由图,1,可看出患者多次服药后，体内药品浓度缓慢趋于极限值,y,c,。,图,1,体内药物浓度的分布图,浓度,y,c,方法,2,方法,1,y,0,0 T 2T 3T 4T 5T 6T,t(,时间,),服药方法,2,假设 患者开始服药，就采用剂量,y,c,(,身体所需量,),且每间隔时间,T,继续服药，使体内药品浓度达到,y,c,，,若药品浓度变化仍遵循,dy/dt=,ky,的规律,那么,t0,，,T),时，体内药品浓度,y(t)=y,c,e,-kT,当,t=T,时，根据假设，需加剂量,y1,，使,y,1,十,y,c,e,-kT,=y,c,所以,y,1,=(1-e,-kT,)y,c,由上式可知，每间隔时间,T,患者服用剂量为,y,1,(,事实上,y,1,=y,0,),，药品在体内分布由,图,1,标出。,图,1,体内药物浓度的分布图,浓度,y,c,方法,2,方法,1,y,0,0 T 2T 3T 4T 5T 6T,t(,时间,),服药方法,2,的特点,优点：使药品在体内浓度从一开始就满足所需水平,缺点：以大的初始剂量作为开始，可能会使身体不适应，产生副作用,通常采用的方法：,初始时刻，令患者服药剂量为,2y,0,，,以后每间隔时间,T,，继续服药,y,0,。,进一步讨论,若患者体内药品浓度的减少遵守规律,dy/dt=,ke,y,（,k,为正常量），,仍以时间间隔,T,为周期服药，在第二次服药前，浓度,y,l,由下式给出,y,l,=,ln(kT,e,-y0,),在第三次服药前,y,2,=,lnkT(l,e,-y0,)+e,-2y0,以此类推，在第,n,十,l,次服药前,:,y,n,=,lnkT(l+e,-y0,+e,-(n-1)y0,)+e,-ny0,则当,y,n,y,c,时,T=e,-yc,(l,e,-y0,)/k,由该式我们也可以确定服药时间间隔。,传染病的传播模型,问题背景,随着卫生设施的改善、医疗水平的提高以及人类文明的不断发展，诸如霍乱、天花等曾经肆虐全球的传染性疾病已经得到有效的控制,但是在世界的某些地区，特别是贫穷的发展中国家，还不时出现传染病流行的情况。,同时，一种更为险恶的传染病,爱滋病，则在全球范围内蔓延。,2000,年底，世界各地的艾滋病患者和病毒携带者高达,3610,万人，全球死于艾滋病的总人数已达到,2180,万。,长期以来，建立数学模型来描述传染病的传播过程，分析受感染人数的变化规律，预报传染病高潮的到来等等，一直是各国关注的课题。,问题提出,被传染的人数与哪些因素有关,?,如何预报传染病高潮的到来了,?,为什么同一地区一种传染病每次流行时，被传染的人数大致不变,?,传染病的传播涉及因素很多，不可能通过一次简单的假设就能建立起完善的数学模型,这里的方法是，先做出最简单的假设，看看会得到什么结果，然后针对不合理或不完善处，逐步修改和增加假设，得到比较满意的模型。,传染病模型一,模型一假设,1),每个病人在单位时间内传染的人数为常数,k,0,2),一人得病后，经久不愈，人在传染期内不 会死亡。,模型一的建立,记时刻,t,的得病人数为,i(t),，开始时有,i,0,个传染病人，则在,t,t+t,时间内增加的病人数为,i(t,十,t),i(t)=k,0,i(t)t,于是得,di(t)/dt=k,0,i(t),i(0)=i,0,（,1,）,其解为,i(t)=i,0,模型一的评注,该结果表明，病人人数将按指数规律无限增加，当,t,时，,i(t),，显然与实际不符。,事实上，一个地区的总人数大致可视为常数,(,不考虑传染病传播时期出生和迁移的人数,),。,在传染病传播期间，一个病人单位时间能传染的人数,k,0,则是在改变的。在初期，,k,0,较大，随着病人的增多，健康者减少，被传染机会也将减少，于是,k,0,就会变小。,所以应该对模型一的假设进行修改。,传染病模型二,记时刻,t,的健康者人数为,s(t),，当总人数不变时，,k,0,应随,s(t),减少而变小。,模型二假设,1),总人数为常数,n,，且,i(t)+s(t)=n,2),单位时间内一个病人能传染的人数与当时健康者人数成正比，比例系数为,k(,传染强度,),。,3),一人得病后，经久不愈，人在传染期内不会死亡。,模型二的建立,据假设,2),，方程,(1),中的,k,0,应变为,ks(t),即,di(t)/dt=ks(t)i(t),i(0)=i,0,(2),将,s(t)=n-i(t),代入上式,得,di(t)/dt=ki(t)n-i(t),i(0)=i,0,(3),其解为,i(t)=n/1+(n/i,0,-1)e,-knt,(4),传染病曲线,i(t),t,曲线见图,1,，这个模型可用来预报传染病较快的疾病前期传染病高峰到来的时间，医学上称,di/dt,t,为传染病曲线。它反映了传染病人的变化率与时间的关系，如图,2,所示。,图,1,i(t),t,曲线,图,2,di/dt,t,传染病曲线,i(t)di/dt,n,n/2,i,0,0 t,1,t 0 t,1,t,传染病高峰时刻,由,i(t)=n/1+(n/i,0,-1)e,-knt,可得,di(t)/dt=kn,2,(n/i,0,-1)e,-knt,/1+(n/i,0,-1)e,-knt,2,令,d,2,i/dt,2,=0,得,di(t)/dt,的极大值点为,t,1,=ln(n/i,0,-1)/kn,由此可见，当传染强度,k,或,n,增加时,t,1,都将变小，即传染病高峰来得快，这与实际吻合。,此处的,k,可由经验和统计数据估计,模型二的评注,由,(4),式,:i(t)=n/1+(n/i,0,-1)e,-knt,当,t,时，,i(t)n,。这意味着最终人人都将被传染，显然这与实际不符,其原因是,假设,3),不合理，应进一步改进。,传染病模型三,有些传染病如伤风、痢疾等愈后免疫力很低，可以假定无免疫性，于是病人被治愈后变成健康者，健康者还可以被感染再变成病人,模型三假设,1),总人数为常数,n,，且,i(t)+s(t)=n,2),单位时间内一个病人能传染的人数与当时健康者人数成正比，比例系数为,k(,传染强度,),。,3),病人每天被治愈的占病人总数的比例为,，,称为日治愈率。病人治愈后成为仍可被感染的健康者。,考虑到条件,3,，模型二的,(2),式应改为,di(t)/dt=ks(t)i(t)-,i(t),i(0)=i,0,(5),将,s(t)=n-i(t),代入上式,得,di(t)/dt=ki(t)n-i(t)-,i(t),i(0)=i,0,(6),其解为,i(t)=(,kt+n/i,0,),-1,(,k=,),和,i(t)=,k/(k-)+(n/i,0,-k/(k-)exp(-(k-)t),-1,(,k,),i,0,i,0,接触数,=1,阈值,感染期内,有效接触感染的健康者人数不超过病人数,i,0,0,t,i,1,1-1/,i,0,t,1,di,/,dt,0,传染病模型四,大多数传染病如天花、流感、肝炎、麻疹等治愈后均有很强的免疫力，所以病愈的人既非健康者，也非病人,(,已感染者,),，他们已经退出传染系统，称为移出者,我们将人员分成,3,类,:,第一类是传染者,(i),类,;,第二类是健康者,(,易受传染者,)(s),类,;,第三类是移出者,(r),类，,(,包括患病死去的人、痊愈后具有长期免疫力的人以及在痊愈以前被隔离起来的人,),令,i(t),，,s(t),和,r(t),分别表示在时刻,t,的第一、二、三类的人数。,模型四的假设,1),总人数为常数,n,，则,i(t),十,s(t),十,r(t)=n (5),2),单位时间内一个病人能传染的人数与当时健康者人数成正比，比例系数为,k(,传染强度,),。,3),单位时间内病愈免疫的人数与当时的病人人数成正比，比例系数为,a(,恢复系数,),。,模型四的建立,由假设,3),，有,dr/dt=ai(t)(6),由于引入了,r(t),，则模型,2,的方程,(2),式应改为,di(t)/dt=ks(t)i(t)-dr/dt (7),上式表示单位时间内病人人数的增加应等于被传染的人数减去病愈的人数。,从,(5)(7),式中消去,di/dt,，并设,r(0)=0,，,s(0)=s,0,。,得,ds/dt=-ks(t)i(t),dr/dt=ai(t),r(0)=0,s(0)=s,0,i(t),十,s(t),十,r(t)=n (8),模型四的求解,记,p=a/k,，称为特征指标，对同一地区，同一种传染病,p,是常数。由方程,(8),可得,ds/dt=-(1/p)s(t)(dr/dt),其解为,s(t)=s,0,e,-(1/p)r(t),s,0,1-(1/p)r(t)+(1/2p,2,)r,2,(t),于是,dr/dt=ai(t)=an-r(t)-s(t),a(n-s,0,)+(s,0,/p-1)r(t)-(s,0,/2p,2,)r,2,(t),再解上式,经过整理简化得,lim,t,r(t)(2p,2,/s,0,)(s,0,/p,1),上式表明,:,当,s,0,很小,(p),时，才存在传染问题。,今设,s,0,=p,十,h(h0,只有这时才有传染问题,),，且,hp,则,lim,t,r(t)2p,2,/(p+h)(h/p)2h,因,lim r(t),表示最终免疫的人数，也即最终痊愈的人数。显然，这也就是传染病传染的人数，于是由上式可得下列结论,模型四的评注,(1),对于同一地区，同一传染病所传染的人数大体上是常数,2h,，该结论与统计结果一致。,(2),当恢复系数,a,变大，传染强度,k,变小时，,p=a/k,变大，相应,h,变小,从而传染病传染的人数,2h,也就减小。,(3),由于,s,0,可视为总人数，所以传染病传染的人数是,p,以上人数,h,的两倍。,