Skip to content

算法Algorithm

Mirror-Prox 方法

Mirror-Prox method · Mirror extragradient

在同一镜像锚点上执行预测与校正,以Bregman三点恒等式证明单调Lipschitz问题的平均间隙率,并完整计算负熵矩阵博弈的双步更新。

概率向量的更新需要保持非负和归一化,欧氏投影并非唯一选择。Mirror-Prox保留外梯度的两次查询,却让每次位移由镜像函数计价。两次子问题共享旧锚点;预测点负责查询和最后的平均,校正点负责把下一轮接起来。

形式陈述 ​

几何、定义域与两次调用 ​

设 C⊆Rd 为非空闭凸集,F:C→Rd 是单调VI的场。在指定范数及其对偶范数下,假定

‖F(x)−F(y)‖∗≤L‖x−y‖,L>0.

沿用Bregman散度的边界扩展:ψ在包含 C 的凸域上有限且凸,在开集 U 内可微;第一槽可属于 C,第二槽必须属于 C∩U。假定 ψ|C关于当前范数为1-强凸,所用子问题均取得唯一解且解留在 C∩U。

对锚点 z∈C∩U,定义

Pz(a)=arg⁡minu∈C{⟨a,u⟩+Dψ(u,z)}.

固定 0<γ≤1/L,从 z0∈C∩U开始,执行

(1)wk=Pzk(γF(zk)),zk+1=Pzk(γF(wk)).

注意第二个下标仍是 zk。假定初始半径 Θ=supu∈CDψ(u,z0) 有限,输出 w¯N=N−1∑k=0N−1wk,则

(2)GM(w¯N)≤ΘγN.

迭代定义及逐个比较点的不等式不要求 C有界;有限的全域间隙预算则额外需要 Θ<∞。例如欧氏全空间可以执行两步法,但非零旋转点的全域间隙仍可能为无穷。

本页认证精确查询下的平均间隙;不从式(2)额外宣布最后锚点有同一速率或一般点收敛。

直觉

一次镜像最优性,接两次不同反馈 ​

子问题的一阶方向条件与Bregman三点恒等式给出:若 p=Pz(a),则对所有 u∈C,

(3)⟨a,p−u⟩≤Dψ(u,z)−Dψ(u,p)−Dψ(p,z).

即使比较点 u在约束边界,式(3)仍可用;可微性只在锚点及实际子问题解处调用。

一轮简记 z,w,q=zk+1。对校正步应用式(3),再对预测步取比较点 q,得到

(4)γ⟨F(w),w−u⟩≤Dψ(u,z)−Dψ(u,q)−Dψ(q,w)−Dψ(w,z)+γ⟨F(w)−F(z),w−q⟩.

对偶范数的基本配对界将交叉项控制为 γL‖w−z‖‖w−q‖。强凸性保证两份被减去的散度至少分别为相应平方范数的一半,故在 γL≤1 时,交叉项全部由它们支付。于是

(5)γ⟨F(w),w−u⟩≤Dψ(u,z)−Dψ(u,q).

固定 u,由单调性可将左边的 F(w)换成 F(u)作为下界。对轮次求和、除以 γN,再对 u取上确界,即得式(2)。这份支付机制利用了场的Lipschitz差,而不是给每次场值本身设一个上界。

例子与边界

两个概率向量的完整预测与校正 ​

考虑行方最小化、列方最大化的矩阵博弈

minp∈Δ2maxq∈Δ2pTMq,M=(1−1−11).

对应场 F(p,q)=(Mq,−MTp)单调:两个交叉内积相消。选乘积范数

‖(a,b)‖=‖a‖12+‖b‖12,‖(c,d)‖∗=‖c‖∞2+‖d‖∞2.

由于 |Mij|≤1,有 ‖Mv‖∞≤‖v‖1,所以可取 L=1。镜像函数为两个负熵之和,散度是两个以自然对数计算的KL之和;由负熵镜像的强凸性,它关于上述乘积范数为1-强凸。

在每个单纯形上,镜像子问题返回

[Pp(γa)]i=pie−γai∑jpje−γaj.

从 p0=(3/4,1/4),q0=(1/2,1/2)出发,取 γ=log⁡2<1。初始反馈为 (0,0;−1/2,1/2),因此预测点为

wp,0=(3/4,1/4),wq,0=(2/3,1/3).

预测反馈变成 (1/3,−1/3;−1/2,1/2)。再从原来的 p0,q0校正,得到

p1=(33+22/3,22/33+22/3)≈(0.6539650591,0.3460349409),q1=(2/3,1/3).

若把第二次的 p0擅自换成 wp,0,第一轮恰好看不出差别,因为两者此时相同;列方若改用 wq,0则会再次乘权,已经不等于上式。必须逐槽检查完整算法,不能靠一个坐标的巧合验证实现。

真正能计算的鞍点间隙 ​

鞍点为均匀策略对,值为零。对任意候选 (p,q),全局鞍点间隙为

maxvpTMv−minuuTMq=|2p1−1|+|2q1−1|.

双线性结构使它也恰好等于该VI的Minty间隙。最坏初始Bregman距离由两个顶点取得:Θ=log⁡4+log⁡2=log⁡8,所以式(2)为 3/N。

第一轮预测点间隙为 5/6,两轮平均预测点间隙约为 0.7141545637,三轮约为 0.5004984598;对应统一上界为 3,1.5,1。实际间隙可以远小于保证。若目标是按此界认证不超过 0.01,需300轮、600次完整场调用和600次乘积镜像子问题。

初始零权重、尺度与随机误差 ​

负熵梯度在零坐标不有限,本页从内部初始化。即使程序强行沿用乘法公式,初始为零的坐标也永远为零。例如行方始终被锁在 (1,0)时,上述博弈间隙至少为一,不能逼近均匀鞍点。

若镜像函数仅为 α-强凸,应相应使用 γL≤α,或先把镜像函数归一化;若场整体放大十倍,Lipschitz常数也放大十倍。仅调整名义步长而不核对范数与尺度,不能保留同一个预算。

紧集本身也不保证任意镜像半径有限,尤其要检查生成函数在比较点处的延拓。带噪反馈、近似子问题或离开可微域,都增加了式(4)之外的误差项;本页不将确定性精确公式冒称随机保证。

推论与应用

两查询换来的保证与计算成本 ​

已有镜像下降对一般有界次梯度使用一次查询,典型平均目标界为 O(N−1/2)。本方法用两次查询利用单调Lipschitz结构,获得 O(N−1)间隙。输入假设和输出量都不同,不能只比较幂次就宣布对任意目标都更快。

若 ψ(z)=‖z‖22/2,两次子问题分别为欧氏投影,得到外梯度法。本页在步长端点仍证明平均界;外梯度页为整列收敛采用严格不等式,两个结论的边界不应混写。

一般每轮有两次场调用、两次镜像子问题和 O(d)状态更新。上述 m×n 矩阵博弈每次场调用计算 Mq及MTp,稠密成本 O(mn),两组归一化成本 O(m+n);保存运行平均另需 O(m+n)空间。指数权重实现可先减去共同最大对数权重,避免上溢;舍入和下溢误差若需要严格认证,须另行纳入数值误差预算。

参考资料
关系图谱14 个相邻概念 · 3 类关系

拖动节点调整位置。

显示关系

显示:依赖

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

上位 / 更一般

暂未标注直接上位概念。

下位 / 直接特例

类型化关系