[论文解读] Variational Principles for Minkowski Type Problems, Discrete Optimal Transport, and Discrete Monge-Ampere Equations

Xianfeng Gu, Feng Luo+2 more|arXiv (Cornell University)|Feb 22, 2013

Geometric Analysis and Curvature Flows参考文献 7被引用 33

一句话总结

本文通过在支撑高上定义凸优化问题，为离散最优传输、闵可夫斯基型问题和离散蒙日-安培方程建立了有限维变分原理。通过使用严格凸的能量泛函，提供了亚历山德罗夫定理的变分证明，其临界点在平移下唯一确定解，并通过海森矩阵上的牛顿法实现高效计算。

ABSTRACT

In this paper, we develop several related finite dimensional variational principles for discrete optimal transport (DOT), Minkowski type problems for convex polytopes and discrete Monge-Ampere equation (DMAE). A link between the discrete optimal transport, discrete Monge-Ampere equation and the power diagram in computational geometry is established.

研究动机与目标

• 提供亚历山德罗夫定理在给定法向量与面面积下，凸多面体存在性与唯一性的有限维变分证明。
• 建立离散最优传输、离散蒙日-安培方程与计算几何中幂图之间联系的变分框架。
• 开发一种利用凸优化与牛顿法高效计算具有有限像的亚历山德罗夫映射的算法。
• 通过能量泛函的严格凸性，证明离散海森映射的无穷小刚性。

提出的方法

• 将离散最优传输问题的解表述为定义在支撑高空间上的凸能量泛函的临界点。
• 将能量泛函定义为超水平集体积与目标面积之差的积分，并通过高度变量的线性项进行调整。
• 利用与体积梯度相关的微分1-形式的闭包性，确保定义良好势函数的存在性。
• 推导能量泛函的海森矩阵，并证明其正定且对角占优，从而保证严格凸性。
• 对能量泛函应用牛顿法，以高效计算最优高度向量的数值解。
• 利用分段线性凸函数与其勒让德变换之间的对偶性，将离散海森矩阵解释为类似沃罗诺伊胞的体积。

实验结果

研究问题

• RQ1 能否为亚历山德罗夫定理（关于给定面面积与法向量的凸多面体存在性与唯一性）构建一个变分证明？
• RQ2 如何将离散最优传输表述为具有几何意义的有限维凸优化问题？
• RQ3 在离散蒙日-安培方程的背景下，能量泛函的海森矩阵具有何种几何解释？
• RQ4 从高度向量到目标体积的梯度映射是否为局部微分同胚？在何种条件下其为全局满射？
• RQ5 能否通过在凸能量泛函上应用牛顿法，高效求解离散蒙日-安培方程的解？

主要发现

• 能量泛函 E(h) 在仿射子空间 H₀ 上严格凸，确保在平移下存在唯一最大化解。
• E(h) 的海森矩阵正定，保证牛顿法的局部收敛性，并表明梯度映射为局部微分同胚。
• 解向量 h 是 E(h) 的唯一最大化解，且对应的分段线性函数 u_h 实现了最小化二次成本的离散最优传输映射。
• 梯度映射 ∇E: H₀ → A 是从 H₀ 到目标体积正卦限的全局微分同胚，证明了解的存在性与唯一性。
• 在每个点 p_i 处，对偶函数的离散海森矩阵等于对应胞的体积 W_i(h)，为离散蒙日-安培方程提供了几何解释。
• 该方法通过牛顿法实现亚历山德罗夫映射的高效数值计算，严格凸性保证了收敛性。