Skip to content

拟 Newton 法

Quasi-Newton method · Variable metric method

由相邻梯度差逐步近似 Hessian 或其逆,从而构造曲率缩放搜索方向的方法族。

条目类型
方法

形式陈述

对无约束可微最小化,拟 Newton 法把一阶方程 f(x)=0 作为目标,却不在每轮计算精确 Hessian。给定 x0 与对称正定矩阵 B0,第 k 轮求解

Bkpk=gk,gk=f(xk),

选择步长 αk>0,令

xk+1=xk+αkpk,sk=xk+1xk,yk=gk+1gk.

随后用 sk,yk 更新 Bk+1,并至少满足割线方程

Bk+1sk=yk.

Bk 近似 Hessian;也可直接维护 HkBk1,以 pk=Hkgk 更新。输入还包括矩阵更新规则、Wolfe 线搜索、容差与预算,输出应带梯度范数、线搜索/曲率状态和迭代次数。方法族由这些接口定义,不自带某个统一收敛率。

直觉

Newton 法用局部二次模型的精确曲率重新缩放梯度,但形成、存储和分解 Hessian 可能昂贵。拟 Newton 法从已经支付的两次梯度中提取曲率:沿实际位移 sk,梯度变化 yk 近似 2f(xk)sk。割线方程要求下一矩阵至少在刚走过的方向上复现这一观测,而对未观测方向作最小、结构化的修改。

正定近似使 gkTpk=gkTBk1gk<0,从而方向可交给下降线搜索;矩阵若变成不定,方向可能上升。拟 Newton 的优势不是“无二阶信息”,而是用梯度差积累近似二阶信息。它比普通梯度下降多存状态并承担矩阵代数成本,却常在局部显著减少迭代数;是否更省时间取决于维数、稀疏性与梯度代价。

例子与边界

取一维

f(x)=14x4+12x2,g(x)=x3+x.

x0=1B0=1 出发,方向 p0=2。若线搜索接受 α0=1/4,则 x1=1/2g1=5/8,所以

s0=12,y0=118,B1=y0s0=114=2.75.

一维割线方程唯一确定新曲率。下一方向

p1=g1B1=5220.2273

比未缩放的 0.625 保守;取完整步得到 x2=3/11,函数值从 f(1/2)=0.140625 降到约 0.03857B1 是区间上平均 Hessian,而非 f(x1)=1.75f(x0)=4,说明割线数据只约束走过方向。

多维中一条割线不足以唯一决定矩阵,必须选择更新准则。若 skTyk0,正定保持可能失败;这在非凸区域、错误梯度或不合适步长下都会出现。线搜索找不到可接受步时,方法应返回失败或采用明确的阻尼策略,不能除以接近零的曲率继续。局部超线性结论还需要解附近 Hessian Lipschitz、极小点 Hessian 正定、步长最终为一及 Dennis–Moré 型逼近条件,不能由“用了梯度差”自动推出。

推论与应用

不同更新由同一接口产生不同方法。BFGS 更新在割线、对称与正定要求下选取特定低秩修正,是本页方法族的一种具体实现;DFP、SR1 则采用不同的矩阵距离或秩结构。把这些公式都称作“拟 Newton”可以共享全局框架,但比较时必须报告维护 B 还是 H、是否跳过曲率对、线搜索条件和有限精度重启。

若目标是严格凸二次型,精确线搜索下经典 BFGS 与共轭方向结构相关,并在精确算术中至多 n 步恢复解;一般非线性问题只在局部近似这个情形。有限内存 L-BFGS 不存 n×n 矩阵,而保留最近若干 (si,yi) 并用双循环计算 Hkgk,适合大规模稠密变量。有限内存是成本选择,不改变曲率条件和线搜索失败边界。

参考资料
  • Jorge Nocedal and Stephen J. Wright, Numerical Optimization, 2nd ed., Springer, 2006,Ch. 6,quasi-Newton methods and local convergence。
  • J. E. Dennis Jr. and Robert B. Schnabel, Numerical Methods for Unconstrained Optimization and Nonlinear Equations, SIAM Classics, 1996,Chs. 6 and 9,secant updates and convergence。
  • Roger Fletcher, Practical Methods of Optimization, 2nd ed., Wiley, 1987,Ch. 3,variable metric methods。
关系图谱15 个相邻概念 · 3 类关系

拖动节点调整位置。

显示关系

显示:依赖

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