Системы неравенств и задачи оптимизации с двусторонними ограничениями на переменные Зоркальцев Валерий Иванович, проф., д.т.н., Заведующий лабораторией «Методов математического моделирования и оптимизации в энергетике» Института систем энергетики им. Л.А. Мелентьева СО РАН, г. Иркутск

2 Исходные положения 1.При сужении класса рассматриваемых задач может расширятся набор их полезных свойств. 2.Разработчик любой математической модели способен оценить априори диапазон возможных значений искомых показателей. 3.Введение ограничений на диапазоны значений искомых показателей может дать новые, полезные свойства рассматриваемой задаче.

3 Общая постановка рассматриваемых задач Найти такой, что: Здесь: X исходное (как-то определяемое) множество реше- ний задачи; g, d заданные (в том числе экспертно вводи- мые) векторы R n, характеризующие диапазоны возможных значений переменных. Вводимое условие (2) означает в частности, что множество допустимых (даже по какой-то части исходных условий) переменных является ограниченным (возможно пустым).

4 I. Введение двусторонних ограничений на переменные может быть полезно для выясне- ния вопроса существования решения у задачи 1.На базе теоремы Вейерштрасса, если исходная задача представляется в виде проблемы оптимизации непрерывной функции на замкнутом множестве. 2.На базе теорем о неподвижной точке (Бауэра, Какутани), если процедура поиска решения задачи представляется в виде непрерывного отображения векторов некоторого замкнутого множества, удовлетворяющих (2), в векторы этого же множества.

5 II. Введение двусторонних ограничений на переменные может быть полезно для идентифи- кации ситуации отсутствия решения у задачи Доказать отсутствия решения у задачи (из-за несовместности ограничений и других причин) часто не просто. Только в редких случаях это сводится к выводу из условий задачи очевидных противоречий. Для систем линейных уравнений и неравенств (и на базе этого для других классов задач) конструктивное доказательство отсутствия решения можно осуществлять на базе теорем об альтернативных системах линейных неравенств.

6 Теоремы об альтернативных системах линейных неравенств Системы линейных неравенств являются непосредственным и очень важным (для моделирования и вычислительной математики) обобщением систем линейных алгебраических уравнений (СЛАУ). Любую СЛАУ можно представить в виде системы линейных неравенств. Обратное не верно. То, что теория систем линейных неравенств не изучается в базовых курсах университетов большое, затянувшееся недоразумение. Общий вид теорем об альтернативах: каждой системе S можно поставить в соответствие систему S *, такую что заранее известно одна и только одна из этих систем совместна.

7 Теорема Фредгольма об альтернативных СЛАУ Либо разрешима система Ax = b, (1) либо разрешима система A T u = 0, b T u=1. (2) Заданы: A матрица m×n, вектор b R m. Алгоритм поиска решения систем (1) или (2): Пусть решение указанной задачи. Тогда: если то решение системы (1), если то решение системы (2).

8 Модификация теоремы Фредгольма при ограничениях на абсолютные значения переменных (теорема Дакса) Пусть задан вектор d R n, d>0. Тогда либо при некотором x R n Ax = b, |x| d,(1) либо при некотором u R m b T u d T |A T u| > 0,(2) где |x|, |A T u| векторы, составленные из абсолютных значений компонент векторов x и A T u. Здесь условие (2) является более слабым, чем в теореме Фредгольма. Если A T u = 0, b T u > 0, то (2) выполняется. Обратное не верно. Для выполнения условия (2) нет необходимости в достижении равенства A T u = 0.

9 Некоторые приложения теорем об альтер- нативных системах линейных неравенств 1.Для выявления ситуации отсутствия решения у задачи: если нашли какое-либо решение для системы S *, то тем самым доказали противоречивость условий исходной системы S. 2.Для выявления избыточных ограничений (условий- следствий из других) избавление от которых не влияет на задачу. 3.Для выявления решений с минимальным набором активных ограничений. 4.Основа теории двойственности линейного (и на базе этого нелинейного) программирования

10 Система линейных неравенств с двусто- ронними ограничениями на переменные (общий случай) Ax = b, g x d. (S) Теорема. Либо разрешима система S, либо при некотором u R m b T u d T (A T u) + g T (A T u) > 0, (S*) где ( ) + (A T u) положительные и отрицательные срезки: для y R m (y + ) j = max{0, y j }, (y ) j = min{0, y j }, j = 1, …, m. Неравенство S* эквивалентно следующей системы линейных неравенств относительно u R m, v R n, w R n b T u d T v g T w > 0, A T u v w = 0, v 0, w 0.

