形式陈述
对无约束可微最小化,拟 Newton 法把一阶方程 公理库 一阶最优性条件 First-order optimality condition · Variational inequality optimality condition 以梯度和所有可行方向的非负内积充要刻画可微凸问题的全局极小点。 ∇ f ( x ) = 0 作为目标,却不在每轮计算精确 Hessian。给定 x 0 与对称正定矩阵 公理库 正定与半正定矩阵 Positive definite matrix · Positive semidefinite matrix · PSD matrix 由二次能量严格为正或非负定义的实对称与复 Hermitian 矩阵。 B 0 ,第 k 轮求解
B k p k = − g k , g k = ∇ f ( x k ) , 选择步长 α k > 0 ,令
x k + 1 = x k + α k p k , s k = x k + 1 − x k , y k = g k + 1 − g k . 随后用 s k , y k 更新 B k + 1 ,并至少满足割线方程
B k + 1 s k = y k . B k 近似 Hessian;也可直接维护 H k ≈ B k − 1 ,以 p k = − H k g k 更新。输入还包括矩阵更新规则、Wolfe 线搜索 公理库 线搜索的 Wolfe 条件 Wolfe conditions · Wolfe line-search conditions · Strong Wolfe conditions 以充分下降与方向导数曲率两项不等式判定线搜索步长是否可接受。 、容差与预算,输出应带梯度范数、线搜索/曲率状态和迭代次数。方法族由这些接口定义,不自带某个统一收敛率。
直觉
Newton 法用局部二次模型的精确曲率重新缩放梯度,但形成、存储和分解 Hessian 可能昂贵。拟 Newton 法从已经支付的两次梯度中提取曲率:沿实际位移 s k ,梯度变化 y k 近似 ∇ 2 f ( x k ) s k 。割线方程要求下一矩阵至少在刚走过的方向上复现这一观测,而对未观测方向作最小、结构化的修改。
正定近似使 g k T p k = − g k T B k − 1 g k < 0 ,从而方向可交给下降线搜索;矩阵若变成不定,方向可能上升。拟 Newton 的优势不是“无二阶信息”,而是用梯度差积累近似二阶信息。它比普通梯度下降多存状态并承担矩阵代数成本,却常在局部显著减少迭代数;是否更省时间取决于维数、稀疏性与梯度代价。
例子与边界
取一维
f ( x ) = 1 4 x 4 + 1 2 x 2 , g ( x ) = x 3 + x . 从 x 0 = 1 、B 0 = 1 出发,方向 p 0 = − 2 。若线搜索接受 α 0 = 1 / 4 ,则 x 1 = 1 / 2 、g 1 = 5 / 8 ,所以
s 0 = − 1 2 , y 0 = − 11 8 , B 1 = y 0 s 0 = 11 4 = 2.75 . 一维割线方程唯一确定新曲率。下一方向
p 1 = − g 1 B 1 = − 5 22 ≈ − 0.2273 比未缩放的 − 0.625 保守;取完整步得到 x 2 = 3 / 11 ,函数值从 f ( 1 / 2 ) = 0.140625 降到约 0.03857 。B 1 是区间上平均 Hessian,而非 f ″ ( x 1 ) = 1.75 或 f ″ ( x 0 ) = 4 ,说明割线数据只约束走过方向。
多维中一条割线不足以唯一决定矩阵,必须选择更新准则。若 s k T y k ≤ 0 ,正定保持可能失败;这在非凸区域、错误梯度或不合适步长下都会出现。线搜索找不到可接受步时,方法应返回失败或采用明确的阻尼策略,不能除以接近零的曲率继续。局部超线性结论还需要解附近 Hessian Lipschitz、极小点 Hessian 正定、步长最终为一及 Dennis–Moré 型逼近条件,不能由“用了梯度差”自动推出。
推论与应用
不同更新由同一接口产生不同方法。BFGS 更新 公理库 BFGS 更新 BFGS update · Broyden-Fletcher-Goldfarb-Shanno update 以满足割线方程的对称秩二修正更新 Hessian 或逆 Hessian 正定近似。 在割线、对称与正定要求下选取特定低秩修正,是本页方法族的一种具体实现;DFP、SR1 则采用不同的矩阵距离或秩结构。把这些公式都称作“拟 Newton”可以共享全局框架,但比较时必须报告维护 B 还是 H 、是否跳过曲率对、线搜索条件和有限精度重启。
若目标是严格凸二次型,精确线搜索下经典 BFGS 与共轭方向结构相关,并在精确算术中至多 n 步恢复解;一般非线性问题只在局部近似这个情形。有限内存 L-BFGS 不存 n × n 矩阵,而保留最近若干 ( s i , y i ) 并用双循环计算 H k g k ,适合大规模稠密变量。有限内存是成本选择,不改变曲率条件和线搜索失败边界。
参考资料
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。