[论文解读] A new optimal transport distance on the space of finite Radon measures

Stanislav Kondratyev, Léonard Monsaingeon+1 more|arXiv (Cornell University)|May 28, 2015

Geometric Analysis and Curvature Flows参考文献 21被引用 46

一句话总结

本文提出了一种针对质量可变的非负有限Radon测度的新最优传输距离，基于使用非守恒连续性方程的改进版Benamou-Brenier公式。该框架支持形式黎曼结构，支持梯度流分析，并证明了在适应度驱动的扩散模型中指数收敛至平衡态，将奥托的微分学推广至非质量守恒的动力学。

ABSTRACT

We introduce a new optimal transport distance between nonnegative finite Radon measures with possibly different masses. The construction is based on non-conservative continuity equations and a corresponding modified Benamou-Brenier formula. We establish various topological and geometrical properties of the resulting metric space, derive some formal Riemannian structure, and develop differential calculus following F. Otto's approach. Finally, we apply these ideas to identify an ideal free distribution model of population dynamics as a gradient flow and obtain new long-time convergence results.

研究动机与目标

• 开发一种适用于有限Radon测度的度量，以容纳质量变化，克服经典Wasserstein距离仅限于等质量测度的局限性。
• 通过改进的Benamou-Brenier公式，将奥托的形式黎曼微积分推广至非守恒动力学。
• 建立所得度量空间的几何与拓扑性质，包括完备性与测地线结构。
• 识别并分析MacCall与Cosner的适应度驱动扩散模型在新形式体系下的梯度流性质。
• 推导具有显式速率的长期收敛结果，解决先前研究中依赖反证法而缺乏显式速率的缺陷。

提出的方法

• 通过最小化动能与势能作用的动态公式提出新距离：$ \int \rho(|\nabla u|^2 + |u|^2) \, dx \, dt $，其中 $ \rho $ 通过非守恒连续性方程 $ \partial_t \rho + \text{div}(\rho \nabla u) = \rho u $ 演化。
• 利用非守恒连续性方程建模质量的产生与消亡，其中 $ u $ 表示驱动扩散与增长的适应度。
• 为有限Radon测度空间赋予类似于奥托框架的形式黎曼结构，支持一阶与二阶微分计算。
• 建立拓扑性质：该度量度量窄收敛，关于弱-*收敛是下半连续的，且该空间为完备测地空间。
• 通过最优速度场解的连续性方程表征该空间中的Lipschitz曲线与测地线。
• 将该框架应用于适应度驱动的扩散模型，证明其在新度量下为梯度流。

实验结果

研究问题

• RQ1 能否构造一种新的最优传输距离，允许质量变化且无需矩或紧性假设？
• RQ2 该新距离是否支持形式黎曼结构，从而在奥托工作的精神下实现微分计算？
• RQ3 MacCall与Cosner的适应度驱动扩散模型是否在此新形式体系中为梯度流？
• RQ4 能否利用该框架建立向理想自由分布的长期收敛性，并获得显式收敛速率？
• RQ5 所得有限Radon测度度量空间的拓扑与几何性质为何？

主要发现

• 新距离度量窄收敛，且关于弱-*收敛是下半连续的。
• 配备该距离的有限Radon测度空间为完备测地空间，且紧支集测度在此空间中稠密。
• 该空间中的Lipschitz曲线与测地线由非守恒连续性方程在最优速度场下的解表征。
• 适应度驱动的扩散模型被识别为该形式体系下的梯度流，从而可推导出其指数收敛至理想自由分布。
• 证明了一类广义Beckner型不等式：$ \Phi(\int \rho) \int |\rho - m|^2 \leq \int \rho |\rho - m|^2 + \int \rho |\nabla(\rho - m)|^2 $，其中 $ \Phi $ 仅依赖于区域与 $ m $，该不等式为收敛结果提供了理论基础。
• 通过该不等式量化了收敛速率，表明在 $ L^2 $-型范数下收敛至平衡态呈指数衰减，解决了文献中此前缺乏显式速率的缺陷。