血管外給藥單室模型的最優(yōu)控制

劉愛(ài)紅, 楊 光

(沈陽(yáng)師范大學(xué) 數(shù)學(xué)與系統(tǒng)科學(xué)學(xué)院, 沈陽(yáng) 110034)

?

血管外給藥單室模型的最優(yōu)控制

劉愛(ài)紅, 楊 光

(沈陽(yáng)師范大學(xué) 數(shù)學(xué)與系統(tǒng)科學(xué)學(xué)院, 沈陽(yáng) 110034)

考慮血管外給藥單室模型的脈沖最優(yōu)控制問(wèn)題。首先,利用脈沖微分方程和常微分方程重新改寫(xiě)傳統(tǒng)的血管外給藥微分方程模型,從而既突出給藥這一脈沖過(guò)程,又將一次給藥模型推廣為周期性多次給藥模型;然后,求解新建立的脈沖微分方程,進(jìn)而確定每次給藥后吸收部位藥量和體內(nèi)藥量的表達(dá)式;接下來(lái),利用一階和二階導(dǎo)數(shù)推導(dǎo)出每次給藥后藥量的變化規(guī)律和藥量與給藥次數(shù)之間的內(nèi)在聯(lián)系,并給出一個(gè)療程內(nèi)體內(nèi)藥量的最高值與最低值;再次,構(gòu)建描述體內(nèi)藥量積蓄程度的目標(biāo)函數(shù),利用最優(yōu)化方法,確定最佳給藥次數(shù)、時(shí)間間隔和劑量,使得藥物在體內(nèi)的積蓄最少,對(duì)人體的傷害最小;最后,Matlab仿真結(jié)果表明該方法有效、可行。

單室; 血管外給藥; 藥量; 給藥方案; 脈沖微分方程; 最優(yōu)控制

藥物動(dòng)力學(xué)是一門(mén)藥效學(xué)與數(shù)學(xué)有效結(jié)合的交叉學(xué)科,模型的建立與給藥方案的設(shè)計(jì)一直是其研究的主要內(nèi)容之一,到目前為止,有大量的相關(guān)文獻(xiàn)。對(duì)于單室血管外給藥模型,在模型建立方面,大多數(shù)都建立了如下的微分方程模型:

其中:Xa為吸收部位可吸收的藥量;X為體內(nèi)的藥量。這個(gè)模型有一定的局限性。首先,這個(gè)模型不能直觀(guān)的反映出周期性給藥方式下的藥量變化;其次,此模型僅體現(xiàn)了給藥后藥量變化的連續(xù)性,卻未能體現(xiàn)給藥這一行為的特點(diǎn)。在給藥方案設(shè)計(jì)方面,少有文獻(xiàn)將最優(yōu)化方法與積蓄藥量對(duì)肝腎等器官的傷害程度結(jié)合在一起共同去確定最佳給藥方案。因此,本文將在模型建立和給藥方案設(shè)計(jì)這兩方面進(jìn)行大膽的嘗試研究。

首先,本文考慮周期性給藥模式,運(yùn)用脈沖微分方程改寫(xiě)傳統(tǒng)的單室血管外給藥微分方程模型;其次,計(jì)算出一個(gè)療程內(nèi)體內(nèi)的最高、最低藥量;然后,運(yùn)用最優(yōu)化理論,確定使得體內(nèi)積蓄藥量最少,對(duì)人體害程度最小的最佳給藥間隔與給藥量;最后,通過(guò)仿真模擬證明本文建立模型的有效性。

1 血管外給藥單室模型的建立

本文給出如下假設(shè):

[H1] 給藥是一個(gè)脈沖過(guò)程,首次給藥量為x0,維持給藥量為I;

[H2] 給藥的時(shí)間間隔恒為τ,1 d內(nèi)的給藥次數(shù)為n,1 d的給藥次數(shù)不會(huì)超過(guò)5,否則太過(guò)繁瑣,即n≤5且n∈N,nτ=24。

[H3] 藥量在人體內(nèi)積蓄程度用體內(nèi)藥量在給藥時(shí)間內(nèi)的積分來(lái)表示;

