形式陈述
设 C ⊆ R n 为非空闭凸集,f : C → R 为凸函数。在选定范数下取定义于 C 的 1 -强凸镜像函数 ψ ,并假定 ψ 在所有迭代锚点可微、下述子问题均取得解且解仍在其可微域;令 D ψ 为其Bregman 散度 公理库 Bregman 散度 Bregman divergence 凸函数值高于其一阶切平面的余量,用来度量与生成函数几何相适配的偏离。 。给定 x 0 ∈ ri C 、次梯度 g k ∈ ∂ f ( x k ) 与步长 η k > 0 ,镜像下降更新
x k + 1 = argmin x ∈ C { η k ⟨ g k , x ⟩ + D ψ ( x , x k ) } . 在无约束且 ∇ ψ 可逆时,一阶方程化为
∇ ψ ( x k + 1 ) = ∇ ψ ( x k ) − η k g k , 即先在对偶坐标做加法步,再映回原空间。约束情形的子问题最优性与 Bregman 三点恒等式给出,对任意 x ∈ C ,
η k ⟨ g k , x k + 1 − x ⟩ ≤ D ψ ( x , x k ) − D ψ ( x , x k + 1 ) − D ψ ( x k + 1 , x k ) . 这条不等式是镜像版的一阶最优性 公理库 一阶最优性条件 First-order optimality condition · Variational inequality optimality condition 以梯度和所有可行方向的非负内积充要刻画可微凸问题的全局极小点。 证书,也是后续望远镜求和的核心。
直觉
欧氏梯度默认各坐标用同一平方距离计价,但许多可行域并不服从这种几何。概率向量靠近边界时,加减一个固定坐标量容易产生负数;稀疏多面体的欧氏半径也可能随维数放大。镜像函数把原空间弯曲成合适的对偶坐标:线性化损失仍提供方向,而 Bregman 散度决定离旧点多远算“昂贵”。
更新式中散度的槽位不能交换。变量位于第一个槽 D ψ ( x , x k ) ,旧点提供切平面锚;一般 Bregman 散度不对称,交换会产生另一算法。强凸性使散度至少控制平方范数,从而把线性损失的对偶范数界转成位移和 regret 界。镜像下降不是把梯度名字换掉,而是在不改变一阶信息的前提下选择更匹配约束的几何。
例子与边界
在概率单纯形 Δ n 的相对内部取关于 ‖ ⋅ ‖ 1 为 1 -强凸的负熵
ψ ( x ) = ∑ i x i log x i , D ψ ( x , y ) = ∑ i x i log x i y i . 用 Lagrange 乘子处理 ∑ i x i = 1 ,更新得到乘法权重
x k + 1 , i = x k , i e − η k g k , i ∑ j x k , j e − η k g k , j . 取 x k = ( 1 / 2 , 1 / 2 ) 、g k = ( 1 , 0 ) 、η k = log 3 ,未归一化权重为 ( 1 / 6 , 1 / 2 ) ,归一化后恰为 ( 1 / 4 , 3 / 4 ) 。两坐标保持非负且和为一,无需事后欧氏投影。若某初始坐标为零,负熵梯度在边界不有限,乘法更新也永远无法恢复该坐标;标准理论因而从相对内部开始,或采用适当广义定义。
镜像函数若不关于目标范数强凸,三点散度不能控制实际位移,标准率会失去常数。次问题若不能精确求解,也需记录误差。对非凸目标,线性化加 Bregman 稳定项仍可运行,但本页的凸 regret/最优性推论不成立。步长过大可能让函数值上升;镜像下降的基础分析常控制加权平均或 regret,而不是保证最后一点单调。
推论与应用
若 ψ 关于 ‖ ⋅ ‖ 为 1 -强凸且 ‖ g k ‖ ∗ ≤ G ,将三点不等式与 Young 不等式结合并求和,得到任意 x ∈ C 的
∑ k = 0 T − 1 η k ⟨ g k , x k − x ⟩ ≤ D ψ ( x , x 0 ) + 1 2 ∑ k = 0 T − 1 η k 2 ‖ g k ‖ ∗ 2 . 若对所选比较器 x 有正数 D 满足 D ψ ( x , x 0 ) ≤ D ,固定 η = 2 D / ( G 2 T ) 时,固定凸目标的平均点误差为 O ( G D / T ) ;在线版本把同一估计解释为平均 regret。在 n 个坐标的单纯形上采用负熵并从 x 0 = ( 1 / n , … , 1 / n ) 均匀初始化时,才有对所有比较器成立的 D ψ ( x , x 0 ) ≤ log n ;任意内点初始化不享有这个统一常数。
若选 ψ ( x ) = 1 2 ‖ x ‖ 2 2 ,D ψ 化为平方欧氏距离,子问题完成平方后得到投影梯度法 公理库 投影梯度法 Projected gradient method · Gradient projection method 在每次梯度步后投影回闭凸可行集以保持可行性的约束一阶算法。 ;因此后者只是在欧氏镜像映射下的严格特殊情形。负熵选择则产生指数权重,适用于在线学习与概率单纯形。实际选择镜像函数时还要同时检查定义域、强凸常数、对偶范数和子问题成本,缺一项都不能仅凭漂亮几何宣称更快。
参考资料
Arkadi Nemirovski and David Yudin, Problem Complexity and Method Efficiency in Optimization , Wiley, 1983,Ch. 3,mirror-type methods and complexity。
Amir Beck and Marc Teboulle, “Mirror Descent and Nonlinear Projected Subgradient Methods for Convex Optimization,” Operations Research Letters 31(3), 2003, 167–175,Theorem 4.1。
Sébastien Bubeck, Convex Optimization: Algorithms and Complexity , Foundations and Trends in Machine Learning 8(3–4), 2015,§§4.2–4.3,mirror descent and the negative-entropy simplex setup。