РАЗРАБОТКА И ЭКСПЕРИМЕНТАЛЬНОЕ ИССЛЕДОВАНИЕ МЕТОДА ВЕТВЕЙ И ГРАНИЦ ДЛЯ РЕШЕНИЯ ЗАДАЧИ ВОГНУТОГО ПРОГРАММИРОВАНИЯ НА ОСНОВЕ ДОБАВЛЕНИЯ ОГРАНИЧЕНИЙ К ДОПУСТИМОМУ МНОЖЕСТВУ
Введение 3
Постановка задачи и алгоритм ее решения 6
Общая схема метода ветвей и границ для задачи вогнутого программирования 9
Постановка задачи и ее программная реализация 11
Алгоритм решения задачи в узле 12
Методы и поля, используемые в программе 20
Эксперименты и выводы 32
Заключение 40
Список литературы
Как известно, зачастую задачи экономики, физики, экологии, молекулярной биологии и медицины сводятся к задачам невыпуклой глобальной оптимизации. Сложность таких задач заключается в том, что они носят многоэкстремальный характер, т.е. предполагают возможность существования большого числа локальных решений, которые могут быть далеки от глобального [9].
Часто функция стоимости в таких задачах может быть представлена в виде разности двух выпуклых функций [2]. Например, в задаче поиска равновесия по Нэшу в биматричной игре [6]. Также к задачам невыпуклой оптимизации относятся двухуровневые задачи [5], которые возникают во многих областях человеческой деятельности, где появляется необходимость введения процедуры отделения множеств (задача определения кредитоспособности клиента банка, прогнозирование и диагностика в области медицины); задачи линейной дополнительности (задача торможения, задача рыночного равновесия). Кроме того, задачи целочисленного линейного программирования могут быть сведены к задачам невыпуклого программирования.
Ввиду распространенности задач невыпуклой оптимизации разрабатываются алгоритмы и методы их решения.
Рассмотрим основные подходы к решению данных задач:
1. Методы, основанные на разбиении допустимого множества [7]. При данном походе множество допустимых решений разбивается на подмножества, которые при объединении образуют исходное множество. Далее решается исходная задача, но уже на каждом из полученных подмножеств, и затем сравниваются результаты, и выбирается наилучший. Так, множество решений сужается для каждой подзадачи, однако стоит отметить, что растет количество задач.
2. Методы отсекающих плоскостей. Данные методы являются основой для создания точных и приближенных алгоритмов решения задач целочисленного программирования [3]. На каждой новой итерации сужается множество допустимых решений путем отсечения некоторой его части, тем самым уменьшая окрестность оптимального решения. И на полученном множестве отыскивается решение исходной задачи. К достоинствам таких методов можно отнести то, что они не требуют неопределенного размера оперативной памяти для хранения дерева решений, однако при каждом отсечении плоскости увеличивается множество ограничений, кроме того, при приближении к решению отсекающие плоскости начинают практически совпадать, что затрудняет решение задачи.
3. Сеточный метод. Сегодня данный подход широко применяется на практике. На множество допустимых решений накладывается сетка, имеющая различную плотность, которая зависит от предположений о положении глобального оптимума на допустимом множестве. В каждом узле рассчитывается значение целевой функции. Полученные значения сравниваются между собой, и определяется наилучшее из них. Это дает возможность отыскания решения задачи наиболее близкого к точному.
4. Метод ветвей и границ. Метод применяется для отыскания глобального решения задачи. На начальном этапе определяется верхняя и нижняя оценки исходной задачи, если они равны, то метод прекращает работу. В противном случае множество переменных разбивается на несколько подмножеств, объединение которых совпадает с исходным множеством. Данные подмножества образуют дерево подзадач, на которых рекурсивно применяется исходный метод. Алгоритм продолжает работу до тех пор, пока не будет найдено решение задачи. Данный метод является одним из наиболее надежных при решении задач, к его недостаткам можно отнести необходимость полностью решать задачи в узлах дерева, что при достаточно большой размерности задачи требует значительных затрат времени [3].
Для повышения эффективности современные алгоритмы часто являются комбинацией вышеизложенных методов решения задач оптимизации, кроме того в некоторых применяются параллельные вычисления [7].
На текущий момент разработаны последовательные и параллельные алгоритмы МВиГ, решающие задачу вогнутой минимизации дифференцируемой квадратичной функции:
min->/(x), (1)
XED
на непустом, замкнутом, выпуклом многограннике следующего вида:
D = [aL
С учетом проведенных экспериментов, можно предположить, что целесообразнее использовать распредененный алгоритм, приведенный в текущей работе.
1. Андрианова А.А., Коннов И.В., Проблемы теоретической кибернетики. // Материалы XVII международной конференции. Казань: Отечество. 2014. С. 23-26.
2. Груздева Т. В., Стрекаловский А. С., Орлов А. В., Дружинина О. В., Негладкие задачи минимизации разности двух выпуклых функций // Выч. мет. Программирование. Том 12. № 4. 2011. С. 384-396.
3. Захарова Е. М., Минашина И. К. Обзор методов многомерной оптимизации // Информационные процессы. Том 14. № 3. 2014. С. 256¬274.
4. Коннов И.В. Нелинейная оптимизация и вариационные неравенства // Казань: Казан. ун-т, 2013. С.188-195
5. Малышев А. В., Стрекаловский А. С., Глобальный поиск гарантированных решений в квадратично-линейных задачах двухуровневой оптимизации // Изв. Иркутского гос. ун-та. Сер. Математика. Том 4. № 1. 2011. С. 73-82
6. Орлов А. В., Стрекаловский А. С., О поиске ситуаций равновесия в биматричных играх // Автомат. и телемех. Выпуск 2. 2004. С. 55-68.
7. Посыпкин М.А. Методы решения задач конечномерной оптимизации в
распределенной вычислительной среде // САИТ-2009: Труды
конференции. М., 2009. С. 729-740.
8. Розинова. Н.С. Условия оптимальности в задаче максимизации разности двух выпуклых функций. // Известия вузов. Математика. №10. 2010. С. 87-91
9. Стрекаловский А. С. Современные методы решения невыпуклых задач оптимизации и оптимального управления // Известия Иркутского государственного университета. Сер. Математика. Том 2. №1. 2009. С. 245-256.
10. Стрекаловский А. С., Кузнецова А. А., Яковлева Т. В. О численном решении задач невыпуклой оптимизации // Сиб. журн. вычисл. матем. Том 4. № 2. 2001. С. 185-199.
11. Сухарев А. Г., Тимохов А. В., Федоров В. В., Курс методов оптимизации. — Учеб. пособие. — 2-е изд. — М.: ФИЗМАТЛИТ, 2005. — C. 147-160.
12. Konnov I.V. Sign reversion approach to concave minimization problems // Optim Lett. — 2010. — No.4. — P. 491-500.