Skip to content

近端梯度法

Proximal gradient method · Forward-backward splitting

对复合目标的光滑项取显式梯度步、对非光滑凸项取隐式近端步的算法。

条目类型
算法

形式陈述

考虑复合问题

minxRnF(x)=g(x)+h(x),

其中 g 可微,h 为 proper、下半连续凸函数且其 prox 可计算。给定 x0、步长 αk>0、容差与预算,近端梯度法迭代

xk+1=proxαkh(xkαkg(xk)).

第一部分是梯度下降的显式前向步,第二部分是近端算子的隐式后向步,因此也称 forward-backward splitting。算法定义只规定更新与输出状态;若要保证凸全局收敛,另需 g$L$-光滑凸函数h 凸、解存在,并选择例如 0<αk1/L

定义近端梯度映射

Gα(x)=1α(xproxαh(xαg(x))).

更新即 x+=xαGα(x)。由 prox 的最优性条件,Gα(x)=0 当且仅当

0g(x)+h(x),

这正是复合凸问题的一阶最优性条件。停止时应报告 Gα(x),不能在 h 不可微时误报不存在的 F(x)

直觉

若直接对 g+h 做梯度步,h 的尖角可能没有梯度;若每轮对整个 F 做 prox,又可能把本来容易求梯度的 g 放进昂贵子问题。近端梯度恰好按计算结构拆分:先用 g 的局部线性信息提出候选,再让 h 的 prox 在降低惩罚与靠近候选之间精确权衡。算法利用的是“一个梯度便宜、另一个 prox 便宜”,而不只是代数上能写成两项。

h=0 时,prox 是恒等映射,算法退化为梯度下降;在 h 是集合指标函数时,prox 变成投影。这个统一视角也解释了软阈值:1 惩罚不会在显式梯度中产生精确零,但其 prox 能把一整段输入直接送到零。步长仍由光滑项 g 的曲率控制,而不是由整个非光滑目标的某个虚构 Hessian 控制。

例子与边界

取一维复合目标

F(x)=12(x3)2+|x|.

g(x)=(x3)2/2L=1h(x)=|x| 的 prox 是软阈值。用 α=1/2x0=0,有

x1=S1/2(012(3))=S1/2(1.5)=1,x2=S1/2(112(2))=S1/2(2)=1.5,x3=S1/2(2.25)=1.75,

其中 Sτ(z)=sign(z)(|z|τ)+。函数值依次为 4.5,3,2.625,2.53125,趋向唯一解 x=2F(x)=2.5。若取紧步长 α=1,无论初值为何,梯度候选都是 3,一次软阈值即到 2;这是一维特殊结构,不代表一般问题一步收敛。

边界来自拆分而非记号。若 h 的 prox 只能由昂贵内层迭代近似,必须控制内层误差;若 g 非凸,固定点通常只表示复合驻点而非全局极小。低估 L 可能破坏充分下降,可用回溯实际检查由下降引理给出的二次模型。把两个都非光滑的项硬指定一个为“光滑项”,或把约束集合非凸时的多值投影当单值 prox,均超出标准理论。

推论与应用

在凸假设与 α1/L 下,更新的最优性不等式和光滑二次上界合并可得

F(xk+1)F(x)xkx2xk+1x22α.

求和后得到 F(xk)F=O(1/k);若复合目标强凸,可进一步得到几何率。加上 Nesterov 权重便形成 FISTA 的 O(1/k2) 率,但主序列可能不单调,证明需要新的势函数。这里控制的是复合函数值或梯度映射,不应偷换成每个坐标恢复正确支持集。

Lasso、稀疏逻辑回归、总变差模型和带简单凸约束的经验风险最小化都可使用这一结构。Moreau 包络也通过 xprox(x) 产生光滑梯度,但它对整个函数取 infimal convolution;近端梯度只对 h 取 prox。二者公式相似而优化对象不同,诊断时应先写清子问题。

参考资料
  • Amir Beck, First-Order Methods in Optimization, SIAM, 2017,Ch. 10,proximal gradient methods。
  • Neal Parikh and Stephen Boyd, Proximal Algorithms, Foundations and Trends in Optimization 1(3), 2014,§4.2,proximal gradient and splitting interpretation。
  • Amir Beck and Marc Teboulle, “A Fast Iterative Shrinkage-Thresholding Algorithm for Linear Inverse Problems,” SIAM Journal on Imaging Sciences 2(1), 2009, 183–202,§§2–4,ISTA and FISTA rates。
关系图谱11 个相邻概念 · 3 类关系

拖动节点调整位置。

显示关系

显示:依赖

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

上位 / 更一般

暂未标注直接上位概念。

下位 / 直接特例

类型化关系