Skip to content

算法Algorithm

近似近端梯度的误差预算

Inexact proximal gradient · Proximal subproblem residual budget

用可核验的近端子问题次梯度残差分配内层精度,并把累计误差传到外层平均目标差。

形式陈述 ​

考虑复合凸目标 F=g+h。g凸且梯度为 L-Lipschitz,L>0;h proper、闭凸,最小点 x∗ 存在。固定 0<α≤1/L。精确梯度已可用,但prox需要内层迭代,本页不要求把它解到无限精度。

第 k 轮从 xk 出发,内层返回 z=xk+1∈domh 和一个可验证的次梯度 sk+1∈∂h(z),并核验

(1)ek+1=∇g(xk)+z−xkα+sk+1,‖ek+1‖≤δk+1.

这里 δk+1是子问题最优性残差的上限,不是内层目标值差。输入还包括整数预算 T≥1、误差序列、初始距离有效上界 D0≥‖x0−x∗‖和已知距离上界 R,要求 ‖xk+1−x∗‖≤R。例如 h的定义域是直径不超过 R 的有界闭凸集,这个条件自动成立;无约束问题则必须另证有界性,不能凭观察到的有限轨道宣称成立。

返回平均点 x¯T=T−1∑k=0T−1xk+1、每轮残差证据、实际内层工作量和完成状态。本文证明

(2)F(x¯T)−F∗≤‖x0−x∗‖22αT+RT∑k=1Tδk.

某轮内层未能提供式(1)的证据,应报告预算不足;一个看起来稳定的点不能代替承诺的残差。实际报告式(2)时,以已知 D02替代未知初始距离平方。

直觉

精确prox会让模型的梯度与惩罚次梯度完全平衡。式(1)允许有一小段“未平衡力”,但把它的大小记录下来。外层每轮为这段力支付至多 Rδk,因此误差不是消失了,而是进入最终账本。

令

Qk(u)=g(xk)+⟨∇g(xk),u−xk⟩+‖u−xk‖22α+h(u).

式(1)正是给出了 ek+1∈∂Qk(z)。强凸子问题还有两个可用推论:若精确解为 pk,则 ‖z−pk‖≤α‖ek+1‖;且 Qk(z)−Qk(pk)≤α‖ek+1‖2/2。第一条由强单调性与Cauchy–Schwarz得到;第二条由强凸下界对比较点最小化得到。因此一个小残差确实控制内层点和内层值,但反方向不能随意使用。

例子与边界

固定容差造成真实偏移 ​

取 F(x)=(x−2)2/2+|x|,最优点 x∗=1,最优值 3/2。令 α=1/2、x0=1,内层总返回正点并取 s=1,使残差恒为 e=1/10。式(1)变成

2xk+1−xk−1=1/10,xk+1=xk/2+11/20.

于是 x1=21/20、x2=43/40,一般地 xk=11/10−(1/10)2−k,趋于 11/10。由于正半轴上 F(x)−F∗=(x−1)2/2,极限目标误差为 1/200。每次内层都满足同一个小容差,仍会留下不消失的外层误差。

若已知初始距离不超过1、R=2、α=1/2,要求式(2)的右侧不超过 0.02,可取 T=100,并选 δk=1/(4k2)。用 ∑k=1∞k−2≤1+∫1∞t−2dt=2,得上界 1/100+1/100=0.02。这是一份可检查的外层预算;是否能经济地实现后期很小的残差,仍要看内层求解器。

很小的目标误差不保证很小的次梯度残差 ​

取 Q(u)=|u|+u2/2,精确最小点为零。令 z=ε>0,则 Q(z)−Q(0)=ε+ε2/2→0;但 ∂Q(z)={1+ε},最小次梯度范数反而趋于1。靠近尖点的目标值可以非常准确,而普通次梯度残差依然大。

若只有 Qk(z)−minQk≤ϵk,强凸性只能直接给 ‖z−pk‖≤2αϵk。要把这个误差传入外层,应使用相应的点扰动或 ϵ-次微分分析,不能把 ϵk直接当作式(1)的 δk。具体说,s∈∂ϵh(z)表示对全部 u有 h(u)≥h(z)+⟨s,u−z⟩−ϵ。若内层给出这一证据,并在式(1)中使用该 s,则下面一步式右侧额外加 ϵ,平均界额外加各轮 ϵ的平均。仅有 Qk的目标误差不等于已经给出了这份关于 h的证据。残差、目标误差和 ϵ-次梯度是三种不同接口。

推论与应用

一步不等式与完整累计 ​

由式(1),有 s=e−∇g(x)+(x−z)/α∈∂h(z)。任取 u∈domh,用其次梯度不等式、g的凸下界与光滑上界相加,得到

(3)F(z)−F(u)≤‖x−u‖2−‖z−u‖22α−(12α−L2)‖z−x‖2+⟨e,z−u⟩.

这里误差内积的正号来自 s=e−∇g(x)+(x−z)/α;若把残差定义改成相反方向,必须同步改符号。取 u=x∗,因 α≤1/L 可丢掉非正位移项,再以距离上界估计最后一项不超过 Rδk+1。

从 k=0 到 T−1 求和,平方距离望远镜消去;去掉最后的非负距离,并由Jensen不等式将平均函数值换成平均点函数值,便得到式(2)。证明没有使用目标单调,因此不能把平均点保证自动替换为最后一点保证。

若 ∑kδk<∞,式(2)保持 O(1/T);若 δk=d/k,得到 O((1+log⁡T)/T);若固定 δk=d,该账本只给额外的 Rd。这里只讨论不加速版本,FISTA的增长权重会改变误差累计要求,不能把同一容差表直接搬过去。

内外层费用与误差位置 ​

若第 k 轮内层为达到 δk用了 Nk 次工作,总费用要报告 T次光滑梯度、∑kNk 次内层工作以及次梯度证据的计算成本。若模型本身的梯度也近似,那个误差需一同计入式(1),并说明是否可确定地上界。

ADMM的两份近端子问题也可能采用内层求解,但其外层势与本页不同。本文的 R∑δk结论只属于这里的近端梯度协议,不是给所有拆分算法通用的许可。旧近端点方法已经给出prox点误差可求和时的点收敛条件,并展示近似点离解很近却处于有效域之外的反例;本页新增的是复合子问题残差、有效域内候选和有限预算目标平均界。

两道短自测 ​

  1. 在式(2)中取初始距离上界1、α=1/2,R=2,T=100、各轮 δk=0.01,可报告多大界?答案:0.01+2⋅0.01=0.03,增加 T只消除第一项。
  2. 内层给 Q(z)−Q∗≤10−8,是否能直接填 δ=10−8?答案:不能;这是目标误差,尖点例说明普通次梯度残差仍可接近1。必须取得式(1)证据,或换用与目标容差匹配的分析。
参考资料
关系图谱14 个相邻概念 · 2 类关系

拖动节点调整位置。

显示关系

显示:依赖

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