Формула Дэвидона – Флетчера – Пауэлла

редактировать

Формула Дэвидона – Флетчера – Пауэлла (или DFP ; названный в честь Уильяма К. Дэвидона, Роджера Флетчера и Майкла Дж.Д. Пауэлла ) находит решение секущего уравнения, которое наиболее близко к текущей оценке и удовлетворяет условие кривизны. Это был первый квазиньютоновский метод, обобщающий метод секущих на многомерную задачу. Это обновление поддерживает симметрию и положительную определенность матрицы Гессе.

. Для функции f (x) {\ displaystyle f (x)}f (x) ее gradient (∇ f {\ displaystyle \ nabla f}\ набла f ) и positive-definite матрица Гессе B {\ displaystyle B}В , ряд Тейлора равен

f (xk + sk) = f (xk) + ∇ f (xk) T sk + 1 2 sk TB sk +…, {\ displaystyle f (x_ {k} + s_ {k}) = f (x_ {k}) + \ nabla f (x_ {k}) ^ {T} s_ {k} + {\ frac {1} {2}} s_ {k} ^ {T} {B} s_ {k} + \ dots,}{\ displaystyle f (x_ {k} + s_ {k}) = f (x_ {k}) + \ набла f (x_ {k}) ^ {T} s_ {k} + {\ frac {1} {2}} s_ {k} ^ {T} {B} s_ {k} + \ dots,}

и ряд Тейлора самого градиента (секущее уравнение)

∇ f (xk + sk) Знак равно ∇ е (xk) + B sk +… {\ displaystyle \ nabla f (x_ {k} + s_ {k}) = \ nabla f (x_ {k}) + Bs_ {k} + \ dots}{\ displaystyle \ nabla f (x_ {k} + s_ {k}) = \ nabla f (x_ {k}) + Bs_ {k} + \ dots}

используется для обновления B {\ displaystyle B}В .

Формула DFP находит решение, которое является симметричным, положительно определенным и наиболее близко к текущему приблизительному значению B k {\ displaystyle B_ {k}}B_ {k} :

В К + 1 знак равно (I - γ kyksk T) B k (I - γ kskyk T) + γ kykyk T, {\ displaystyle B_ {k + 1} = (I- \ gamma _ {k} y_ {k} s_ {k} ^ {T}) B_ {k} (I- \ gamma _ {k} s_ {k} y_ {k} ^ {T}) + \ gamma _ {k} y_ {k} y_ {k} ^ {T},}{\ displaystyle B_ {k + 1} = (I- \ gamma _ {k} y_ {k} s_ {k} ^ {T}) B_ {k} (I- \ gamma _ { k} s_ {k} y_ {k} ^ {T}) + \ gamma _ {k} y_ {k} y_ {k} ^ {T},}

где

yk = ∇ f (xk + sk) - ∇ f (xk), {\ displaystyle y_ {k} = \ nabla f (x_ {k} + s_ {k}) - \ nabla f (x_ {k}),}{\ displaystyle y_ {k} = \ nabla f (x_ {k} + s_ {k}) - \ nabla f (x_ {k}),}
γ k = 1 yk T sk, {\ displaystyle \ gamma _ { k} = {\ frac {1} {y_ {k} ^ {T} s_ {k}}},}{\ displaystyle \ gamma _ {k} = {\ frac {1} {y_ {k} ^ {T} s_ {k}}},}

и B k {\ displaystyle B_ {k}}B_ {k} симметричная и положительно определенная матрица.

Соответствующее обновление обратного приближения Гессе H k = B k - 1 {\ displaystyle H_ {k} = B_ {k} ^ {- 1}}{\ displaystyle H_ {k} = B_ {k} ^ {- 1}} определяется как

H k + 1 = H k - H kykyk TH kyk TH kyk + sksk T yk T sk. {\ displaystyle H_ {k + 1} = H_ {k} - {\ frac {H_ {k} y_ {k} y_ {k} ^ {T} H_ {k}} {y_ {k} ^ {T} H_ {k} y_ {k}}} + {\ frac {s_ {k} s_ {k} ^ {T}} {y_ {k} ^ {T} s_ {k}}}.}{\ displaystyle H_ {k + 1} = H_ {k} - {\ frac {H_ {k} y_ {k} y_ {k} ^ {T} H_ {k}} {y_ {k} ^ {T} H_ {k} y_ {k}}} + {\ frac {s_ {k} s_ {k} ^ {T}} {y_ {k} ^ {T} s_ {k}}}.}

B {\ displaystyle B}В считается положительно определенным, а векторы sk T {\ displaystyle s_ {k} ^ {T}}s_ {k} ^ {T} и y {\ displaystyle y}y должен удовлетворять условию кривизны

sk T yk = sk TB sk>0. {\ displaystyle s_ {k} ^ {T} y_ {k} = s_ {k} ^ {T} Bs_ {k}>0.}{\displaystyle s_{k}^{T}y_{k}=s_{k}^{T}Bs_{k}>0.}

Формула DFP довольно эффективна, но вскоре была заменена формулой Формула Бройдена – Флетчера – Гольдфарба – Шанно, которая является его двойным (меняя роли y и s).

См. Также
Ссылки
Дополнительная литература
Последняя правка сделана 2021-05-17 03:53:50
Содержание доступно по лицензии CC BY-SA 3.0 (если не указано иное).
Обратная связь: support@alphapedia.ru
Соглашение
О проекте