Skip to content

方法Method

近端点方法

Proximal point method · Proximal point algorithm

反复精确求解带二次稳定项的子问题以逼近闭凸函数最小点的方法。

形式陈述 ​

考虑在 Hilbert 空间上最小化 proper、下半连续凸函数 f(允许取 +∞,但不恒为 +∞,也不取 −∞),并假设 argminf≠∅。给定 x0 与步长 λk>0,经典近端点方法迭代

xk+1=proxλkf(xk)=argminx{f(x)+12λk‖x−xk‖2}.

由近端算子的最优性条件,

xk−xk+1λk∈∂f(xk+1),

所以它是隐式次梯度步。任意最小点 x∗ 都是 prox 的固定点;反过来,固定点满足 0∈∂f(x∗),因而是全局最小点。若步长有统一下界 λk≥λ―>0,则精确迭代在有限维中收敛到某个最小点;一般 Hilbert 空间中标准结论是弱收敛。若每步仅近似求解,一个常用充分条件是绝对误差可求和,即 ‖xk+1−proxλkf(xk)‖≤εk 且 ∑kεk<∞;在同一步长下界条件下,仍有相应的点收敛结论。

精确迭代的收敛论证 ​

取任意最小点 z,近端牢固非扩张性给出

‖xk+1−z‖2+‖xk+1−xk‖2≤‖xk−z‖2.

因此轨道有界,且相邻位移趋于零。步长下界使 gk+1=(xk−xk+1)/λk→0;由次梯度不等式,

0≤f(xk+1)−f(z)≤⟨gk+1,xk+1−z⟩⟶0.

有限维中,用Bolzano–Weierstrass 定理抽取收敛子列。下半连续性使其极限仍为最小点;到这个点的距离单调不增,并有子列趋于零,所以整列范数收敛。

一般 Hilbert 空间中改用有界序列的弱子列。闭凸子水平集也是弱闭的,所以 f 弱下半连续,每个弱聚点都是最小点。若 u,v 是两个弱聚点,到它们的距离都存在极限,而恒等式

‖xk−u‖2−‖xk−v‖2=2⟨xk,v−u⟩+‖u‖2−‖v‖2

沿两条弱子列取极限后,迫使 u=v。任意子列还可抽出弱收敛子列,因此唯一弱聚点就是整列的弱极限。这证明了所述精确迭代结论;带可求和误差的版本另用扰动估计,不能直接把上面的精确等式原样套入。

直觉

每一步最小化“原目标加距离惩罚”。原目标鼓励向更低的函数值移动,二次项则提高远离当前点的代价;步长 λk 越大,距离惩罚越弱。

最优性条件在新点 xk+1 计算次梯度,所以更新是隐式的。若 f(x)=|x|,子问题可直接用软阈值求解;若变量之间通过复杂损失耦合,计算一次 prox 本身就是一个优化问题。因此算法的每轮成本取决于近端子问题的结构。

在显式存储的 n 维实现中,每个外层步骤调用一次 prox,另需 O(n) 的向量操作与状态存储;prox 自身的时间、工作空间以及近似求解所需的内层精度工作另计。下述各向同性二次函数的缩放更新,以及逐坐标绝对值之和对应的软阈值更新,都可用 O(n) 次算术操作计算。

例子与边界

对 f(x)=μ2‖x‖2(μ>0),近端子问题可直接求导,得到

xk+1=xk1+λkμ.

固定步长 λ 下,‖xk‖=(1+λμ)−k‖x0‖,序列以几何速度收缩到唯一最小点 0。这个计算也说明它不同于梯度步 xk+1=(1−λμ)xk:前者的分母形式对任意正步长都保持收缩,后者则受稳定步长范围限制。

同一个目标也能看出两种近端算法的差别。在耦合 Lasso 算例中,取 f=P、x0=0、λ0=1/3,整个目标的近端点求解 minx{P(x)+(3/2)‖x‖2}。猜测两坐标均为正,最优性方程为

(5115)x=(21/2),x=(19/48,1/48).

两坐标确实为正,而且 AT(b−Ax)−(1,1)=(19/16,1/16)=3x,故满足完整的凸最优性条件。这不同于只对 ℓ1 取 prox 的近端梯度第一步 (2/3,1/6);两者处理同一目标,却解不同的子问题。

步长的累积决定了二次例子能推进多远。若 λk=2−k−1,则 ∑kλk=1,而

|xk|=|x0|∏j<k(1+μλj)≥e−μ|x0|.

当 x0≠0 时,迭代始终与最小点保持正距离。统一正下界排除了这种总步长有限的情形;最小点存在则为固定点论证提供目标。

可求和的近端点误差控制的是迭代点,不自动保证近似点的函数值也收敛。例如令 f=δ{0}、λk=1,每次精确近端点都是零;却可取近似序列 xk=2−k,其逐步近端误差 2−k−1 可求和。序列确实趋于最小点零,但每个 f(xk)=+∞。因此要对非精确迭代报告目标值,还需另外保证候选在有效域内并控制函数值误差。

推论与应用

有限维时,近端点方法把凸优化问题写成近端算子的固定点迭代;其形式陈述还允许一般 Hilbert 空间。对精确迭代,近端算子的牢固非扩张性给出

‖xk+1−x∗‖2≤‖xk−x∗‖2−‖xk+1−xk‖2.

每一步都使到任意最小点的距离不增,并控制相邻位移的平方和。若 f−μ‖⋅‖2/2 凸(μ>0),两条带二次余量的次梯度不等式相加,得到 ⟨u−v,p−q⟩≥μ‖p−q‖2。代入近端关系 u=(x−p)/λ、v=(y−q)/λ,再用Cauchy–Schwarz 不等式,可得 ‖p−q‖≤‖x−y‖/(1+λμ)。因此固定步长 prox 的收缩因子为 1/(1+λμ),Banach 不动点定理进一步给出唯一最小点与几何收敛。

对固定步长的有限维情形,也可直接应用单调算子页的预解迭代收敛定理,取 A=∂f。Rockafellar 将这一迭代推广到极大单调算子的 resolvent,从而覆盖变分不等式与鞍点系统。该条目的二维旋转算例完整展示了一个不来自凸次微分的隐式迭代:显式步使长度增长,预解算子却使轨道旋转并收缩到零点。Douglas–Rachford、ADMM 等分裂方法则把相关问题组织成若干较易计算的子步骤,利用各部分的结构降低一次完整近端求解的成本。交替方向乘子法具体推导了约束 x=z 下的两次 prox 与乘子更新,并用双残差及 Fejér 不等式证明有限维点收敛;其中的耦合 Lasso 例子可以与本页对整个目标求 prox 的第一步直接比较。

复合凸问题的原对偶系统进一步把次微分与斜对称耦合放入乘积空间,在同一个新点求解两行隐式关系。它给出每步可解和有限维整列收敛的完整论证,并把一个具体预解算子消元为区间投影,得到 (xk,yk)=(1−2−k,1)(k≥1)的轨道。那是整个原对偶系统的 resolvent,其计算成本与对单个函数求 prox 仍需分别判断。

参考资料
  • R. Tyrrell Rockafellar, “Monotone Operators and the Proximal Point Algorithm,” SIAM Journal on Control and Optimization 14(5), 1976, 877–898。
  • Neal Parikh and Stephen Boyd, Proximal Algorithms, Foundations and Trends in Optimization 1(3), 2014,§4.1, proximal minimization。
关系图谱11 个相邻概念 · 2 类关系

拖动节点调整位置。

显示关系

显示:依赖

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