CC BY
26
УДК 519.85 Вестник СПбГУ. Математика. Механика. Астрономия. 2018. Т. 5 (63). Вып. 1 МБС 65К10, 15А80, 90С48, 65Р50
О.
решении одной многомерной задачи тропической оптимизации с использованием разрежения матриц*
Н. К. Кривулин, В. Н. Сорокин
Санкт-Петербургский государственный университет,
Российская Федерация, 199034, Санкт-Петербург, Университетская наб., 7—9
Для цитирования: Кривулин Н. К., Сорокин В. Н. О решении одной многомерной задачи тропической оптимизации с использованием разрежения матриц // Вестник СПбГУ. Математика. Механика. Астрономия. 2018. Т. 5(63). Вып. 1. С. 91-104. https://doi.org /10.21638/11701^рЬи01.2018.110
В работе предлагается полное решение задачи минимизации функции, заданной на векторах с элементами из тропического (идемпотентного) полуполя. Рассматриваемая задача тропической оптимизации возникает, например, когда требуется найти наилучшее, в смысле метрики Чебышёва, приближенное решение тропических векторных уравнений, и встречается в различных приложениях, включая задачи планирования, размещения и принятия решений. Для решения задачи сначала находится минимальное значение целевой функции, предлагается описание множества решений в форме системы неравенств и приводится одно из решений. Далее с помощью разрежения матрицы задачи строится расширенное множество решений, а затем полное решение в виде некоторого семейства подмножеств. Описываются процедуры, позволяющие сократить число подмножеств, которые необходимо исследовать при построении полного решения. Показано, как полное решение задачи может быть записано в параметрическом виде в компактной векторной форме. Полученное решение обобщает известные результаты, которые обычно сводятся к получению одного из решений и не позволяют найти все множество решений задачи. Для иллюстрации основных результатов работы приводится пример численного решения задачи на множестве трехмерных векторов.
Ключевые слова: идемпотентное полуполе, тропическая оптимизация, чебышёвская аппроксимация, полное решение, разрежение матриц.
1. Введение. Тропическая (идемпотентная) алгебра представляет собой раздел математики, который изучает свойства полуколец и полуполей с идемпотентным сложением и их приложения. Исследованиям в этой области посвящено немало работ, включая [1-7], а также детальный обзор литературы в [8]. Использование моделей и методов тропической алгебры позволяет ряд нелинейных в обычном смысле задач приводить к линейной форме в терминах некоторого идемпотентного полукольца или полуполя, что во многих случаях облегчает решение и упрощает интерпретацию результатов. Такой подход находит все более широкое применение при решении практических задач в различных областях, включая задачи планирования, размещения, принятия решений и др.
Одним из направлений развития тропической математики является разработка методов решения задач оптимизации, которые могут быть сформулированы и решены
* Работа выполнена при финансовой поддержке РГНФ в рамках научного проекта № 16-02-00059. (¡5 Санкт-Петербургский государственный университет, 2018
в терминах идемпотентной математики (задач тропической оптимизации). Имеется целый ряд практических задач (см., например, [9-13]), которые сводятся к наилучшему приближенному решению в смысле метрики Чебышёва векторного уравнения Ax = p, где A и p обозначают заданные матрицу и вектор, x — неизвестный вектор, а произведение матрицы на вектор понимается в смысле тропической математики.
Рассмотрим, например, задачу оптимального планирования времени начала работ некоторого проекта с целью обеспечения заданных директивных сроков завершения работ. В этой задаче берется идемпотентное полуполе с операцией max в роли сложения и арифметического сложения в роли умножения, координаты вектора x имеют смысл времени начала выполнения для каждой работы, которое требуется определить, вектор p задает директивные сроки завершения работ, матрица A — ограничения снизу на величину интервалов времени между началом и завершением работ.
Проблема чебышёвской аппроксимации для решения рассматриваемого уравнения сводится к нахождению векторов x, на которых достигается минимум в задаче
min (Ax)-p ф p- Ax,
где матричные и векторные операции понимаются в смысле идемпотентной алгебры.
Исследованию задачи посвящен ряд работ, опубликованных в различное время, включая [3, 10-13]. Представленные в этих работах результаты обычно описывают получение одного из решений и не позволяют найти все множество решений задачи.
Для нахождения решения задачи в статьях [14-17] предлагается подход, при котором вводится дополнительный параметр для обозначения минимума целевой функции, а затем задача сводится к решению параметризованных неравенств. С помощью такого подхода были получены решения в компактной векторной форме для рассматриваемой задачи, а также некоторых ее вариантов, включая задачи с ограничениями.
Цель настоящей работы состоит в исследовании и решении задачи тропической оптимизации с целевой функцией более общего вида. На основе применения методов решения с использованием разреженных матриц, описанного в [18], находится полное решение задачи и его представление в компактной векторной форме.
Работа построена следующим образом. Сначала находится минимальное значение целевой функции задачи, предлагается описание множества решений в форме системы неравенств и приводится одно из решений. Далее с помощью разрежения матрицы задачи находится расширенное множество решений, а затем полное решение в виде некоторого семейства подмножеств. Предлагаются процедуры, позволяющие сократить число подмножеств, которые необходимо исследовать при построении полного решения. Показано, как полное решение задачи может быть записано в параметрическом виде в компактной векторной форме. Для иллюстрации полученных результатов приводится пример численного решения задачи на множестве трехмерных векторов.
2. Основные понятия тропической алгебры. Рассмотрим необходимые в дальнейшем основные понятия и предварительные результаты тропической (идемпотентной) алгебры [1-7]. Идемпотентное полуполе определяется как набор (X, ф, ©, 0,1), где X обозначает непустое множество с заданными на нем ассоциативными и коммутативными операциями сложения ф с нейтральным элементом 0 (ноль) и умножения © с нейтральным элементом 1 (единица).
Сложение обладает свойством идемпотентности: для любого элемента x G X выполняется равенство x ф x = x. Умножение дистрибутивно относительно сложения и обратимо: для любого x = 0 существует x-1 такой, что x © x-1 = 1.
Идемпотентность сложения задает на множестве X частичный порядок: x < y тогда и только тогда, когда x ф y = y. Отсюда следует, что неравенство x ф y < z эквивалентно двум неравенствам x < z и y < z. Кроме того, операции сложения и умножения являются монотонными, а операция обращения — антитонной. Будем предполагать, что указанный частичный порядок продолжен до линейного на множестве X. Ниже операции max и min понимаются в смысле указанного линейного порядка.
Целые степени определяются обычным образом: x0 = 1, xp = x © xp-1, x-p = (x-1)p, 0p = 0 для всех x = 0 и натуральных p. В дальнейшем будем считать полуполе алгебраически полным в том смысле, что введенная операция возведения в степень может быть распространена на случай степеней с рациональным показателем. Далее символ © для упрощения записи будем опускать.
В прикладных задачах часто встречаются следующие вещественные идемпотент-ные полуполя: Rmax,+ = (RU{—to}, max, +, —то, 0), Rmin,+ = (RU{+to}, min, +, +то, 0), Rmax,x = (R+ U {0}, max, x, 0,1) и Rmin,x = (R+ U {+то}, min, x, +то, 1). Здесь R+ = {x G R|x > 0} обозначает множество положительных вещественных чисел.
Рассмотрим полуполе Rmax,+ , которое обычно называют (max, +)-алгеброй. В нем множество X определено как множество действительных чисел R, расширенное путем добавления числа —то. Роль тропического сложения играет операция взятия максимума, а в качестве умножения берется арифметическое сложение. Для любого x = 0 существует обратный по умножению элемент x-1, равный противоположному числу —x в обычной арифметике. Степень xy определена для всех x, y G R и соответствует арифметическому произведению xy. Порядок, индуцированный идемпотентным сложением, совпадает с обычным линейным порядком на R.
Обозначим через Xmxn множество матриц размера m x n над полуполем X.
Нулевая матрица имеет все элементы равными 0. Матрицу, в которой отсутствуют нулевые строки (столбцы), назовем регулярной по строкам (столбцам). Матрица без нулевых строк и нулевых столбцов называется регулярной.
Операции сложения ф и умножения © матриц выполняются по стандартным правилам с заменой обычных скалярных операций сложения и умножения на тропические. Свойства монотонности скалярных операций распространяются на операции над матрицами, для которых отношение < вводится как покомпонентное.
Будем называть мультипликативно сопряженным транспонированием ненулевой матрицы A G Xmxn операцию преобразования в матрицу A- G Xnxm, элементы которой определяются по правилу: а- = а--1, если aji = 0, и а- = 0 в противном случае.
Рассмотрим множество Xnxn квадратных матриц порядка n. Матрица является диагональной, если все ее недиагональные элементы равны нулю. Единичной называется диагональная матрица I, у которой все элементы на диагонали равны единице.
Для регулярных по строкам матриц A выполняется неравенство AA- > I.
Матрица, состоящая из одного столбца (строки), образует вектор-столбец (вектор-строку). Далее все векторы, если не указано иначе, считаются вектор-
столбцами. Множество вектор-столбцов размерности n обозначим через Xn. Нулевой вектор имеет все компоненты равными 0. Вектор является регулярным, если у него нет нулевых компонент. Вектор, все компоненты которого равны 1, обозначается как 1 = (1,..., 1)T.
Для ненулевого вектора x G Xn справедливо равенство x-x = 1. Если вектор x является регулярным, то выполняется неравенство xx- > I.
Выпуклой линейной комбинацией векторов xi,..., x^ называется выражение вида uixi ф • • • ф ukxk, где числа u1,... ,uk удовлетворяют условию u1 ф • • • ф uk = 1. Это условие означает, что u < 1 для всех индексов i = 1,... ,k, и по крайней мере для одного i выполняется равенство uj = 1.
Пусть заданы матрица A G Xmxn, а также вектор d G Xm, и требуется найти все векторы x G Xn, удовлетворяющие неравенству
Ax < d. (1)
Решение этой задачи обеспечивается следующим утверждением, полное доказательство которого приводится, например, в работах [7, 19].
Лемма 1. Для любой регулярной по столбцам матрицы A и регулярного вектора d 'решение неравенства (1) имеет вид
x < (d-A)-.
3. Задачи тропической оптимизации. Задачи тропической оптимизации обычно состоят в минимизации или максимизации некоторой целевой функции, заданной на векторах над идемпотентным полуполем. Такие задачи возникают, например, при исследовании уравнения Ax = p, для которого требуется найти точное или приближенное решение [3, 10-13].
Сначала предположим, что заданы матрица A G Xnxn и вектор p G Xn. Пусть требуется найти регулярные векторы x G Xn, которые решают задачу
min (Ax)-p ф p-Ax. (2)
Решение этой задачи обеспечивает наилучшее приближенное решение уравнения Ax = p в смысле чебышёвской метрики. Исследование задачи было проведено, например, в работах [14, 16, 17], где приводится частичное решение в следующем виде.
Лемма 2. Для любых регулярных матрицы A и вектора p минимум в задаче (2) равен Д = ((A(p-A)-)-p)1/2 и достигается при x = A(p-A)-.
Обобщением задачи (2) с дополнительным вектором q G Xn является задача нахождения регулярных векторов x, которые обеспечивают
min (Ax)-p ф q-Ax. (3)
Заметим, что при A = I получается задача минимизации функции x-p ф q-x, полное решение которой было получено в [15, 17].
В настоящей работе задача (3) рассматривается в следующей форме. Пусть заданы матрица A G Xnxn и векторы p, q G Xn. Требуется найти все регулярные векторы x G Xn, на которых достигается
min (Ax)-p ф q-x. (4)
В работах [15, 17] было получено частичное решение этой задачи.
Ниже для решения задачи (4) сначала так же, как в работах [15, 17], будет определен минимум целевой функции и получено одно из решений. Затем для полного решения задачи строится система неравенств, которая определяет множество всех решений. На основе использования разреженных матриц находится более широкое множество решений, а затем полное решение задачи в форме семейства подмножеств решений.
4. Предварительный анализ задачи. Цель этого раздела состоит в том, чтобы найти минимум целевой функции, охарактеризовать множество решений и описать некоторые свойства такого множества. Для этого будем использовать подход [15-17, 19, 20], при котором вводится параметр для обозначения минимума целевой функции, а затем задача оптимизации сводится к решению параметризованного неравенства. Докажем следующее утверждение.
Лемма 3. Пусть А —регулярная по строкам матрица, а р и д —регулярные векторы. Тогда минимум в задаче (4) 'равен
В частности, минимум достигается при х = Дд.
Доказательство. Обозначим минимум целевой функции через в. Рассмотрим равенство в = (Ах)-р ф д-х и заметим, что для всех регулярных векторов х справедливо в > 0. Так как в — минимум, то можно заменить равенство на неравенство в > (Ах)-р ф д-х, которое эквивалентно двум неравенствам в > (Ах)-р и в > д-х.
Решая с помощью леммы 1 первое неравенство относительно р, а второе относительно х, приходим к системе из двух неравенств р < вАх и х < вд.
После подстановки второго неравенства системы в первое, получаем р < в2 Ад, затем в > ((Ад)-р)1/2 = Д > 0, и таким образом находим нижнюю границу для в. Проверим, что эта граница достигается при х = Дд. Действительно, после подстановки в целевую функцию имеем (Ах)-р ф д-х = Д-1(Ад)-р ф = Д.
Следовательно, нижняя граница Д для в оказывается точной и поэтому определяет минимум целевой функции.
Заменив в системе в на Д, а затем, умножив обе части первого неравенства на Д-1, получаем, что все регулярные решения задачи (4) должны удовлетворять системе (6). Учитывая эквивалентность преобразований, верно и обратное утверждение: любой регулярный вектор х, удовлетворяющий системе (6), является решением задачи (4). □
Следствие 1. Множество решений задачи (4) вместе с любыми решениями содержит их всевозможные выпуклые линейные комбинации.
Доказательство . Рассмотрим случай выпуклой линейной комбинации двух решений. Пусть х1 и х2 —решения задачи (4), а и и и —такие числа, что и ф и = 1.
Так как векторы х1 и х2 являются решениями задачи (4), то они определяются системой (6), а их выпуклая комбинации удовлетворяет соотношениям
Д=((Ад)-р)1/2,
(5)
а все регулярные решения х определяются системой неравенств
Ах > Д-1р, х < Дд.
(6)
А(и1х1 ф и2х2) > и1Д 1р ф и2Д 1р = (и1 ф и2)Д 1р = Д 1р, и1х1 ф и2х2 < и1 Дд ф и2Дд = (и1 ф и2)Дд = Дд.
Следовательно, вектор uixi©М2Ж2 также является решением задачи (4). Полученный результат легко обобщается на случай с произвольным количеством решений. □
5. Разрежение матрицы задачи. Для расширения множества решений задачи (4) применим метод с использованием разрежения матриц, предложенный в [18]. Сначала преобразуем матрицу задачи A = (a«) в матрицу A = (a«), приравнивая нулю все элементы, которые строго меньше порогового значения по следующему правилу:
Ja«, если а™ > Д-2р4а-1;
a« = s «
I 0, в противном случае.
Далее матрицу A, полученную при помощи такого преобразования, будем называть разреженной матрицей задачи. Свойства этой матрицы отражает следующий результат.
Лемма 4. Замена матрицы A на A не меняет множество 'решений задачи (4).
Доказательство. Сначала проверим, что переход к разреженной матрице сохраняет минимальное значение Д целевой функции, полученное в лемме 3. Для упрощения выкладок рассмотрим величину Д2. Определим индексы k и s следующим образом:
k = arg max (ajiqi © • • • © ainq„)-1pj, s = arg max afc«q«,
1<j<n 1<j<n
и запишем выражение для Д2 в виде
n
Д2 = (Aq)-p = 0(ai1 q1 © • • • © ajnqn)-1pj = i=1
= (afc1 q1 © • • • © afcnqn)-1Pfc = (afcsqs)-1pfc.
Регулярность по строкам матрицы A и регулярность вектора q гарантируют, справедливость неравенства a^q! © • • • © ajnqn > 0 для всех i = 1,..., n. Кроме того, так как вектор p регулярный, то Д > 0, а значит и afcs > 0.
Рассмотрим произвольную строку i матрицы A. Из формулы для Д2 видно, что выполняется неравенство Д2 > (a^ q1 ©• • •©ajnqn)-1 pj. Оно эквивалентно неравенству ai1q1 © • • • © ajnqn > Д-2р^ которое верно только при условии, что выполняется a«q« > Д-2pj для некоторого индекса j. Отсюда заключаем, что в каждой строке i матрицы A существует по меньшей мере один элемент a«, который удовлетворяет условию
ajj > Д-2piq—1. (7)
Рассмотрим строку k матрицы A, чтобы убедиться в справедливости ay < Д—2pkq—1 для всех j в этой строке. Если ay = 0, то это неравенство, очевидно, выполняется. При условии ak« > 0 имеем (ak«q«)-1pk > (ak1 q1 © • • • © aknqn)-1pk = Д2, из чего следует, что неравенство также выполняется. Так как Д2 = (afcsqs)-1pfc, то в строке k есть элементы, которые превращают неравенство (7) в равенство, но нет ни одного элемента, для которого (7) становится строгим неравенством.
Предположим теперь, что (7) не выполняется для некоторых i и j. Учитывая, что pj > 0, можно записать a« < Д-2р^-1 < (aj1 q1 ©• • •©ajnqn)q—1, откуда получаем
а, < ац ф • • • ф а^п^п. В этом случае уменьшение а, путем понижения значения а, вплоть до 0 не влияет на величину ал ф • • • ф и, следовательно, на значение
Д2 > (ал91 ф • • • ф а*пд„)-1р*.
Проверим, что все элементы а,, которые не удовлетворяют неравенству (7), можно обнулить, оставляя без изменений множество регулярных решений задачи (4).
Рассмотрим систему (6). Записав первое неравенство системы в скалярной форме, получим, что для каждого i выполняется а^ж ф • • • ф а4„ж„ > Д-1р^
Учитывая регулярность векторов х и д, умножая обе части второго неравенства системы на регулярный вектор д- слева, можно привести это неравенство к эквивалентному неравенству д-х < Д. В скалярной форме получаем ф- • •фд- 1жп < Д.
Предположим, что в матрице А имеется элемент, скажем а,, удовлетворяющий условию а, < Д-2р^-1 и тем самым нарушающий неравенство (7). Тогда будем иметь а,ж, < Д-2р^д^ж, < Д-2(д-1Ж1 ф — ф< Д-1р < алЖ1 ф — фа.гпжп. В этом случае слагаемое а,ж, не вносит никакого вклада в значение суммы алЖ1 ф • • • ф Следовательно, можно положить а, = 0, не изменяя множество решений
системы.
Осталось понять, что обнуление элементов а,, которые не удовлетворяют неравенству (7), равносильно замене матрицы А на А. □
Заметим, что ненулевые элементы матрицы А отвечают условию а, > Д-2^^-1. Тогда для элементов матрицы А- = (а,) справедливо соотношение а, < Д21. В матричной форме можем записать неравенство А- < Д2др-, которое будет использовано ниже.
6. Полное решение задачи. В этом разделе будет представлено полное решение задачи (4) в виде семейства решений, которое задается множеством матриц, полученных из разреженной матрицы задачи, продолжая процесс обнуления ее элементов.
Начнем с того, что расширим решение х = Дд, полученное в лемме 3, до некоторого множества с помощью следующего утверждения.
Лемма 5. Пусть А —регулярная по строкам разреженная матрица задачи (4), где р и д —регулярные векторы, и Д = ((Ад) р)1/2. Тогда минимум в задаче (4) равен Д и достигается на любом векторе х, который удовлетворяет условию
Д-1А-р < х < Дд. (8)
Доказательство. Заметим, что определяемое неравенством (8) множество не пусто. Действительно, из того, что матрица А является разреженной, следует неравенство А- < Д2др-, после умножения которого на Д-1р справа получаем Д-1А-р < Дд.
Покажем, что все векторы х, которые удовлетворяют двойному неравенству (8), являются решениями системы (6). В силу того, что правое неравенство в (8) совпадает со вторым в системе (6), достаточно проверить выполнение Ах > Д-1р.
Так как матрица А регулярна по строкам, то I < АА-. Умножив левое неравенство из (8) на А слева, получаем Д-1р < Д-1АА-р < Ах, откуда следует требуемое неравенство системы (6). □
Теперь сформулируем лемму, которая обеспечивает полное решение задачи (4).
Лемма 6. Пусть А — регулярная по строкам разреженная матрица задачи (4), где р и д — регулярные векторы, и Д = ((Ад) р)1/2. Обозначим через А множество матриц, полученных из А посредством сохранения по одному ненулевому элементу в каждой строке и обнулением остальных.
Тогда минимум в задаче (4) 'равен Д, а все регулярные решения х образуют семейство решений, каждое из которых определяется условием
Д-1А-р < х < Дд, А1 е А. (9)
Доказательство. Как было показано выше, все решения задачи (4) задаются системой (6). Следовательно, для доказательства леммы требуется проверить, что каждому решению системы (6) соответствует некоторая матрица А1 е А и наоборот.
Пусть вектор х является решением системы (6). Рассмотрим векторное неравенство Ах > Д-1р и исследуем каждое из соответствующих скалярных неравенств, чтобы определить максимальное слагаемое в левой части. Несложно видеть, что найдется матрица А1 е А с ненулевыми элементами в каждой строке, которые соответствуют этим максимальным слагаемым.
Теперь предположим, что х — решение неравенства (9) для некоторой матрицы А1. Так как матрица А1 регулярна по строкам и А1 < А, то выполняется I < А1А- < АА-. Умножая левую часть неравенства (9) на А слева, получаем Ах > Д-1АА-р > Д-1 р. Из этого следует, что х является решением системы (6), а значит и задачи (4). □
Отметим, что различные множества решений из семейства, описанного в лемме 6, могут иметь пересекающиеся подмножества.
Для получения решения в параметрическом виде докажем следующую теорему.
Теорема. Пусть А — регулярная по строкам разреженная матрица задачи (4), где р и д — регулярные векторы, и Д = ((Ад) р)1/2. Обозначим через А множество матриц, полученных из А посредством сохранения по одному ненулевому элементу в каждой строке и обнулением остальных, а через В — матрицу, столбцами которой являются векторы Ь1 = Д-1А-р для всех матриц А1 е А.
Тогда минимум в задаче (4) равен Д, а все регулярные решения х имеют вид
х = Ви ф V, V < Дд, 1ти = 1. (10)
Доказательство. Пусть множество А содержит к матриц. Для каждого г = 1,..., к возьмем матрицу А^ е А и определим левую часть 6 = Д-1 А—р неравенства (9). Учитывая, что правая часть неравенства всегда равна Дд, семейство решений описывается следующим набором двойных неравенств:
Ь1 < х1 < Дд,..., < хй < Дд.
Рассмотрим произвольные векторы х1,..., х^, для которых соответствующие двойные неравенства выполняются. По следствию 1 множество решений задачи (4) замкнуто относительно операции взятия выпуклой линейной комбинации. Следовательно, если м ф • • • ф = 1, то вектор -щх! ф • • • ф х^ также является решением.
Покажем, что общее решение задачи можно записать в виде
«161 ф — ф < х < Дд, М1 ф---ф = 1. (11)
Если вектор х —решение задачи (4), то он принадлежит семейству решений, определяемому нижней границей Ь = Д-1А-р, которая задается некоторой матрицей € А. Чтобы представить вектор х в форме (11), положим коэффициент м равным 1, а все остальные коэффициенты обнулим.
Пусть вектор х удовлетворяет системе (11). Тогда существует такой индекс ^ для которого выполняются соотношения м = 1 и Д-1А-р = < х < Дд, откуда следует, что х принадлежит семейству решений, образованному матрицей
Построим матрицу В = (&1,..., ), используя векторы Ь в качестве столбцов, и введем вектор и = (м1,..., )т. Теперь двойное неравенство и равенство системы (11) можно записать в векторной форме как Ви < х < Дд и 1ти = 1. С помощью вектора V = («1,..., «п)т запишем двойное неравенство в эквивалентной форме в виде системы, состоящей из равенства х = Ви ф V и неравенства V < Дд. Объединение этой системы с условием \ти = 1 приводит к решению (10). □
В качестве следствия полученного результата представим решение задачи (2), которая является частным случаем задачи (4).
Следствие 2. Пусть А — регулярная по строкам разреженная матрица задачи (2), где р —регулярный вектор, и Д = ((А(р-А)-)-р)1/2. Обозначим через А множество матриц, полученных из А посредством сохранения по одному ненулевому элементу в каждой строке и обнулением остальных, а через В — матрицу, столбцами которой являются векторы Ь1 = Д-1А-р для всех матриц А1 € А.
Тогда минимум в задаче (4) равен Д, а все регулярные решения х имеют вид
х = Ви ф V, V < Д(р-А)-, 1ти = 1.
7. Процедуры построения полного решения. Перебор всевозможных матриц А1 € А для построения подмножеств семейства решений задачи в соответствии с результатом сформулированной теоремы может представлять определенные трудности. Кроме того, некоторые семейства решений могут содержаться в других семействах и поэтому при записи общего решения могут быть отброшены. Ниже описываются процедуры, позволяющие во многих случаях сократить число подмножеств, которые необходимо учесть при построении общего решения.
Рассмотрим следующую процедуру построения нижних границ в неравенстве (9) для семейства решений, которая заключается в последовательном выборе по одному ненулевому элементу в строках разреженной матрицы А.
Предположим, что в матрице Аа зафиксированы элементы в некоторых строках, в результате чего получена матрица А = (а,). Пусть выбран ненулевой элемент в одной из оставшихся строк, например элемент аг, в строке г и столбце значение которого фиксируется, а остальные элементы в этой строке замещаются нулями.
Всякий вектор х = (ж$), который является решением задачи (4), удовлетворяет системе (6), в частности, ее первому неравенству в форме Ах > Д-1 р. Скалярное неравенство для строки г, где все элементы кроме аг, равны 0, записывается в виде аг,ж, > Д-1рг, что эквивалентно неравенству ж, > Д-1а-,1рг.
В столбце ] возьмем элемент а,, который расположен в одной из еще не рассмотренных строк е. При выполнении условия а,р-1 > аг,р-1, которое эквивалентно условию а, > аг,р-1ря, имеем а,ж, > аг,р-1ряД-1а-,1рг > Д-1ря. Тогда неравенство аЯ1 Ж1 ф • • • ф ая„ж„ > Д-1ря в системе (6) выполняется вне зависимости от значений ж; для всех I = В этом случае дальнейшее исследование ненулевых элементов
ая; в строке в не может дать новых решений. Эти элементы можно заменить нулями, не нарушая неравенство, что уменьшает число альтернатив, которые требуется рассмотреть.
Заметим, что для проверки условия р-1 > arjр-1 удобно предварительно умножить каждую строку г матрицы А на элемент р—1, а затем исследовать элементы полученной матрицы, которую обозначим через А' = (а^).
Если в какой-то строке матрицы А' преобладают элементы, являющиеся минимальными ненулевыми элементами в своих столбцах, то, вероятно, стоит начать перебор с подобной строки. Это позволит поочередно фиксировать такие элементы и не рассматривать строки, содержащие остальные ненулевые элементы соответствующих столбцов, а при отсутствии в таком столбце нулевых элементов сразу получать одну из границ.
Отсюда видно, что выбор строки в матрице может сыграть существенную роль для уменьшения количества рассматриваемых подмножеств семейства решений.
Некоторые нижние границы подмножеств, полученные с помощью процедуры построения полного семейства решений, могут быть несущественными в том смысле, что их удаление не повлияет на полученное множество решений. Нетрудно видеть, например, что если для каких-то двух границ Ь и Ь выполняется неравенство Ь < Ь, то граница Ь может быть удалена из списка без потери решений.
Сформулируем критерий для отбрасывания несущественных границ.
Предложение. Пусть В — матрица, столбцы которой определяют нижние границы для некоторого набора подмножеств семейства решений (9), а Ь1 = Д-1А1 р, где А1 е А, — еще одна граница. Тогда граница Ь1 является несущественной, если выполняется неравенство
1Т(Ь- В)- > 1. (12)
Доказательство. Подмножество с новой нижней границей Ь1 уже принадлежит семейству решений и может не учитываться, если найдется выпуклая комбинация нижних границ подмножеств семейства, которая не меньше, чем 61. Учитывая, что имеющиеся нижние границы подмножеств образуют столбцы матрицы В, это условие эквивалентно выполнению неравенства Ви < Ь1 для некоторого вектора и такого, что 1ти = 1.
Решение неравенства относительно и с помощью леммы 1 дает и < (Ь-В)-. Для выполнения условия 1ти = 1 необходимо и достаточно, чтобы хотя бы одна координата вектора (Ь-В)- была не меньше 1. Это условие записывается в виде
неравенства 1т(Ь-В)- > 1, которое является условием отбрасывания границы Ь1.
□
8. Численный пример. Рассмотрим задачу (4), которая формулируется для полуполя Ктах,+ при п = 3. Предположим, что имеются матрица А и векторы р и д следующего вида:
Найдем все регулярные векторы х, которые решают задачу. Сначала по формуле (5) определим минимум Д, для чего последовательно вычислим
М = I 0 ^ , Д2 = (Ад)-р = 4, Д = 2. Для нахождения разреженной матрицы задачи построим пороговую матрицу
Д-2рд-
-5 -6 -7 -1 -2 -3 -2 -3 -4
Заменяя нулем 0 = -то те элементы матрицы А, которые строго меньше элементов матрицы Д-2рд-, формируем матрицу А, а умножая каждую строку А на соответствующий элемент вектора р-, находим матрицу А'. В результате получим
5 5 2 1 0 3
5
2
21
0
-5 0 -7 520
Применим процедуру построения полного решения. Заметим, что в матрице А' во второй и в третьей строках все ненулевые элементы являются наименьшими ненулевыми элементами в своих столбцах. В силу этого следует начать с выбора соответствующих элементов в матрице А, что позволит сократить перебор разреженных матриц. Для определенности начнем со второй строки.
Сначала зафиксируем элемент Й21. Учитывая, что в матрице А' элемент ¡0)1 не превосходит элементов а'ц и а^1, можно сохранить без изменений соответствующие элементы ац и аз1 матрицы А, а остальные элементы этой матрицы заменить нулями.
В результате получаем матрицу А1, нижнюю границу Ь1 = Д-1А-р, а также матрицу границ В1 = (Ь1 ) в виде
А1
500 -1 0 0 200
3 0 0
В1
3 0 0
Теперь зафиксируем элемент а2з во второй строке матрицы А. В силу того, что для элементов матрицы А' выполняется условие а2з < а1з, в матрице А можно оставить элемент Й13, а все другие элементы первой строки заменить нулями.
В оставшейся третьей строке следует по очереди фиксировать каждый ее ненулевой элемент. Получим следующие матрицы и соответствующие им нижние границы:
А
2=
0 0 2 0 0 3
20
0
0 0 2 А3 = I 0 0 -3
Ь2 = Д-1 А-р
Ьз = Д-1 А3р:
01
0
Выясним, являются ли полученные границы существенными. Проверка условия (12) приводит к следующим результатам:
1Т(Ь-В1)- =0 = 1, 1т(Ь-В1)- = 0 < 1.
Ь
1
Для границы условие оказывается выполненным, а потому ее можно не учитывать. Заметим, что непосредственная проверка выполнения неравенства bi < 62 подтверждает этот результат. Однако для границы Ьз условие не выполняется и ее следует добавить к имеющейся границе bi. Окончательно получим матрицу границ B = (bi, Ьз).
Общее решение (10) задачи (4) записывается в форме
/ 3 0 \ / 3 \
x = I 0 0 I u ф v, v < I 4 I , 1Tu = 1.
v° ч \4
Представление решения в виде (9) дает два подмножества решений:
xi =3, Х2 < 4, хз < 5; х1 < 3, 0 < х2 < 4, х3 = 5.
В заключение авторы благодарят рецензентов за ряд важных замечаний и предложений, которые были учтены при подготовке окончательного текста статьи.
Литература
1. Baccelli F., Cohen G., Olsder G. J., Quadrat J.-P. Synchronization and Linearity. In: Wiley Series in Probability and Statistics. Chichester: Wiley, 1993. 514 p.
2. Маслов В. П., Колокольцов В. Н. Идемпотентный анализ и его применение в оптимальном управлении. М.: Физматлит, 1994. 144 с.
3. Cuninghame-Green R. A. Minimax algebra and applications. In Ser.: Advances in Imaging and Electron Physics / ed. by P. W. Hawkes. San Diego: Academic Press, 1994. Vol. 90. P. 1-121. https: //doi.org/ 10.1016/S 1076-5670 (08)70083-1
4. Golan J. S. Semirings and Affine Equations over Them. In Ser.: Mathematics and Its Applications. Springer, 2003. Vol. 556. 241 p. https://doi.org/10.1007/978-94-017-0383-3
5. Heidergott B, Olsder G. J., van der Woude J. Max Plus at Work. In: Princeton Series in Applied Mathematics. Princeton: Princeton Univ. Press, 2006. 226 p.
6. Gondran M., Minoux M. Graphs, Dioids and Semirings. In Ser.: Operations Research / Computer Science Interfaces. New York: Springer, 2008. Vol. 41. 383 p. https://doi.org/10.1007/978-0-387-75450-5
7. Кривулин Н. К. Методы идемпотентной алгебры в задачах моделирования и анализа сложных систем. СПб.: Изд-во С.-Петерб. ун-та, 2009. 256 с.
8. Glazek K. A Guide to the Literature on Semirings and Their Applications in Mathematics and Information Sciences. Springer, 2002. 392 p. https://doi.org/10.1007/978-94-015-9964-1
9. Воробьев Н. Н. Экстремальная алгебра положительных матриц // Elektron. Informationsverarb. Kybernet. 1967. Bd. 3, N 1. S. 39-72.
10. Cuninghame-Green R. A. Projections in minimax algebra // Math. Program. 1976. Vol. 10. P. 111-123. https://doi.org/10.1007/BF01580656
11. Zimmermann K. Some optimization problems with extremal operations // Mathematical Programming at Oberwolfach II. Vol. 22 of Mathematical Programming Studies / ed. by B. Korte, K. Ritter. Berlin: Springer, 1984. P. 237-251. https://doi.org/10.1007/BFb0121020
12. Cechlarova K., Cuninghame-Green R. A. Soluble approximation of linear systems in max-plus algebra // Kybernetika. 2003. Vol. 39, N 2. P. 137-141.
13. Butkovic P., Tam K. P. On some properties of the image set of a max-linear mapping // Tropical and Idempotent Mathematics. Vol. 495 of Contemporary Mathematics / ed. by G. L. Litvinov, S. N. Sergeev. Providence: AMS, 2009. P. 115-126. https://doi.org/10.1090/conm/495/09694
14. Кривулин Н. К. О решении линейных векторных уравнений в идемпотентной алгебре // Математические модели. Теория и приложения. Вып. 5: Сб. науч. статей / под ред. М. К. Чиркова. СПб.: ВВМ, 2004. С. 105-113.
15. Krivulin N. A new algebraic solution to multidimensional minimax location problems with Chebyshev distance // WSEAS Trans. Math. 2012. Vol. 11, N 7. P. 605-614.
16. Krivulin N. Solution of linear equations and inequalities in idempotent vector spaces // Int. J. Appl. Math. Inform. 2013. Vol. 7, N 1. P. 14-23.
17. Krivulin N., Zimmermann K. Direct solutions to tropical optimization problems with nonlinear objective functions and boundary constraints // Mathematical Methods and Optimization Techniques in Engineering / ed. by D. Biolek, and H. Walter, I. Utu, C. von Lucken. WSEAS Press, 2013. P. 86-91.
18. Krivulin N. Algebraic solution of tropical optimization problems via matrix sparsification with application to scheduling // J. Log. Algebr. Methods Program. 2017. Vol. 89. P. 150-170. https://doi.org/10.1016/jjlamp.2017.03.004
19. Krivulin N. Extremal properties of tropical eigenvalues and solutions to tropical optimization problems // Linear Algebra Appl. 2015. Vol. 468. P. 211-232. https://doi.org/10.1016/j.laa.2014.06.044
20. Krivulin N. K., Sorokin V. N. Solution of a tropical optimization problem with linear constraints // Vestnik St. Petersburg Univ. Math. 2015. Vol. 48, N 4. P. 224-232. https://doi.org/10.3103/S1063454115040081
Статья поступила в редакцию 7 июля 2017 г.; рекомендована в печать 21 сентября 2017 г. Контактная информация:
Кривулин Николай Кимович —д-р физ.-мат. наук, доц.; [email protected] Сорокин Владимир Николаевич — [email protected]
On solution of a multidimensional tropical optimization problem using matrix sparsification
N. K. Krivulin, V. N. Sorokin
St. Petersburg State University, Universitetskaya nab., 7-9, St. Petersburg, 199034, Russian Federation
For citation: Krivulin N. K., Sorokin V. N. On solution of a multidimensional tropical optimization problem using matrix sparsification. Vestnik SPbSU. Mathematics. Mechanics. Astronomy, 2018, vol. 5 (63), issue 1, pp. 91-104. https://doi.org/10.21638/11701/spbu01.2018.110
In this work a complete solution for a problem of minimizing a function defined on vectors with elements from a tropical (idempotent) semifield is proposed. The tropical optimization problem under consideration arises, for instance, when one needs to find the best, in the sense of Chebyshev metric, approximate solution for tropical vector equations, and occurs in various applications, including scheduling, location and decision-making problems. To solve the problem, first, the minimum value of the objective function is obtained, a characterization of the solution set in the form of a system of inequalities is proposed, and one of the solutions is presented. Next, by using a sparsification of the matrix in the problem, an extended set of solutions, and then a complete solution in the form of a family of subsets are derived. Procedures are described that allow reducing the number of subsets, which one needs to examine when constructing the complete solution. It is shown how the complete solution can be represented in parametric way in a compact vector form. The solution obtained extends known results, which are usually confined to the derivation of one solution and do not allow to find the entire solution set. To illustrate the main results of the work numerical examples are given for the solution of a problem in the set of three-dimensional vectors.
Keywords: idempotent semifield, tropical optimization, Chebyshyov approximation, complete solution, matrix sparsification.
References
1. Baccelli F., Cohen G., Olsder G. J., Quadrat J.-P., ¡Synchronization and Linearity. In: Wiley Series in Probability and Statistics (Wiley, Chichester, 1993, 514 p.).
2. Maslov V. P., Kolokoltsov V. N., Idempotent Analysis and Its Applications to Optimal Control Theory (Nauka, Moscow, 1994, 144 p.) [in Russian].
3. Cuninghame-Green R. A., Minimax algebra and applications. In Ser. Advances in Imaging and Electron Physics 90, 1-121 (ed. by P. W. Hawkes, Academic Press, San Diego, 1994). https://doi.org/10.1016/S1076-5670(08)70083-1
4. Golan J. S., Semirings and Affine Equations over Them. In Ser. Mathematics and Its Applications 556 (Springer, 2003, 241 p.). https://doi.org/10.1007/978-94-017-0383-3
5. Heidergott B, Olsder G. J., van der Woude J., Max Plus at Work. In: Princeton Series in Applied Mathematics (Princeton Univ. Press, Princeton, 2006, 226 p.).
6. Gondran M., Minoux M., Graphs, Dioids and Semirings. In Ser. Operations Research / Computer Science Interfaces 41 (Springer, New York, 2008, 383 p.). https://doi.org/10.1007/978-0-387-75450-5
7. Krivulin N. K., Methods of Idempotent Algebra for Problems in Modeling and Analysis of Complex Systems (Saint Petersburg Univ. Press, St. Petersburg, 2009, 256 p.) [in Russian].
8. Glazek K., A Guide to the Literature on Semirings and Their Applications in Mathematics and Information Sciences (Springer, 2002, 392 p.). https://doi.org/10.1007/978-94-015-9964-1
9. Vorob'ev N. N., "The extremal algebra of positive matrices", Elektron. Informationsverarb. Kybernet 3(1), 39-72 (1967) [in Russian].
10. Cuninghame-Green R. A., "Projections in minimax algebra", Math. Program 10, 111-123 (1976). https://doi.org/10.1007/BF01580656
11. Zimmermann K., "Some optimization problems with extremal operations", Mathematical Programming at Oberwolfach II. Vol. 22 of Mathematical Programming ¡Studies, 237-251 (eds B. Korte, K. Ritter; Berlin, Springer, 1984). https://doi.org/10.1007/BFb0121020
12. Cechlarova K., Cuninghame-Green R. A., "Soluble approximation of linear systems in max-plus algebra", Kybernetika 39(2), 137-141 (2003).
13. Butkovic P., Tam K. P., "On some properties of the image set of a max-linear mapping', Tropical and Idempotent Mathematics. Vol. 495 of Contemporary Mathematics, 115-126 (ed. by G. L. Litvinov, S. N. Sergeev. Providence, AMS, 2009). https://doi.org/10.1090/conm/495/09694
14. Krivulin N. K., "On solution of linear vector equations in idempotent algebra", Mathematical Models. Theory and Applications, Issue 5, 105-113 (ed. by M. K. Chirkov. St. Petersburg, VVM, 2004) [in Russian].
15. Krivulin N., "A new algebraic solution to multidimensional minimax location problems with Chebyshev distance", WSEAS Trans. Math. 11(7), 605-614 (2012).
16. Krivulin N., "Solution of linear equations and inequalities in idempotent vector spaces", Int. J. Appl. Math. Inform. 7(1), 14-23 (2013).
17. Krivulin N., Zimmermann K., "Direct solutions to tropical optimization problems with nonlinear objective functions and boundary constraints", Mathematical Methods and Optimization Techniques in Engineering, 86-91 (ed. by D. Biolek, and H. Walter, I. Utu, C. von Lucken; WSEAS Press, 2013).
18. Krivulin N., "Algebraic solution of tropical optimization problems via matrix sparsification with application to scheduling', J. Log. Algebr. Methods Program 89, 150-170 (2017). https://doi.org/10.1016 /j.jlamp.2017.03.004
19. Krivulin N., "Extremal properties of tropical eigenvalues and solutions to tropical optimization problems", Linear Algebra Appl. 468, 211-232 (2015). https://doi.org/10.1016/j.laa.2014.06.044
20. Krivulin N. K., Sorokin V. N., "Solution of a tropical optimization problem with linear constraints", Vestnik St. Petersburg Univ. Math. 48(4), 224-232 (2015). https://doi.org/10.3103 /S1063454115040081
A u t h o r's i n f o r m a t i o n:
Krivulin Nikolai K. — [email protected]
Sorokin Vladimir N. — [email protected]
