[论文解读] Optimal transport maps in Monge-Kantorovich problem

Luigi Ambrosio|ArXiv.org|Apr 24, 2003

Geometric Analysis and Curvature Flows参考文献 1被引用 34

一句话总结

本文通过 Γ-收敛的奇异扰动方法，建立了在 Monge-Kantorovich 问题中最优传输映射的存在性，证明了标准代价 $|x-y|$ 的最优映射可作为正则化问题中 $c_ ho = \|x-y\|^{1+\epsilon}$ 解的极限获得。关键结果表明，当初始测度 $\mu$ 绝对连续时，极限映射在扰动下既最优又稳定，从而在凸性和正则性条件下，解决了欧几里得和黎曼几何设定下的长期存在性问题。

ABSTRACT

In the first part of the paper we briefly decribe the classical problem, raised by Monge in 1781, of optimal transportation of mass. We discuss also Kantorovich's weak solution of the problem, which leads to general existence results, to a dual formulation, and to necessary and sufficient optimality conditions. In the second part we describe some recent progress on the problem of the existence of optimal transport maps. We show that in several cases optimal transport maps can be obtained by a singular perturbation technique based on the theory of $Γ$-convergence, which yields as a byproduct existence and stability results for classical Monge solutions.

研究动机与目标

• 解决当代价为欧几里得距离时，Monge 问题中最优传输映射存在的长期悬而未决问题。
• 发展一种变分正则化技术，从 Kantorovich 对偶公式中选择由传输映射诱导的极值计划。
• 建立代价函数扰动下最优映射的稳定性，特别是在 $\epsilon \to 0$ 时 $c_\epsilon = \|x-y\|^{1+\epsilon}$ 的极限情形。
• 将存在性结果扩展至严格凸范数之外，包括晶体范数，使用高阶变分校正。
• 阐明传输密度及其唯一性在选择最优计划中的作用，特别是当 $\mu$ 绝对连续时。

提出的方法

• 使用 Γ-收敛分析正则化代价 $c_\epsilon(x,y) = \|x-y\|^{1+\epsilon}$ 在 $\epsilon \downarrow 0$ 时最优传输映射的极限，证明其收敛于原始代价 $c(x,y) = \|x-y\|$ 的解。
• 应用二级变分原理：极限计划在所有最优 Kantorovich 计划中最小化 $\int \|x-y\|\ln\|x-y\|\,d\gamma$，从而确保唯一性与映射诱导。
• 采用 Sudakov 类型的射线分解，利用代价函数的几何结构与测度的分解。
• 在局部坐标下利用面积公式与 Fubini 型定理处理分解下的绝对连续性。
• 建立传输密度 $\sigma$ 与质量优化问题解之间的 1-1 对应关系，其中 $\sigma(B) = \int \mathcal{H}^1(B \cap [x,y])\,d\gamma(x,y)$。
• 应用对偶公式 $\min(\text{MK}) = \sup\left\{ \int h\,d\mu + \int k\,d\nu : h(x) + k(y) \leq c(x,y) \right\}$ 推导最优性条件。

实验结果

研究问题

• RQ1 在什么条件下，Monge 问题中代价为 $c(x,y) = \|x-y\|$ 的最优传输映射存在？
• RQ2 最优传输映射能否作为正则化问题 $c_\epsilon = \|x-y\|^{1+\epsilon}$ 解的极限获得？
• RQ3 在所有最优 Kantorovich 计划中，对 $\int \|x-y\|\ln\|x-y\|\,d\gamma$ 的二次最小化是否能选择唯一一个由传输映射诱导的计划？
• RQ4 初始测度 $\mu$ 的绝对连续性是否为欧几里得设定下最优传输映射存在的必要条件？
• RQ5 Γ-收敛方法能否扩展至非欧几里得范数（如晶体范数）？

主要发现

• 对于紧支集测度 $\mu, \nu$ 且满足 $\mu \ll \mathcal{L}^n$ 的情形，当 $\mu$ 绝对连续时，$c(x,y) = \|x-y\|$ 的最优传输映射存在。
• 最优映射作为求解正则化问题 $c_\epsilon(x,y) = \|x-y\|^{1+\epsilon}$ 的映射 $\psi_\epsilon$ 的极限获得，且该极限在 $\epsilon \downarrow 0$ 时稳定。
• 极限计划 $\gamma_0$ 不仅对原始 Kantorovich 问题最优，还最小化二次泛函 $\int \|x-y\|\ln\|x-y\|\,d\gamma$，从而确保唯一性与映射诱导。
• 当 $\mu$ 或 $\nu$ 绝对连续时，由 $\sigma(B) = \int \mathcal{H}^1(B \cap [x,y])\,d\gamma(x,y)$ 定义的传输密度 $\sigma$ 是唯一的。
• 作为 $\psi_\epsilon$ 极限构造的映射 $\psi$ 与早期工作（如 Caffarelli 等人）中的结果一致，但新证明通过变分选择建立了稳定性与唯一性。
• 对于晶体范数，采用高阶正则化 $c_\epsilon(x,y) = \|x-y\| + \epsilon|x-y| + \epsilon^2|x-y|\ln|x-y|$ 可实现向最优传输映射的收敛，三重变分问题选择出唯一一个诱导映射的计划。