Skip to content

算法Algorithm

外梯度法

Extragradient method · Korpelevich method

在同一原锚点上先预测再校正,用两次场求值和投影处理单调旋转,证明有限维点收敛及平均预测点的有界域间隙率。

直接沿反馈场走一步,可能绕着平衡点越转越远。外梯度先问“按当前方向走到的预测点,会给出什么反馈”,再用那里的反馈从原点重新出发。第二次不是在预测点上继续走;这个锚点安排正是收敛证明能够抵消旋转误差的原因。

形式陈述 ​

两次投影,第二步仍从原点出发 ​

设 C⊆Rd 非空、闭且凸,F:C→Rd 单调且 L-Lipschitz,L>0。假定变分不等式有解。固定 0<γ<1/L,从 x0∈C 开始,精确执行

(1)yk=PC(xk−γF(xk)),xk+1=PC(xk−γF(yk)).

PC采用闭凸最近点投影。每个预测点和主点都在 C 中,所以仅需在 C 上能计算 F。有限维中,xk及yk收敛到同一个VI解。

若另有有限半径

R2=supu∈C‖x0−u‖2<∞,

输出预测点平均 y¯N=N−1∑k=0N−1yk,则

(2)GM(y¯N)≤R22γN.

这项速率控制平均预测点的Minty间隙,不是声称最后一点的距离为 O(1/N)。集合有界用于统一间隙常数;点收敛部分不要求 C 有界。

直觉

两条最近点不等式怎样支付场的变化 ​

一轮中简记 x=xk,y=yk,p=xk+1。第二次投影给任意 u∈C 的

2γ⟨F(y),p−u⟩≤‖x−u‖2−‖p−u‖2−‖x−p‖2.

第一次投影则给 γ⟨F(x),y−p⟩≤⟨x−y,y−p⟩。把目标内积拆成 ⟨F(y),y−u⟩,相加并展开平方,得到

2γ⟨F(y),y−u⟩≤‖x−u‖2−‖p−u‖2−‖x−y‖2−‖y−p‖2+2γ⟨F(y)−F(x),y−p⟩.

Lipschitz条件和 2ab≤a2+b2 将最后一项控制为 γL(‖x−y‖2+‖y−p‖2)。因此关键估计是

(3)2γ⟨F(y),y−u⟩≤‖x−u‖2−‖p−u‖2−(1−γL)(‖x−y‖2+‖y−p‖2).

两个短位移支付两处反馈差。若第二步改从 y 出发,上面第一行不再使用同一个平方距离中心,不能原样沿用式(3)。

为什么得到整列收敛 ​

取任一解 u=x∗。单调性与VI条件给 ⟨F(y),y−x∗⟩≥0。由于 1−γL>0,式(3)表明到每个解的距离不增,且 ∑k‖xk−yk‖2<∞。于是主点有界,预测位移趋于零。

由有限维Bolzano–Weierstrass定理取 xkj→x^。投影和 F 连续,ykj−xkj→0,故 x^=PC(x^−γF(x^)),它是VI解。再把距离不增结论中的解选成 x^:这份距离沿子列趋零,因此整列也趋零。预测点与主点相差趋零,亦有同一极限。

例子与边界

旋转场的全部线性轨道 ​

先取 C=R2、F=J,其中 J(x1,x2)=(−x2,x1),L=1。投影是恒等映射,式(1)化成

y=(I−γJ)x,p=((1−γ2)I−γJ)x.

因为 JT=−J、JTJ=I,有

‖p‖2=(1−γ2+γ4)‖x‖2.

当 0<γ<1 时系数小于一,故收敛到唯一零点。取 γ=1/2,x0=(1,0),前两轮为

y0=(1,−1/2),x1=(3/4,−1/2),y1=(1/2,−7/8),x2=(5/16,−3/4).

主点平方长度每轮乘 13/16;只有一次前向更新时则每轮乘 5/4。算法多出的查询不是更精细地沿旧方向走,而是修正方向本身。

图中三条轨道使用同一初值和步长。右图比较主更新轮数,不表示两次场求值与一次完整预解的运算成本相同;这里的隐式轨道沿用单调算子页的精确矩阵。

有界盒中的平均输出证书 ​

改取 C=[−1,1]2、x0=(1/2,0),步长仍为 1/2。主点长度不断缩小,预测点长度至多为 5/4/2<1,所以两类点都在盒内,投影不改变上面的线性公式。

此时 R2=(3/2)2+1=13/4,而盒旋转的 GM(x)=|x1|+|x2|。因而任何 N≥1 都能同时报告实际平均间隙和保守上界 13/(4N)。例如要按这个上界保证间隙不超过 0.01,取 N≥325 即可;需要650次场求值和650次投影调用,盒投影在本例虽不改变值,验证其不活跃仍是模型分析的一部分。

步长端点、解存在与误差来源 ​

旋转例中若 γ=1,主点更新为 −Jx,非零初值形成四周期;若 γ>1,平方放大系数大于一。这给出不能把严格步长条件任意放宽的具体反例,但没有说所有问题在端点都失败。

若没有解,距离证明没有可比较的目标。例如全空间常值场 F(x)=1 是单调且满足任意正Lipschitz上界,却无零点,迭代持续平移。非单调场、带偏查询或近似投影也需要新的分析;名称相同不能保留当前常数。

推论与应用

平均间隙的望远镜证明 ​

对固定 u∈C,单调性给 ⟨F(u),yk−u⟩≤⟨F(yk),yk−u⟩。式(3)求和、舍去非正项,有

2γ∑k=0N−1⟨F(u),yk−u⟩≤‖x0−u‖2.

除以 2γN,再对 u取上确界,便得式(2)。不能把平均 y¯N 换成 x¯N 或最后一点而仍称使用同一证明;求和里的实际评估点是 yk。

每轮有两次 F 与两次投影,另需 O(d) 向量运算和 O(d) 工作存储;运行平均无需保存全部历史。若 F(x)=Mx+b,稠密矩阵每次求值为 O(d2),稀疏存储 s 个条目时为 O(s+d);投影本身的子问题成本另计。

在Mirror-Prox中选平方欧氏生成函数,两次Bregman子问题就严格化为式(1),故本法是其欧氏特例。一般范数可以改善有界域常数,但要重新计算强凸系数与对偶Lipschitz界,不能只替换图中的距离名称。

Tseng的前向–后向–前向分裂在法锥情形只做一次投影,再加显式场差校正;校正状态可能离开可行集。本页的两次投影始终保持可行,因此只要求场在 C上有定义。两种更新在无约束零算子情形重合,一般约束下并不相同。

参考资料
  • G. M. Korpelevich, “The Extragradient Method for Finding Saddle Points and Other Problems,” Ekonomika i Matematicheskie Metody12(4),747–756,1976,方法原始出处;书目与下列原文参考文献相核。
  • Arkadi Nemirovski, Prox-Method with Rate of Convergence O(1/t),2004,§§2–3,VI误差与双步近端机制。本文欧氏式(3)、点收敛和旋转轨道直接推导。
  • Ernest K. Ryu and Stephen Boyd, A Primer on Monotone Operator Methods,2016,§7.2,区别全空间外梯度与一般Tseng校正。
关系图谱13 个相邻概念 · 4 类关系

拖动节点调整位置。

显示关系

显示:依赖

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

上位 / 更一般

下位 / 直接特例

暂未标注直接特例。

类型化关系