用吸收马尔可夫链估算药物在体内的平均转运时间.pdfVIP

用吸收马尔可夫链估算药物在体内的平均 转运时间 丁勇 南京医科大学数学教研室，江苏南京 （210029） 摘 要：本文提出了导出矩阵的概念，解决了药物转运速率常数转化为转移概率的问 题，从而可用吸收马尔可夫链对药物在体内的转运过程进行研究。推导了用于估算药 物在体内各部位平均转运时间的基本矩阵。 关键词：导出矩阵，吸收马尔可夫链，基本矩阵，药物动力学 中图分类号：O211.6 Q332 0 引言 [1~2] 近年来，随机过程在生命科学中的研究得到了越来越多的关注 。本文对吸收 马尔可夫链在药动学房室模型中的应用进行了探讨，提出了导出矩阵的概念，解决了 药物转运速率常数转化为转移概率的问题，为估算药物在体内转运过程的吸收、分布 和消除的平均时间提供了新的方法。 1 导出矩阵 [ ] 定义 设有矩阵A aij m ×n 和任一常数δ ，对整数1≤l ≤m 构造一个新矩阵 ~ ~ ~ [ ] A aij m ×n ，其中： al ,j δ al ,j j 1, 2, L, n; j ≠l n n n n ~ ~ a a − a a −δ a l l ∑ l j ∑ l j ∑ l j ∑ l j , , , , , j 1 j 1 j 1 j 1 j ≠l j ≠l ~ i 1, 2, L, l −1, l +1, L, n j 1, 2, L, n ai ,j ai ,j ； ~ ~ 称A 是由A 、δ和l 导出的矩阵，记为A A[δ | l ] ；称此变换为导出变换。 ~ 导出矩阵A 的元素除第l 行外，都与A 相同；第l 行元素除主对角线元素外，与 A 只差一个常数因子δ ；各行元素之和与A 相同。 对 多 个 δ 和 l ， 1≤i ≤t ， 可 对 A 进 行 多 次 导 出 ， 记 为 i i ~ A A[δ | l ][δ | l ]L[δ | l ] 。 1 1 2 2 t t 容易证明，导出矩阵有如下的一些简单性质： −1 （1） 对于任意矩阵A ，有A[1| l ] A ；对δ ≠0 ，有A[δ | l ][δ | l ] A 。 （2 ） 对于任意对角阵A ，有A[δ | l ] A （3 ）对于任意n 阶方阵A ，有A[0 | 1][0 | 2]L[0 | n] 为对角阵，对角线元素为 -1- n ∑aij i 1,2, L, n 。 j 1 （4