Материал из MachineLearning.
(Различия между версиями)
|
|
| Строка 22: |
Строка 22: |
| | Следущее (k+1)-e приближение будем искать в виде | | Следущее (k+1)-e приближение будем искать в виде |
| | | | |
| - | u_{k+1}=u_k+alpha_k(overline{u_k}-u_k), 0<alpha_k<1. | + | <tex>u_{k+1}=u_k+alpha_k(overline{u_k}-u_k), 0<alpha_k<1. </tex> |
Версия 17:39, 23 ноября 2008
Постановка задачи
Общие понятния
Если минимизируемая функция дважд непрерывно дифференцируема и производные , J''(u))
просто вычисляются, то можно применять методы минимизации второго порядка, которые используют квадратичную часть разложения функции в ряд Тейлора. Поскольку квадратичная часть разложения аппроксимирует функцию гораздо точнее, чем линейная, то естесвенно ожидать, что методв второго порядка сходятся быстрее, чем методы первого. Метод Ньютона, имеющий квадратичную скорость сходимости на классе сильно выпуклых функций.
Говорят, что последовательность 
сходитcz к 
с линейной скоростью или со скоростью геометрической прогресси (со знаменателем q), если начиная с некоторго номера, выполняется неравенство ;)
при выполнении неравенства 
, где 
, говорят о сверхлинейной скорости сходимости последованояти к , а если здесь 
, т. е. 
, то говорят о скорости сходимсоти порядка s. При s=2, говорят о квадратичной скорости сходимости.
Метод Ньютона
Рассмотирим метод Ньютона для задачи
-> inf; u U,)
где  C^2(U))
, U - выпуклое замкнутое множество из E^n. Пусть 
∈U - некоторое начальное приближение. Если известно k-е приближение 
, то приращение функции J(u)∈ )
в точек можно представить в виде
-J(u_k)=<J'(u_k),u-u_k>+1/2*<J''(u_k)(u-u_k),u-u_k>+o(|u-u_k|^2))
Возьмем квадратичную часть этого приращения
=<J'(u_k),u-u_k>+1/2*<J''(u_k)(u-u_k),u-u_k> )
и определим вспомогательное приближение из условий
=inf J_k(u).)
Следущее (k+1)-e приближение будем искать в виде
, 0<alpha_k<1. )
