资源描述：
《基于matlab的药代动力学实验数据处理及参数估计 》由会员上传分享，免费在线阅读，更多相关内容在学术论文-天天文库。
1、基于MATLAB的药代动力学实验数据处理及参数估计【关键词】药代动力学；,,实验数据处理；,,MATLAB；,,评价,,,, 摘要：介绍了MATLAB在药代动力学中实验数据处理的基本方法并对利用MATLAB进行双指数模型拟合与常用EXCEL软件进行双指数模型拟合以及参数求解的方便性、准确性作出了评价。 关键词：药代动力学；实验数据处理；MATLAB；评价 药代动力学是临床药理学、药剂学、毒理学等科学的专业基础，对于新药设计、剂型改革、合理用药等方面都有直接指导意义。如何对实验数据进行正确的处理以及选择正确的处理方式是准确表达药物动力学过程的关键，以便对临床用药和药物设计等方
2、面提供可靠的依据。目前大多数的药物的药代动力学过程基本符合二室模型，其原理是将人体分为两个房室即中央室和周边室，以中央室的消除为主不考虑周边室的消除，并假设药物的吸收交换和消除都是一级过程。 1模型及参数 二室模型的药物动力学〔1〕方程为：C=Ae-αt+Be-βt其中C表示血药浓度，t表示时间，A、B、α、β表示各项常系数。二室模型的其它参数计算公式分别为：①药物的分布相半衰期：t12α=0.693α②药物的消除半衰期：t12β=0.693β③药物从周边室向中央室的分布速率常数：k21=Aβ+BαA+B④药物中央室消除速率常数：k10=αβk21⑤药物的中央室向周边室的分布
3、速率常数：k12=α+β-k10-k21⑥药物的中央室分布容积：V1=D0A+B⑦药物的消除率：CL=V1k10以上是药代动力学过程中的各个参数的计算公式，求得模型中的各项常系数即可取得各个参数值。 2MATLAB处理数据〔2〕 以磺胺噻唑(D0=250mg/kg)二室模型动力学参数测定为例，利用MATLAB对实验数据浓度时间作非线性曲线拟合。实验测得的数据为：时间(min)251020306090120150180210浓度(mg/L)766659470300239180160145132123112MATLAB中引入了两个命令来进行曲线的拟合定义非线性方程的命令(inli
4、ne)和拟合曲线方程命令(nlinfit)，其进行曲线拟合的程序为：%输入数据t=〔251020306090120150180210〕；c=〔766659470300239180160145132123112〕；%x0为给A，alpha,B,beta初值x0=〔0000〕；%定义模型函数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)通过该程序得到了各个
5、系数的值为：A=683.1323α=0.0908B=208.0097β=0.0030将各系数带入方程就直接得到了磺胺噻唑药代动力学模型方程为：C=683.1323e-0.0908t+208.0068e-0.003t将其模型函数与实验点利用MATLAB作图，见图1，其调用的命令为plot与ezplot。从图形看，拟合曲线和实际数据点吻合很好，此外引入一个可决系数〔3〕R2值作进一步说明，R2一般大于0.99说明效果较好，越接近于1就越好。其计算公式为：R2=1-∑ni=1(Ci-i)／∑ni=1(Ci-)通过MATLAB按照此计算公式得出可决系数R2为099877，拟合效果较好。
6、 图1磺胺噻唑药代动力学 略 3MATLAB与EXCEL处理比较 在用EXCEL处理二室模型的数据时，常把模型分为分布相和消除相，采用分段拟合来进行，由于指数可以转化为对数而成直线，以此分两段直线进行拟合而得各个系数以及各个参数值，具体处理过程略。由于其不是把整个模型连续起来且前后选点的个数均受到实验、自身条件等因素不能正确选择，处理时可能带来误差、计算的麻烦及处理的不准确性，另外EXCEL对整个拟合不能做出很好的评价。用EXCEL处理得到的磺胺噻唑药代动力学模型方程为：C=593.3711e-0.077752t+206.7486e-0.00291t利用EXCEL拟合曲线与实
7、验点的可决系数R2为0.98860.99877，说明拟合效果不是很好，较MATLAB差，不能更好的表达磺胺噻唑的药代动力学过程，因此EXCEL处理的其他参数也不能更准确的计算。将利用MATLAB和EXCEL处理的参数分别带入前面的参数计算公式得到各个参数的值列于下表：参数t12αt12βk21k10k12V1CLMATLAB7.6332231.000.02350.01160.05870.28053.5EXCEL8.9130238.140.02220.01020.04830.312
