Skip to content

Bregman 散度

Bregman divergence

凸函数值高于其一阶切平面的余量,用来度量与生成函数几何相适配的偏离。

定义与非负性

C 是内积空间中的凸集,ψ:CR 在所讨论点可微且凸。对 x,y(其中 ψ(y) 存在),定义

Dψ(x,y)=ψ(x)ψ(y)ψ(y),xy.

凸函数的一阶条件给出

ψ(x)ψ(y)+ψ(y),xy,

所以 Dψ(x,y)0。它描绘函数图像在 x 处比从 y 出发的切平面高出多少;改变生成函数 ψ,也就改变了“偏离”的几何。

ψ 关于范数 α-强凸的,则

Dψ(x,y)α2xy2.

严格凸通常保证 xy 时散度为正,但这并未使它成为距离。

三点恒等式

直接展开定义可得

Dψ(x,z)Dψ(x,y)Dψ(y,z)=ψ(y)ψ(z),xy.

欧氏平方距离的勾股关系在非欧几里得优化中由这个恒等式取代。镜像下降证明会把一次更新的线性损失,与两个相邻 Bregman 散度之差相连;对时间求和后,大部分项望远镜消去。

两种几何

ψ(x)=12x22,便有

Dψ(x,y)=12xy22.

因此欧氏投影梯度可以看成 Bregman 几何的特例。对概率单纯形内部取负熵 ψ(p)=ipilogpi,则

Dψ(p,q)=ipilogpiqi=DKL(pq).

这说明 KL 散度是特定方向、特定定义域上的 Bregman 散度。负熵的镜像更新自然产生乘法权重,而非欧氏加法步;这不是换一种记号,而是几何与约束集合相匹配的结果。

为什么它不是度量

一般有 Dψ(x,y)Dψ(y,x)。例如在正实数上取 ψ(u)=uloguDψ(1,2)=1log2,而 Dψ(2,1)=2log21。它也通常不满足三角不等式;即使平方欧氏例子也不满足普通三角不等式。因而“散度”只承诺非负及同点为零,不承诺度量空间结构。

另一个边界来自定义域。负熵在单纯形边界可连续延拓函数值,但梯度含 logpi,在 pi=0 处并非有限。书写更新与三点恒等式时,应把第二个参数限制在相对内部,或显式采用次梯度与广义 Bregman 散度;不能默认任意凸函数处处可微。

Bregman 投影也揭示参数方向的重要性。给定凸集 C,常定义 ΠCψ(y)=argminxCDψ(x,y):优化变量在第一个槽位,梯度锚点是旧点 y。交换成 Dψ(y,x) 一般得到另一问题;只有平方欧氏生成函数的对称特例看不出差别。

取一维 ψ(u)=12u2,三点恒等式退化为平方展开;取负熵并在单纯形上最小化“线性损失 ηg,pDKL(pq)”,一阶条件给 piqieηgi。这条乘法更新具体展示了生成函数如何把同一线性信号转成不同几何中的一步,而不是把 KL 仅当作评价预测分布的统计量。

参考资料
  • Lev M. Bregman, “The Relaxation Method of Finding the Common Point of Convex Sets,” USSR Computational Mathematics and Mathematical Physics, 1967.
  • Amir Beck, Marc Teboulle, “Mirror Descent and Nonlinear Projected Subgradient Methods,” Operations Research Letters, 2003.