查看“︁沃瑟斯坦度量”︁的源代码

'''沃瑟斯坦度量'''（{{lang-en|Wasserstein metric}}）又称为'''沃瑟斯坦距离'''（{{lang|en|Wasserstein distance}}）、'''坎托罗维奇–鲁宾施泰因度量'''（{{lang|en|Kantorovich–Rubinstein metric}}），在数学中是指某一给定[[度量空间]]<math>M</math>上[[概率分布]]之间的[[度量空间|距离函数]]，其名称源于俄裔美国数学家{{le|列昂尼德·沃瑟斯坦|Leonid Vaserstein}}。

直观而言，可以将每个概率分布都看作''<math>M</math>''上堆积的单位数量的土堆，则沃瑟斯坦度量便是指将一个土堆变成另一土堆的最小代价，可定义为需要移动的土壤量乘以需要移动的平均距离。这一问题于1781年由法国数学家[[加斯帕尔·蒙日]]首次正式提出。基于上述土堆类比，该度量在计算机科学中也被称为[[推土机距离]]。

“沃瑟斯坦距离”这一名称是苏联数学家{{le|罗兰·多布鲁申|Roland Dobrushin}}于1970年提出的，他在沃瑟斯坦关于描述大型自动机系统的马尔可夫过程的工作中了解到了这一概念。<ref>{{Cite journal |vauthors=Vaserstein LN |year=1969 |title=Markov processes over denumerable products of spaces, describing large systems of automata |url=http://www.mathnet.ru/links/b90e28ab49fa878b215cbbdfa0235264/ppi1811.pdf |journal=Problemy Peredači Informacii |volume=5 |issue=3 |page=64–72}}</ref>不过早在1939年，另一位苏联数学家[[列昂尼德·坎托罗维奇]]就在研究货物的最优运输问题时最早定义了该度量。<ref>{{Cite journal |vauthors=Kantorovich LV |year=1939 |title=Mathematical Methods of Organizing and Planning Production |journal=Management Science |volume=6 |issue=4 |page=366–422 |doi=10.1287/mnsc.6.4.366 |jstor=2627082}}</ref>因此一些学者认为使用“坎托罗维奇度量”或“坎托罗维奇距离”来称呼此度量更为恰当。

== 定义 ==
对于一个[[度量空间]]<math>(M,d)</math>，假设''<math>M</math>''上每个[[博雷尔测度|博雷尔概率测度]]都是[[拉东测度]]（即该度量空间为拉东空间）。对有限<math>p\ge1</math>，<math>P_p(M)</math>表示''<math>M</math>''上所有有<math>p</math>阶[[矩 (數學)|矩]]的概率测度<math>\mu</math>的集合，即''<math>M</math>''中存在<math>x_0</math>满足

: <math>\int_M d(x, x_0)^{p} \, \mathrm{d} \mu (x) < \infty.</math>

于是，<math>P_p(M)</math>中两个概率测度''<math>\mu</math>''和<math>\nu</math>之间的'''<math>p</math>'''阶'''沃瑟斯坦距离'''可定义为

: <math>W_p (\mu, \nu):=\left( \inf_{\gamma \in \Gamma (\mu, \nu)} \int_{M \times M} d(x, y)^p \, \mathrm{d} \gamma (x, y) \right)^{1/p},</math>

其中<math>\Gamma(\mu,\nu)</math>表示<math>M \times M</math>上所有测度的集合，即''<math>\mu</math>''与''<math>\nu</math>''的所有{{le|耦合 (概率论)|Coupling (probability)|耦合}}组成的集合。

此外，沃瑟斯坦度量也可等效地定义为

: <math>W_{p} (\mu, \nu) = \left( \inf \operatorname{\mathbf{E}} \big[ d( X, Y )^p \big] \right)^{1/p},</math>

其中<math>\mathbf{E}[Z]</math>表示[[随机变量]]<math>Z</math>的[[期望值]]，[[下确界]]则由随机变量<math>X</math>和<math>Y</math>的所有联合分布确定，它们分别对应边缘分布''<math>\mu</math>''和<math>\nu</math>。

== 与最优运输问题的联系 ==
[[File:Transport-plan.svg|link=//upload.wikimedia.org/wikipedia/commons/thumb/7/74/Transport-plan.svg/220px-Transport-plan.svg.png|thumb|x轴与y轴上分别绘有两个一维分布<math>\mu</math>和<math>\nu</math>，中间是两者一种可能的联合分布（不唯一），而这一联合分布即定义了<math>\mu</math>和<math>\nu</math>之间的一个运输方案。]]
理解上述定义的一种方法是将其与[[最优运输]]问题联系起来。对于空间<math>X</math>上的某一质量分布<math>\mu(x)</math>，我们希望以某种方式运输其质量，使其转化同一个空间中的另一分布<math>\nu(x)</math>，即将“土堆”<math>\mu</math>转换为“土堆”<math>\nu</math>。两者的总质量必须相同，因而可以不失一般性地假设<math>\mu</math>和<math>\nu</math>是总质量为1的概率分布。此外还需假设代价函数

: <math>c(x,y) \mapsto [0,\infty),</math>

表示从点<math>x</math>运输质量到点<math>y</math>的代价。一个从<math>\mu</math>到<math>\nu</math>的运输方案可以用函数<math>\gamma(x,y)</math>来描述，该函数表明从<math>x</math>移动到<math>y</math>的质量。一个运输方案<math>\gamma(x,y)</math>必须满足以下性质

: <math>
\begin{align}
\int \gamma(x,y) \,\mathrm{d} y = \mu(x),\\
\int \gamma(x,y) \,\mathrm{d} x = \nu(y),
\end{align}
</math>

前者表示从某一点<math>x</math>移到其他所有点的土堆总质量必须等于最初该点<math>x</math>上的土堆质量，后者则表示从所有点移到某一点<math>y</math>的土堆总质量必须等于最终该点<math>y</math>上的土堆质量。

这一性质相当于要求<math>\gamma</math>是边缘分布<math>\mu</math>与<math>\nu</math>对应的[[联合分布|联合概率分布]]。于是可得到运输方案<math>\gamma</math>的总代价

: <math>
\iint c(x,y) \gamma(x,y) \, \mathrm{d} x \, \mathrm{d} y = \int c(x,y) \, \mathrm{d} \gamma(x,y).
</math>

方案<math>\gamma</math>并不是唯一的，所有可能的运输方案中代价最低的方案即为最优运输方案。最优运输方案的代价为

: <math>
C = \inf_{\gamma \in \Gamma(\mu, \nu)} \int c(x,y) \, \mathrm{d} \gamma(x,y).
</math>

如果两点之间移动的代价等于两点之间的距离，那么最小代价则与<math>W_1</math>距离等价。

== 参见 ==
* [[動土距離]]
* [[传输理论]]
* [[全变差距离]]
* [[柯尔莫哥洛夫-斯米尔诺夫检验]]

== 参考文献 ==
{{Reflist}}

{{Authority control}}
[[Category:統計距離]]
[[Category:概率分布理论]]
[[Category:度量几何]]
[[Category:测度论]]