形式陈述
Monge 映射可能因为原子不能拆分而无解。若只规定每个出发地供给多少、每个目的地接收多少,搬运怎样成为一个可行的优化问题?
对Polish 空间公理库Polish 空间Polish space可分且可完全度量化的拓扑空间;用换度量、可数乘积和子空间例子区分拓扑性质与指定度量的完备性。 上的 Borel 概率分布公理库概率分布Probability distribution · Law可测空间上总质量为一的测度;随机变量的律是由样本概率推出的一类分布。 ,令 为所有边缘分别为 的联合分布公理库联合分布Joint distribution · 联合概率分布多个随机元素组成的向量所推出的概率测度,完整记录边缘与依赖结构。。给定有有限常数下界的 Borel 成本 ,Kantorovich 问题为
目标采用Lebesgue 积分公理库Lebesgue 积分Lebesgue integral从简单函数积分出发,按单调逼近定义非负、扩展值与可积函数的积分。;成本有下界,所以不存在正负无穷抵消。 称为输运计划, 表示从 送往 的质量。由乘积测度公理库乘积测度Product measure在乘积 σ-代数上把可测矩形的测度规定为边测度乘积的测度。构成的独立计划 总是可行,所以边缘为概率时可行集合非空。若 还下半连续,且至少有一个有限成本计划,则存在有限成本的最优计划。减去成本的常数下界后,就归约到非负成本版本;所有计划的总质量均为一,平移不改变最优计划。
有限分布 、 在所有 有限时,将问题化为线性规划:
直觉
输运矩阵的每一行是一份源质量的去向清单,每一列是一处目标收到的总量。固定边缘后,能够选择的是依赖结构;独立计划只是众多选择之一,通常并不最省成本。
存在性来自可行计划集合的紧性,而非成本一定严格凸。两组固定边缘在 Polish 空间上都是紧的。给定 ,分别取遗漏质量小于 的紧集 ,则任意耦合都满足 。Prokhorov 定理公理库Prokhorov 定理Prokhorov theorem · Prokhorov's theorem · 普罗霍罗夫定理在 Polish 空间上把概率测度族的一致紧性等价为弱拓扑中的相对紧性。因此给成本极小化序列提供弱收敛子列。对任意有界连续 , 在乘积上仍是有界连续探针,故极限保留第一边缘;第二边缘同理。最后,对平移后的非负下半连续成本使用Portmanteau 的积分下界公理库Portmanteau 定理Portmanteau theorem · Portmanteau lemma · 移植定理把概率测度的弱收敛等价改写为开闭集不等式、半连续函数不等式与连续集上的测度收敛。,得到极限计划的成本不超过下确界。成本线性,最优解可能不唯一。
例子与边界
取源分布 ,目标 ,成本为平方距离。源按 排列,目标按 排列,则
由四个边缘方程,所有可行计划可写成
成本为 ,所以取 得
这个计划把位于零点的半单位质量拆成两份,分别送到 和 。本例没有可行Monge 映射公理库Monge 输运问题Monge transport problem用确定映射搬运概率质量,先检查推前约束是否可行,再讨论最小成本。,因为两个源原子的质量均为 ,确定映射无法在目标 处恰好产生 质量。松弛真正扩大了可行集合。
若目标总质量不等于源总质量,以上行列和约束就不相容;非平衡输运需要另外允许创造、销毁或惩罚边缘误差。负成本若无下界,或缺少下半连续性,也可能破坏上述直接法条件。
推论与应用
Kantorovich 对偶公理库Kantorovich 对偶性Kantorovich duality用满足两侧势之和不超过成本的对偶函数,构造可核验的输运最优性证书。用两侧势函数为成本提供下界,能够证实计划最优而不枚举所有方案。若成本取距离的 次幂,再开 次方,就得到Wasserstein 距离公理库Wasserstein 距离Wasserstein distance用最小平均搬运距离定义概率测度的度量,并用矩条件说明它比弱收敛多控制什么。;不是每个任意成本的最小值都构成度量。
参考资料