Skip to content

算法Algorithm

镜像下降法

Mirror descent method · Bregman gradient method

用强凸镜像映射生成的 Bregman 几何执行线性化损失更新的约束一阶算法。

形式陈述 ​

设 D⊆Rn 为凸集,U⊆D 为开凸集,ψ:D→R 为凸函数并在 U 可微。约束集 C⊆D 非空、闭且凸,目标 f:C→R 凸,并且 ψ|C 关于选定范数为 1-强凸。使用Bregman 散度的边界扩展:第一参数可在 D 中,作切平面的第二参数必须在 U 中。

假定下述子问题均取得解,且解仍在 C∩U。给定 x0∈C∩U、有限次梯度 gk∈∂f(xk) 与步长 ηk>0,镜像下降更新如下;这里以 f 在 C 外取 +∞ 的扩展解释其次微分。

xk+1=argminx∈C{ηk⟨gk,x⟩+Dψ(x,xk)}.

在无约束且 ∇ψ 可逆时,一阶方程化为

∇ψ(xk+1)=∇ψ(xk)−ηkgk,

即先在对偶坐标做加法步,再映回原空间。约束情形的子问题最优性与 Bregman 三点恒等式给出,对任意 x∈C,

ηk⟨gk,xk+1−x⟩≤Dψ(x,xk)−Dψ(x,xk+1)−Dψ(xk+1,xk).

这条不等式是镜像版的一阶最优性证书,也是后续望远镜求和的核心。

直觉

欧氏梯度默认各坐标用同一平方距离计价,但许多可行域并不服从这种几何。概率向量靠近边界时,加减一个固定坐标量容易产生负数;稀疏多面体的欧氏半径也可能随维数放大。镜像函数把原空间弯曲成合适的对偶坐标:线性化损失仍提供方向,而 Bregman 散度决定离旧点多远算“昂贵”。

全空间情形的镜像更新

图示取 C=D=U=Rn,且 ψ 可微并强凸;此时凸共轭的梯度满足 ∇ψ∗=(∇ψ)−1。一般约束集上的更新仍须求解前面的 Bregman 极小化子问题,不能只用这条全空间逆映射代替。

更新式中散度的槽位不能交换。变量位于第一个槽 Dψ(x,xk),旧点提供切平面锚;一般 Bregman 散度不对称,交换会产生另一算法。强凸性使散度至少控制平方范数,从而把线性损失的对偶范数界转成位移和 regret 界。镜像下降不是把梯度名字换掉,而是在不改变一阶信息的前提下选择更匹配约束的几何。

例子与边界

本例的 log 均为自然对数。取 D=[0,∞)n、U=(0,∞)n,负熵在零坐标处按 0log⁡0=0 延拓。约束集为闭概率单纯形 C=Δn,其与 U 的交正是相对内部;负熵在 C 上关于 ‖⋅‖1 为 1-强凸

ψ(x)=∑ixilog⁡xi,Dψ(x,y)=∑ixilog⁡xiyi.

用 Lagrange 乘子处理 ∑ixi=1,更新得到乘法权重

xk+1,i=xk,ie−ηkgk,i∑jxk,je−ηkgk,j.

取 xk=(1/2,1/2)、gk=(1,0)、ηk=log⁡3,未归一化权重为 (1/6,1/2),归一化后恰为 (1/4,3/4)。两坐标保持非负且和为一,无需事后欧氏投影。若某初始坐标为零,负熵梯度在边界不有限,乘法更新也永远无法恢复该坐标;标准理论因而从相对内部开始,或采用适当广义定义。

镜像函数若不关于目标范数强凸,三点散度不能控制实际位移,标准率会失去常数。次问题若不能精确求解,也需记录误差。对非凸目标,线性化加 Bregman 稳定项仍可运行,但本页的凸 regret/最优性推论不成立。步长过大可能让函数值上升;镜像下降的基础分析常控制加权平均或 regret,而不是保证最后一点单调。

推论与应用

若 ψ 关于 ‖⋅‖ 为 1-强凸且 ‖gk‖∗≤G,取正整数 T,将三点不等式与 Young 不等式结合并求和,得到任意 x∈C 的

∑k=0T−1ηk⟨gk,xk−x⟩≤Dψ(x,x0)+12∑k=0T−1ηk2‖gk‖∗2.

若对所选比较器 x 有正数 B 满足 Dψ(x,x0)≤B,取正的有效上界 G>0,固定 η=2B/(G2T)。对固定凸目标,令 x¯T=T−1∑k=0T−1xk,凸性与上述估计给出目标值差

f(x¯T)−f(x)≤1T∑k=0T−1⟨gk,xk−x⟩≤BηT+ηG22=G2BT.

这控制平均点的目标值,不直接控制 ‖x¯T−x‖。在线版本若 gk∈∂fk(xk),同一估计控制 T−1∑k(fk(xk)−fk(x)),即平均 regret。在 n 个坐标的单纯形上采用负熵并从 x0=(1/n,…,1/n) 均匀初始化时,才有对所有比较器成立的 Dψ(x,x0)≤log⁡n;任意内点初始化不享有这个统一常数。

若选 ψ(x)=12‖x‖22,Dψ 化为平方欧氏距离,在 f 可微且选 gk=∇f(xk) 时,子问题完成平方后得到投影梯度法;对不可微凸 f 则是投影次梯度更新;因此后者只是在欧氏镜像映射下的严格特殊情形。欧氏情形可取 D=U=Rn,所以迭代点位于 C 的边界也不妨碍镜像函数可微。例如 C=[0,∞)、f(x)=2x、x0=1、η0=1 时,平方欧氏镜像更新得到 x1=0;它在约束边界,却仍在 U 中。负熵选择则产生指数权重,适用于在线学习与概率单纯形。实际选择镜像函数时还要同时检查定义域、强凸常数、对偶范数和子问题成本,缺一项都不能仅凭漂亮几何宣称更快。

随机复合镜像下降保留这里的几何,再把条件无偏次梯度与新点结构项放在不同槽位。它完整保留结构项首尾差,证明更新前平均的期望界,并计算目标熵正则与镜像熵同时出现时的幂次更新;只有镜像几何并不表示目标已被正则化。

参考资料
  • 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。
关系图谱13 个相邻概念 · 3 类关系

拖动节点调整位置。

显示关系

显示:依赖

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

上位 / 更一般

暂未标注直接上位概念。

下位 / 直接特例

类型化关系