RUS  ENG
Полная версия
ЖУРНАЛЫ // Компьютерные исследования и моделирование // Архив

Компьютерные исследования и моделирование, 2022, том 14, выпуск 2, страницы 417–444 (Mi crm976)

Эта публикация цитируется в 2 статьях

МАТЕМАТИЧЕСКИЕ ОСНОВЫ И ЧИСЛЕННЫЕ МЕТОДЫ МОДЕЛИРОВАНИЯ

Применение градиентных методов оптимизации для решения задачи Коши для уравнения Гельмгольца

Н. В. Плетневa, П. Е. Двуреченскийb, А. В. Гасниковacd

a Московский физико-технический институт, Россия, 141701, г. Долгопрудный, Институтский пер., д. 9
b Институт прикладного анализа и стохастики им. Вейерштрасса, Германия, 10117, г. Берлин, Моренштрассе, д. 39
c Кавказский математический центр, Адыгейский государственный университет, Россия, 385000, Республика Адыгея, г. Майкоп, ул. Первомайская, д. 208
d Институт проблем передачи информации им. А. А. Харкевича Российской академии наук, Россия, 127051, г. Москва, Большой Каретный переулок, д. 19, стр. 1

Аннотация: Статья посвящена изучению применения методов выпуклой оптимизации для решения задачи Коши для уравнения Гельмгольца, которая является некорректной, поскольку уравнение относится к эллиптическому типу. Задача Коши формулируется как обратная задача и сводится к задаче выпуклой оптимизации в гильбертовом пространстве. Оптимизируемый функционал и его градиент вычисляются с помощью решения краевых задач, которые, в свою очередь, корректны и могут быть приближенно решены стандартными численными методами, такими как конечно-разностные схемы и разложения в ряды Фурье. Экспериментально исследуются сходимость применяемого быстрого градиентного метода и качество получаемого таким образом решения. Эксперимент показывает, что ускоренный градиентный метод — метод подобных треугольников — сходится быстрее, чем неускоренный метод. Сформулированы и доказаны теоремы о вычислительной сложности полученных алгоритмов. Установлено, что разложения в ряды Фурье превосходят конечно-разностные схемы по скорости вычислений и улучшают качество получаемого решения. Сделана попытка использовать рестарты метода подобных треугольников после уменьшения невязки функционала вдвое. В этом случае сходимость не улучшается, что подтверждает отсутствие сильной выпуклости. Эксперименты показывают, что неточность вычислений более адекватно описывается аддитивной концепцией шума в оракуле первого порядка. Этот фактор ограничивает достижимое качество решения, но ошибка не накапливается. Полученные результаты показывают, что использование ускоренных градиентных методов оптимизации позволяет эффективно решать обратные задачи.

Ключевые слова: обратные задачи, выпуклая оптимизация, оптимизация в гильбертовом пространстве, методы первого порядка, быстрые градиентные методы, неточный оракул.

УДК: 519.85

Поступила в редакцию: 13.02.2022
Принята в печать: 13.02.2022

DOI: 10.20537/2076-7633-2022-14-2-417-444



© МИАН, 2024