МАТЕМАТИЧЕСКОЕ МОДЕЛИРОВАНИЕ И ОБРАБОТКА ДАННЫХ
Краткое сообщение УДК 519.87
Б01: 10.18101/2304-5728-2023-2-49-54
ОПТИМИЗАЦИОННАЯ ФОРМАЛИЗАЦИЯ ПРОЦЕДУРЫ ЛЕЧЕНИЯ ДИАБЕТА
© Аксенюшкина Елена Владимировна
кандидат физико-математических наук, доцент, Байкальский государственный университет Россия, 664003, г. Иркутск, ул. Ленина, 11 [email protected]
Аннотация. Приводится формализация задачи, связанной с терапией диабета на основе последовательных инъекций инсулина с целью минимизации квадра -тичного отклонения текущей концентрации глюкозы от нормального уровня. Моделирование этой процедуры проводится в непрерывно-дискретном варианте в терминологии оптимального управления. Фазовая система описывается дифференциальными уравнениями, допустимые управления, фазовые ограничения и целевой функционал формализуются в дискретном времени на заданной сетке значений. В результате построена специальная задача квадратичного программирования, которая допускает эффективное решение за конечное число итераций.
Ключевые слова: процедура лечения диабета, оптимизация в непрерывно-дискретном времени, задача квадратичного программирования.
Для цитирования
Аксенюшкина Е. В. Оптимизационная формализация процедуры лечения диабета // Вестник Бурятского государственного университета. Математика, информатика. 2023. № 2. С. 49-54.
Введение
Возьмем за основу модель лечения сахарного диабета, представленную в [1, 2] в форме линейно-квадратичной задачи оптимального управления, но без необходимых ограничений на управление и фазовые переменные.
В статье в рамках подхода из [3] проводится построение более адекватной задачи, связанной с терапией диабета через последовательные инъекции инсулина с целью минимизации среднеквадратичного отклонения текущей концентрации глюкозы в крови пациента от нормального уровня на заданном промежутке времени.
Моделирование этой ситуации реализуется в непрерывно-дискретном варианте: фазовая система описывается дифференциальными уравнениями (непрерывное время), допустимые управления, фазовые ограничения и целевой функционал формализуются в дискретном времени по заданной сетке значений (ступенчатые функции, условия неотрицательности в узловых точках и сумма квадратов отклонений). В результате построена
специальная задача квадратичного программирования, которая допускает численное решение за конечное число итераций.
1 Допустимые управления и соответствующие траектории
Пусть в момент времени V е [0, Т]:
и (V) — инъекция инсулина;
х(У) — концентрация глюкозы в крови, причем х - нормальный уровень глюкозы;
хо(/) — концентрация гормонов.
Укажем естественные ограничения на переменные
и^) е [0, и+], х(0 > 0, х > 0, х0(0 > 0.
Дифференциальная связь между функциями и(1), Х0(0, х(?) для V е [0, Т] описывается линейной системой [1, 2]:
Х0(0 = - СХ0(0 + и(0, Х0(0) = 0, (1)
х^) = - ах($) - Ъх0^), х(0) = х. с положительными коэффициентами а, Ъ, с и начальным состоянием х > 0.
Проведем дискретизацию отрезка времени [0, Т] с помощью заданного набора точек (моментов инъекций)
0 = {0 < Ч < ... < -1 < _ Т .
Выделим промежутки Tj _ [Vу-1, tj), ] _ 1, т и определим соответствующие характеристические функции
_Г 1, V е Т], С _{ 0, V е [0, Т]\^.
Введем набор значений yj, j _ 1, т (объем инъекций в момент tj-1) с условием yj е [0, и+], у > 0 и сформируем вектор у _ (>1,..., ут) (суммарная инъекция за период терапии).
Образуем допустимое управление для системы (1) (стратегия терапии)
т
и^,у) Уjа (V), V е [0, Т].
j _1
Это кусочно-постоянная (ступенчатая) функция со значениями и(^-1, >) _ У]-, j _ 1, т .
Рассмотрим решение х0^,у), х^,у) системы (1), соответствующее управлению и (V, у), V е [0, Т]. Для j _ 1, т определим опорные траектории х0 j (), х j (V) как решения задач Коши:
=-сх0 +Х] ((), х0(0) = 0, х = -ах -Ьх0, х(0) = 0. Тогда справедливо следующее представление для решения системы (1)
т
х0(Ь,У) = 2х0-(Ь)У- , (2)
-=1
т
-аЬ
х(Ь, У) = 2 X] (Ь)У] + хе-а, (3)
]=1
которое проверяется непосредственным дифференцированием.
2 Фазовое ограничение и целевой функционал
Конкретизируем ситуацию с опорными траекториями. Функция
Х0 - (Ь) представляется по формуле:
(
0, Ь е [0, Ь--1),
Х0 ] (Ь) =
1е~сг{есг -е^-1), Ь е [г^, -
С
1 е~С (еС ] - еС ]-1), Ье [ Ь,, Т].
СЬ,-
е (е -
V с
Отсюда Х0- > 0, Ь е (Ь--1, Т]. Следовательно, на основании (2) выполняется условие положительности для всех траекторий х$(Ь, у), у е [0, и+], У1 > 0, Ь е (0, Т].
Опорная траектория х- (Ь) определяется уравнением (задача Коши):
X = - ах - Ьх0- (Ь), х(0) = 0 .
Отсюда получаем нулевой участок: х-(Ь) = 0, Ь е [0, Ь -, ] = 1, т . В
частности, для узловых точек , к = 1, т -1 имеем х- (^) = 0, - > к . Кроме того, по формуле Коши:
Ьк _
х- (Ьк) = -Ье~а1к |еаЬх0- (Ь)<#, к = 1, т . 0
Это значит, что для - < к х- ) < 0 .
Таким образом, согласно представлению (3) справедлива формула:
к _
х(к, У) = 2 х- (Ьк)У- + хе~а1к , к = 1, т, -=1
причем х- (Ьк) < 0 .
Это выражение допускает возможность отрицательных значений для , У) (за счет суммы), что противоречит содержательному смыслу этих значений (концентрация глюкозы).
Ограничения неотрицательности траектории х(^, у), ? е [0, Т] будем выдерживать в узловых точках с помощью переменных у^,..., ут:
х^к, у) > 0 , к = 1, т . Это система линейных неравенств с нижней треугольной матрицей:
(Х1(^) ^
Х1С2) Х1 (^2 )
ух1(^ ) х2(^т ) ... хт ((т ) 0
Перейдем на уровень целевого требования минимизировать отклоне -ние «глюкозной» траектории х(^, у), ^ е [0, Т] от желаемого норматива х . Реализуем эту установку на сетке узловых точек с формализацией откло -нения в среднеквадратическом смысле. В результате приходим к целевой функции:
т
1
+
Р(У) = 2 2 (х((к, у) - х) , 2 к=1
которая подлежит минимизации.
В развернутой записи с обозначением Гк = хе а*к - х получаем
1 т к 1 т к
р(у)=2 2(2 х] ({к) У]+Гк)2=-г 2(2 х] ({к) у^)2
2 к=1 ]=1 2 к=1 ]=1
т к 1 т
+2(2х] ({к)у]) гк + 2 2 Гк.
к=1 ]=1 2 к=1
Подчеркнем, что (р(у) - выпуклая квадратичная функция, которая имеет следующую структуру р( у) = р( у{) + р2( у1, у2) + ... + + Рт(у1, ..., ут). Укажем формулы для частных производных:
— ( ) т к т
=2 (2 х] ({к) у]) х5 (гк) + 2 х, ('к) Гк =
^5 к=, ]=1 к=,
т к _
= 2 (2 х] ((к )у] + гк )х5 (^к ), 5 =1 т . к=5 ] =1
В результате задача оптимизации процесса лечения диабета (в рамках метода наименьших квадратов) на основе инъекций в объеме у] в избранные моменты времени 1к формулируется следующим образом:
j(y) ® min, yj е[0,u+], j = 1, m , (4)
k _
^ xj (tk)yj + xe - atk > 0, k = 1, m . j=1
Это задача квадратичного программирования, которая допускает эффективное решение за конечное число итераций хорошо известными процедурами (метод особых точек, метод сопряженных градиентов) [4, 5].
По части корректности поставленной задачи (4) отметим, что ограничения на переменные yj , j = 1,m совместны: для достаточно малых положительных значений yj ограничения-неравенства выполняются, поскольку xe atk > 0, k = 1, m . Следовательно, задача (4) имеет решение * * *
y = (y1,..., ym), которое реализует оптимальную стратегию терапии.
Заключение
Проведена формализация процедуры лечения диабета на основе последовательных инъекций инсулина. В результате построена непрерывно-дискретная модель этого процесса с целевой установкой на нормализацию уровня глюкозы в крови. В оптимизационной классификации получена задача квадратичного программирования, которая допускает эффективное численное решение.
Дальнейшее исследование в этом направлении связано с проведением вычислительного эксперимента на основе реальных числовых данных из медицинской практики.
Литература
1. Lenhart S., Workman J. T. Optimal Control Applied to Biological Models. Chapman & Hall/CRC Mathematical and Computational Biology Series, 2007. 257 p.
2. Martin Eisen. Mathematical Methods and Models in the Biological Sciences. Prentice Hall, Englewood Cliffs, New Jersey, 1988.
3. Srochko V. A., Aksenyushkina E. V. On Resolution of an Extremum Norm Problem for the Terminal State of a Linear System // Известия Иркутского государственного университета. Сер. Математика. 2020. Т. 34. С. 3-17. Текст: непосредственный.
4. Васильев Ф. П. Методы оптимизации. Москва: МЦНМО, 2011. 620 с. Текст: непосредственный.
5. Измаилов А. Ф., Солодов М. В. Численные методы оптимизации. Москва: Физматлит, 2005. 304 с. Текст: непосредственный.
Статья поступила в редакцию 13.06.2023; одобрена после рецензирования 14.06.2023; принята к публикации 23.06.2023.
OPTIMIZATION FORMALIZATION OF DIABETES
TREATMENT PROCEDURE
Elena V. Aksenyushkina
Cand. Sci. (Phys. and Math.), A/Prof,
Baikal State University
11 Lenin Street, Irkutsk, 664003, Russia
Abstract. A formalization of the problem associated with diabetes therapy based on successive injections of insulin in order to minimize the quadratic deviation of the current glucose concentration from the normal level is given. Modeling of this procedure is carried out in a continuous-discrete version in the terminology of optimal control. The phase system is described by differential equations, permissible controls, phase constraints and the target functional are formalized in discrete time on a given grid of values. As a result, a special quadratic programming problem is constructed, which allows an effective solution in a finite number of iterations.
Keywords: diabetes treatment procedure, optimization in continuous-discrete time, quadratic programming problem.
For citation
Aksenyushkina E. V. Optimization Formalization of Diabetes Treatment Procedure // Bulletin of Buryat State University. Mathematics, Informatics. 2023. N. 2. P. 49-54.
The article was submitted 13.06.2023; approved after reviewing 14.06.2023; accepted for publication 23.06.2023.