Квазиньютоновский метод - Quasi-Newton method
Было высказано предположение, что Квазиньютоновский метод наименьших квадратов и Квазиньютоновский метод обратных наименьших квадратов быть слился в эту статью. (Обсуждать) Предлагается с июля 2020 года. |
Квазиньютоновские методы - это методы, используемые для поиска нулей или локальных максимумов и минимумов функций в качестве альтернативы методу Ньютона. Их можно использовать, если Якобиан или же Гессен недоступен или слишком дорог для вычислений на каждой итерации. Полный" Метод Ньютона требует якобиан для поиска нулей или гессиан для поиска экстремумов.
Поиск нулей: поиск корней
Метод Ньютона найти нули функции нескольких переменных определяется как , куда это левый обратный из Матрица якобиана из оценивается для .
Строго говоря, любой метод, заменяющий точный якобиан с приближением - это квазиньютоновский метод.[1] Например, хордовый метод (где заменяется на для всех итераций) - простой пример. Приведенные ниже методы для оптимизация относятся к важному подклассу квазиньютоновских методов, методов секущих.[2]
Использование методов, разработанных для поиска экстремумов, для поиска нулей не всегда является хорошей идеей, поскольку большинство методов, используемых для поиска экстремумов, требуют, чтобы используемая матрица была симметричной. Хотя это верно в контексте поиска экстремумов, это редко выполняется при поиске нулей. «Хорошие» и «плохие» методы Бройдена два метода, обычно используемых для поиска экстремумов, которые также могут применяться для поиска нулей. Другие методы, которые можно использовать: метод обновления столбца, то обратный метод обновления столбца, то квазиньютоновский метод наименьших квадратов и квазиньютоновский метод обратных наименьших квадратов.
Совсем недавно квазиньютоновские методы стали применяться для поиска решения множественных связанных систем уравнений (например, задач взаимодействия жидкость-структура или задач взаимодействия в физике). Они позволяют найти решение, решая каждую составляющую систему отдельно (что проще, чем глобальная система) циклическим, итеративным способом, пока не будет найдено решение глобальной системы.[2][3]
Поиск экстремумов: оптимизация
Отметив, что поиск минимума или максимума скалярнозначной функции есть не что иное, как поиск нулей градиент этой функции квазиньютоновские методы могут быть легко применены для поиска экстремумов функции. Другими словами, если это градиент , затем поиск нулей вектор-функции соответствует поиску экстремумов скалярнозначной функции ; якобиан теперь становится гессенской . Главное отличие в том, что матрица Гессе - это симметричная матрица, в отличие от якобиана, когда поиск нулей. Большинство квазиньютоновских методов, используемых при оптимизации, используют это свойство.
В оптимизация, квазиньютоновские методы (частный случай методы с переменной метрикой) являются алгоритмами поиска локальных максимумы и минимумы из функции. Квазиньютоновские методы основаны на Метод Ньютона найти стационарная точка функции, где градиент равен 0. Метод Ньютона предполагает, что функцию можно локально аппроксимировать как квадратичный в области вокруг оптимума и использует первую и вторую производные для нахождения стационарной точки. В более высоких измерениях метод Ньютона использует градиент и Матрица Гессе второй производные минимизируемой функции.
В квазиньютоновских методах матрица Гессе не требует вычисления. Гессен обновляется путем анализа последовательных векторов градиента. Квазиньютоновские методы являются обобщением секущий метод найти корень первой производной для многомерных задач. В нескольких измерениях уравнение секущей недоопределенный, и квазиньютоновские методы различаются тем, как они ограничивают решение, обычно путем добавления простого обновления низкого ранга к текущей оценке гессиана.
Первый квазиньютоновский алгоритм был предложен Уильям К. Дэвидон, физик, работающий на Аргоннская национальная лаборатория. Он разработал первый квазиньютоновский алгоритм в 1959 году: Формула обновления DFP, который позже был популяризирован Флетчером и Пауэллом в 1963 году, но сегодня используется редко. Наиболее распространенными квазиньютоновскими алгоритмами в настоящее время являются Формула SR1 (для «симметричного ранга один») BHHH метод, широко распространенный Метод BFGS (независимо предложено Бройденом, Флетчером, Гольдфарбом и Шенно в 1970 году) и его расширение с низким объемом памяти L-BFGS. Класс Бройдена представляет собой линейную комбинацию методов DFP и BFGS.
Формула SR1 не гарантирует, что матрица обновления будет поддерживать положительная определенность и может использоваться для неопределенных задач. В Метод Бройдена не требует, чтобы матрица обновления была симметричной, и используется для нахождения корня общей системы уравнений (а не градиента) путем обновления Якобиан (а не гессенское).
Одно из главных преимуществ квазиньютоновских методов перед Метод Ньютона это то Матрица Гессе (или, в случае квазиньютоновских методов, его приближение) переворачивать не нужно. Метод Ньютона и его производные, такие как методы внутренней точки, требуют инвертирования гессиана, что обычно реализуется путем решения система линейных уравнений и часто бывает довольно дорого. Напротив, квазиньютоновские методы обычно дают оценку прямо.
Как в Метод Ньютона, используется приближение второго порядка для нахождения минимума функции . В Серия Тейлор из вокруг итерации
где () это градиент, и приближение к Матрица Гессе[4]. Градиент этого приближения (по ) является
и установка этого градиента на ноль (что является целью оптимизации) обеспечивает шаг Ньютона:
Гессенское приближение выбран, чтобы удовлетворить
который называется секущее уравнение (серия Тейлора самого градиента). Более чем в одном измерении является недоопределенный. В одном измерении решение для и применение шага Ньютона с обновленным значением эквивалентно секущий метод. Различные квазиньютоновские методы различаются выбором решения секущего уравнения (в одном измерении все варианты эквивалентны). Большинство методов (но с исключениями, такими как Метод Бройдена ) искать симметричное решение (); кроме того, варианты, перечисленные ниже, могут быть мотивированы поиском обновления это как можно ближе к в некоторых норма; то есть, , куда есть некоторые положительно определенная матрица что определяет норму. Примерное начальное значение часто бывает достаточно для достижения быстрой сходимости, хотя нет общей стратегии выбора [5]. Обратите внимание, что должно быть положительно-определенным. Неизвестный обновляется с применением шага Ньютона, рассчитанного с использованием текущей приближенной матрицы Гессе :
- , с выбран для удовлетворения Условия Вульфа;
- ;
- Градиент, вычисленный в новой точке , и
используется для обновления приблизительного гессенского , или прямо обратное с использованием Формула Шермана – Моррисона.
- Ключевым свойством обновлений BFGS и DFP является то, что если положительно определен, и выбирается так, чтобы выполнялись условия Вульфа, то также положительно определен.
Наиболее популярные формулы обновления:
Другие методы - это метод Пирсона, метод Маккормика, симметричный метод Пауэлла Бройдена (PSB) и метод Гринштадта.[2]
Связь с обращением матрицы
Когда - выпуклая квадратичная функция с положительно определенным гессианом , можно было бы ожидать, что матрицы квазиньютоновским методом сходиться к обратному гессиану . Это действительно так для класса квазиньютоновских методов, основанных на обновлениях с наименьшими изменениями.[6]
Известные реализации
Реализации квазиньютоновских методов доступны на многих языках программирования. Известные реализации включают:
- GNU Octave использует форму BFGS в своих
fsolve
функция, с регион доверия расширения. - Mathematica включает квазиньютоновские решатели.[7]
- В Библиотека NAG содержит несколько процедур[8] для минимизации или максимизации функции[9] которые используют квазиньютоновские алгоритмы.
- В MATLAB Панель инструментов оптимизации, то
fminunc
функция использует (среди других методов) BFGS квазиньютоновский метод.[10] Многие из методов с ограничениями из набора инструментов оптимизации используют BFGS и вариант L-BFGS.[11] - р с
оптим
универсальная программа оптимизатора использует BFGS метод с использованиемmethod = "BFGS"
.[12] - Scipy.optimize имеет fmin_bfgs. в SciPy расширение на Python, то
scipy.optimize.minimize
функция включает, среди других методов, BFGS выполнение.[13]
Смотрите также
- Метод BFGS
- Метод Бройдена
- Формула обновления DFP
- Метод Ньютона
- Метод Ньютона в оптимизации
- Квазиньютоновский метод наименьших квадратов
- Формула SR1
Рекомендации
- ^ Бройден, К. Г. (1972). «Квазиньютоновские методы». В Мюррей, W. (ред.). Численные методы безусловной оптимизации. Лондон: Academic Press. С. 87–106. ISBN 0-12-512250-0.
- ^ а б c Хельтерман, Роб (2009). «Аналитическое исследование квазиньютоновского метода наименьших квадратов для задач взаимодействия». Докторская диссертация, Гентский университет. Получено 2014-08-14.
- ^ Роб Хэлтерман, Дирк Ван Эстер, Даан Верлейен (2015). «Ускорение решения физической модели внутри токамака с помощью (обратного) метода обновления столбца». Журнал вычислительной и прикладной математики. 279: 133–144. Дои:10.1016 / j.cam.2014.11.005.CS1 maint: использует параметр авторов (ссылка на сайт)
- ^ https://mathinsight.org/taylors_theorem_multivariable_introduction
- ^ Нокедаль, Хорхе; Райт, Стивен Дж. (2006). Численная оптимизация. Нью-Йорк: Спрингер. стр.142. ISBN 0-387-98793-2.
- ^ Роберт Мансел Гауэр; Питер Рихтарик (2015). «Рандомизированные квазиньютоновские обновления представляют собой алгоритмы обращения линейно сходящейся матрицы». arXiv:1602.01768 [math.NA ].
- ^ http://reference.wolfram.com/mathematica/tutorial/UnconstrainedOptimizationQuasiNewtonMethods.html
- ^ Группа численных алгоритмов. "Указатель ключевых слов: квазиньютон". Руководство библиотеки NAG, Mark 23. Получено 2012-02-09.
- ^ Группа численных алгоритмов. «E04 - Минимизация или максимизация функции» (PDF). Руководство библиотеки NAG, Mark 23. Получено 2012-02-09.
- ^ http://www.mathworks.com/help/toolbox/optim/ug/fminunc.html
- ^ http://www.mathworks.com/help/toolbox/optim/ug/brnoxzl.html
- ^ [1]
- ^ http://docs.scipy.org/doc/scipy/reference/generated/scipy.optimize.minimize.html
дальнейшее чтение
- Bonnans, J. F .; Gilbert, J. Ch .; Лемарешаль, К.; Сагастизабал, К.А. (2006). Численная оптимизация: теоретические и численные аспекты (Второе изд.). Springer. ISBN 3-540-35445-X.
- Флетчер, Роджер (1987), Практические методы оптимизации (2-е изд.), Нью-Йорк: Джон Уайли и сыновья, ISBN 978-0-471-91547-8.
- Нокедаль, Хорхе; Райт, Стивен Дж. (1999). «Квазиньютоновские методы». Численная оптимизация. Нью-Йорк: Спрингер. С. 192–221. ISBN 0-387-98793-2.
- Press, W. H .; Теукольский, С. А .; Vetterling, W. T .; Фланнери, Б. П. (2007). «Раздел 10.9. Квазиньютон или методы переменной метрики в многомерности». Числовые рецепты: искусство научных вычислений (3-е изд.). Нью-Йорк: Издательство Кембриджского университета. ISBN 978-0-521-88068-8.
- Весы, Л. Э. (1985). Введение в нелинейную оптимизацию. Нью-Йорк: Макмиллан. С. 84–106. ISBN 0-333-32552-4.