11 Представление альтернативной системы в виде задачи линейного программирования b T u d T v g T w max (1) A T u v w = 0, v 0, w 0. (2) Задача всегда имеет допустимые условиям (2) решение. При любом u R m векторы v = (A T u) + + e, w = (A T u) e (здесь е – вектор из единиц) удовлетворяют (2) в строгой форме (могут служить стартовой точкой для алгоритмов внутренних точек) Если задача (1), (2) не имеет решения (целевая функция не ограничена сверху при условиях (2)), то имеет решение неравенство S *. Если задача (1), (2) имеет оптимальное решения, то вектор множителей Лагранжа ограничений-равенств составляет реше- ние исходной системы S.

12 III. Модель расчета допустимых режимов электроэнергетических систем Исходная задача: найти вектор y (1) где квадратичная вектор-функция. Линеаризованная, упрощенная (в результате применения к (1) метода приведенного градиента) задача: найти вектор x (2)

13 Сопоставление вариантов алгоритмов внутрен- них точек на задачах определения допустимых режимов электроэнергетических систем АлгоритмЧисло итераций на задачах несовместныхсовместных 6*740*802*719*19201*201 A B C D11558 E A, B – прямые алгоритмы C, D – двойственные алгоритмы Е – самосопряженный алгоритм

14 Сопоставление вариантов алгоритмов внутрен- них точек на задачах определения допустимых режимов электроэнергетических систем Задачи/АлгоритмABCD Совместные (30x80)–(41х80) 9.5(3–13)13.4(7–17)10.2(6–13)5.9(3–8) Несовместные (30x80)–(41х80) 1.6(1–5)1.6(1–7)1(все 1) Тестовый пример (19х19) Тестовый пример (201х201)

15 Некоторые выводы из экспериментов 1. Высокая эффективность алгоритмов внутренних точек, уси- ливающаяся двусторонними ограничениями на переменные. 2. Высокая эффективность введенного критерия несовместности. 3. Большая польза от сочетания экспериментальных и теоретиче- ских исследований: новый вариант алгоритмов внутренних точек (D) гораздо эффективнее, чем affine scaling method (A) и dual affine scaling method (C), в т.ч. из-за того, что множители Лагранжа сходятся быстрее, чем исходные переменные

16 Необходимо вводить в учебные курсы алгоритмы внутренних точек 1. Теперь общепризнано, что они высокоэффективны даже при решении задач линейного программирования. 2. Проще и изящнее: алгоритмы симплекс-метода на их фоне подобны квадратному колесу на фоне круглого. 3. Имеются наилучшие полиномиальные оценки объемов вычислений от размерности задачи. 4. Легко переносятся на многие классы задач нелинейной оптимизации, естественно сочетаются с другими методами. 5. Были созданы и активно используются с 70-х годов в СССР (ИМ СО РАН, СЭИ СО РАН, ВЦ РАН).

17 IV. Задача линейного программирования с двусторонними ограничениями Двойственная задача Имеется только один случай отсутствия решения – противоречи- вость ограничений исходной задачи, неограниченность целевой функции двойственной задачи. Легко сформировать стартовую точку для двойственного алгоритма внутренних точек: при любом

18 V. Задача выпуклого программирования Исходная задача (P) Двойственная задача (D) Здесь сопряженные по Лежандру-Фенхелю строго выпуклые функции вещественного аргумента.

19 Задача выпуклого программирования Задачи (P) и (D) представляют случай симметричной двойствен- ности в оптимизации – когда двойственная к двойственной задаче оптимизации совпадает с исходной задачей. Приложения. 1. Регуляризация в линейной оптимизации (ИММ УрО РАН, Еремин И.И.). 2. Решение большеразмерных систем линейных неравенств (ВЦ РАН, Голиков А.И., Евтушенко Ю.Г.). 3. Поиск равновесия Нэша. 4. Поиск термодинамического равновесия. 5. Модели потокораспределения (нелинейные транспортные задачи, электрические и гидравлические цепи).

20 Задача определения "узких" мест в Единой системе газоснабжения с позиций обеспечения её живучести Были выполнены расчеты для двух реальных сетей: –Укрупненная сеть ЕСГ (21 узел, 28 ветвей) –Подробная сеть ЕСГ (337 узла, 589 ветвей) Цель расчетов: 1)определить возможности системы по удовлетворению потребности узлов в ресурсе, когда возможен «нагруженный» режим передачи по ветвям сети 2)определить «узкие» места ЕСГ – ветви, где передача осуществляется в «нагруженном» режиме 3)проранжировать «узкие» места

21 Расчет подробной сети ЕСГ Количество узлов: 337 Количество ветвей: 589 Количество итераций метода внутренних точек: 51 Время счета: 2 сек (Процессор - Intel Core i-5 650, 3200МГц)

22 СПАСИБО ЗА ВНИМАНИЕ!