Компьютерные инструменты в образовании, 2018 № 1:31-50 УДК: 519.87
http://ipo.spb.ru/journal
РЕШЕНИЕ МИНИМАКСНОЙ ЗАДАЧИ РАЗМЕЩЕНИЯ В ТРЕХМЕРНОМ ПРОСТРАНСТВЕ С ПРЯМОУГОЛЬНОЙ МЕТРИКОЙ*
Плотников П. В.1, Кривулин Н. К.1 1 Санкт-Петербургский государственный университет, Санкт-Петербург, Россия
Аннотация
Рассматривается минимаксная задача размещения точечного объекта в трехмерном пространстве с прямоугольной метрикой (/^метрикой) и предлагается ее прямое аналитическое решение при помощи методов тропической (идем-потентной) математики. Сначала задача записывается в терминах тропической математики как задача тропической оптимизации, вводится параметр для обозначения минимума целевой функции, и задача сводится к решению параметризованной системы неравенств. Эта система решается относительно одной из переменных, а условия существования решений используются для нахождения оптимальных значений второй переменной с помощью вспомогательной задачи оптимизации. Затем вспомогательная задача решается аналогичным образом, и находится значение третьей переменной. Полученное общее решение преобразуется в набор прямых решений, записанных в компактной форме для различных случаев соотношений между исходными параметрами задачи.
Ключевые слова: задача 1-центра, трехмерное пространство, прямоугольная метрика, идемпотентное полуполе, тропическая оптимизация, полное решение.
Цитирование: Плотников П. В., Кривулин Н. К. Решение минимаксной задачи размещения в трехмерном пространстве с прямоугольной метрикой // Компьютерные инструменты в образовании. 2018. № 1. С. 31 -50.
1. ВВЕДЕНИЕ
Задачи размещения составляют важный класс задач оптимизации, который имеет многочисленные приложения в различных областях науки и техники [1-3]. Например, при решении задач пространственного размещения объектов информационной системы часто требуется определить оптимальное положение управляющего объекта относительно системы управляемых объектов. Такого рода научно-технические задачи встречаются при планировании размещения центров управления, хранения и обработки данных, собранных с видеокамер системы видеонаблюдения в здании [4].
Рассматриваемая система видеонаблюдения состоит из трех главных компонент: 1) видеокамеры, производящие входной поток видеоинформации, 2) система передачи
*Работа выполнена при финансовой поддержке Российского фонда фундаментальных исследований (грант № 18-010-00723).
сигнала от камер до центра управления, хранения и обработки данных и 3) центр управления. В качестве объектов, относительно которых необходимо решать задачу размещения, выступают видеокамеры внутреннего и наружного наблюдения. Суть задачи размещения (в общем задачу называют задачей 1-центра) состоит в поиске оптимального положения центра управления системой. При этом необходимо минимизировать функцию расстояния до самой дальней камеры, чтобы снизить влияние шумов и повысить качество сигнала, благодаря уменьшению длины кабеля. Внутри одного этажа кабели чаще всего прокладывают вдоль линий разделения пола, стен и потолка. Междуэтажные перекрытия проходятся по вертикальным шахтам. Поэтому можно считать, что все изгибы кабеля осуществляются под прямыми углами, что позволяет измерять его длину при помощи прямоугольной метрики (/1-метрики).
Существуют различные подходы к решению задачи 1-центра в прямоугольной метрике [5-8]. В качестве примера можно рассматривать методы и алгоритмы линейного и смешанного целочисленного линейного программирования. Алгоритмический подход обычно обеспечивает численное нахождение одного из решений и требует применения специальных программных средств. Такой подход не позволяет получить все решения аналитически в виде явных формульных зависимостей, в удобном для дальнейшего анализа и непосредственных расчетов виде. Возможностью получать решения в аналитическом виде обладает подход на основе применения методов тропической математики.
Тропическая математика — это перспективная и быстро развивающаяся область прикладной математики и алгебраической информатики, которая занимается изучением теории и приложений полуколец с идемпотентным сложением [9-19]. Основное преимущество использования тропического подхода состоит в том, что для целого ряда задач он дает возможность получать полные решения в виде явных формул, удобных для представления, анализа и интерпретации получаемых результатов. Это позволяет эффективно решать задачу размещения в случаях, когда алгоритмическое решение оказывается по тем или иным причинам невозможным или нецелесообразным. Кроме того, процедуры, построенные на основе методов тропической математики, часто имеют меньшую вычислительную сложность в сравнении с известными итерационными алгоритмами.
Подход на основе методов тропической математики находит свое применение при решении задач в различных областях науки и техники, включая задачи принятия решений при анализе результатов оценки альтернатив на основе парных сравнений [20], и задачи планирования сроков выполнения проектов [21]. На основе этого подхода в работах [22-26] получены решения минимаксных задач размещения точечного объекта на плоскости с прямоугольной и чебышевской метрикой без ограничений и с ограничениями на допустимую область размещения.
Настоящая работа направлена на дальнейшее развитие и применение методов тропической математики для решения задач размещения. В статье рассмотрена минимаксная задача размещения точечного объекта в трехмерном пространстве с прямоугольной метрикой. Решение основано на подходе, развитом в работах [27-29], который заключается в замене исходной задачи оптимизации, записанной в терминах тропической математики, на параметризованную систему неравенств с последующим ее решением. На основе этого подхода найдено прямое полное решение задачи, представленное в терминах тропической математики, а также в обычной форме.
2. ЗАДАЧА РАЗМЕЩЕНИЯ В ТРЕХМЕРНОМ ПРОСТРАНСТВЕ
В этом разделе рассматривается минимаксная задача размещения точечного объекта в трехмерном пространстве с прямоугольной метрикой. Она состоит в поиске на основе уже имеющегося набора объектов новой точки, расстояние от которой до самого дальнего объекта из набора было бы минимальным.
Пусть заданы две точки х = (xi,x2,x3)T и y = (yi,y2,y3)T в пространстве R3. Расстояние между ними в прямоугольной метрике вычисляется по формуле
p(x, y) = |xi - yi | + |Х2 - У2 | + |X3 - Уз (1)
Рассмотрим набор точек vi = (гц, r2i, r3i)T е R3 и чисел ¡i е R для всех i = 1,..m. Числа ¡i могут отражать дополнительные затраты, связанные с перемещением между объектами. Зададим функцию ф(х), которая определяет максимальное расстояние от точки х = (xi, x2,x3)T до набора точек vi следующим образом:
ф(х) = max (p(vi,х) + ¡f) = max(|rii - xi | + |г2i - x2l + |r3i - x3| + ¡f) = ф^,x2,x3). (2)
i<i <m i<i <m
Тогда рассматриваемая минимаксная задача размещения точечного объекта в трехмерном пространстве записывается так:
min ф(х). (3)
xeR3 (3)
Для решения задачи и получения полного решения в явном аналитическом виде ниже будет использован аппарат тропической математики.
3. ЭЛЕМЕНТЫ ТРОПИЧЕСКОЙ МАТЕМАТИКИ
В этом разделе приведем обзор основных понятий тропической математики [9-19], на которые будет опираться дальнейшее изложение. Рассмотрим числовое множество X с заданными на нем ассоциативными и коммутативными операциями сложения Ф и умножения ®, которые образуют на X коммутативное полукольцо с нулем O и единицей 1. Такое полукольцо обозначим через (X, 0,1, Ф, ®). Сложение будем считать идемпотент-ным (то есть для любого числа x е X выполняется x Ф x = x), а умножение - обратимым (то есть для каждого x ф O существует обратный элемент x-i такой, что x ® x-i = 1). В силу того, что (X \ {0}, 1, ®) образует коммутативную группу по умножению, описанную структуру (X, 0,1, Ф, ®) обычно называют идемпотентным полуполем.
Операция возведения в степень с целым показателем вводится стандартным образом. Для любого ненулевого числа x е X и натурального n определим x0 = 1, xn = x ® xn-i, x-n = (x-i)n и 0n = O. Будем считать, что операция возведения в целую степень x ф O может быть распространена на случай рационального показателя степени. Далее для упрощения математических выкладок знак умножения ® в алгебраических выражениях, как обычно, будет опускаться.
В силу того, что сложение идемпотентно, можно определить отношение частичного порядка < так: x < y тогда и только тогда, когда x Ф y = y. Из этого определения следует пара неравенств x < x Ф y и y < x Ф y, а также равносильность неравенства x Ф y < z и системы неравенств x < z и y < z для любых x, y, z e X. Нетрудно проверить, что выполняются свойства монотонности операций сложения и умножения. Так из условия x < y
для любого z вытекают неравенства x Ф z < y Ф z и xz < yz. Кроме того, для любой пары x, y Ф O из неравенства x < y следует x-1 > y-1.
Далее будем дополнительно предполагать, что введенный частичный порядок явля-
1/2 1/2
ется линеиным. Несложно проверить неравенство x1/2y1/2 < x Ф y, которое является тропическим аналогом неравенства между геометрическим и арифметическим средними.
В качестве примера алгебраической структуры рассматриваемого типа можно привести вещественное идемпотентное полуполе Rmax,+ = (R и {-те}, -те, 0, max, +), которое обычно называют (max, +)-алгеброй. Нулевым элементом в этом полуполе является -те, единичным — число 0. Любому числу x е R можно сопоставить обратный x-1, который совпадает с противоположным числом -x в обычной алгебре. Для любой пары чисел x, y е R определена степень xy, значение которой равно арифметическому произведению xy. Заметим, что порядок, порожденный операцией идемпотентного сложения, совпадает с обычным линейным порядком на множестве R.
4. ЗАДАЧИ ТРОПИЧЕСКОЙ ОПТИМИЗАЦИИ С ОДНОЙ И ДВУМЯ ПЕРЕМЕННЫМИ
В этом разделе рассматриваются задачи оптимизации, которые формулируются и решаются в терминах тропической математики (задачи тропической оптимизации). Краткий обзор таких задач можно найти, например, в работе [27]. Ниже предлагается прямое решение задач оптимизации с одной и двумя переменными, заданных в терминах общего идемпотентного полуполя.
Предположим, что необходимо найти ненулевые решения u е X задачи
min au-1/2 Ф bu112 Ф cu-1 Ф du Ф e, (4)
ueX v '
где a, b, c, d, e — числа из X, удовлетворяющие условию a, b, c, d, e > O.
Для получения аналитического решения задачи в явном виде сначала вводится параметр, обозначающий минимум целевой функции. Затем задача сводится к нахождению всех решений параметризованной системы неравенств, условие существования которых используется для определения величины введенного параметра. Все решения рассматриваемой задачи описываются следующим утверждением.
Лемма 1. Минимум в задаче (4) равен л = a1/2b1/2 Ф a2/3d1/3 Ф b2/3c1/3 Ф c1/2d1/2 Ф e и достигается тогда и только тогда, когда
u = (л-2a2 Ф л^с)1-"^2^ Ф n~1d)-a, 0 < a < 1. (5)
При этом справедливы следующие утверждения:
1) если л = a1/2b1/2, то u = ab-1;
2) если л = a2/3d1/3, то u = a2/3d-2/3;
3) если л = b2/3c1/3, то u = b-2/3c2/3;
4) если л = с 1/2d1/2, то u = с 1/2d-1/2;
5) если л = e, то u = (a2e-2 Ф ce-1)1-a(b2e-2 Ф de-1)-a, для любого 0 < a < 1.
Доказательство. Обозначим минимум целевой функции в задаче (4) через л. Тогда все решения задачи определяются уравнением
au-1/2 Ф bu1/2 Ф cu-1 Ф du Ф e = л-
В силу того, что л — минимальное значение целевой функции, равенство можно заменить на неравенство
аи-1/2 Ф Ьи1/2 Ф си-1 Ф йиФ е < л-
Полученное неравенство эквивалентно системе
аы-1/2 < Ьи112 < си-1 < йи < л, е < ¡л. (6)
Заметим, что перемножение соответствующих частей первого и второго неравенств дает неравенство аЬ < л2, откуда с учетом условия леммы следует, что л > а1/2 Ь1/2 > О.
Решение относительно переменной и первых четырех неравенств (6) дает
2 2 2 2 1 1 и > л а , и < л Ь , и > л с, и < ¡й .
Объединим первое и третье неравенство в одно и получим и > л-2 а2 Ф л'1 с. Запишем второе и третье неравенства как неравенства относительно и-1, чтобы получить и-1 > л-2Ь2 и и-1 > л-1й. Затем объединим их в одно неравенство и-1 > л-2Ь2 Ф л-1й и перейдем к эквивалентному неравенству и < (л-2Ь2 Ф л-1й)-1.
Теперь для переменной и можно записать двойное неравенство
л-2а2 Фл'1 с < и < (л-2Ь2 Ф л-1й)-1. (7)
Множество значений и, удовлетворяющих этому неравенству, не пусто, если выполняется неравенство
(л-2 а2 Ф л-1 с) < (л-2Ь2 Ф л-1й)-1, которое после умножения обеих частей на л-2Ь2 Ф л-1й принимает форму
(л-2а2Ф л-1 с)(л-2Ь2Ф л-1й) < 1.
Раскроем скобки в левой части, а затем заменим полученное неравенство эквивалентной системой неравенств
л-4а2Ь2 < 1, л-3а2й < 1, л-3Ь2с < 1, л-2сй < 1.
Решая неравенства этой системы относительно л, получим
л > а1/2Ь1/2, л > а2/3й1/3, л > Ь2/3с1/3, л > с 1/2й1/2.
Полученная система, с учетом последнего неравенства из (6), равносильна неравенству
л > а1/2Ь1/2 Ф а2/3й1/3 Ф Ь2/3с1/3 Ф с 1/2й1/2 Ф е.
Заметим, что переход от исходной задачи к этому неравенству, в котором л обозначает минимум целевой функции, включал только эквивалентные преобразования. Из этого следует, что полученное неравенство задает точную нижнюю границу для л, а значит, знак неравенства можно заменить на знак равенства:
л = а1/2 Ь1/2 Ф а2/3 й1/3 Ф Ь2/3с1/3 Ф с1/2 й1/2 Ф е.
Двойное неравенство (7) для и можно записать с помощью параметра 0 < а < 1, представив его в виде (5).
Рассмотрим возможные значения л. Пусть л = а
1/21,1/2
После подстановки в формулу
(7) и применения свойств монотонности идемпотентного сложения получим
аЪ-1 < а-1Ъ-1 а2 Ф а-1/2Ъ-1/2с < и < (а-1Ъ-1Ъ2 Ф а-1/2Ъ-1/2й)-1 < аЬ-1,
откуда можно заключить, что и = аЪ-1.
Аналогичным образом можно проверить, что при л = а2/3й1/3 выполняется равенство и = а2/3й-2/3, при л = Ъ2/3с1/3 — равенство и = Ъ-2/3с2/3, а при л = с 1/2й1/2 — равенство и = с1/2й-1/2.
Приусловии л = е из (5) получаем, что и = (а2е-2 Фсе-1)1-а(Ъ2е-2 Фйе ~1)~а для любого 0 < а < 1. □
Теперь рассмотрим задачу тропической оптимизации с двумя переменными. Пусть заданы числа а, Ъ, с, й, е, /, g, Н, к > О. Требуется найти ненулевые решения и, V е X задачи
Ш1П аи
и^ еХ
1 Ф Ъv 1 Ф си Ф (IV Ф еи 1V 1 Ф/и 1V Ф guv 1 Ф Нuv Ф к.
(8)
Разделим решение задачи (8) на два этапа. Сначала решим ее относительно одной переменной, сводя к задаче (4). Затем рассмотрим несколько частных случаев и найдем значение второй переменной. Все решения задачи (8) сформулируем в виде теоремы.
Теорема 1. Введем обозначения:
а1 = Ъ1/2/1/2 Ф й1/2е1/2,
Ъ1 = Ъ1/2Н1/2 Ф й1/2g1/2,
С1 = е1/2/1/2 Ф а,
й1 = g1/2 Н1/2 Ф с,
е 1 = Ъ1/2й1/2 Ф /1/2g1/2 Ф е 1/2Н1/2 Ф к.
(9)
Тогда минимум в задаче (8) равен л = аЦ/2Ъ1/2 Ф а2/3й^/3 Ф Ъ2/3с^/3 ф с^12й^12 Ф е1, и справедливы следующие утверждения:
1) если л = а1/2Ъ1/2, то и = а1Ъ1 1 и
1
V =
|Ъ3/4 / - 1/4 ъ - 1/2, 13/4 е1/4 Ь1/2,
если а1 = Ъ1/2/1/2, если а1 = й 1/2е1/2;
2) если л = а2/3й! , то и = а2/3й- 2/3 и
Ъ2/3 ^-1/3 ^-1/3
V =
Ъ2/3 / -
й
1
-2/3е1/3й1/3
если а1 = Ъ1/2/1/2, если а1 = й 1/2е1/2;
3) если л = Ъ2/3с^/3, то и = Ъ-
,-2/3с2/3 и
V =
Ъ2/3Н-1/3с-1/3, если Ъ1 = Ъ1/2Н1/2, если Ъ1 = й1/2 g1/2;
й-2/3 g 1/3 с^3
4) если л = с1/2й^/2, то и = с^2^
1/2 1/2
V =
' е 1/2 /-1/2 g 1/2Н-1/2,
если с1 = е1/2 /1/2, если й1 = g1/2 Н1/2,
1/2с 1/2ЪФ а 1 еФ с 1 g)1 а®
Г1/2с-1/2й Ф а-1 /Ф с-1Н)-а, если с1 = а,й1 = с;
и
5) если л = e1t то u = (a^e-2 Ф c1e-1)1 a(b2e-2 Ф d^-1) a и
b1/2d-1/2, если e1 = b1/2d1/2,
f-1/2g 1/2u, если e1 = f1/2g1/2,
v = ^ e 1/2h-1/2u-1, если e1 = e 1/2 h1/2, (bk-1 Ф k~1et~1 Ф k-1 gt)1_e<8>
®(k-1d Ф k-1 ft-1 Ф k~1ht)~e, если e1 = k;
где а, в — вещественные числа, удовлетворяющие условиям 0 < a < 1, 0 < в < 1.
Доказательство. Пусть л — минимум целевой функции в задаче (8). Тогда все решения задачи определяются неравенством
au-1 Ф bv-1 Ф cuФ dvФ eu-1 v-1 Ф f u-1 vФ guv-1 Ф huv Ф k < л.
Это неравенство эквивалентно системе
au-1 < л, bv-1 < л, cu < л, dv < л, 1 i 1 1 (10) eu v < л, fu v < л, guv < huv < л, k <
Перемножение соответствующих частей первого и третьего неравенств дает ac < л2, откуда следует, что л > a1/2c1/2 > O.
Решим систему относительно переменной v, считая u параметром. Рассмотрим неравенства, в которых имеется переменная v, и представим их в виде
v > л~1Ь, v < лd-1, v >^^eu-1, v < лf"1u, v > л"1 gu, v < ^h^u-1.
Объединив вместе соответствующие неравенства, получим двойное неравенство
л-1(Ь Ф eu-1 Ф gu) < v < л^ Ф fu-1 Ф hu)-1. (11)
Для того, чтобы множество значений v было не пусто, необходимо выполнение условия
(b Ф eu-1 Ф gu)(d Ф fu-1 Ф hu) < л2.
Раскроем скобки слева и извлечем квадратный корень из обеих частей неравенства. После добавления тех неравенств из системы (10), которые не зависят от v, получим
л > (b1/2 f1/2 Ф d1/2e 1/2)u-1/2 Ф (b1/2h1/2 Ф d1/2g 1/2)u 1/2Ф
Ф (e1/2 f1/2 Ф a)u-1 Ф (g 1/2h1/2 Ф c)uФ b 1/2d1/2 Ф f1/2g1/2 Ф e 1/2h1/2 Ф k.
В силу того, что л является минимумом целевой функции, а переход от исходной задачи к последнему неравенству содержал только эквивалентные преобразования, полученное неравенство задает точную нижнюю границу для л, которая зависит от u. Рассмотрим задачу оптимизации
min (b1/2 f1/2 Ф d1/2e1/2) u-1/2 Ф (b1/2h1/2 Ф d1/2g1/2) u1/2 Ф (e1/2 f1/2 Ф a) u-1Ф
ueX
Ф^ 1/2h1/2 Ф c)uФ b 1/2d1/2 Ф f1/2g1/2 Ф e 1/2h1/2 Ф k.
(12)
С учетом обозначений (9) задача (12) принимает следующий вид:
min alu~1/2 Ф blu112 Ф c\u~l Ф dluФ ei.
ueX
Применение леммы 1 для решения задачи дает минимальное значение л = ai/2 bi/2 Ф af/3dl/3 Ф bf'cl13 Ф c^d^2 Ф eb причем выполняются условия:
1) если л = al/2bl'2, то u = aib-i;
2) если л = a2/3dl'3, то u = a2/3d-2/3;
3) если л = b2/3c l/3, то u = b-2/3c2/3;
4) если л = cl12dl/2, то u = cl/2d-l/2;
5) если л = el, то u = (a2e-2 Ф cle- l)l ~a(b2xe-2 Ф dle- l)-a, 0 < а < l .
Рассмотрим различные значения, которые может принимать величина л, и уточним полученное решение.
а .. - а 1/2 ь 1/2
4.1. Случай ц = aj/2bj
Предположим, что л = а{ Ъ1 . В этом случае и = а1Ъ-1 , а двойное неравенство (11) для V можно записать в виде
а-1 /2Ъ- 1/2(ЪФ а- 1 Ъ]_еФ а1Ъ- 1 g) < V < а^/2Ъ1 /2(й Ф а- 1Ъ1/Ф а1Ъ- 1Н)- 1.
Рассмотрим все значения, которые может принимать а1, определенное по формулам (9). Найдем соответствующие представления для V и л.
Пусть а1 = Ъ1/2/1/2. Двойное неравенство для V в этом случае преобразуется к виду
Ъ3/4/-1/4ъ-1/2 < Ъ3/4/-1/4Ъ-1/2 Ф Ъ-3/4/-3/4Ъ^/2еФ Ъ1/4/ 1/4Ъ-3/2g < V <
< (Ъ-1/4/-1/4Ъ-1/2й Ф Ъ-3/4/1/4Ь1/2 Ф Ъ1/4/ 1/4Ъ-3/2Н)-1 < Ъ3/4/-1/4Ъ-1/2, 3/4 1/4 1/2
откуда следует, что V = Ъ / 1/4Ъ-1/2.
Применяя такие же рассуждения, находим, что при а1 = й 1/2е1/2 выполняется равенство V = й-3/4 е1/4 Ъ^/2.
4.2. Случай ц = я;273^173
Допустим, что л = а2/3й^/3, тогда и = а2/3й-2/3, а двойное неравенство (11) для V принимает следующий вид:
а-2^^Ф а-2/3й2/3еФ а2/3йГ2/3^ < V < а2/3<3(йФ а-2/3й2/3/Ф а2^-2^)-1.
Рассмотрим различные значения, которые может иметь а1. Если а1 = Ъ1/2/1/2, то из двойного неравенства для V следует, что
2/3 1/3 1/3 2/3 1/3 1/3 2/3 2/3 1/3
b2/3 f-l/3di"l/3 < b2/3 f-l/3di"l/3 Ф b-2/3 f ~2/3dl/3e Ф d-1g < v <
< (b-l/3 f -l/3d-l/3d Ф b-2/3f l/3dl/3 Ф d-lh)i1 < b2/3f-l/3d-l/3,
или что v = b2/3f 1/3d-1/3.
Повторяя рассуждения для случая al = dl/2e1/2, получим равенство v = d 2/3el/3dl/3
4.3. Случай = Ъ1'3с113
Перейдем к рассмотрению случая л = Ь2/3 с^/3. При таком условии имеем и = Ь-2/3с2/3, а двойное неравенство (11) принимает вид
Ь-2/3с-1/3(ЬФ Ь2/3с-2/3еФ Ь-2/3с2/3я) < V < Ь2/3с^/3(й Ф Ь2/3с-2/3/Ф Ь-2/3с2/3Н)-1.
Все значения, которые может принимать Ь1, определяются формулами (9). Рассмотрим соответствующие представления для V и л.
Пусть Ь1 = Ь1/2 Н1/2, тогда двойное неравенство для V записывается так:
Ь2/3й-1/3с-1/3 < Ь2/3Н-1/3с-1/3 Ф с-1е Ф Ь-2/3Н-2/3с1/3я < V <
< (Ь-1/3^-1/3с-1/3йФ с-1 /Ф Ь-2/3Н 1/3с1/3)-1 < Ь2/3к-1/3с-1/3,
откуда следует, что V = Ь2/3^-1/3с-1/3.
Аналогично, при Ь1 = й1/2я1/2 получаем V = й-2/3я 1/3с^/3.
4.4. Случай ц = с}/24}/2
Предположим, что л = с^й^2. В этом случае и = с^й-1/2, а двойное неравенство (11) для V принимает вид
сГ1/2й-1/2(ЬФ сГ1/2й1/2еФ с1/2йГ1/2я) < V < с1/2й1/2(й Ф сГ1/2й1/2/Ф с^й-1^)-1.
Рассмотрим все значения, которые могут принимать с1 и й1, исходя из (9). Пусть с1 = е1/2/1/2. Двойное неравенство для V в этом случае преобразуется к виду
е1/2/-1/2 < е-1/4/-1/4й-1/2ЬФ е1/2/-1/2 Ф й--1 я < V <
< (е-1/4/-1/4й-1/2йФ е-1/2/1/2 Ф й-1«-1 < е1/2/-1/2,
откуда следует, что V = е1/2/-1/2.
В случае й1 = я 1/2Н1/2 аналогичные рассуждения дают V = я1/2 Н-1/2. Пусть с1 = а и й1 = с. Тогда двойному неравенству для V можно придать параметрическую форму
V = (а-1/2с-1/2ЬФ а-1еФ с-1 я)1-а(а-1/2с-1/2й Ф а-1 /Ф с-1Н)-а 0 < а < 1. 4.5. Случай ц = в\
Сначала предположим, что л = Ь1/2й1/2. Тогда двойное неравенство (11) для V можно записать следующим образом:
Ь-1/2й-1/2(ЬФ еи-1 Ф яи) < V < Ь1/2й 1/2(йФ /и-1 Ф Ни)-1.
После дальнейших преобразований получим
Ь1/2й-1/2 < Ь-1/2й-1/2ЬФ Ь-1/2й-1/2еи-1 Ф Ь-1/2й-1/2яи < V <
< (Ь-1/2й-1/2йФ Ь-1/2й-1/2/и-1 Ф Ь-1/2й-1/2Ни)-1 < Ь1/2й-1/2,
откуда вытекает равенство V = Ь1/2й-1/2.
Если л = /1/2g1/2, то, проводя аналогичные рассуждения, получим V = / 1/2g1/2 и, где
и = (а2/-1 g-1 Ф с1 /-1/2g-1/2)1-а(Ъ2/-1 g-1 Ф й1 /-1/2g-1/2)-а, 0 < а < 1.
При л = е1/2 Н1/2 приходим к равенству V = е 1/2Н-1/2и-1, где
и = (а2е-1Н-1 Ф с1е-1/2Н-1/2)1-а(Ъ2е-1Н-1 Ф й1е-1/2Н-1/2)-а, 0 < а < 1.
Пусть л = к. Тогда двойные неравенства (11) для х и (5) для и можно записать с использованием неравенств 0 < а < 1 и 0 < в < 1 в параметрической форме:
и = (а2к-2 Ф с1к-1)1-а(Ъ2к-2 Ф й1к-1)-а, V = (Ък-1 Ф к-1еи-1 Ф к-1 gu)1-в(k-1й Ф к-1 /и-1 Ф к-1Ни)-в.
Исследование этого случая завершает доказательство теоремы. □
Сформулируем полученный результат в компактной форме.
Следствие 1. Введем обозначения
а1 = Ъ1/2/1/2 Ф й1/2е1/2, Ъ1 = Ъ1/2Н1/2 Ф й1/2g1/2, с1 = е1/2/1/2 Ф а, й1 = g1/2Н1/2 Ф с, е 1 = Ъ1/2й1/2 Ф /1/2g1/2 Ф е 1/2Н1/2 Ф к.
Тогда минимум в задаче (8) равен л = а1/2Ъ^/2 Ф а2/3й^/3 Ф Ъ2/3с^/3 Ф с^12й^12 Ф е1 и достигается тогда и только тогда, когда
и = (л-2а2 Ф л-1с)1-а(л-2ЪФ л-1й)-а,
V = (л-1Ъ Ф л-1еи-1 Ф л-1 gu)1-а(л-1d Ф л-1 /и-1 Ф л-1Ни)-а, где а — вещественное число, удовлетворяющее условию 0 < а < 1.
Доказательство. При доказательстве теоремы 1 было показано, что множество решений (и, V) задачи (8) описывается при помощи двойного неравенства
л-1 ЪФ л-1еи-1 Ф л-1 gu < V < (л-1й Ф л-1 /и-1 Ф л-1Ни)-1,
где величина и является решением задачи (4).
Решение задачи (4) с помощью леммы 1 дает результат в параметрической форме (5), откуда следует утверждение следствия. □
5. ЗАДАЧА ТРОПИЧЕСКОЙ ОПТИМИЗАЦИИ С ТРЕМЯ ПЕРЕМЕННЫМИ
В этом разделе рассматривается задача тропической оптимизации с тремя переменными, для которой находится аналитическое решение в явном виде. Сначала вводится параметр для обозначения минимума целевой функции, а затем задача сводится к решению параметризованной системы неравенств относительно одной из неизвестных. Условия существования решений системы позволяют найти оптимальное значение для двух других неизвестных с помощью вспомогательной задачи, которая имеет вид (8). Решение такой задачи с использованием теоремы 1 дает величину введенного параметра. Полученное общее решение задачи преобразуется в набор частных прямых решений для различных соотношений между исходными параметрами.
Пусть заданы числа a, b, c, d, e, f, g, h e X и требуется найти все ненулевые решения u, v, w e X задачи
-1 -1 -1 I -1 -1 -1 -1 7 -1
min au v w Ф bu v w Ф cu vw Ф du vwФ
u,v,w eX
Ф euv- 1 w- 1 Ф fuv- 1 w Ф guvw- 1 Ф huvw. (13)
Следующий результат описывает все решения этой задачи.
Теорема 2. Введем обозначения:
a1 = a 1/2d1/2 Ф b1/2c1/2, b1 = a1/2 f1/2 Ф b1/2e1/2, C1 = e 1/2h1/2 Ф f1/2g1/2,
dl = c 1/2h1/2 Ф d1/2g1/2, e 1 = a1/2b1/2, f1 = c1/2d1/2, g1 = e1/2f1/2, (14)
h1 = g1/2h1/2, k1 = a1/2h1/2 Ф b1/2f1/2 Ф c1/2 f1/2 Ф d1/2e1/2,
a2 = b{12f1/2 Ф dl12e1/2, b2 = b^h^2 Ф d^2gl1/2, C2 = e}/2f1/2 Ф ab d2 = g\/2h 1 /2 Ф C1, e2 = b1 /2d1 /2 Ф f1 /2g1/2 Ф e1 /2h1 /2 Ф ^.
(15)
Тогда минимум в задаче (13) равен л = a1/2 b1/2 Ф a^l3d^13 фЬ^/3c^13 Ф c^l2d^12 Фe2, и справедливы следующие утверждения:
1) если л = a1/2b1/2, то u = a2b- 1,
_ i b3/4 f--1/4b- 1/2, если a2 = b1/2f1/2, v I df^e1^2, если a2 = d\/2e12,
r -1 -1 -1 -1 -1 -1 -1 -1 л1-а
w = (л 1 au 1 v 1 Ф л cu 1 v Ф л euv 1 Ф л 1 guv)1 a® ®(л - 1bu -1 v - 1 Ф л - 1du -1 v Ф л -1 fuv - 1 Ф л - 1huv)- a;
2) если л = a2/3di/3, то u = a?/3dr2/3
2 1^2 * — 2 2 *
_ i b2/3/f1/3d2-1/3, если a2 — b1/2f1/2,
V = 1 л-2/^1/^1/3 „„„,, „ _ J1/2„1/2
-1 -1 -1 -1 -1 -1 -1 -1 л1— a
w — (л 1 au 1 v 1 Ф л cu 1 v Ф л 1euv 1 Ф л 1 guv)1 ®(Л - 1bu -1 v - 1 Ф л - 1du -1 v Ф л -1 fuv 2 1 Ф л - 1huv)- a;
3) если л — b?/3c21/3; mo u — b2 2/3c2/3
2 ^2 * и*-*-' I* — >-'2 c'2 *
_ i b2/3h-1/3c2"1/3, если b2 — b1/2h1/2, v I d-2/3g1/3c21/3, если b2 — d1/2g1/2,
-1 -1 -1 -1 -1 -1 -1 -1 л1-а
w — (л 1 au 1 v 1 Ф л cu 1 v Ф л euv 1 Ф л 1 guv)1 a® ®(л - 1bu 21 v 2 1 Ф л - 1du 21 v Ф л -1 fuv 2 1 Ф л - 1huv)2 a;
1/2 1/2 1/2 1/2
4) если л — c1/2d1/2, mo u — c2d2
e1/2ff1/2, если C2 — e1/2f1/2,
1/2 1/2 1/2 1/2
v—
g^'2h^1/2, еслии d2 — g^h1
(a2_1/2c21/2b1 Ф a^1e1 Ф c^g1)12a®
2"1/2c21/2d1 Ф a2"1 f1 Ф c2_1h1)2a, если c2 — a1; d2 — c1;
w =
-1 -1 -1 -1 -1 -1 -1 -1 ч1-в
w = (л au v Ф л cu v Ф л euv Ф л guv)
®(Л~1Ьи-1 v-1 Ф ¡~ldu~l vФ ¡л'1 fuv-1 Ф лi~lhuv)~e;
5) если ц = e2, то u = (a^e-2 Ф c2e-1)1_a(b:2e-2 Ф d2e-1)-a,
(b\/2 d-1/2, если e2 = b1 /2 d1/2, ff1/2g\/2u, если e2 = f1/2gj/2, eJ/2h-1/2u-1, если e2 = e 1/2hj/2,
r -1 -1 -1 -1 -1 -1 -1 -1 ч1-в
w = (л au v Ф л cu v Ф л euv Ф л guv)
®(л_1 bu-1 v-1 Ф n~ldu~l v Ф л-1 fuv-1 Ф ¡л~1huv)~в;
6) если л = fc1, то u = (a^k-2 Ф c2 k-1)1-a(b2k-2 Ф d2 k-1)-a,
v = (b1fc-1 Ф k-1e 1 u-1 Ф k-1 g1u)1-e(fc-1d1 Ф k-1 f1u-1 Ф fc-1h1u)"e,
' a1/2h-1/2u-1 v-1, если k1 = a1/2h1/2,
c1/2f-1/2 u-1 v, если k1 = c1/2f1/2,
d-1/2e1/2uv-1, если k1 = d 1/2e1/2,
b-1/2g1/2 u v, если k1 = b1/2g1/2,
где a, в — вещественные числа, удовлетворяющие условиям 0 < a < 1, 0 < в < 1.
Доказательство. Применяя те же рассуждения, что и при решении задач с одной и двумя переменными, обозначим минимум целевой функции в задаче (13) через ц. Тогда все решения задачи определяются неравенством
1 1 1 1 1 1 1 1 au v w Ф bu v w Ф cu vw Ф du vwФ
1 1 1 1 Ф euv w Ф fuv w Ф guvw Ф huvw < л,
которое эквивалентно системе неравенств
1 1 1 1 1 1 1 1 au v w < л, bu v w < л, cu vw < л, du vw < л,
euv-1 w-1 < л, fuv-1 w < л, guvw-1 < л, huvw < ц.
Перемножение соответствующих частей первого и восьмого неравенств дает неравенство ah < л2, откуда с учетом условия леммы следует, что л > a1/2h1/2 > O. Решим систему относительно w, считая u и v параметрами:
-1 -1 -1 I -1 -1 -1 7-1 -1
w > л au v , w < лЬ uv, w > л cu v, w < лd uv ,
-1 -1 J—1 -1 -1 7-1-1-1
w > л euv , w < л1 u v, w > л guv, w < лh u v . Записывая правую и левую части неравенства для w, получим
Л~1(au-1 v-1 Ф cu-1 v Ф euv-1 Ф guv) < w < лЬ^1 v-1 Ф du-1 v Ф fuv-1 Ф huv)-1. (16) Множество значений w, которые удовлетворяют неравенству (16), не пусто, если (au-1 v-1 Ф cu-1 v Ф euv-1 Ф guv)(bu-1 v-1 Ф du-1 v Ф fuv-1 Ф huv) < л2-
Раскрывая скобки слева и извлекая квадратный корень, получим неравенство
л > (a1/2d1/2 Ф bl/2c1/2) u-1 Ф (al/2 f1/2 Ф bl/2e1/2) v-1 Ф (e1/2hl/2 Ф f1/2g1/2) uФ
Ф (cl/2h1/2 Ф d1/2gl/2)vФ al/2bl/2u-1 v-1 Ф cl/2dl/2u-1 vФ e1/2f l/2uv-1Ф
Ф g1/2hl/2uvФ al/2h1/2 Ф c1/2fl/2 Ф d1/2e1/2 Ф b1/2g1/2.
Так как л является минимумом целевой функции задачи, то как и при доказательстве теоремы 1, можно заключить, что полученное неравенство задает точную нижнюю границу для л- Следовательно, решение исходной задачи сводится к решению задачи
min (al/2d1/2 Ф bl/2cl/2)u-1 Ф (a1/2 f1/2 Ф bl/2e1/2) v-1 Ф (el/2hl/2 Ф f1/2g1/2)uФ u,v eX
Ф(cl/2hl/2 Ф dl/2gl/2)vФ al/2bl/2u-1 v-1 Ф cl/2dl/2u-1 vФ e1/2fl/2uv-^ (17) Фgl/2hl/2uvФ al/2hl/2 Ф c1/2fl/2 Ф dl/2e1/2 Ф bl/2g1/2.
С учетом обозначений (14) задача (17) может быть записана в более короткой форме
min alu-1 Ф blv-1 Ф cluФ dlvФ elu-1 v-1 Фflu-1 vФ gluv-1 Ф hluvФ kl.
u , ve X
Решение такой задачи предложено в теореме 1. С учетом обозначений (15), минимум в задаче (13) равен л = a2/2b2/^a2/3d2/^b2/3c;i/^c;i/2d2/^e2, и справедливы следующие утверждения:
1) если л = ai/2bl/2, то u = a2b-1 и
,3/4 -1/4 b-i/2
v
b3/4 fi
d-3/4e 1 14 b^12,
ест a2 = bl fi ,
e /
/2 /2 ест a2 = d i ei
2) если л = a?/3dl /3, то u = a2/3d.
2/3 2/3
и
v
_ i b2/3f2- 1 /3d.- 1 /3 1d-
2/3e 1/3 d 1/3 el d2 ,
3) если л = b?/3ci/3, то u = b.
2/3 2/3
2 c2 и
v
rb2/3hl-l/3c2-l/3, d-2/3gl/3cl/3,
d1 g1 c2 ,
4) если л = c2/2d^/2, то u = c2/2d-l/2 и 1/2 -1/2
если a2 = bi/2f2l/2, если a2 = d^l2e^12;
если b2 = b11/2h11/2,
если b2 = dl'2gl /2;
v=
ei/2 fi gl/2 h -1 /2,
(a-1/2 c- 1 /2b2 Ф a- 1 e2 Ф c- 1 g2)l-a<s ®(a-l/2Cl- 1 /2di Ф a- 1 fi Ф c- 1 hi)-a
5) если л = e2, то u = (a2,e22 Ф c2e2 2)2 a(b;2e22 Ф d2e2 2) b^d- 1 /2,
fi- 1 /2 gl/2u, v = <| e-/2h -l/2u- 1,
(blk- 1 Ф k- 1 e2u- 1 Ф k- 1 g2 u) »(kj'd Ф k-1 f2u-1 Ф k-1 h2u)-^
1 -a
если с2 = е^/2/11/2, если й2 = g^l2Н^l2,
если с2 = а1, й2 = с1; и
если е2 = Ъ1/2й1/2, если е2 = /11/2gl1/2, если е2 = еЦ/2НЦ/2,
если е2 = к1;
где а, в — вещественные числа, удовлетворяющие условиям 0 < а < 1 и 0 < в < 1 •
Во всех случаях, кроме последнего, формулу для w в виде двойного неравенства (16) запишем с использованием параметра 0 < в < 1 следующим образом:
/ -1 -1 -1 -1 -1 -1 -1 -1 ч1-в
w = (л au v Ф ц cu v Ф ц euv Ф ц guv)
® (¡~1bu~1 v-1 Ф ¡-1du-1 v Ф ¡-1 fuv-1 Ф ¡Г1Н^)~в.
Покажем, что в этом случае w может быть записано без использования дополнительного параметра. Действительно, пусть ц = k1 = a1/2h1/2. Тогда двойные неравенства (16) для w можно записать следующим образом:
a 1/2h-1/2u-1 v-1 < a-1/2h-1/2au-1 v-1 Ф a-1/2h-1/2cu-1 vФ
Ф a~ll2h~ll2euv-1 Ф a-1/2h-1/2guv < w < (a-1/2h-1/2bu-1 v-1 Ф a-1/2h-1/2du-1 vФ
Ф a-1/2h-1/2 fuv-1 Ф a~ll2h~ll2huv )-1 < a 1/2h-1/2u-1 v-1,
что означает w = a1/2h-1/2u-1 v-1.
Аналогичным путем устанавливаем, что при ц = k1 = c1/2f1/2, выполняется равенство w = c1/2f-1/2u-1 v, при л = k1 = d1/2e1/2 получаем равенство w = d-1/2e 1/2uv-1, при Л = k1 = b1/2g1/2 - равенство w = b-1/2g1/2 uv.
Рассмотрение этих случаев завершает доказательство. □
6. ПРИЛОЖЕНИЕ К ЗАДАЧЕ РАЗМЕЩЕНИЯ
Рассмотрим задачу размещения (3) и представим ее в терминах (max, +)-алгебры. Расстояние между двумя векторами в прямоугольной метрике (1) с использованием операций идемпотентного полуполя Rmax,+ записывается так:
p(x, у) = (х-1 У1 Ф у-1 Х1)(Х-1 У2 Ф У-1 Х2)(Х-1 уэ Ф у-1 Х3).
Для заданного набора точек тг- = (г1г-, г2г-, Г3г-)T е R3 и чисел ¡i е R, где i = 1,..., m, целевая функция (2) принимает вид
m
ф(Х1, Х2, Хэ) = 0 ¡i (Х-1 Гц Ф Г- Х1)(Х-1 r*2i Ф Г- Х2)(Х-1^ Ф Г- Х3). i = 1
Раскроем скобки и введем дополнительные обозначения
m m m m
a = 0 ¡irЦГ2г'Г3i, b = 0 ¡iГ1 iГ2iГ-, c = 0 ¡1Т11Т2-1Г3,-, d = 0 Ьгцг^г-1,
i=1 i=1 i=1 i=1
m m m m
(18)
Ф^' 1 / Г2 i Г3 i, f = 0 ¡i^-1 Г2гГ3г1, g = 0 ¡i Г11Г21Г3 i, h = 0 ¡i Г1 / Г21Г3/' i 1 i 1 i 1 i 1
Перегруппировав слагаемые, запишем целевую функцию так:
ф(х\, х2,х3) = ах- 1 х- 1 х- 1 Ф Ьх- 1 х-1 х3 Ф сх- 1 х2х- 1 Ф dx- 1 х2х3Ф
Ф ех1 х- 1 х- 1 Ф fx1 х- 1 х3 Ф gx1 х2х-1 Ф кх1 х2х3.
Задача (3) принимает форму задачи (8), где и = х1, V = х2 и ш = х3. Применение теоремы 2 дает следующий результат.
Теорема 3. В дополнение к (18) введем обозначения:
а1 = а1/2й1/2 Ф Ъ1/2с1/2, Ь1 = а1/2/1/2 Ф Ъ1/2е1/2, с1 = е 1/2Н1/2 Ф /1/2g1/2, й1 = с 1/2Н1/2 Ф й1/2g1/2, е 1 = а1/2Ъ1/2, /1 = с1/2й1/2, gl = е1/2/1/2, Н1 = g1/2й1/2, к1 = а112Н112 Ф Ъ1/2/1/2 Ф с1/2/1/2 Ф й1/2е1/2, а2 = Ъ1/2/1/2 Ф й1/2е1/2, Ъ2 = Ь^2^2 Ф й1/2gl1/2, с2 = е1/2/1/2 Ф а1, й2 = gl1/2h1/2 Ф с1, е2 = Ъ1/2<2 Ф /1/2gl1/2 Ф е1/2й1/2 Ф кь
Тогда минимум в задаче (3) равен л = а^/2Ъ^/2 Ф а^13й^13 Ф Ъ|/3с^/3 Ф с^/2й^12 Ф е2 и справедливы следующие утверждения:
если л а1/2Ъ1/2, то
1) если л = а1/2Ъ1/2, то х1 = а2Ъ-1,
_ ] Ъ3/4/1-1/4Ъ2-1/2, если а2 = Ъ1/2/1/2, Х2 | й-3^1/4Ъ1/2, если а2 = ¿1/2е1/2;
-1 -1 -1 -1 -1 -1 -1 -1 л1-а
х3 = (л ах1 х2 Ф л сх1 х2 Ф л ех1 х2 Ф л gx1 х2) ® ф^^х-1х-1 Ф л^йх-1 х2 Фл-1 /х1 х-1 Ф л~1Ьх1 х2)-а;
2) если л = а?/3й2/3, то х1 = а2/3й2 2/3,
_ ] Ъ2/3/1-1/3й2-1/3, если а2 = Ъ1/2/1/2, х2 |dl-2/3el1/3d21/3, если а2 = й\/2е\/2;
г -1 -1 -1 -1 -1 -1 -1 -1 л1-а
х3 = (л ах1 х2 Ф л сх1 х2 Ф л ех1 х2 Ф л gx1 х2) ® ф^^х-1х-1 Ф л-1йх-1 х2 Ф л-1 /х1 х-1 Ф л~1Ьх1 х2)-а;
3) если л = Ъ2/3с21/3, то х1 = Ъ- 2/3с2/3,
ÍЪ2/3hl-1/3c2-1/3, если Ъ2 = Ъ1/2й1/2, I й^3gl1/3c21/3, если Ъ2 = й1/2gl1/2;
-1 -1 -1 -1 -1 -1 -1 -1 л1-а
х3 = (л ах1 х2 Ф л сх1 х2 Ф л ех1 х2 Ф л gx1 х2) ® ф^^х-1х-1 Ф л-1йх-1 х2 Ф л-1 /х1 х-1 Ф л~1Ьх1 х2)-а;
4) если л = с^/2й^/2, то х1 = с^/2й2 1/2,
х2 =
е1/2/-1/2, если с2 = е1/2/1/2,
g1/2 Й -1/2, если й2 = gl1/2 й1/2, (а1"1/2c2l/2Ъ1 Ф а-1 е1 Ф с-1 g1)1-а<8>
®(a1"1/2с-1''2^ Ф а-1 /1 Ф с-1^)^, если с2 = а1, й2 = с1;
-1 -1 -1 -1 -1 -1 -1 -1 л1-в
х3 = (л ах1 х2 Фл сх1 х2 Ф л ех1 х2 Ф л gx1 х2)
»(л^Ъх-1 х2-1 Ф л-1^-1 х2 Ф л-1 /х1 х2-1 Ф л~1Ьх1 х2)-в;
5) если л = e2, то Х1 = (a|e2 2 Ф c2e2 1)1 a(b|e2 2 Ф d2e2 1) а,
г -1 -1 -1 -1 -1 -1 -1 -1 ч1-в
Х3 = (л ax1 Х2 Ф л cx1 Х2 Ф л ex1 Х2 Ф л gx1 Х2)
®(Л 1bx- 1 x2 1 Ф л 1dx- 1 x2 Ф л 1 fx1 x2 1 Ф л 1hx1 x2) в;
6) если л = e2 = k1; то x1 = (a|k- 2 Ф c2k- 1)1 a(b|k- 2 Ф d2k- 1) а,
x2 = (b1k- 1 Ф k- 1e1 Х- 1 Ф k- 1 g1 x1)1 1 d1 Ф k- 1 f1 x- 1 Ф k- 1h1 x1) в,
Х3 =
Х3 | d-1/2e1/2x1 x-1, если k1 = d 1/2e1/2, b-1/2g1/2x1 x2, если k1 = b1/2g1/2;
• a1/2h-1/2x-1 x-1, если k1 = a 1/2h1/2, c1/2f-1/2x-1 x2, если k1 = c1/2f1/2,
где а, в — вещественные числа, удовлетворяющие условиям 0 < а < 1, 0 < в < 1-
Сформулируем результат в терминах обычной математики в виде следствия. Следствие 2. Введем обозначения:
a = max (¡i + Гц + Г2{ + Г3i), b = max (¡i + Гц + Г2г- - Г3г-), c = max (¡i + Г1г- - Г2г- + Г3г-),
1<i <m 1<i <m 1<i<m
d = max (¡i + Г1г- - Г2г- - Г3г-), e = max (¡i - Гц + Г2г- + Г3г-), f = max (¡i - Гц + Г2г- - Г3{),
1<i<m 1<i<m 1<i<m
g = max (¡i - Г 1i - Г2i + Г3г-), h = max (¡i - Гц - Г2г- - Г3г-),
1<i<m 1<i<m
a1 = max(a + d, b + c)/2, b1 = max(a + f, b + e)/2, c1 = max(e + h, f + g)/2, d1 = max(c + h, d + g)/2, e1 = (a + b)/2, f1 = (c + d)/2, g1 = (e + f)/2, h1 = (g + h)/2, k1 = max(a + h, b + g, c + f, d + e)/2, a2 = max(b1 + f1; d1 + e1)/2, b2 = max(b1 + h1t d1 + g1)/2, c2 = max((e 1 + f1)/2, a1), d2 = max((g1 + h1)/2, c1), e2 = max(b1 + d\, f1 + gb e 1 + h1,2k1)/2.
Тогда минимум в задаче (3) равен л = max((a2 + b2)/2,(2a2 + d2)/3,(2b2 + c2)/3, (c2 + d2)/2, e 2) и справедливы следующие утверждения:
1) если л = (a2 + b2)/2, то
x3 = (1 - а) max(^ + a - x1 - x2, -л + c - x1 + x2,
- л + e + x1 - x2, -л + g + x1 + x2) - а max(-л + b - x1 - x2,
- л + d - x1 + x2, -л + f + x1 - x2, -л + h + x1 + x2);
Х1 = a2 - b2,
Х2 =
(3b1 - f1)/4 - b2/2, если a2 = (b1 + f1)/2, (-3d1 + e1)/4 + b2/2, если a2 = (d1 + e1)/2;
2) если ¡1 = (2а2 + d2)/3, то
х1 = (2а2 - 2d2)/3,
[ (2Ь1 - ^ - d2)/3, если а2 = (Ь1 + /0/2,
х2 =
[(-2d1 + е1 + d2)/3, если а2 = (d1 + е1)/2;
х3 = (1 - а) шах(-^ + а - х1 - х2, + с - х1 + х2,
- ¡1 + е + х1 - х2, + g + х1 + х2) - ашах(-^ + Ь - х1 - х2,
- ¡1 + d - х1 + х2, + f + х1 - х2, + Н + х1 + х2);
3) если ¡1 = (2Ь2 + с2)/3, то
х1 = (-2Ь2 + С2)/3,
[ (2Ь1 - Н1 - с2)/3, если Ь2 = Ь1/2Н1/2,
х2 = 1 1
[ (-2d1 + g1 + с2)/3, если Ь2 = (d1 + g1)/2;
х3 = (1 - а) шах(-^ + а - х1 - х2, + с - х1 + х2,
- ¡1 + е + х1 - х2, + g + х1 + х2) - ашах(-^ + Ь - х1 - х2,
- ¡1 + d - х1 + х2, + f + х1 - х2, + Н + х1 + х2);
4) если ¡1 = (с2 + d2)/2, то
х1 = (С2 - d2)/2,
(е1 -/!)/2, если с2 = (е 1 + ^/2,
- Н1)/2, если d2 = (g1 + Н1)/2,
(1 - а) шах((-а1 - с1)/2 + Ь1; -а1 + е 1,-С1 + g 1)--а шах((-а1 - с1)/2 + d1,
-а1 + /1, -С1 + Н1), если С2 = а1, d2 = С1;
х3 = (1 - в) шах(-^ + а - х1 - х2, + с - х1 + х2,
- ¡1 + е + х1 - х2, + g + х1 + х2) - вшах(-^ + Ь - х1 - х2,
- ¡1 + d - х1 + х2, + f + х1 - х2, + Н + х1 + х2);
5) если ¡1 = е2, то
х1 = (1 - а)шах(2а2 - 2е2, с2 - е2) - ашах(2Ь2 - 2е2, d2 - е2),
{(Ь1 - d1)/2, если е2 = (Ь1 + d1)/2,
(-Л + gl)/2 + хь если е2 = (Л + gl)/2, (е1 - Н1)/2 - х1, если е2 = (е1 + Н1)/2;
х3 = (1 - в) шах(-^ + а - х1 - х2, + с - х1 + х2,
- ¡1 + е + х1 - х2, + g + х1 + х2) - вшах(-^ + Ь - х1 - х2,
- ¡1 + d - х1 + х2, + а + х1 - х2, + Н + х1 + х2);
х2 =
6) если j = e2 = ki, то
xi = (1 - a)max(2a2 - 2ki, c2 - k1) - amax(2b2 - 2ki, d2 - k1), x2 = (1 - ^)max(b1 - ki, -k1 + e1 - xi, -k1 + g1 + x1)--вmax(-k1 + di, -k1 + f1 - xi, -k1 + h1 + x1), (a - h)/2 - xi - x2, если ki = (a + h)/2, (c -f)/2 -xi + x2, если ki = (c + f)/2,
x3 =
(-d + e)/2 + xi -x2, если ki = (d + e)/2, (-b + g)/2 + xi + x2, если ki = (b + g)/2;
где a, в — вещественные числа, удовлетворяющие условиям 0 < a < i, 0 < в < i.
Оценим вычислительную сложность полученного решения задачи размещения. Ясно, что трудоемкость расчетных формул следствия 2 определяется сложностью вычисления параметров a, b, c, d, e, f, g и h. Для вычисления каждого такого параметра необходимо выполнение m операций определения максимума величин, полученных за фиксированное число сложений, что дает линейную по m сложность. Применение других формул требует конечного числа операций, которое не зависит от числа точек m.
7. ЗАКЛЮЧЕНИЕ
В работе рассмотрена минимаксная задача размещения точечного объекта в трехмерном пространстве с прямоугольной метрикой. Такая задача встречается на практике, например, при размещении центра управления системой камер видеонаблюдения в здании. На основе применения методов тропической оптимизации получено полное прямое решение задачи, представленное в форме, которая является удобной для дальнейшего анализа решений и пригодной для непосредственных расчетов с невысокой вычислительной сложностью. Дальнейшие исследования могут включать решение задач размещения в трехмерном пространстве с дополнительными ограничениями на допустимую область размещения, а также рассмотрение задач с произвольной размерностью пространства размещения.
Список литературы
1. Sule D. R. Logistics of Facility Location and Allocation / New York: Marcel Dekker, 2001.
2. Farahani R. Z., Hekmatfar M. Facility Location / Contributions to Management Science. Heidelberg: Physica-Verlag, 2009.
3. Eiselt H. A., Marianov V. Foundations of Location Analysis / New York: Springer, 2011. Vol. 155 of International Series in Operations Research and Management Science.
4. Krivulin N. Using tropical optimization to solve constrained minimax single-facility location problems with rectilinear distance // Computational Management Science, 2017. Vol. 1, № 4. P. 493-518.
5. Francis R. L. A geometrical solution procedure for a rectilinear distance minimax location problem // AIIE Transactions, 1972. Vol. 4. № 4. P. 328-332.
6. Hansen P., Peeters D., Thisse J. F. Constrained location and the Weber-Rawls problem // North-Holland Mathematics Studies, 1981. Vol. 59. P. 147-166.
7. Zabudskii G. G. A minimax planar location problem with forbidden zones: its solution algorithm // Automation and Remote Control, 2004. Vol. 65. № 2. P. 241-247.
8. Nobakhtian S., Dehkordi A. R. A fast algorithm for the rectilinear distance location problem // Mathematical Methods of Operations Research, 2018. Vol. 87. № 1. P. 1-18.
9. Cuninghame-Green R. A. Minimax algebra / Berlin: Springer-Verlag, 1979.
10. Baccelli F., Cohen G., Olsder G. J., Quadrat J. P. Synchronization and Linearity. Wiley Series in Probability and Statistics / Chichester: Wiley, 1993.
11. Маслов В. П., Колокольцов В. Н. Идемпотентный анализ и его применение в оптимальном управлении / М.: Физматлит, 1994.
12. Cuninghame-Green R. A. Minimax algebra and applications // Advances in Imaging and Electron Physics, Vol. 90 / Ed. by P. W. Hawkes. San Diego: Academic Press, 1994. P. 1-121.
13. Golan J. S. Semirings and Affine Equations Over Them / New York: Springer, 2003. Vol. 556 of Mathematics and Its Applications.
14. Heidergott B, Olsder G. J., van der Woude J. Max Plus at Work / Princeton Series in Applied Mathematics. Princeton: Princeton Univ. Press, 2006.
15. McEneaney W. M. Max-Plus Methods for Nonlinear Control and Estimation / Systems and Control: Foundations and Applications. Boston: Birkhauser, 2006.
16. ItenbergI., Mikhalkin G., ShustinE. I. Tropical Algebraic Geometry / Basel: Birkhauser, 2007. Vol. 35 of Oberwolfach Seminars. doi: 10.1007/978-3-7643-8310-7.
17. Gondran M., Minoux M. Graphs, Dioids and Semirings / New York: Springer, 2008. Vol. 41 of Operations Research/ Computer Science Interfaces.
18. Кривулин Н. К. Методы идемпотентной алгебры в задачах моделирования и анализа сложных систем / СПб.: Издательство Санкт-Петербургского государственного университета, 2009. 255 с.
19. Butkovic P. Max-linear Systems. Springer Monographs in Mathematics / London: Springer, 2010.
20. Кривулин Н. К., Агеев В. А., Гладких И. В. Применение методов тропической математики для оценки альтернатив на основе парных сравнений // Вестник Санкт-Петербургского университета. Прикладная математика. Информатика. Процессы управления, 2017. № 1. С. 27-41.
21. KrivulinN. K. Tropical optimization problems in time-constrained project scheduling // Optimization, 2017. Vol. 66. № 2. P. 205-224.
22. Krivulin N. K. An extremal property of the eigenvalue of irreducible matrices in idempotent algebra and solution of the Rawls location problem // Vestnik St. Petersburg Univ. Math, 2011. Vol. 44. № 4. P. 272-281.
23. Krivulin N. K. A new algebraic solution to multidimensional minimax location problems with Chebyshev distance // WSEAS Transactions on Mathematics, 2012. Vol. 11. № 7. P. 605-614.
24. Krivulin N. K. Complete solution of a constrained tropical optimization problem with application to location analysis // Relational and Algebraic Methods in Computer Science / Cham: Springer, 2014. Vol. 8428 of Lecture Notes in Computer Science. P. 362-378.
25. KrivulinN. K., Plotnikov P. V. On an algebraic solution of the Rawls location problem in the plane with rectilinear metric // Vestnik St. Petersburg University: Mathematics, 2015. Vol. 48. № 2. P. 75—81.
26. Krivulin N.K., Plotnikov P.V. Using tropical optimization to solve minimax location problems with a rectilinear metric on the line // Vestnik St. Petersburg University: Mathematics, 2016. Vol. 49. № 4. P. 340-349.
27. Krivulin N. A multidimensional tropical optimization problem with nonlinear objective function and linear constraints // Optimization, 2015. Vol. 64. № 5. P. 1107-1129.
28. Krivulin N. Extremal properties of tropical eigenvalues and solutions to tropical optimization problems // Linear Algebra and its Applications, 2015. Vol. 468. P. 211-232.
29. Krivulin N. Direct solution to constrained tropical optimization problems with application to project scheduling// Computational Management Science, 2017. Vol. 14. № 1. P. 91-113.
Поступила в редакцию 19.01.2018, окончательный вариант — 22.02.2018.
Computer tools in education, 2018 № 1:31-50
http://ipo.spb.ru/journal
DIRECT SOLUTION OF MINIMAX LOCATION PROBLEM IN THREE-DIMENSIONAL SPACE WITH RECTILINEAR METRIC
Plotnikov P. V.1, Krivulin N. K.1 1 Saint Petersburg State University, Saint Petersburg, Russia
Abstract
The minimax problem of placing a point object in a three-dimensional space with a rectangular metric (/^metric) is considered and its direct analytical solution is proposed using the methods of tropical (idempotent) mathematics. First, the problem is written in terms of tropical mathematics as a problem of tropical optimization, a parameter is introduced to denote the minimum of the objective function and the problem reduces to solving a parametrized system of inequalities. This system is solved with respect to one of the variables, and the conditions for the existence of solutions are used to find the optimal values of the second variable using an auxiliary optimization problem. Then the auxiliary problem is solved in a similar way and the value of the third variable is found. The obtained general solution is transformed into a set of direct solutions written in a compact form for different cases of relationships between the initial parameters of the problem.
Keywords: 1-center problem, three-dimensional space, rectilinear metric, idempotent semifield, tropical optimization, complete solution.
Citation: P. V. Plotnikov, N. K. Krivulin, "Direct Solution of Minimax Location Problem in Three-Dimensional Space with Rectilinear Metric," Computer tools in education, no. 1, pp. 31-50, 2018 (in Russian).
Received 19.01.2018, The final version — 22.02.2018.
Pavel V. Plotnikov, postgraduate student, Department of Statistical Modeling, SPbGU, [email protected]
Nikolai K. Krivulin, Dr. of physical and mathematical sciences, professor, SPbGU; 199034, Russian Federation, Saint Petersburg, Universitetskaya nab. 7-9, [email protected]
© Наши авторы, 2018. Our authors, 2018.
Плотников Павел Владимирович, аспирант кафедры статистического моделирования СПбГУ, [email protected]
Кривулин Николай Кимович, доктор физико-математических наук, профессор кафедры статистического моделирования СПбГУ; 199034, Россия, Санкт-Петербург, Университетская наб. 7-9, [email protected]