Skip to content

BFGS 更新

BFGS update · Broyden-Fletcher-Goldfarb-Shanno update

以满足割线方程的对称秩二修正更新 Hessian 或逆 Hessian 正定近似。

条目类型
算法

形式陈述

拟 Newton 法中,给定位移与梯度差

sk=xk+1xk,yk=f(xk+1)f(xk),

并假设曲率条件 skTyk>0。若 Bk 是 Hessian 的对称正定近似,BFGS 更新为

Bk+1=BkBkskskTBkskTBksk+ykykTykTsk.

它满足割线方程 Bk+1sk=yk。若维护逆近似 Hk=Bk1,令 ρk=(ykTsk)1,则等价的稳定形式为

Hk+1=(IρkskykT)Hk(IρkykskT)+ρkskskT.

输入必须包含 BkHk、有效曲率对及数值容差;输出除新矩阵外还应表明更新、阻尼或拒绝状态。更新公式本身不选择 xk+1,方向和步长仍由外层优化算法负责。

直觉

第一项先从旧矩阵中移除沿 sk 的旧曲率预测,第二项再植入观测到的新关系 skyk;其余方向尽量保留旧信息。修正是秩二的,所以比重算 Hessian 便宜。对称性使它仍可解释为二次模型,割线条件让最近一次梯度变化被精确复现,正定性则保证下一搜索方向下降。

正定保持可以直接看出。对任意非零 z

zTBk+1z=zTBkz(zTBksk)2skTBksk+(zTyk)2skTyk.

前两项由 Bk-内积下的 Cauchy–Schwarz 非负;若恰为零,则 zsk 共线,而最后一项因 skTyk>0 严格为正。因此 Bk+1 仍是正定矩阵。曲率条件不是避免除零的技术细节,而是方向性质的核心。

例子与边界

B0=I,s=(1,0)T,y=(2,1)T,

sTy=2>0,且

B1=IssT+yyT2=(2113/2).

直接相乘得 B1s=(2,1)T=y,割线条件成立;首个顺序主子式为 2、行列式为 2,故矩阵正定。逆公式给

H1=(3/41/21/21),

B1H1=I,可逐项复核两种更新确实一致。这不是换数字展示同一标量公式:非对角项具体显示一个方向的曲率观测如何改变耦合方向。

若改成 y=(1,0)T,则 sTy=1,更新的最后分母为负,正定证明完全失效。非凸目标上负曲率对很常见;常用选择是跳过更新、Powell 阻尼或改用容许不定近似的 SR1,但每种策略都改变实际算法,必须记录。若 sTy 虽正却接近舍入误差,两个大秩一项相减会损害对称性和条件数。线搜索失败、梯度不一致或变量尺度悬殊同样会污染曲率对。

推论与应用

pk 是下降方向且步长满足弱 Wolfe 曲率条件,则

skTyk=αk(ϕ(αk)ϕ(0))αk(c21)ϕ(0)>0,

所以线搜索为 BFGS 正定更新提供可验证的曲率证书。对光滑强凸目标,在标准有界水平集与线搜索条件下可建立全局收敛;在极小点附近再加 Hessian 正定及 Lipschitz 等条件,完整 BFGS 通常具有超线性局部收敛。函数值下降、矩阵正定与超线性速度是三层不同结论。

大规模问题使用 L-BFGS 保存最近 m(si,yi),通过双循环递推计算 Hkgk,把存储从 O(n2) 降到 O(mn)。它复用了 BFGS 的曲率对,却不显式形成本页矩阵;有限记忆会丢弃旧方向信息,因此不能直接声称与完整 BFGS 有相同有限终止性质。实现还应定期检查 sTy、对称误差和搜索方向内积,发现失效时给出可审计状态。

参考资料
  • Charles G. Broyden, “The Convergence of a Class of Double-rank Minimization Algorithms,” Journal of the Institute of Mathematics and Its Applications 6, 1970, 76–90;Roger Fletcher, “A New Approach to Variable Metric Algorithms,” Computer Journal 13(3), 1970, 317–322。
  • Donald Goldfarb, “A Family of Variable-Metric Methods Derived by Variational Means,” Mathematics of Computation 24(109), 1970, 23–26;David F. Shanno, “Conditioning of Quasi-Newton Methods for Function Minimization,” Mathematics of Computation 24(111), 1970, 647–656。
  • Jorge Nocedal and Stephen J. Wright, Numerical Optimization, 2nd ed., Springer, 2006,§§6.1–6.2 and 7.2,BFGS, positive definiteness, and L-BFGS。
关系图谱8 个相邻概念 · 3 类关系

拖动节点调整位置。

显示关系

显示:依赖

  1. 前置三跳
  2. 前置二跳
  3. 前置一跳
  4. 当前条目
  5. 后续一跳
  6. 后续二跳
  7. 后续三跳
文字版关系按与当前条目的最短距离分组
类型化关系