[H4] 吸收部位可吸收藥量的變化為τ-周期解,F為生物利用度,從第2次給藥開(kāi)始,每次給藥后吸收部位藥量均為首次給藥量,則

Xa((i-1)τ-0)+FI=Fx0,i=2,3,…。

[H5] 中毒臨界值hmax與藥效臨界值hmin均是已知的常數(shù)。

[H6] 本文考慮短期給藥且一個(gè)療程為6 d,即144 h。

依據(jù)上述假設(shè),可得新模型:

由式(2)前4項(xiàng)可得

由[H4]和式(3)可得

注1 根據(jù)[H4]及式(3),易得第i次給藥后吸收部位可吸收藥量的表達(dá)式為

2 最佳給藥策略的確定

最佳給藥方案使體內(nèi)藥量在一個(gè)療程內(nèi)積蓄最小,并且從第2次給藥時(shí)刻起體內(nèi)藥量均在藥效臨界值與中毒臨界值之間。首先確定每次給藥后體內(nèi)藥量的表達(dá)式,然后根據(jù)最優(yōu)化方法給出最佳給藥方案的最優(yōu)化模型。

定理1 對(duì)于血管外給藥模型(2),在一個(gè)療程,即在[0,144]內(nèi),第i次給藥后體內(nèi)藥量的表達(dá)式為

其中:(i-1)τ≤t≤iτ,i=1,2,…,6n;nτ=24。

證明 第1次給藥后體內(nèi)藥量的表達(dá)式(4)式已給出。第2次給藥后,體內(nèi)藥量的初值為

等價(jià)于求解如下微分方程:

即

同理可得第i次給藥后體內(nèi)藥量的表達(dá)式為式(7)。

定理2 在整個(gè)療程內(nèi),體內(nèi)最高藥量為最后一次給藥期間的最高藥量,且為

證明 令θ=t-(i-1)τ, (i-1)τ≤t≤iτ,i=1,2,…,6n,nτ=24, 0≤θ≤τ。

則第i次給藥后體內(nèi)藥量的表達(dá)式(7)變?yōu)?/p>

由于

因此,?θ, 0≤θ≤τ,i次給藥后體內(nèi)藥量隨著i增大而增大。

θ*

由于

所以當(dāng)θ滿(mǎn)足式(11)時(shí),此時(shí)刻藥量為第i次給藥后體內(nèi)最高藥量。綜上可知,最后一次給藥后體內(nèi)藥量X6n(t)的最高值即為是整個(gè)療程體內(nèi)最高藥量,即為(9)。

注2 為了使藥物達(dá)到藥效且不對(duì)人體造成中毒,要求第2次給藥在τ時(shí)刻的藥量與體內(nèi)最高藥量均在藥效臨界值與中毒臨界值之間,即滿(mǎn)足下列條件即可

hmin

定理3 在一個(gè)療程內(nèi),體內(nèi)藥物的積蓄f為

證明 由[H3]、式(10)易得。

且最佳給藥間隔τ*及維持給藥量I*為

3 仿真與討論

設(shè)有某種藥物ka=1h-1,k=0.1h-1,hmin=5 000 μg,hmax=50 000μg,F=0.8,將數(shù)據(jù)帶入上述最優(yōu)化模型中,運(yùn)用Matlab軟件求解可得

此時(shí),一個(gè)療程內(nèi)吸收部位與體內(nèi)藥量的變化圖像如圖1、圖2所示。圖1說(shuō)明給藥確實(shí)是一個(gè)脈沖過(guò)程,且吸收部位藥量的變化是周期的。此外,如此的給藥方式保證了體內(nèi)藥量積蓄最小,且每次給藥后的藥量都在藥效之上,中毒臨界值之下,這符合最佳給藥策略的期望,模型有效可行。

圖1 吸收部位藥量變化Fig.1 Variety of dose in absorption site

圖2 體內(nèi)藥量變化Fig.2 Variety of dose in vivo

4 結(jié) 論

