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

기술적 요약: SO(n) 및 SE(n) 위의 확률 미분 방정식에 대한 기하학적 밀레틴 scheme

문제 제기
본 논문은 리 군 $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$는 접공간에서의 드리프트 및 확산 장 (field) 이며, $\Sigma_N$은 이토 보정에서 비롯된 법선 공간의 "고정 - 드리프트 (pinning-drift)"항이다.

핵심적인 과제는 **1/2 보다 큰 강수렴 속도 (SCR)**를 달성하는 **기하학 보존 수치 scheme(GPNS)**을 개발하는 것이다. 기존 방법들은 다음과 같은 트레이드오프에 직면해 있다:

1. **Euler-Maruyama 변형 (예: G-EM, TaSP-EM)**은 기하학을 보존하지만 SCR 이 1/2(또는 $1/2 - \epsilon$) 로 제한된다.
2. **고차 방법 (예: Magnus 전개 기반의 S-RKMK)**은 더 높은 차수를 달성할 수 있지만, 종종 비선형 확산에 실패하는 Stratonovich 형식에 의존하거나 계산적으로 처리 불가능한 무한 급수 절단을 수반한다.
3. 표준 유클리드 밀레틴 scheme은 SCR 1 을 달성하지만, 드리프트 및 확산 증분 내의 법선 성분으로 인해 리 군에서 벗어나는 궤적을 초래하여 다양체 제약을 보존하지 못한다.

방법론: 접공간 매개변수화 - 보정 밀레틴 (TaSP-CM)
저자들은 리 대수 내의 보정된 밀레틴 증분과 기하학적 재구성 단계를 결합하여 유클리드 밀레틴 방법을 리 군으로 확장한 TaSP-CM scheme 을 제안한다.

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. 보정된 밀레틴 증분:
   이 scheme 은 리 대수 내에서 증분 $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$) 의 경우, 이 scheme 은 이중 이토 적분의 근사가 필요하다. 저자들은 계산 비용을 감당할 수 없게 하지 않으면서 필요한 정확도를 유지하기 위해 가우스 변수의 절단된 푸리에 급수를 포함하는 알려진 근사 방법 (알고리즘 2) 을 활용한다.

4. SE(n) 으로 확장:
   반직접곱 구조 $SE(n) = SO(n) \ltimes \mathbb{R}^n$을 활용하여 방법이 자연스럽게 확장된다. 회전 성분은 $SO(n)$ 동역학을 따르는 반면, 병진 성분은 회전 행렬을 통해 결합된 표준 유클리드 SDE 로 진화한다.

주요 기여

• 계산적으로 실현 가능한 최초의 SCR 1 scheme: 본 논문은 가환 및 비가환 잡음 하에서 SCR 1 을 달성하는 $SO(n)$및$SE(n)$ 위의 비선형 SDE 에 대한 최초의 계산적으로 실현 가능한 GPNS 를 소개한다.
• 수렴성의 이론적 증명: 저자들은 TaSP-CM scheme 이 SCR 1 을 달성함을 보여주는 엄밀한 증명 (정리 1 및 2) 을 제공한다. 분석은 오차를 TaSP-CM 궤적과 "고정된 장 (frozen-field)"유클리드 밀레틴 급수 간의 차이로 분해하여, 기하학적 보정 항이 평균 제곱에서 고차 오차 ($O(\delta^2)$) 만을 도입함을 보인다.
• SE(n) 으로 일반화: 프레임워크는 $SE(n)$으로 성공적으로 확장되어 (정리 4), 강체 포즈 추정 문제에 대해 SCR 1 보장을 유지한다.
• 알고리즘적 구현: 논문은 제곱근 조정 및 이중 적분 근사 처리를 포함한 구현을 위한 구체적인 알고리즘 (알고리즘 1 및 3) 을 상세히 기술한다.

결과
가환 및 비가환 잡음 설정 하에서 $SO(n)$및$SE(3)$에 대한 수치 실험이 수행되었다:

• 수렴 속도: TaSP-CM scheme 은 이론적 예측과 일치하는 일관된 경험적 SCR 1 을 보였다. 반면, TaSP-EM 및 S-RKMK 방법은 SCR 1/2 을 보였으며, 표준 유클리드 밀레틴 (Eu-Mil) scheme 은 기하학적 드리프트로 인해 수렴에 실패했다.
• 기하학 보존: 이 scheme 은 $SO(n)$의 경우 직교성을 유지하여$10^{-11}$에서$10^{-13}$ 수준의 기하학적 편차를 보였으며, 이는 다른 기하학 보존 방법과 비교 가능한 수준이다. 반면 Eu-Mil 은 다양체에서 상당히 벗어났다.
• 효율성: 비가환 잡음에 대한 계산 비용은 $O(\delta^{-3/2})$ 가우스 호출로 스케일링되어 유클리드 밀레틴 방법과 일치하며, 유사한 정확도를 달성하기 위해 저차 기하학 방법이 요구하는 $O(\delta^{-2})$ 비용보다 상당한 개선을 제공한다.

의의
본 논문은 리 군 위의 SDE 수치 분석에서 중요한 간극을 해소했다고 주장한다. 본 연구 이전에는 비선형 SDE 에 대해 기하학적 제약을 동시에 보존하면서도 1/2 보다 큰 강수렴 차수를 달성하는 계산적으로 실현 가능한 방법이 존재하지 않았다. TaSP-CM 방법은 밀레틴 scheme 을 곡면 다양체로 성공적으로 확장함으로써, 위성 제어, 로봇 포즈 추정, 객체 추적, 그리고 다양체 위의 확률적 제어와 같이 높은 수치 정확도와 엄격한 기하학적 충실도를 요구하는 응용 분야를 위한 실용적인 도구를 제공한다. 저자들은 이 접근 방식이 스테이펠 (Stiefel), 아핀 군과 같은 다른 행렬 다양체 및 리만 확산 모델의 응용으로의 미래 확장에 문을 열었다고 언급한다.