Geometric Milstein Scheme for Stochastic Differential Equations on SO(n) and SE(n)

技術的概要：SO(n) および SE(n) 上の確率微分方程式に対する幾何学的ミルステイン法

問題定義
本論文は、リー群 $SO(n)$（特殊直交群）および$SE(n)$（特殊ユークリッド群）上で発展する非線形自律型伊藤確率微分方程式（SDE）の数値積分を取り扱います。これらの方程式は以下の形をとります：
$$dX(t) = X B_I(X(t))dt + \frac{1}{2}X \Sigma_N(X)dt + \sum_{j=1}^d B_j(X(t))dW_j(t)$$
ここで、$B_I$ および $B_j$ は接空間におけるドリフト場および拡散場であり、$\Sigma_N$ は伊藤補正に起因する法線空間における「ピンニング・ドリフト」項です。

核心的な課題は、**1/2 を超える強収束率（SCR）を達成する幾何学保存数値スキーム（GPNS）**を開発することです。既存の手法は以下のトレードオフに直面しています：

1. オイラー・マルヤマ変種（例：G-EM、TaSP-EM）は幾何学を保存しますが、SCR は 1/2（または 1/2 - $\epsilon$）に制限されます。
2. 高次精度手法（例：マグナス展開に基づく S-RKMK）はより高い次数を達成できますが、しばしばストラトノビッチ形式（非線形拡散に対しては失敗する）に依存するか、計算的に扱いにくい無限級数の切断を伴います。
3. 標準的なユークリッド空間のミルステイン法は SCR 1 を達成しますが、多様体制約を保存せず、ドリフトおよび拡散増分における法線成分により、軌道がリー群から逸脱してしまいます。

手法：接空間パラメータ化–修正ミルステイン法（TaSP-CM）
著者らは、リー代数における修正されたミルステイン増分と幾何学的再構成ステップを組み合わせることで、ユークリッド空間のミルステイン法をリー群に拡張するTaSP-CMスキームを提案します。

1. 接空間パラメータ化（TaSP）：
   多様体上で直接発展させるのではなく、$SO(n)$上の局所解$R(t)$ を以下のように表現します：
   $$R(t) = R_0 + R_0 Z(t) + R_0 C(Z(t))$$
   ここで、$Z(t) \in \mathfrak{so}(n)$ はリー代数における確率的変位であり、$C(Z) = \sqrt{I - Z^\top Z} - I$ は $R(t)^\top R(t) = I$ を保証する法線空間の調整項です。これにより、幾何学保存と確率積分が分離されます。

2. 修正ミルステイン増分：
   このスキームは、リー代数における増分 $Z(t)$ をミルステイン型の更新式を用いて計算します。決定的な点は、標準的なユークリッド空間の微分 $\nabla B_j(B_{j'})$ を内在的共変微分（CovD） $D^O_{B_j}(B_{j'})$ に置き換えることです。
   $$D^O_{B_j}(B_{j'}) = \nabla_{B_j} B_{j'} + \frac{1}{2}[B_j, B_{j'}]$$
   この補正は、リー代数の非可換構造（リー括弧積を通じて）を考慮し、増分が接空間内に留まることを保証します。これにより、多様体制約を破る偽の法線成分の導入を回避します。

3. 非可換ノイズの処理：
   非可換ノイズ（$[B_j, B_{j'}] \neq 0$）の場合、このスキームは二重伊藤積分の近似を必要とします。著者らは、必要な精度を維持しつつ計算コストの増大を招かないよう、ガウス変数の切断フーリエ級数を用いた既知の近似手法（アルゴリズム 2）を利用します。

4. SE(n) への拡張：
   半直積構造 $SE(n) = SO(n) \ltimes \mathbb{R}^n$ を活用し、この手法は自然に拡張されます。回転成分は $SO(n)$ のダイナミクスに従い、並進成分は回転行列を介して結合された標準的なユークリッド空間の SDE として発展します。

主要な貢献

• 計算的に実行可能な最初の SCR 1 スキーム： 本論文は、可換および非可換ノイズの両方において強収束率 1 を達成する、$SO(n)$および$SE(n)$ 上の非線形 SDE に対する計算的に実行可能な最初の GPNS を導入します。
• 収束の理論的証明： 著者らは、TaSP-CM スキームが SCR 1 を達成することを示す厳密な証明（定理 1 および 2）を提供します。この解析では、誤差を TaSP-CM 軌道と「凍結場」ユークリッド空間のミルステイン級数の差に分解し、幾何学的補正項が二次以上の誤差（平均二乗で $O(\delta^2)$）のみを導入することを示しています。
• SE(n) への一般化： この枠組みは $SE(n)$ へ成功裡に拡張され（定理 4）、剛体姿勢推定問題に対して SCR 1 の保証を維持します。
• アルゴリズム的実装： 論文は、平方根の調整および二重積分近似の処理を含む、実装のための具体的なアルゴリズム（アルゴリズム 1 および 3）を詳述しています。

結果
可換および非可換ノイズの設定下で、$SO(n)$および$SE(3)$ 上で数値実験が行われました：

• 収束率： TaSP-CM スキームは、理論的予測と一致する経験的 SCR 1 を一貫して示しました。対照的に、TaSP-EM および S-RKMK 手法は SCR 1/2 を示し、標準的なユークリッド空間のミルステイン法（Eu-Mil）は幾何学的逸脱により収束に失敗しました。
• 幾何学保存： このスキームは、他の幾何学保存手法と同程度の幾何学的偏差（$10^{-11}$ から $10^{-13}$ のオーダー）で直交性（$SO(n)$ の場合）を維持しました。一方、Eu-Mil は多様体から著しく逸脱しました。
• 効率性： 非可換ノイズにおける計算コストは、ガウス関数呼び出しで $O(\delta^{-3/2})$ とスケーリングし、ユークリッド空間のミルステイン法と一致します。これは、同程度の精度を達成するために低次幾何学手法が要求する $O(\delta^{-2})$ のコストに対して、大幅な改善です。

意義
本論文は、リー群上の SDE に関する数値解析における重要なギャップを埋めると主張しています。本研究以前には、非線形 SDE に対して、基礎となる幾何学的制約を保存しつつ、強収束次数を 1/2 超に達成する計算的に実行可能な手法は存在しませんでした。曲がった多様体へのミルステイン法の成功した拡張により、TaSP-CM 手法は、衛星制御、ロボット姿勢推定、物体追跡、および多様体上の確率制御など、高い数値精度と厳密な幾何学的忠実度を必要とする応用分野に対する実用的なツールを提供します。著者らは、このアプローチが、スタイフェル多様体やアフィン群などの他の行列多様体、およびリーマン拡散モデルへの応用への将来の拡張への扉を開くと指摘しています。