本文在周期性給藥的前提下,結(jié)合脈沖微分方程和常微分方程建立了新的血管外給藥單室模型,突出了給藥的脈沖過(guò)程,將一次給藥模型拓展為周期性多次給藥模型。另外,考慮了藥量在體內(nèi)的積蓄程度,運(yùn)用最優(yōu)化理論,構(gòu)造了給藥方案設(shè)計(jì)的最優(yōu)控制模型,確定的最佳給藥間隔、給藥量使得藥物在體內(nèi)的積蓄最小,對(duì)人體的傷還最小。Matlab仿真結(jié)果證明了方案的科學(xué)可行性。

[ 1 ]梁文權(quán). 生物藥劑學(xué)與藥物動(dòng)力學(xué)[M]. 北京:人民衛(wèi)生出版社, 2007.

[ 2 ]楊光,劉愛(ài)紅,王爽. 單室模型的最佳給藥方案[J]. 沈陽(yáng)師范大學(xué)學(xué)報(bào)(自然科學(xué)版), 2015,33(3):337-340.

[ 3 ]楊光. 害蟲(chóng)綜合治理模型與最優(yōu)控制策略[J]. 沈陽(yáng)師范大學(xué)學(xué)報(bào)(自然科學(xué)版), 2011,29(4):482-485.

[ 4 ]薛定宇,陳陽(yáng)泉. 基于MATLAB/Simulink的系統(tǒng)仿真技術(shù)與應(yīng)用[M]. 北京:清華大學(xué)出版社, 2002.

[ 5 ]段春林. 臨床給藥方案探討[J]. 中國(guó)現(xiàn)代藥物應(yīng)用, 2010,4(2):163.

[ 6 ]汝玲. 臨床藥師參與個(gè)體化給藥方案設(shè)計(jì)的臨床實(shí)踐與案例分析[J]. 中國(guó)醫(yī)院用藥評(píng)價(jià)與分析, 2013,13(3):273-276.

[ 7 ]楊進(jìn)波. 創(chuàng)新性藥物臨床試驗(yàn)劑量和給藥方案的探索與確定[J]. 中國(guó)臨床藥理學(xué)與治療學(xué), 2008,13(8):841-846.

[ 8 ]劉昌孝. 我國(guó)藥物動(dòng)力學(xué)研究50年發(fā)展概述[J]. 天津中醫(yī)藥大學(xué)學(xué)報(bào), 2008,27(3):127-134.

[ 9 ]李琳. 淺談給藥方案設(shè)計(jì)[J]. 中國(guó)繼續(xù)醫(yī)學(xué)教育, 2015,7(3):218.

[10]楊慶之. 最優(yōu)化方法[M]. 北京:科學(xué)出版社, 2015.

Optimal control for one-compartment extravascular administration model

LIU Aihong, YANG Guang

(College of Mathematics and Systems Science, Shenyang Normal University, Shenyang 110034, China)

This paper considers the issue of impulsive optimal control for one-compartment extravascular administration model. Initially, impulsive differential equations and ordinary differential equations are used to rewrite the traditional extravascular administration differential equation model so as to highlight the impulsive process of dosing and promote the once dosing model to periodic dosing model. Furthemore, expressions of doses in vivo and absorption site after each dosing are determined by solving the fresh impulsive differential equation. The varying pattern of dose and the inner relationship between dose and number of dosing are obtained by the first and second order derivates, and the highest and lowest doses of drugs during a course of treatment are given. Additionally, a target function which describes the accumulation extent of drugs in vivo is constructed, the best dosing numbers, time interval and doses are obtained by optimization methods, which makes the accumulation extent of drugs in vivo minimum as well as brings the least harm to human body. Finally, the simulation result operated by Matlab shows the validity and feasibility of the proposed method.

one-compartment; extravascular administration; dose; dosing strategy; impulsive differential equations; optimal control

2016-03-02。

遼寧省科技廳自然科學(xué)基金資助項(xiàng)目(2014020120)。

劉愛(ài)紅(1990-),女,遼寧北票人,沈陽(yáng)師范大學(xué)碩士研究生; 楊 光(1964-),女,遼寧撫順人,沈陽(yáng)師范大學(xué)教授,博士。

1673-5862(2016)03-0287-04

O232

A

10.3969/ j.issn.1673-5862.2016.03.007