形式陈述
设 f : R n → R ∪ { + ∞ } 为 proper、下半连续凸函数,λ > 0 。它的 Moreau 包络定义为
e λ f ( x ) = inf y { f ( y ) + 1 2 λ ‖ y − x ‖ 2 2 } . 二次项使关于 y 的子问题强凸,故极小点唯一,并等于近端算子 公理库 近端算子 Proximal operator · Proximity operator 在降低凸函数值与保持靠近输入点之间取得精确平衡的单值算子。 p ( x ) = prox λ f ( x ) 。包络处处有限且连续可微,梯度为
∇ e λ f ( x ) = 1 λ ( x − p ( x ) ) . 由于 prox firmly nonexpansive,I − p 也 firmly nonexpansive,因此包络具有常数 1 / λ 的Lipschitz 梯度 公理库 Lipschitz 梯度 Lipschitz gradient · Lipschitz continuous gradient 梯度映射的变化量由点间距离乘统一常数控制的正则性条件。 。此外
e λ f ( x ) ≤ f ( x ) , inf e λ f = inf f , argmin e λ f = argmin f . 当 λ ↓ 0 时,二次惩罚迫使 p ( x ) 靠近 x ,包络从下方点态逼近 f ;λ 越大则平滑越强、下偏越明显。
直觉
在每个查询点 x ,包络不直接读取 f ( x ) ,而是在所有候选 y 中寻找“函数值低且离 x 不太远”的最佳折中。最优候选 p ( x ) 给包络值,位移 x − p ( x ) 则给斜率。即使原函数有尖角,二次惩罚使最优候选随 x 稳定移动,于是尖角被圆滑成可微过渡区。
λ 同时控制两件相反的事:小 λ 保真但梯度常数 1 / λ 大,数值步长更受限;大 λ 更平滑却把函数压得更低。包络并不是对数据做局部平均,而是一次精确的最小化卷积,所以它能保留全局极小点。计算梯度也并非免费:公式要求能够计算 prox;若近端子问题本身困难,写出光滑包络并没有消除计算成本。
例子与边界
取一维 f ( y ) = | y | 。近端点是软阈值
p ( x ) = sign ( x ) max { | x | − λ , 0 } , 代回可得 Huber 型包络
e λ f ( x ) = { x 2 2 λ , | x | ≤ λ , | x | − λ 2 , | x | > λ . 其导数在内区间为 x / λ ,外区间为 sign ( x ) ,在 x = ± λ 连续;斜率变化的最大比值正是 1 / λ 。例如 λ = 2 , x = 3 时,p ( x ) = 1 ,包络值 1 + ( 3 − 1 ) 2 / 4 = 2 ,也等于 3 − 1 ,可直接复算。
若 f = δ C 是非空闭凸集的指标函数,则 e λ f ( x ) = dist ( x , C ) 2 / ( 2 λ ) ,把硬约束变成平方距离。凸性不可随意删除:非凸 f 的近端点可能多值,包络也可能不可微。即使在凸情形,e λ f ( x ) ≤ f ( x ) 不表示它在每一点都近似良好;靠近陡峭跳变或取大 λ 时偏差可以显著。包络保留极小点,但不保留所有函数值、Hessian 或统计解释。
推论与应用
梯度公式把不可微最优性变成光滑残差:
∇ e λ f ( x ) = 0 ⟺ x = prox λ f ( x ) ⟺ 0 ∈ ∂ f ( x ) . 因此可用 ‖ x − p ( x ) ‖ / λ 衡量离凸极小点的程度。对集合指标函数,这就是点到其投影的缩放向量;对 | x | ,它在尖角附近线性连续地替代不连续的符号次梯度。Moreau 分解还把 f 与其凸共轭的 prox 连接起来,形成原—对偶残差分解。
需要区分包络梯度下降与近端梯度法 公理库 近端梯度法 Proximal gradient method · Forward-backward splitting 对复合目标的光滑项取显式梯度步、对非光滑凸项取隐式近端步的算法。 。前者对整个 f 计算 prox 后沿 e λ f 的梯度移动;后者针对 g + h ,只对 h 取 prox、对光滑的 g 取显式梯度。二者都出现 x − prox ( x ) ,但子问题、固定点和收敛常数不同。Moreau 平滑还用于单调算子的 Yosida 逼近和增广 Lagrangian 分析,不过必须把 prox 求解误差纳入最终证书。
参考资料
Jean-Jacques Moreau, “Proximité et dualité dans un espace hilbertien,” Bulletin de la Société Mathématique de France 93, 1965, 273–299。
Heinz H. Bauschke and Patrick L. Combettes, Convex Analysis and Monotone Operator Theory in Hilbert Spaces , 2nd ed., Springer, 2017,Ch. 12 and Proposition 12.30,Moreau envelopes and gradients。
Neal Parikh and Stephen Boyd, Proximal Algorithms , Foundations and Trends in Optimization 1(3), 2014,§§3.1–3.2,Moreau envelope and decomposition。