基于MATLAB的药代动力学实验数据处理及参数估计.docxVIP

基于MATLAB的药代动力学实验数据处理及参数估计 【关键词】 药代动力学；,,实验数据处理；,,MATLAB；,,评价,,,,  摘 要： 介绍了MATLAB在药代动力学中实验数据处理的基本 方法 并对利用MATLAB进行双指数模型拟合与常用EXCEL软件进行双指数模型拟合以及参数求解的方便性、准确性作出了评价。  关键词： 药代动力学； 实验数据处理； MATLAB； 评价  药代动力学是临床药 理学 、药剂学、毒理学等 科学 的专业基础，对于新药设计、剂型改革、合理用药等方面都有直接指导意义。如何对实验数据进行正确的处理以及选择正确的处理方式是准确表达药物动力学过程的关键，以便对临床用药和药物设计等方面提供可靠的依据。 目前 大多数的药物的药代动力学过程基本符合二室模型，其原理是将人体分为两个房室即中央室和周边室，以中央室的消除为主不考虑周边室的消除，并假设药物的吸收交换和消除都是一级过程。  　1 模型及参数  二室模型的药物动力学［1］方程为：C=Ae-αt+Be-βt其中C表示血药浓度，t表示时间，A、B、α、β表示各项常系数。二室模型的其它参数 计算 公式分别为：① 药物的分布相半衰期：t1 α= α② 药物的消除半衰期：t1 β= β③ 药物从周边室向中央室的分布速率常数： k21=Aβ+Bα A+B④ 药物中央室消除速率常数： k10=αβ k21⑤ 药物的中央室向周边室的分布速率常数： k12=α+β-k10-k21⑥ 药物的中央室分布容积： V1=D0 A+B⑦ 药物的消除率： CL=V1k10以上是药代动力学过程中的各个参数的计算公式，求得模型中的各项常系数即可取得各个参数值。  　2 MATLAB处理数据［2］  以磺胺噻唑(D0=250mg/kg)二室模型动力学参数测定为例，利用MATLAB对实验数据浓度时间作非线性曲线拟合。实验测得的数据为： 时间(min) 10 0 0 0 0 120 150 180 10浓度(mg/L) 6 59 70 00 39 180 160 145 132 123 11 MATLAB中引入了两个命令来进行曲线的拟合定义非线性方程的命令(inline)和拟合曲线方程命令(nlinfit)，其进行曲线拟合的程序为：% 输入数据 t=［ 120 150 18010］； c=［80 160 ］；% x0为给A，alpha, B, beta初值 x0=［0 0 0 0］；% 定义模型函数f=inline(’x(1)*exp(-x(2)*t)+x(3)*exp(-x(4)*t)’，’x’，’t’)；% 进行非线性曲线拟合 ［b,r,j］=nlinfit(t,c,f,x0)；b=b’;% 得到各个系数值 A=b(1)，alpha=b(2),B=b(3),beta=b(4)通过该程序得到了各个系数的值为： A= α= B= β=将各系数带入方程就直接得到了磺胺噻唑药代动力学模型方程为：C=+将其模型函数与实验点利用MATLAB作图，见图1，其调用的命令为plot与ezplot。从图形看，拟合曲线和实际数据点吻合很好，此外引入一个可决系数［3］R2值作进一步说明，R2一般大于说明效果较好，越接近于1就越好。其计算公式为：R2=1-∑n i=1(Ci-i)／∑n i=1(Ci-)通过MATLAB按照此计算公式得出可决系数R2为，拟合效果较好。 图1 磺胺噻唑药代动力学　略  　3 MATLAB与EXCEL处理比较  在用EXCEL处理二室模型的数据时，常把模型分为分布相和消除相，采用分段拟合来进行，由于指数可以转化为对数而成直线，以此分两段直线进行拟合而得各个系数以及各个参数值，具体处理过程略。由于其不是把整个模型连续起来且前后选点的个数均受到实验、自身条件等因素不能正确选择，处理时可能带来误差、 计算 的麻烦及处理的不准确性，另外EXCEL对整个拟合不能做出很好的评价。用EXCEL处理得到的磺胺噻唑药代动力学模型方程为：C=+利用EXCEL拟合曲线与实验点的可决系数R2为 图2 略  饮酒动力学模型其拟合的模型方程为：C=+可决系数R2为09932，效果较好。此外，MATLAB还介绍了更多的非线性拟合方法如命令(curvefit)、(lsqcurvefit)、(lsqnonlin)等。  　4 结论  本 研究 对利用MATLAB处理数据和EXCEL处理数据进行了