[论文解读] Measure-valued spline curves: an optimal transport viewpoint

Yongxin Chen, Giovanni Conforti+1 more|arXiv (Cornell University)|Jan 9, 2018

Advanced Numerical Analysis Techniques参考文献 13被引用 5

一句话总结

本文通过最优传输理论，将经典样条曲线推广至概率测度空间，提出测度值样条曲线。通过在Wasserstein空间中最小化均方加速度，构建了经验测度或高斯测度的光滑插值曲线；当边际分布为高斯分布时，可通过半定规划（SDP）获得显式解，确保方差演化过程的正定性与光滑性。

ABSTRACT

The aim of this article is to introduce and address the problem to smoothly interpolate (empirical) probability measures. To this end, we lift the concept of a spline curve from the setting of points in a Euclidean space that that of probability measures, using the framework of optimal transport.

研究动机与目标

• 开发一种系统性框架，用于对经验测度或随时间演化的概率测度进行光滑插值。
• 将经典样条曲线从欧氏空间中的点推广至Wasserstein空间中的概率测度。
• 确保插值曲线保持关键性质，如方差的正定性与演化过程的光滑性。
• 提供计算上可行的公式化方法，特别是针对高斯测度，利用半定规划求解。
• 建立与多边际最优传输及流体动力学公式的联系。

提出的方法

• 基于在Wasserstein空间中最小化均方加速度的变分样条公式。
• 通过最优传输几何将经典样条概念推广至测度值曲线，将测度视为黎曼流形类空间中的点。
• 将问题表述为具有连续时间测度联合分布约束的多边际最优传输问题。
• 通过构造条件独立的高斯向量图模型，确保时间序列上均值与协方差的一致性。
• 通过半定规划（SDP）求解所得优化问题，其复杂度在时间点数量上线性增长。
• 利用涉及Φ(t,0)和M(t,0)的矩阵函数，推导出插值协方差矩阵的显式参数形式，确保演化过程的光滑性与正定性。

实验结果

研究问题

• RQ1 如何将经典样条插值从R^d中的点推广至Wasserstein空间中的概率测度？
• RQ2 何种最优传输理论公式可使测度值样条曲线保持光滑性与方差正定性？
• RQ3 插值问题能否被简化为可计算的优化问题，特别是针对高斯测度？
• RQ4 测度值样条曲线与多边际最优传输之间存在何种结构性关系？
• RQ5 Benamou-Brenier的流体动力学公式如何推广至样条设置？

主要发现

• 测度值样条曲线集中在底层空间中的C²样条曲线上，确保插值轨迹的光滑性。
• 对于高斯测度，插值曲线为具有随时间平滑演化之均值与协方差的高斯流，可通过半定规划计算。
• 协方差演化最优解的秩至多为2d，意味着相邻时间边际之间存在确定性的线性映射。
• 公式(48)为复杂度在N上线性增长的半定规划，相较于原始公式(45)的N⁶复杂度有显著改进。
• 插值方差曲线始终为正且光滑，而标准R上的三次样条可能产生负方差。
• 显式公式(50)利用涉及Φ(t,0)和M(t,0)的矩阵函数，提供了时间可变协方差矩阵的闭式表达。