Павлов В. А., Носовець О. К., Давидько О. Б., Дяк А. І.,Шапошник О. Г. Змішано-цілочисельна задача оптимізації лікувальної стратегії // Міжнародний науковий журнал "Інтернаука". — 2018. — №21.
Технічні науки
УДК 004.02 +616.1
Павлов Володимир Анатолійович
кандидат технічних наук, доцент,
доцент кафедри біомедичної кібернетики
Національний технічний університет України
«Київський політехнічний інститут імені Ігоря Сікорського»
Павлов Владимир Анатолиевич
кандидат технических наук, доцент,
доцент кафедры биомедицинской кибернетики
Национальный технический университет Украины
«Киевский политехнический институт имени Игоря Сикорского»
Pavlov Vladimir
Candidate of Technical Sciences, Docent,
Associate Professor of the Department of Biomedical Cybernetics
National Technical University of Ukraine
«Igor Sikorsky Kyiv Polytechnic Institute»
Носовець Олена Костянтинівна
кандидат технічних наук,
доцент кафедри біомедичної кібернетики
Національний технічний університет України «Київський
політехнічний інститут імені Ігоря Сікорського»
Носовец Елена Константиновна
кандидат технических наук,
доцент кафедры биомедицинской кибернетики
Национальный технический университет Украины
«Киевский политехнический институт имени Игоря Сикорского»
Nosovets Olena
Candidate of Technical Sciences,
Associate Professor of the Department of Biomedical Cybernetics
National Technical University of Ukraine
«Igor Sikorsky Kyiv Polytechnic Institute»
Давидько Олександр Богданович
магістрант
Національного технічного університету України
«Київський політехнічний інститут імені Ігоря Сікорського»
Давидько Александр Богданович
магистрант
Национального технического университета Украины
«Киевский политехнический институт имени Игоря Сикорского»
Davydko Oleksandr
Master Student of the
National Technical University of Ukraine
«Igor Sikorsky Kyiv Polytechnic Institute»
Дяк Андрій Ігорович
студент
Національного технічного університету України
«Київський політехнічний інститут імені Ігоря Сікорського»
Дяк Андрей Игоревич
студент
Национального технического университета Украины
«Киевский политехнический институт имени Игоря Сикорского»
Diak Andrii
Student of the
National Technical University of Ukraine
«Igor Sikorsky Kyiv Polytechnic Institute»
Шапошник Ольга Геннадіївна
студент
Національного технічного університету України
«Київський політехнічний інститут імені Ігоря Сікорського»
Шапошник Ольга Геннадьевна
студент
Национального технического университета Украины
«Киевский политехнический институт имени Игоря Сикорского»
Shaposhnyk Olha
Student of the
National Technical University of Ukraine
«Igor Sikorsky Kyiv Polytechnic Institute»
ЗМІШАНО-ЦІЛОЧИСЕЛЬНА ЗАДАЧА ОПТИМІЗАЦІЇ ЛІКУВАЛЬНОЇ СТРАТЕГІЇ
СМЕШАННО-ЦЕЛОЧИСЛЕННАЯ ЗАДАЧА ОПТИМИЗАЦИИ ЛЕЧЕБНОЙ СТРАТЕГИИ
MIXED-INTEGER OPTIMIZATION PROBLEM OF TREATMENT STRATEGY
Анотація. В роботі розглянуто персоніфікований розрахунок лікувальної стратегії. Модель прогнозу стану об'єкту допускає нелінійну залежність від його початкового стану та параметрів об'єкту. Задача оптимізації після підстановки початкового стану пацієнта та його параметрів приймає вид лінійної змішано-цілочисельної задачі. Розглянуто приклад розрахунку оптимальної стратегії.
Ключові слова: оптимізація, змішано-цілочисельне лінійне програмування, моделювання, прогнозування, початковий стан об'єкту, лікувальна стратегія.
Аннотация. В работе рассмотрен персонифицированный расчет линейной лечебной стратегии. Модель прогноза состояния объекта допускает нелинейную зависимость от его начального состояния. Задача оптимизации после подстановки начального состояния пациента и его параметров принимает вид линейной смешано-целочисленной задачи. Рассмотрен пример расчета оптимальной стратегии.
Ключевые слова: оптимизация, смешанно-целочисленное линейное программирование, моделирование, прогнозирование, начальное состояние объекта, лечебная стратегия.
Summary. A personalized calculation of linear treatment strategy is considered. The forecast model admits a nonlinear dependence on initial state and parameters of object. After the substitution of the initial state and parameters of the patient the optimization task takes the form of a linear mixed-integer problem. An example of an optimal treatment strategy calculation is considered.
Key words: optimization, mixed integer linear programming, modeling, forecasting, the initial state of object, treatment strategy.
Постановка проблеми. Подання задач оптимального управління і задач прийняття рішень в еквівалентному вигляді задачами математичного програмування досить відомо. Зазначимо, монографію [1], де розглянуто загальний підхід для такого представлення. Однак, в кожній конкретній задачі управління або прийняття рішень з'являється свій привід для розгляду окремих оптимізаційних моделей, що дозволяють знаходити компроміс між обчислювальними ресурсами що використовуються, точністю рішення, часом розрахунку. Такі моделі дозволяють виявити можливості для поліпшення результатів розрахунків, що робить доцільним їх розгляд для практичного застосування.
Аналіз останніх досліджень і публікацій. Дослідження складають праці таких фахівців, як О.І. Шелудько [3], В.В Ванін, О.В. Павлов [4] , A.В. Павлов [5].
Мета дослідження: розробити програмний комплекс вирішення проблеми переходу від суттєво-математичного формалізму практичної задачі оптимізації до вхідного формалізму програмного інструменту вирішення задач оптимізації змішано-цілочисельного лінійного класу.
Виклад основного матеріалу. Припустимо, що змістовна постановка задачі прийняття рішення описує статистичними даними стан об'єкту до застосування управлінського рішення - підматрицею, початкове значення критеріальної змінної - вектором, відомі управлінські рішення, що застосовано - підматрицею, кінцевий стан об'єкту отриманий, внаслідок застосування управлінських рішень - підматрицею та якість стану об'єкту, що було досягнуто - вектором, то зазначена статистика може бути подана блочною матрицею об'єкт-властивості Х:
(1)
(2)
Будемо далі формувати функціонал і обмеження оптимізаційної задачі таким чином, щоб не вийти за межі класу лінійних змішано-цілочисельних задач. При цьому, використовуючи (1), ми можемо моделювати для кожної змінної кінцевого стану та критеріальної змінної, нелінійні співвідношення-моделі вигляду:
Тоді можливо записати оптимізаційну задачу у наступному вигляді:
(3)
Тут принципово важливо те, що в оптимізаційній задачі для конкретного стану передісторії об'єкту ми не виходимо за рамки ЗЦЛП задачі.де деякі з управлінь u можуть приймати тільки цілочисельні значення.
Дійсно, оскільки рядки підматриць відомі нам, як стан об'єкту до оптимізації, то підставляючи конкретні її значення у моделі, що отримано, ми замість нелінійних членів маємо відповідні константи, що налаштовують моделі на його передісторію. Сама ж задача оптимізації (3) після такої підстановки має лінійний вигляд (4).
(4)
Задача розрахунку оптимального керуючого впливу може бути приведена до вигляду (4) не тільки при наявність матриці (1) спостережень за одним і тим же об'єктом, але і при наявності статистики для множини досить однорідних об'єктів: статистика подається матрицею об'єкт-властивості:Задачі (3,4) можуть бути використані для прийняття оптимальних рішень в разі достатнього рівня адекватності подання реакцій об'єкту статистичними моделями виду (2).
(5)
(6)
де, та де кожний рядок підматриці містить g характерних параметрів відповідного об'єкту.Однорідність розуміється в сенсі можливості подання адекватними статистичними моделями співвідношень (2) за даними матриці об'єкт-властивості (5), де відповідні рядки матриці Х відносяться вже не до різних варіантів переходу одного і того ж об'єкта з різних початкових станів у відповідні кінцеві, а описують перехід з деякого початкового стану в кінцевий для різних об'єктів. У такому випадку для врахування особливостей кожного об'єкту матриця об'єкт-властивості доцільно розширюється за рахунок врахування характерних параметрів об'єктів, що досліджуються. Блочна матриця вихідних даних завдання тоді має вигляд:
Вид задачі оптимізації при цьому практично не змінюється, а в співвідношення (2) додаються члени моделі, що пов'язані зі характерними параметрами об'єкту.
(7)
Відповідна задача оптимізації при цьому набуває вигляду:
(8)
При необхідності оптимізаційного розрахунку керуючих впливів для деякого об'єкту, в наявні моделі виду (7) підставляються значення його передісторії і параметри, тим самим ми підстроюємо систему обмежень на передісторію і параметри об'єкта, для якого будемо шукати оптимальне значення керуючих змінних. Задача (8) при цьому приймає лінійний вид щодо змінних, що розраховуються, аналогічний (4). В результаті, підхід, що ми розглянули, дозволяє формувати частково нелінійні моделі для опису процесів, що досліджуються, тим самим збільшуючи точність прогнозу наших моделей і в той же час не ми виходимо на етапі оптимізації з класу ЗЦЛП задач.
Інтерес до прикладних задач того виду, що було розглянуто, викликаний тим, що в разі можливості подання реакції біологічного об'єкта у вигляді такої дискретної моделі виникає можливість не тільки оптимізувати процес клінічних випробувань в процесі самого випробування, але і відкривається можливість індивідуальних налаштувань лікарських впливів для конкретного об'єкту з урахуванням його індивідуальних параметрів і стану перед лікувальним періодом.
В якості прикладу було обрано пошук оптимальної стратегії проведення операції та лікування хворого у післяопераційний період, після виконання операції аортокоронарного шунтування. База даних налічує 129 записів, кожен з яких відповідає окремому пацієнтові та 271 змінну. Для даного дослідження з них було відібрано 21 (табл. 1).
Таблиця 1
Застосовані у дослідженні змінні
Позначення |
Атрибут |
Короткий опис |
X1 |
Возраст |
вік пацієнта |
X2 |
Время наблюдения |
час, який пацієнт провів під лікарським наглядом |
X3 |
NYHA |
функціональний клас хворого на серцеву недостатність |
X4 |
КСО |
кінцевий систолічний об’єм при госпіталізації пацієнта |
X5 |
КСР |
кінцевий систолічний розмір при госпіталізації пацієнта |
X6 |
К-во КА |
кількість коарктованих артерій |
X7 |
Жизнеспособный миокард |
відсоток життєздатних тканин міокарду |
X8 |
Лактат |
частка солей молочної кислоти у крові пацієнта при госпіталізації |
X9 |
ИСЛС |
систолічний тиск при госпіталізації |
X10 |
САД1 |
сатурація гемоглобіну змішаної венозної крові киснем при госпіталізації, |
X11 |
SvO2 |
коефіцієнт утилізації кисню при госпіталізації |
X12 |
КУО2 |
коефіцієнт утилізації кисню при госпіталізації |
X13 |
САД4 |
систолічний тиск після операції та лікування. |
X14 |
SvO2(4) |
сатурація гемоглобіну змішаної венозної крові киснем після операції та подальшого лікування. |
X15 |
КУО2(4) |
коефіцієнт утилізації кисню після операції та подальшого лікування. |
U1 |
К-во шунтов |
кількість аортокоронарних шунтів, що були імплантовані пацієнтові під час операції |
U2 |
Длительность |
час прийому препаратів після операції |
U3 |
Эфедрин |
доза препарату «Ефідрин» |
U4 |
Фенамин |
доза фенаміну |
U5 |
Тирамин |
доза тираміну |
I |
Продолжительность жизни |
Тривалість життя пацієнта після операції |
Дані подано у вигляді матриці Х об'єкт-властивості, де кожен об'єкт (рядок) описано набором властивостей (стовпці). Матриця містить інформацію про стан об'єктів до та після застосування лікувальної дії. Таким чином, можемо її подати у вигляді блокової матриці, подібної (6), де - підматриця параметрів, - підматриця станів до застосування управління, - підматриця станів після застосування управління, - вектор значень критеріальної змінної після лікування,- підматриця дозувань лікувальних препаратів, що застосовано.
Задачу розрахунку сформулюємо наступним чином: розрахувати оптимальне поєднання доз лікувальних препаратів та кількості шунтів для конкретного пацієнта виходячи з отримання максимальної тривалості життя після операції при заданих обмеженнях на значення інших змінних стану.
Використовуючи дані блокової матриці Х, отримаємо моделі прогнозу кінцевого стану об'єкту від вхідних станів та управлінь у вигляді (2). Для цього використаємо версію модифікованого алгоритму МГУА з комбінаторної селекцією і ортогоналізації змінних (МАКСО) [3-5] в який введена функція протекції змінним в лінійній формі.
В результаті розрахунку, було одержано наступні моделі:
Запишемо формальну постановку оптимізаційної задачі, ґрунтуючись на моделях, що отримано. Для формування обмежень визначимо границі змінних з матриці (табл. 2). Також треба зазначити, що дана задача відноситься до класу ЗЦЛП, через те, що змінна U1 обмежена лише цілими значеннями.
Таблиця 2
Обмеження змінних з матриці
Змінна |
Нижня границя |
Верхня границя |
САД4 (Х13) |
18 |
81 |
SvO2(4) (X14) |
0.35 |
0.99 |
КуО2(4) (Х15) |
21 |
50 |
З урахуванням виразів знайдених моделей сформуємо оптимізаційну задачу :
На даному етапі задача сформована у вигляді (13). Наступним кроком перетворимо (13) до її часткового виду, налаштованого на конкретний об'єкт. Для цього було обрано об’єкт з такими показниками Хp та Хb: Возраст = 73, NYHA = 4, Время наблюдения = 36, КСО госпит = 151.2, КСР госпит = 5.5, К-во КА = 4, Жизнеспособный миокард = 85, Лактат = 1.2, ИСЛС = 420, САД1 = 78, SvO2(1) = 0.62, КУО2(1) = 37.3737374, К-во шунтов 4, Длительность = 70, Эфедрин = 0, Фенамин = 0, Тирамин = 240, Продолжительность жизни = 26, САД4 = 85, SvO2(4) = 0.72, КУО2(4) = 28.847465. Після налаштування, вирази моделей набудуть такого вигляду:
Для зручності роботи з задачами того типу, що було розглянуто, була розроблена спеціалізована програмна система. Система дозволяє в інтерактивному режимі вводити повні вирази моделей, що визначають цільову функцію та кінцеві стани об'єкту в зручних для оператора позначеннях, формувати обмеження, підставляти в них значення передісторії, параметрів об'єкту, розраховувати вхідні параметри для ЗЦЛП задачі. Оператор формує оптимізаційні задач в зручних для себе позначеннях, паралельно довизначає їх системою змінних для структури даних програмного інструменту оптимізації. Після завершення введення формалізмів оптимізаційної задачі проводиться підстановка в неї параметрів і початкових умов стану об'єкту. Система формує підсумкову лінійну форму та визначає параметри для структури даних на вході програмного інструменту вирішення ЗЦЛП задачі. Потім відбувається розрахунок оптимального рішення.
Проілюструємо вищесказане: вигляд інтерфейсу системи у браузері зображено на рис. 1.
Рис. 1. Початок роботи з системою
Джерело: розробка автора
Вводимо в систему інформацію про змінні, що відносяться до підматриці Хb і U. Формуємо цільову функцію I та набір обмежень. Моделі кінцевих станів формують набір обмежень для кінцевого стану об'єкту - помістимо вирази відповідних моделей у подвійні нерівності, що обмежують їх значення (рис. 2):
Рис. 2. Сформована задача оптимізації
Джерело: розробка автора
Тепер здійснимо налаштування вільного члену моделей на передісторію конкретного об'єкту. Для цього використаємо значення атрибутів Хb для, наприклад, першого рядка і підставимо їх у моделі. В результаті будемо мати:
Рис. 3. Налаштування моделей на об'єкт
Джерело: розробка автора
Після виконання розрахунків, сервер повертає інформацію: статус, оптимальні значення змінних управління і кінцеві значення змінних стану (рис. 4).
Рис. 4. Результати розрахунку
Джерело: розробка автора
Результати, що отримано задовольняють встановленим обмеженням, розраховано оптимальні значення для змінної критерію і змінних управління. У таблиці 3 наведена порівняльна характеристика розрахованих даних та реальних значень:
Таблиця 3
Порівняння розрахункових та реальних даних
Змінна |
Розрахункове значення |
Реальне значення |
Продолжительность жизни |
40.73715335381506 |
26 |
Длительность |
241.7351 |
70 |
Количествошунтов |
5 |
4 |
Тирамин |
240 |
240 |
Количествошунтов |
5 |
4 |
Тирамин |
240 |
240 |
Фенамин |
0 |
0 |
Эфедрин |
0 |
0 |
САД4 |
81.00000019143623 |
85 |
SvO2(4) |
0.5992951327303244 |
0.72 |
KyO2(4) |
30.881212201036025 |
28.847465 |
Висновки. Розрахунок лікувальної стратегії представлено, як задача оптимізації прогнозу стану об'єкту, що подається дискретними, лінійними за управліннями, моделями. Моделі кінцевих станів описуються, як нелінійні за початковими станами об'єкту та лінійно-цілочисельні за управліннями, що дозволяє підвищити точність прогнозу моделей і не вийти з класу ЗЦЛП задач в задачі оптимізації. Підхід, що запропоновано дозволяє здійснити індивідуальний розрахунок лікувальних впливів для конкретного об'єкту з урахуванням його індивідуальних параметрів і стану перед лікувальним процесом. Розроблено програмний комплекс, що дозволяє у інтерактивному режимі вводити вирази моделей в зручних для оператора позначеннях, формувати обмеження, підставляти в них значення передісторії, параметрів об'єкту, розраховувати вхідні параметри для структури даних на вході програмного інструменту вирішення ЛП задачі. Розглянуто приклад розрахунку оптимальної лікувальної стратегії.
Література
References