Skip to content

算法Algorithm

信赖域法的接受与半径更新

Trust-region method

用实际下降与预测下降的比值接受或拒绝模型步,并调整信赖半径,使局部二次模型服务于真实目标下降。

形式陈述 ​

考虑无约束光滑优化问题 minf(x)。在 xk 处计算梯度 gk,选择对称模型矩阵 Bk,建立

mk(p)=f(xk)+gkTp+12pTBkp,‖p‖≤Δk.

Bk 可以是真实Hessian,也可以是曲率近似。内层信赖域二次子问题提供候选 pk。若 gk≠0,要求候选至少取得固定比例的 Cauchy 下降:存在与迭代无关的 κ∈(0,1],使

predk=mk(0)−mk(pk)≥κ2‖gk‖min{Δk,‖gk‖1+‖Bk‖}>0.

Cauchy 步是沿 −gk 在球内最小化模型的步。上式保证内层确实利用了当前的一阶下降信息,不要求每轮精确求出二次子问题的全局解。

再实际计算 f(xk+pk),令

aredk=f(xk)−f(xk+pk),ρk=aredkpredk.

下面固定一套可执行规则。输入 0<Δ0≤Δmax,接受阈值 η=0.1:

  1. 若 ρk≥η,取 xk+1=xk+pk;否则保留 xk+1=xk。
  2. 若 ρk<1/4,取 Δk+1=Δk/2。
  3. 若 ρk>3/4 且步到达球面,取 Δk+1=min(2Δk,Δmax);其余情形保持半径。

接受点与更新半径是两个判断。例如 ρk=0.2 时,步被接受,但下一轮半径会缩小。若候选预测下降非正、求值非有限或超出函数定义域,按失败候选处理并缩小半径,不计算一个没有含义的比值。

直觉

局部模型像一张附有使用范围的地图。内层负责在地图允许的区域里找到好方向,外层负责检查地图预测与真实地形是否一致。一次坏候选不会强迫算法移动;位置可以不变,而地图的使用半径变小。

线搜索通常先给定下降方向,再选择这条射线上的长度。信赖域先给定允许范围,再从模型选方向;半径改变后,方向也可能改变。只有一维或某些特殊模型中,缩半径才恰好等同于沿原方向缩步。

先拒绝,再缩小半径
例子与边界

模型预测下降十,真实目标却上升六 ​

取 f(x)=x4−2x2+x,从 x0=0 出发。此时 f(0)=0,g0=1,B0=−4,二次模型为

m0(p)=p−2p2.

半径 Δ0=2 时,模型在区间 [−2,2] 内的最小点是 p=−2。预测下降为 0−(−10)=10,但真实函数值 f(−2)=6,实际下降为 −6。所以 ρ=−0.6,拒绝候选,仍留在 x=0,半径缩成一。

重解同一位置的新子问题,得到 p=−1。这次预测下降为三,真实函数值 f(−1)=−2,实际下降为二,所以 ρ=2/3。候选被接受,下一点是 x1=−1,半径保持一。拒绝之后要重新解较小范围内的模型;这里恰好方向没变,是这个一维例子的特点。

下一轮的 g1=1,B1=8,无约束模型最优步 p=−1/8 位于球内。预测下降为 1/16,实际下降为 223/4096,故

ρ=223256≈0.8711.

它被接受,但因为没有碰到球面,半径仍不扩大。最终是否达到要求,要继续检查新点梯度;高比值只表示当前模型预测准确。

比值好,不等于已最优 ​

对 f(x)=x2/2 使用精确二次模型,任何非零下降步都有 ρ=1。若在 x=10 只走到 9.9,模型仍完全准确,却离最优点很远。因此终止量应包括 ‖∇f(x)‖,而不是只检查 ρ。

反过来,若 gk=0,Cauchy 下降下界退化为零。取 f(x)=−x2、x=0,梯度为零却是局部极大点。只检验一阶驻性的版本可以在此结束;要求排除负曲率时,还需内层检测 Bk 的负方向,并使用相应二阶停止标准。

推论与应用

缩小半径为什么能够修复失真 ​

设 f 在所有试探线段的邻域内有 L-Lipschitz 梯度,且 ‖Bk‖≤M。Taylor 的一阶余项给出

|f(xk+p)−mk(p)|≤12(L+M)‖p‖2.

若 ‖gk‖≥ε>0 且半径足够小,Cauchy 下界为 predk≥(κ/2)εΔk。于是

|1−ρk|≤(L+M)Δkκε.

误差是半径的二次量,保证下降是一阶量,所以半径充分小时,比值趋近一。只要梯度还显著非零,算法就不能永远以“模型不可信”为由拒绝所有步。

进一步假设 f 在迭代水平集上有下界,求值精确,模型矩阵一致有界,而且每轮均满足上述下降要求。则固定规则产生的无限序列满足

lim infk→∞‖gk‖=0.

证明用反证法:若尾部梯度始终至少为 ε,足够小的半径就不会再被缩小。因为每次只减半,所有尾部半径存在一个正下界。无限拒绝又会迫使半径减到必能接受的范围,所以成功步有无限多个。每个成功步通过 Cauchy 下界使 f 至少下降一个固定正数,最终违反下有界性。这个结论保证任意小的一阶残差会出现;它没有把非凸驻点升级成全局最优点。

工作量与实际停止 ​

一次试探需要模型求解以及一次真实目标求值;接受后通常重算梯度,拒绝时可复用旧点的梯度与模型矩阵。若显式组装并分解稠密 Bk,主要线性代数成本可达 O(n3);使用截断 CG 时,应报告 Hessian–向量积次数及目标、梯度求值次数,而不是只报外层轮数。

程序可在梯度达到目标容差时返回“一阶条件达到”,在达到预算、半径小到无法产生可分辨位移或内层失败时返回对应状态。小步长本身可能来自尺度失衡、舍入或坏模型,不构成驻性证书。若需要约束优化,必须把约束线性化、可行性恢复或罚函数机制加入外层;仅把普通 Newton 步放进一个球,并没有处理业务约束。

参考资料
  • Jorge Nocedal and Stephen J. Wright, Numerical Optimization, 2nd ed., 2006,§§4.1–4.2;作者目录。
  • Jorge J. Moré and D. C. Sorensen, Computing a Trust Region Step, 1981,§4 的信赖域外层与收敛分析。
  • David Bindel, Trust regions, Cornell CS 4220, 2026-04-20,“Adapting the trust region”;本页明确采用减半规则并逐项推导自造四次算例。
关系图谱14 个相邻概念 · 3 类关系

拖动节点调整位置。

显示关系

显示:依赖

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