[論文レビュー] 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空間における平均二乗加速度の最小化により、経験的またはガウス測度の滑らかな補間を構築し、周辺分布がガウス分布の場合には半正定値計画法を用いて明示的な解が得られ、時間経過における分散の正値性と滑らかさが保証される。

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) を含む行列関数を用いて、時間に依存する共分散行列の閉形式表現が得られる。