[論文レビュー] Entropic approximation for multi-dimensional martingale optimal transport

Hadrien De March|arXiv (Cornell University)|Dec 28, 2018

Spectral Theory in Mathematical Physics参考文献 1被引用数 12

ひとこと要約

本稿では、多変数マルティングル-最適輸送問題を解くために、ニュートン法と組み合わせた新しいエントロピー正則化アルゴリズムを提案し、収束速度の向上と数値的性能の向上を達成した。理論的収束レート、エントロピー正則化の普遍的誤差境界、および非凸順序のマージナルを扱う実用的解決策を提供する。

ABSTRACT

We study the existing algorithms that solve the multidimensional martingale optimal transport. Then we provide a new algorithm based on entropic regularization and Newton's method. Then we provide theoretical convergence rate results and we check that this algorithm performs better through numerical experiments. We also give a simple way to deal with the absence of convex ordering among the marginals. Furthermore, we provide a new universal bound on the error linked to entropy.

研究の動機と目的

• 有限マージナルを伴う多変数マルティングル-最適輸送問題を解く際の計算的課題に対処すること。
• 従来のBregman射影法やSinkhorn法よりも優れた性能を示す数値的に効率的なアルゴリズムの開発。
• 提案されたエントロピー正則化アプローチの理論的収束レートの提供。
• マージナルが凸順序を満たさない場合に対処する実用的な手法の導入。
• マルティングル-最適輸送におけるエントロピー近似の誤差に対する普遍的境界の確立。

提案手法

• マルティングル-最適輸送の双対問題を解くために、エントロピー正則化とニュートン法に基づく新しいアルゴリズムを提案。
• 初期段階でBregman射影ステップを適用し、収束速度を向上させるために、切り替えをトリガとするハイブリッド戦略を採用。
• 最適化問題を滑らかにすることで、反復解法による効率的計算を可能にするためにエントロピー正則化を適用。
• エントロピー解のギブス構造を活用し、双対最適化子の明示的公式を導出。
• 誤差低減の度合いや反復回数に基づき、Bregman射影ステップとニュートンステップの間で動的切り替え基準を実装。
• 有限の一次モーメントを持つコスト関数 $ c $ とターゲット測度 $ \mu, \nu $ を用いて、離散グリッド上での輸送問題を定義。

実験結果

リサーチクエスチョン

• RQ1 エントロピー正則化とニュートン法を組み合わせることで、多変数マルティングル-最適輸送ソルバーの収束速度が著しく向上するか？
• RQ2 提案されたエントロピー正則化アルゴリズムの理論的収束レートは何か？
• RQ3 ターゲット測度 $ \nu $ が $ \mu $ に対して凸順序を満たさない場合、アルゴリズムはどのように処理できるか？
• RQ4 この文脈におけるエントロピー正則化によって生じる誤差の普遍的境界は何か？
• RQ5 多次元設定において、最適カップリングが観察されたように、少数の写像に集中するのか？

主な発見

• Bregman射影と切断ニュートンステップを組み合わせたハイブリッドアルゴリズムは、1次元および2次元の両方のテストケースにおいて、純粋なBregman法やニュートン法よりも計算時間が短い。
• $ 160 \times 160 $ グリッドと $ \varepsilon = 7.4 \times 10^{-3} $ の条件下で、2次元問題においてもハイブリッドアルゴリズムが最も優れた性能を示した。
• エントロピー近似誤差は普遍的に有界であり、特定のコスト関数やマージナルに依存しない理論的保証を提供する。
• 数値実験により、最適カップリングが1次元では最大2つの写像に集中することが確認され、これはすべての $ \mu, \nu $ に対して成り立つと推測される。
• 2次元では、条件付きカーネルが条件値に応じて1〜4点に集中しており、単純な普遍的構造は存在しないことが示された。
• 非凸順序マージナルに対しても、簡単で実用的な修正により効果的に対処でき、より広範な適用可能性が得られた。