Участник:Margarita/Градиентный спуск

Материал из Алговики
Перейти к навигации Перейти к поиску

автор: М.В.Зайцева

Постановка задачи

Рассмотрим задачу поиска минимума функции [math]\displaystyle{ f(x)\colon \mathbb{R}^n \to \mathbb{R} }[/math]

[math]\displaystyle{ \begin{align} f(x) \rightarrow \min_{x \in \mathbb{R}^n} \end{align} }[/math]

Будем считать, что у функции [math]\displaystyle{ f(x) }[/math] существуют частные производные.

Если в задаче требуется найти максимум, вместо [math]\displaystyle{ f(x) }[/math] следует брать [math]\displaystyle{ -f(x) }[/math].

Общее описание алгоритма

Градиентный спуск — это метод нахождения локального экстремума функции с помощью движения вдоль направления градиента. На каждой итерации для минимизации функции в направлении градиента используются методы одномерной оптимизации.

Градиентный спуск активно применяется в области машинного обучения, так как часть процесса машинного обучения — это поиск максимальной точности или минимизация частоты ошибок. Метод градиентного спуска используется для поиска минимальной ошибки путем минимизации функции цены.

Математическое описание алгоритма

Основная идея [1] метода состоит в том, чтобы идти в направлении антиградиента [math]\displaystyle{ -\nabla f }[/math] (отсюда название метода), то есть в направлении, в котором функция быстрее всего убывает при бесконечно малом движении из данной точки:

[math]\displaystyle{ \begin{align} x_{k+1} = x_{k} - \alpha_{k}\nabla f(x_{k}), \end{align} }[/math]

где [math]\displaystyle{ \alpha_{k} }[/math] выбирается по одному из следующих правил:

  • постоянной, в таком случае метод может не сходиться;
  • дробным шагом, то есть во время итерации длина шага делится на некоторое число;
  • наискорейшим спуском[2], где [math]\displaystyle{ \alpha_{k} = \operatorname{arg}\min_{\alpha} f(x_{k} - \alpha\nabla f(x_{k})) }[/math].

Вычислительное ядро алгоритма

Макроструктура алгоритма

Схема реализации последовательного алгоритма

Алгоритм:

  1. Выбираем начальную точку [math]\displaystyle{ x_{0}, \varepsilon }[/math].
  2. Вычисляем [math]\displaystyle{ x_{k+1} = x_{k} - \alpha_{k}\nabla f(x_{k}) }[/math], где [math]\displaystyle{ \alpha_{k} }[/math] выбирается одним из описанных выше способом.
  3. Если выполнено условие останова, то возвращаем [math]\displaystyle{ x_{k+1} }[/math], иначе переходим к шагу 2.

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

  • [math]\displaystyle{ \| x_{k+1} - x_{k} \| \lt \varepsilon }[/math];
  • [math]\displaystyle{ \| f(x_{k+1}) - f(x_{k}) \| \lt \varepsilon }[/math];
  • [math]\displaystyle{ \| \nabla f(x_{k+1}) \| \lt \varepsilon }[/math].

Последовательная сложность алгоритма

Информационный граф

Ресурс параллелизма алгоритма

Входные и выходные данные алгоритма

Свойства алгоритма

Программная реализация алгоритма

Особенности реализации последовательного алгоритма

Локальность данных и вычислений

Возможные способы и особенности параллельной реализации алгоритма

Масштабируемость алгоритма и его реализации

Динамические характеристики и эффективность реализации алгоритма

Выводы для классов архитектур

Существующие реализации алгоритма

Литература

  1. ↑ , Н. Н. Калиткин “Численные методы.” Москва «Наука», 1978
  2. ↑ A. D. Belegundu and T. R. Chandrupatla. Optimization Concepts and Applications in Engineering, chapter 3. Prentice Hall, 1999