- •ВВЕДЕНИЕ
- •1.1. Формулировка задачи линейного программирования
- •1.2. Распределение специализированных бригад скорой помощи по категориям больных
- •1.3. Некоторые другие задачи на распределение
- •1.4. Разработка комплексной лекарственной терапии
- •1.5. Выработка оптимального плана массового лечения
- •1.6. Определение линейных разделяющих функций
- •2.1. Метод динамического программирования
- •2.4. Выравнивание символьных последовательностей
- •3. ПРИМЕНЕНИЕ ТЕОРИИ ИГР ДЛЯ ОПТИМИЗАЦИИ КЛИНИЧЕСКИХ РЕШЕНИЙ В ХИРУРГИИ
- •3.1. Игры и методы их решения
- •3.2. Элементы теории статистических решений
- •3.3. Критерии принятия решений в условиях неопределенности
- •3.4. Принятие решений при острых хирургических заболеваниях органов брюшной полости
- •3.5. Минимизация риска хирургического вмешательства в онкологии
- •Список литературы
- •ПРИЛОЖЕНИЯ
1.3. Некоторые другие задачи на распределение
Большое число задач рационального распределения сводится к транспортной задаче ЛП. В данном разделе приведены примеры таких задач, возникающих в здравоохранении.
А. Распределение численности людей при массовом обследовании. При периодических профосмотрах, массовом обследовании населения, массовых профилактических лечебных мероприятиях возникает следующая задача. Имеется m предприятий П1, П2, , Пm с числом работающих на них людей
N1, N2, , Nm соответственно. Имеется также n медпунктов M1, M2, , Mn , расположенных на некотором расстоянии от предприятий. Затраты на транспортировку одного человека из предприятия Пj в медпункт Mi известны и составляют cij . Всех людей на данных предприятиях необходимо обследовать в течение Т дней. Дневная нагрузка на все медпункты составляет
(N1 + N2 + + Nm ) = N1 + N2 + + Nm = b1 +b2 + +bm .
T T T T
Она должна соответствовать суммарной пропускной способности медпунктов. Если для каждого из них пропускная способность равна Q1, Q2, , Qn , то должно выполняться
n |
m |
∑ai = ∑bj . |
|
i=1 |
j=1 |
Требуется так распределить число людей xij , направляемых от предприятия
|
n m |
Пj в медпункт |
Mi , чтобы суммарные транспортные затраты L = ∑ ∑cij xij |
|
i=1 j=1 |
были минимальными.
Ясно, что данная задача так же, как и предыдущая, сводится к составлению ограничений (1.7) и (1.8), т. е. приводит к основной задаче ЛП.
Далее сформулирована одна и та же (с математической точки зрения) задача ЛП в трех различных содержательных постановках.
Б. Распределение врачей разной специальности по разным группам больных в условиях массового заболевания. В условиях массового заболева-
ния имеется 2 категории врачей, A1 и A2 , которые обслуживают 3 группы больных, B1, B2, B3 (различающихся по стадиям заболевания). Дневное чис-
15
ло вызовов по группе больных B j составляет b j . Число врачей категории Ai равно ki . Так как каждый врач в день может обслужить N больных, то все врачи категории Ai в день могут обслужить ai больных, где ai = Nki . Считается, что общее число вызовов точно равно общему числу выездов, т. е. a1 +a2 = b1 +b2 +b3. Пусть xij – число больных из группы B j , которых обслуживает врач из категории Ai . Пусть также качество обслуживания cij
больного из B j врачом из Ai определяется матрицей cij . Нужно рассчитать
оптимальное (наилучшее с точки зрения получения наибольшего суммарного качества обслуживания) распределение врачей по группам больных, т. е. оп-
тимальные значения элементов матрицы xij .
В. Распределение лекарственных препаратов по различным группам больных. Для лечения трех групп больных B1, B2, B3 применяются два медикаментозных препарата A1 и A2 . Так как общее число доз этих препаратов равно общему числу больных, то каждому больному может быть выдана только одна доза какого-то из этих двух лекарств. Число больных в группе B j равно bj . Число доз препарата Ai равно ai . Эффективность лечения боль-
ного типа B j препаратом Ai равна cij . Пусть xij – число больных группы B j , получающих препарат Ai . Нужно распределить дозы препарата по больным так, чтобы суммарная эффективность лечения была максимальной, т. е. найти оптимальные значения элементов матрицы xij .
Г. Распределение операторов по рабочим местам. Операторов готовят для управления сложными объектами трех видов, B1, B2, B3 . По результатам психофизиологического тестирования все операторы были разделены на две группы, A1 и A2 . Число объектов вида B j составляет b j . Число операторов в
группе Ai равно ai . Общее число объектов управления равно общему числу операторов. Известно, что эффективность работы оператора из группы Ai при управлении объектом B j определяется как элемент cij матрицы сij . Пусть xij – число объектов вида B j , на которые предполагается направить операторов из группы Ai . Нужно распределить всех операторов по объектам
16
так, чтобы суммарная эффективность их работы была максимальной, т. е. найти оптимальные значения элементов матрицы xij .
Данная задача ЛП, изложенная в трех различных постановках, решается аналогично той, которая дана в примере 1.1. В случае, если a1 = 60, a2 = 40,
b1 = 30, b2 = 20, b3 = 50,
|
|
|
|
1 |
0,2 |
0,1 |
|
= |
|
|
|
. |
|||
сij |
|
0,7 |
1 |
||||
|
|
|
0,3 |
|
|||
Эта задача имеет решение |
|
|
|
|
|
|
|
|
|
|
|
30 |
20 |
10 |
|
xij |
|
= |
|
0 |
|
, |
|
|
|
|
|
0 |
40 |
||
полученное также геометрическим способом.
1.4. Разработка комплексной лекарственной терапии
Пусть при лечении некоторого больного его комплексный диагноз состоит из заболеваний D1, D2, , Dm , проявляющихся одновременно. Для ле-
чения этих заболеваний у врача имеется n лекарственных |
препаратов |
A1, A2, , An . Эффективность применения единицы препарата |
Ai для лече- |
ния заболевания D j равна сij (эта величина не обязана быть положитель-
ной). Кроме полезного (в целом) эффекта, каждый препарат обладает некоторой токсичностью. Токсичность единицы препарата Ai равна qi . Так как, принимая все n препаратов в количествах x1, x2, , xn соответственно, мы
рассчитываем на получение положительного лечебного эффекта для всех m заболеваний, то эффективность применения всех препаратов для заболевания D j не должна быть меньше некоторой положительной величины bj . Это
условие дает систему m неравенств
n |
|
∑сij xi ≥bj , |
j =1, , m. |
i=1 |
|
Условия ограничения токсичности всех принимаемых лекарств имеют вид
n
∑qi xi ≤Q ,
i=1
17
где Q – некоторая постоянная величина. Задача формулируется как максимизация суммарного эффекта воздействия принимаемых n препаратов на все m заболеваний, равного
m n
Lс=x∑ ∑ ij i , j=1i=1
при сформулированных ранее ограничениях.
Данная задача приводится к основной задаче ЛП заменой ограниченийнеравенств ограничениями-равенствами и введением (n +1) добавочных не-
отрицательных переменных. Сходная по содержательной постановке задача рассмотрена также в [10]. Вместо лекарственных препаратов в данной задаче можно рассматривать другие лечебные воздействия, физиотерапию, бальнеотерапию, физические нагрузки и т. д., суммарная доза которых также должна быть ограничена.
1.5. Выработка оптимального плана массового лечения
Данная задача взята из [11]. Пусть в результате массовой эпидемии имеется большой контингент больных в количествеN человек, нуждающихся в медицинской помощи. Эти больные находятся в различных состояниях B1, B2, , Bm , соответствующих различной степени тяжести заболевания (рис.1.3). Число
больных, находящихся в состоянии |
|
m |
|
, их относи- |
B j , равно N j |
∑ N j = N |
|||
|
|
|
|
|
|
j=1 |
|
|
|
тельное
18
N1 |
|
Nj |
Nm |
|||
Больные |
|
Больные |
|
Больные |
|
|
B1 |
|
Bj |
|
Bm |
|
|
xi1, Ai1, γi1 |
|
|
xij , Aij , γij |
xim, Aim, γim |
||
|
|
|||||
|
|
|
|
|
||
План |
|
План |
|
План |
|
|
лечения |
|
лечения |
|
лечения |
|
|
P1 |
|
Pi |
|
Pr |
|
|
|
|
Рис. 1.3 |
|
|
||
число сj = N j
N . В нашем распоряжении имеется r различных планов лечения этих больных P1, P2, , Pr , каждый из которых требует использования n медикаментов. Наличный запас медикаментов ограничен и равен M1, M2, , Mn соответственно. Пусть xij – относительное число больных, находящихся в состоянии B j , к которым применяется лечение по плану Pi ; абсолютное число таких больных равно xij N j . Пусть также Aij – эффектив-
ность такого лечения, выражаемая в доле выздоровевших больных, при этом абсолютное число выздоровевших (из класса B j , к которым применено лече-
ние Pi ) равно Aij xij N j . Допустим, |
что требуемые количества медикаментов |
||
для лечения одного больного типа |
B j |
по плану Pi |
заданы в виде вектора |
γij =(γ(ij1), γij(2), , γij(n)), причем γij(1) |
– количество |
первого медикамента, |
|
γij(2) – второго, и т. д.
Примем за критерий качества L выбранной системы лечения (т. е. вы-
бранных значений всех переменных |
xij , i =1, , r; |
j =1, , m ) отношение |
|||||
выздоровевших больных к общему числу больных. Тогда |
|
||||||
|
1 |
m r |
m r |
N j |
|
m r |
|
L = |
|
∑ ∑Aij xij N j = ∑ ∑ |
|
Aij xij |
= ∑ ∑qij xij , |
(1.19) |
|
|
N |
||||||
|
N j=1i=1 |
j=1i=1 |
|
j=1i=1 |
|
||
19
где qij = c j Aij . Так как xij – относительные величины, то по группам больных B1, B2, , Bm должны выполняться равенства
|
|
|
|
|
r |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
∑xij = |
1, |
j =1, , m. |
(1.20) |
||||||
|
|
|
|
|
i=1 |
|
|
|
|
|
|
|
|
|
Из-за ограниченности запасов медикаментов должны также выполняться |
||||||||||||||
неравенства |
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
m |
r |
γ(k )x N |
|
≤ M |
|
, |
k =1, , n . |
|
|||
|
|
|
∑ ∑ |
j |
k |
|
||||||||
|
|
|
j=1i=1 |
ij |
ij |
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
||
Разделив обе части этих неравенств на N , получим |
|
|||||||||||||
|
|
|
|
m |
r |
|
|
|
≤ b , |
k =1, , n , |
(1.21) |
|||
|
|
|
|
∑ ∑ f (k )x |
|
|||||||||
|
|
|
|
j=1i=1 |
ij |
ij |
|
k |
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
||
где f (k ) = |
N j |
γ(k ), |
b = |
Mk . |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
ij |
Nk ij |
k |
N |
|
|
|
|
|
|
|
|
|
|
|
Таким образом, математически задача формулируется так: найти mr неотрицательных переменных xij , которые максимизируют L (1.19) при огра-
ничениях (1.20) и (1.21).
Пример 1.2. Пусть имеется 100 больных, находящихся в трех состояниях, B1, B2, B3 . Число больных каждого вида равно 60, 30 и 10 соответственно. Имеется три плана лечения, P1, P2, P3 , которые предполагают использовать три вида медикаментов. Запасы медикаментов в условных единицах составляют: M1 =1500, M2 =800, M3 =1900. Эффективность имеющихся планов лечения задается матрицей А (табл. 1.2), а требуемые количества лекарств – матрицей γ (табл. 1.3). Нулевые элементы матрицы А, соответствующие прочеркам в матрице γ, связаны с отсутствием планирования применения данных планов лечения к соответствующим группам больных, поскольку неэффективность такого применения очевидна.
|
|
|
Таблица 1.2 |
|
|
|
Таблица 1.3 |
||
|
|
|
|
|
|
|
|
|
|
Pi |
|
Bj |
|
|
|
Pi |
|
Bj |
|
|
|
|
|
|
|
|
|
||
B1 |
B2 |
|
B3 |
|
B1 |
B2 |
B3 |
||
|
|
|
|
||||||
20
P1 |
0,8 |
0 |
1 |
P2 |
0 |
0,8 |
0 |
P3 |
1 |
1 |
0 |
P1 |
(0, 0, |
30) |
– |
(100, 0, 0) |
P2 |
– |
|
(10, 5, 5) |
– |
P3 |
(0, 10, 10) |
(20, 10, 0) |
– |
|
Используя исходные данные, вычислим
c |
= |
|
60 |
= 0,6; |
c |
= |
|
30 |
= 0,3; |
c |
= |
|
10 |
= 0,1. |
|
100 |
100 |
100 |
|||||||||||||
1 |
|
|
2 |
|
|
3 |
|
|
|||||||
Целевая функция (1.19) равна
L = 0,48x11 +0,1x13 +0,24x22 +0,6x31 +0,3x32 .
Ограничения (1.20) и (1.21) для данного примера имеют вид
x |
+ x =1; |
|
|
|
11 |
31 |
|
x22 + x32 =1; |
(1.22) |
||
x |
=1; |
|
|
|
13 |
|
|
3x |
+10x |
+ 6x |
≤15; |
|
|
|||||
|
|
22 |
13 |
|
|
32 |
|
|
|
|
1,5x22 + 6x31 +3x32 ≤8; |
|
|
||||||||
|
|
|
|
|
|
+ 6x31 ≤19. |
|
|||
18x11 +1,5x22 |
|
|||||||||
Последнюю систему неравенств преобразуем в систему равенств, введя |
||||||||||
новые неотрицательные переменные y1, |
y2, y3 . |
|
|
|||||||
y |
|
=15 −3x |
|
−10x |
|
−6x ; |
|
|||
1 |
|
22 |
|
13 |
|
32 |
|
|||
y2 =8 −1,5x22 −6x31 −3x32; |
(1.23) |
|||||||||
y |
|
=19 −18x |
|
|
−1,5x |
− |
6x . |
|
||
3 |
|
11 |
|
22 |
31 |
|
||||
Таким образом, (1.22) и (1.23) образуют систему из 6 уравнений с 8 не- |
||||||||||
известными: x11, x31, x22, |
x32, |
x13, y1, |
|
y2, y3 . Применим геометрический |
||||||
способ решения. Выберем 2 свободных переменных x31 |
и x32 . Остальные 6 |
|||||||||
будут базисными. Выразим базисные переменные через свободные:
x11 =1− x31; |
y1 = 2 −3x32; |
x22 =1− x32; |
y2 = 6,5 −6x31 −1,5x32; |
x13 =1; |
y3 = −0,5 +12x31 +1,5x32 . |
Из условия неотрицательности базисных переменных получаем x31 ≤1; x32 ≤ −4x31 + 4,33;
21
x32 ≤1; x32 ≥ 0,3 −8x31; x32 ≤ 0,66 .
Графически эти условия изображены на рис. 1.4.
Функция цели, выраженная через свободные переменные, равна
L = 0,82 +0,12x31 +0,06x32 ;
L′= L −0,82 = 0,12x31 +0,06x32 ,
и основная прямая L' = 0 имеет вид x32 = −2x31.
Из рис. 1.4 видно, что максимум L′ |
и L достигается в точке F . Таким |
||||
образом, решение задачи: |
|
|
|||
|
|
|
x31 = 0,92; x22 = 0,34; |
||
x32 = 0,66; x13 =1; |
x11 = 0,08; |
||||
L = 0,82 + 0,12 0,92 + 0,06 0,66 = 0,97 . |
|||||
x32 |
|
|
x32 = –4x31 + 4,33 |
||
|
|||||
|
1,0 |
|
|
|
|
0,66 |
|
|
F |
||
L' = max
ОДР
L' = 0
0 |
0,92 |
1,0 |
|
|
x31 |
||||
x32 = –8x31 + 0,5
Рис. 1.4
Следовательно, при выбранной оптимальной схеме лечения относительное число выздоровевших больных максимально и составляет 97 %.
22
