形式陈述
设 C ⊆ R n 为开凸集,f : C → R 可微。若 f 是凸函数 公理库 凸函数 Convex function 函数在任意凸组合处不超过相同权重下函数值的凸组合。 ,且存在 L > 0 使其梯度满足 Lipschitz 条件 公理库 Lipschitz 梯度 Lipschitz gradient · Lipschitz continuous gradient 梯度映射的变化量由点间距离乘统一常数控制的正则性条件。
‖ ∇ f ( x ) − ∇ f ( y ) ‖ 2 ≤ L ‖ x − y ‖ 2 , x , y ∈ C , 则称 f 为 L -光滑凸函数,常记 f ∈ F L 。这两个条件分别给出一阶模型的下界与上界:
0 ≤ f ( y ) − f ( x ) − ⟨ ∇ f ( x ) , y − x ⟩ ≤ L 2 ‖ y − x ‖ 2 2 . 在 C = R n 且 f 二次连续可微时,等价的逐点 Hessian 条件是
0 ⪯ ∇ 2 f ( x ) ⪯ L I . 对任意 x , y ,凸性还把光滑性加强为梯度的余单调性:
⟨ ∇ f ( x ) − ∇ f ( y ) , x − y ⟩ ≥ 1 L ‖ ∇ f ( x ) − ∇ f ( y ) ‖ 2 2 . 这里的 L 是一个有效上界,不必是最小常数;最小者刻画该函数在所选范数下的最坏一阶曲率。
直觉
凸性说切平面永远位于图像下方,L -光滑性则说图像离开切平面的速度至多像曲率为 L 的抛物面。二者合起来,把未知函数在当前点附近夹在“支撑平面”和“支撑平面加二次帽”之间。因此只知道 f ( x ) 与 ∇ f ( x ) ,也能预先控制走一步后的最坏函数值,而不必猜测整张图像。
“光滑”在这里不是“有任意多阶导数”的口语。它要求梯度在整个讨论域上有统一变化上界;C ∞ 函数也可能没有全局有限的 L 。反过来,梯度 Lipschitz 只控制曲率上界,并不强迫曲率为正;凸性是另一个独立条件。缩放目标也会缩放常数:若 f ∈ F L 且 a > 0 ,则 a f ∈ F a L ,所以步长不能脱离目标尺度讨论。
余单调不等式还有一个算法直觉:梯度差不会主要横着摆动,而是必须在位移方向上支付足够内积。这正是映射 I − ∇ f / L 非扩张的来源。它比单纯的 ‖ ∇ f ( x ) − ∇ f ( y ) ‖ ≤ L ‖ x − y ‖ 多用了凸性,不能移植到一般非凸光滑函数。
例子与边界
二次函数
f ( x ) = 1 2 x ⊤ Q x + b ⊤ x + c , Q = Q ⊤ ⪰ 0 , 满足 ∇ f ( x ) = Q x + b ,最小光滑常数是 L = λ max ( Q ) 。取 Q = diag ( 1 , 4 ) 、b = ( − 2 , 0 ) ⊤ ,则 L = 4 ,唯一极小点为 ( 2 , 0 ) 。沿 e 2 方向有
‖ ∇ f ( x + t e 2 ) − ∇ f ( x ) ‖ 2 ‖ t e 2 ‖ 2 = 4 , 所以常数 4 不是松估计;沿 e 1 的比值仅为 1 。在 x = ( 0 , 1 ) 处,切平面余量对增量 h = ( 2 , − 1 ) 恰为 1 2 h ⊤ Q h = 4 ,而二次上界给 L 2 ‖ h ‖ 2 = 10 ,可直接核对夹逼。
边界要分三种。f ( x ) = x 4 凸且无限可微,但 f ″ ( x ) = 12 x 2 无全局上界,因此在 R 上不属于任何 F L ;限制到 [ − R , R ] 后才可取 L = 12 R 2 。f ( x ) = | x | 凸却在零点不可微,也不属于本页函数类。f ( x ) = sin x 的导数是 1 -Lipschitz,却不是凸函数,说明 Lipschitz 梯度本身不能推出全局最优结构。若 L 被低估,基于二次上界选择的固定步长就可能失去下降保证;若高估很多,结论仍正确但迭代会不必要地保守。
推论与应用
对 L -光滑凸函数,取梯度步 x + = x − ∇ f ( x ) / L ,二次上界立刻给
f ( x + ) ≤ f ( x ) − 1 2 L ‖ ∇ f ( x ) ‖ 2 2 . 这条可观测下降把函数值差、梯度范数和迭代复杂度连接起来,是后续梯度下降收敛定理 公理库 梯度下降收敛定理 Gradient descent convergence theorem · Convergence rate of gradient descent 分别给出光滑凸与光滑强凸目标上梯度下降的次线性和几何函数值收敛率。 的基本接口。若再有强凸下界,曲率被 μ I 与 L I 双向夹住,条件数 κ = L / μ 才出现;只有光滑凸性时,极小点可以不唯一,不能宣称迭代点按几何速度收敛。
余单调性还说明 ∇ f / L 是 firmly nonexpansive,从而可用固定点几何分析梯度映射。复合优化中,光滑部分适合显式梯度步,非光滑但 prox 易算的部分适合隐式近端步;加速方法则利用同一个二次上模型,却以额外序列组织历史信息。所有这些算法都需要分别说明输入、步长和停止准则,函数类定义本身不包含某一种实现。
参考资料
Yurii Nesterov, Introductory Lectures on Convex Optimization: A Basic Course , Kluwer, 2004,§§2.1.2–2.1.3,smooth convex functions and gradient inequalities。
Sébastien Bubeck, Convex Optimization: Algorithms and Complexity , Foundations and Trends in Machine Learning 8(3–4), 2015,§3.2,smooth convex optimization。
Heinz H. Bauschke and Patrick L. Combettes, Convex Analysis and Monotone Operator Theory in Hilbert Spaces , 2nd ed., Springer, 2017,Corollary 18.17,Baillon–Haddad theorem。