形式陈述
设 D ⊆ R n 为凸集,U ⊆ D 为开凸集,ψ : D → R 为凸函数并在 U 可微。约束集 C ⊆ D 非空、闭且凸,目标 f : C → R 凸,并且 ψ | C 关于选定范数为 1 -强凸。使用Bregman 散度的边界扩展 理路 Bregman 散度 Bregman divergence 凸函数值高于其一阶切平面的余量,用来度量与生成函数几何相适配的偏离。 :第一参数可在 D 中,作切平面的第二参数必须在 U 中。
假定下述子问题均取得解,且解仍在 C ∩ U 。给定 x 0 ∈ C ∩ U 、有限次梯度 理路 次梯度与次微分 Subgradient · Subdifferential 以全局仿射下界刻画凸函数在不可微点的支撑斜率集合。 g k ∈ ∂ f ( x k ) 与步长 η k > 0 ,镜像下降更新如下;这里以 f 在 C 外取 + ∞ 的扩展解释其次微分。
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 散度决定离旧点多远算“昂贵”。
图片加载失败 全空间情形的镜像更新 图示取 C = D = U = R n ,且 ψ 可微并强凸;此时凸共轭的梯度满足 ∇ ψ ∗ = ( ∇ ψ ) − 1 。一般约束集上的更新仍须求解前面的 Bregman 极小化子问题,不能只用这条全空间逆映射代替。
更新式中散度的槽位不能交换。变量位于第一个槽 D ψ ( x , x k ) ,旧点提供切平面锚;一般 Bregman 散度不对称,交换会产生另一算法。强凸性使散度至少控制平方范数,从而把线性损失的对偶范数界转成位移和 regret 界。镜像下降不是把梯度名字换掉,而是在不改变一阶信息的前提下选择更匹配约束的几何。
例子与边界
本例的 log 均为自然对数。取 D = [ 0 , ∞ ) n 、U = ( 0 , ∞ ) n ,负熵在零坐标处按 0 log 0 = 0 延拓。约束集为闭概率单纯形 C = Δ n ,其与 U 的交正是相对内部;负熵在 C 上关于 ‖ ⋅ ‖ 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 ,取正整数 T ,将三点不等式与 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 有正数 B 满足 D ψ ( x , x 0 ) ≤ B ,取正的有效上界 G > 0 ,固定 η = 2 B / ( G 2 T ) 。对固定凸目标,令 x ¯ T = T − 1 ∑ k = 0 T − 1 x k ,凸性与上述估计给出目标值差
f ( x ¯ T ) − f ( x ) ≤ 1 T ∑ k = 0 T − 1 ⟨ g k , x k − x ⟩ ≤ B η T + η G 2 2 = G 2 B T . 这控制平均点的目标值,不直接控制 ‖ x ¯ T − x ‖ 。在线版本若 g k ∈ ∂ f k ( x k ) ,同一估计控制 T − 1 ∑ k ( f k ( x k ) − f k ( x ) ) ,即平均 regret。在 n 个坐标的单纯形上采用负熵并从 x 0 = ( 1 / n , … , 1 / n ) 均匀初始化时,才有对所有比较器成立的 D ψ ( x , x 0 ) ≤ log n ;任意内点初始化不享有这个统一常数。
若选 ψ ( x ) = 1 2 ‖ x ‖ 2 2 ,D ψ 化为平方欧氏距离,在 f 可微且选 g k = ∇ f ( x k ) 时,子问题完成平方后得到投影梯度法 理路 投影梯度法 Projected gradient method · Gradient projection method 在每次梯度步后投影回闭凸可行集以保持可行性的约束一阶算法。 ;对不可微凸 f 则是投影次梯度更新;因此后者只是在欧氏镜像映射下的严格特殊情形。欧氏情形可取 D = U = R n ,所以迭代点位于 C 的边界也不妨碍镜像函数可微。例如 C = [ 0 , ∞ ) 、f ( x ) = 2 x 、x 0 = 1 、η 0 = 1 时,平方欧氏镜像更新得到 x 1 = 0 ;它在约束边界,却仍在 U 中。负熵选择则产生指数权重,适用于在线学习与概率单纯形。实际选择镜像函数时还要同时检查定义域、强凸常数、对偶范数和子问题成本,缺一项都不能仅凭漂亮几何宣称更快。
随机复合镜像下降 理路 随机复合镜像下降 Stochastic composite mirror descent · 随机复合 Bregman 更新 在非欧氏几何中分开随机次梯度和精确结构项,证明更新前平均的有限预算界,并计算熵正则单纯形的真实更新。 保留这里的几何,再把条件无偏次梯度与新点结构项放在不同槽位。它完整保留结构项首尾差,证明更新前平均的期望界,并计算目标熵正则与镜像熵同时出现时的幂次更新;只有镜像几何并不表示目标已被正则化。
参考资料
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。