模拟退火算法在药物动力学参数反演中的应用.pdfVIP

第24卷第2期 四川理工学院学报(自然科学版) V01．24No．2 of of Science JournalSichuan Edition) Apr．201l 2011年4月 UniversityScienceEngineering(Natural 文章编号：1673—1549(2011)02-0175-03 模拟退火算法在药物动力学参数反演中的应用 钟太勇1，彭先萌2 (1．郧阳师范高等专科学校数学系，湖北十堰442000；2．十堰职业技术学院公共课部，湖北十堰442000) 摘要I介绍了模拟退火算法的原理和求解方法，并将其用于药物动力学参数的反演，在利用残敷 法求得解的基础上确定了各参数的变化区间，在Matlab语言环境下实现了该算法，并与残数法、IGS—M 法的结果相比较，数值结果表明，模拟退火算法简单、准确，能更好地实现最优拟合。 关键词：模拟退火算法；残数法；反演；IGS—M；Matlab 中图分类号：嘲Ol 文献标识码：A 药物动力学研究内容之一是以建立数学模型的方 间作回归直线，最终依次求出所有的参数。该方法又称 法来研究机体内药物的吸收、分布和消除等过程随时间 变化的规律，这种规律我们可以通过函数来表示，其中 ping)等。 关键的问题是要求出药物动力学参数，求参数的方法就 2模拟退火算法 是驰豫曲线拟合的问题，当机体给药后，由于吸收、分 布、生物转化、排泄等动态速率过程，体内的药物浓度时 2．1算法原理 刻处于动态变化之中，其规律是药物浓度C随时间t的 模拟退火算法¨朋1的基本思想是从一给定解开始， 延长而减少，把C，t的关系式描绘成图就是驰豫曲线。 从邻域中随机产生另一个解，接受Metropolis准则允许 在药动学房室模型静脉快速注射给药后主要采用 目标函数在有限范围内变坏，它由一控制参数t决定，其 三种形式的负指数模型来描述，即单室模型C=Ae一、作用类似于物理过程中的温度丁，对于控制参数的每一 取值，算法持续进行“产生一判断一接受或舍去”的迭代 双室模型C=Ael+母和三室模型C=Ael+踟谓 过程，对应着固体在某一恒定温度下的趋于热平衡的过 +R～。笔者主要讨论第三房室模型C=Ae“+＆母+ 程，当控制参数逐渐减小并趋于0时，系统越来越趋于 Pe巾，其中c为t时刻的血药浓度，t为时间，A、a、B、 平衡态，最后系统状态对应于优化问题的全局最优解， JB、P、P为复杂参数。结合残数法，在Matlab语言环境 下采用具有最优保存策略的退火算法对药物动力学参 该过程也称为冷却过程，由于固体退火必须缓慢降温， 数进行反演。 才能使得固体在每一温度下都达到热平衡，最终趋于平 衡状态，因此控制参数t经缓慢衰减，才能确保模拟退火 1残数法 算法最终优化问题的整体最优解。 残数法是把一条多项式指数曲线分解成若干指数 2．2算法具体步骤 项成分的常用药物动力学方法，即当时间足够大时，其 (1)给定模型每一个参数变化范围，在这个范围内 中速率常数较小的指数项t趋近于零，从而略去，以至剩 随机选择一个初始模型m。，并计算相应的目标函数值 下速率常数较小的单项指数式，经过对数转化后做线性