形式陈述
考虑无约束光滑优化问题公理库优化问题Optimization problem在可行解集合上最小化或最大化目标函数的计算问题。 。在 处计算梯度公理库梯度Gradient标量函数微分在内积下对应的向量。 ,选择对称模型矩阵 ,建立
可以是真实Hessian公理库Hessian 矩阵Hessian matrix · Hessian标量函数二阶 Fréchet 导数在坐标中的矩阵表示,用于描述局部曲率与二次近似。,也可以是曲率近似。内层信赖域二次子问题公理库信赖域二次子问题Trust-region subproblem在有限半径内最小化可能不定的二次模型,以半正定移位和互补关系认证全局解,并处理奇异的困难情形。提供候选 。若 ,要求候选至少取得固定比例的 Cauchy 下降:存在与迭代无关的 ,使
Cauchy 步是沿 在球内最小化模型的步。上式保证内层确实利用了当前的一阶下降信息,不要求每轮精确求出二次子问题的全局解。
再实际计算 ,令
下面固定一套可执行规则。输入 ,接受阈值 :
- 若 ,取 ;否则保留 。
- 若 ,取 。
- 若 且步到达球面,取 ;其余情形保持半径。
接受点与更新半径是两个判断。例如 时,步被接受,但下一轮半径会缩小。若候选预测下降非正、求值非有限或超出函数定义域,按失败候选处理并缩小半径,不计算一个没有含义的比值。
直觉
局部模型像一张附有使用范围的地图。内层负责在地图允许的区域里找到好方向,外层负责检查地图预测与真实地形是否一致。一次坏候选不会强迫算法移动;位置可以不变,而地图的使用半径变小。
线搜索公理库线搜索的 Wolfe 条件Wolfe conditions · Wolfe line-search conditions · Strong Wolfe conditions以充分下降与方向导数曲率两项不等式判定线搜索步长是否可接受。通常先给定下降方向,再选择这条射线上的长度。信赖域先给定允许范围,再从模型选方向;半径改变后,方向也可能改变。只有一维或某些特殊模型中,缩半径才恰好等同于沿原方向缩步。
先拒绝,再缩小半径
例子与边界
模型预测下降十,真实目标却上升六
取 ,从 出发。此时 ,二次模型为
半径 时,模型在区间 内的最小点是 。预测下降为 ,但真实函数值 ,实际下降为 。所以 ,拒绝候选,仍留在 ,半径缩成一。
重解同一位置的新子问题,得到 。这次预测下降为三,真实函数值 ,实际下降为二,所以 。候选被接受,下一点是 ,半径保持一。拒绝之后要重新解较小范围内的模型;这里恰好方向没变,是这个一维例子的特点。
下一轮的 ,无约束模型最优步 位于球内。预测下降为 ,实际下降为 ,故
它被接受,但因为没有碰到球面,半径仍不扩大。最终是否达到要求,要继续检查新点梯度;高比值只表示当前模型预测准确。
比值好,不等于已最优
对 使用精确二次模型,任何非零下降步都有 。若在 只走到 ,模型仍完全准确,却离最优点很远。因此终止量应包括 ,而不是只检查 。
反过来,若 ,Cauchy 下降下界退化为零。取 、,梯度为零却是局部极大点。只检验一阶驻性的版本可以在此结束;要求排除负曲率时,还需内层检测 的负方向,并使用相应二阶停止标准。
推论与应用
缩小半径为什么能够修复失真
设 在所有试探线段的邻域内有 -Lipschitz 梯度,且 。Taylor 的一阶余项给出
若 且半径足够小,Cauchy 下界为 。于是
误差是半径的二次量,保证下降是一阶量,所以半径充分小时,比值趋近一。只要梯度还显著非零,算法就不能永远以“模型不可信”为由拒绝所有步。
进一步假设 在迭代水平集上有下界,求值精确,模型矩阵一致有界,而且每轮均满足上述下降要求。则固定规则产生的无限序列满足
证明用反证法:若尾部梯度始终至少为 ,足够小的半径就不会再被缩小。因为每次只减半,所有尾部半径存在一个正下界。无限拒绝又会迫使半径减到必能接受的范围,所以成功步有无限多个。每个成功步通过 Cauchy 下界使 至少下降一个固定正数,最终违反下有界性。这个结论保证任意小的一阶残差会出现;它没有把非凸驻点升级成全局最优点。
工作量与实际停止
一次试探需要模型求解以及一次真实目标求值;接受后通常重算梯度,拒绝时可复用旧点的梯度与模型矩阵。若显式组装并分解稠密 ,主要线性代数成本可达 ;使用截断 CG 时,应报告 Hessian–向量积次数及目标、梯度求值次数,而不是只报外层轮数。
程序可在梯度达到目标容差时返回“一阶条件达到”,在达到预算、半径小到无法产生可分辨位移或内层失败时返回对应状态。小步长本身可能来自尺度失衡、舍入或坏模型,不构成驻性证书。若需要约束优化,必须把约束线性化、可行性恢复或罚函数机制加入外层;仅把普通 Newton 步放进一个球,并没有处理业务约束。
参考资料