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

技术摘要：SO(n) 和 SE(n) 上随机微分方程的几何 Milstein 格式

问题陈述
本文针对在李群 $SO(n)$（特殊正交群）和$SE(n)$（特殊欧几里得群）上演化的非线性自治 Itô 随机微分方程（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$ 是源于 Itô 修正的正交空间中的“钉扎漂移”项。

核心挑战在于开发能够实现强收敛阶（SCR）大于 $1/2$ 的几何保持数值格式（GPNS）。现有方法面临权衡：

1. Euler-Maruyama 变体（如 G-EM、TaSP-EM）保持几何结构，但强收敛阶仅限于 $1/2$（或 $1/2 - \epsilon$）。
2. 高阶方法（如基于 Magnus 展开的 S-RKMK）可以达到更高阶，但通常依赖于 Stratonovich 表述（对非线性扩散失效），或涉及计算上不可行的无穷级数截断。
3. 标准欧几里得 Milstein 格式虽然能达到 SCR 1，但无法保持流形约束，由于漂移和扩散增量中的法向分量，导致轨迹偏离李群。

方法论：切空间参数化 - 修正 Milstein（TaSP-CM）
作者提出了 TaSP-CM 格式，通过将李代数中的修正 Milstein 增量与几何重构步骤相结合，将欧几里得 Milstein 方法推广到李群上。

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. 修正 Milstein 增量：
   该格式利用 Milstein 型更新计算李代数中的增量 $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$），该格式需要近似双重 Itô 积分。作者利用了一种已知近似方法（算法 2），涉及高斯变量的截断傅里叶级数，以在避免不可行计算成本的同时保持所需的精度。

4. 推广至 SE(n)：
   利用半直积结构 $SE(n) = SO(n) \ltimes \mathbb{R}^n$，该方法自然扩展。旋转分量遵循 $SO(n)$ 动力学，而平移分量作为标准欧几里得 SDE 演化，并通过旋转矩阵耦合。

主要贡献

• 首个计算可行的 SCR 1 格式：本文引入了首个针对 $SO(n)$和$SE(n)$ 上非线性 SDE 的计算可行 GPNS，在交换和非交换噪声下均实现了强收敛阶 1。
• 收敛性理论证明：作者提供了严格的证明（定理 1 和定理 2），表明 TaSP-CM 格式达到了 SCR 1。分析将误差分解为 TaSP-CM 轨迹与“冻结场”欧几里得 Milstein 级数之间的差异，表明几何修正项仅引入高阶误差（均方意义下为 $O(\delta^2)$）。
• 推广至 SE(n)：该框架成功扩展至 $SE(n)$（定理 4），为刚体位姿估计问题保持了 SCR 1 保证。
• 算法实现：文章详细说明了具体实现算法（算法 1 和算法 3），包括平方根调整和双重积分近似的处理。

结果
在交换和非交换噪声设置下，对 $SO(n)$和$SE(3)$ 进行了数值实验：

• 收敛阶：TaSP-CM 格式一致表现出 1 的经验强收敛阶，与理论预测相符。相比之下，TaSP-EM 和 S-RKMK 方法表现出 1/2 的强收敛阶，而标准欧几里得 Milstein（Eu-Mil）格式因几何漂移而未能收敛。
• 几何保持：该格式保持了正交性（针对 $SO(n)$），几何偏差量级为$10^{-11}$至$10^{-13}$，与其他几何保持方法相当，而 Eu-Mil 则显著偏离流形。
• 效率：对于非交换噪声，计算成本按 $O(\delta^{-3/2})$ 次高斯调用缩放，与欧几里得 Milstein 方法相当，并显著优于为达到类似精度而需要 $O(\delta^{-2})$ 成本的低阶几何方法。

意义
本文声称填补了李群上 SDE 数值分析中的一个关键空白。在此之前，没有任何计算可行的方法能同时保持底层几何约束并为非线性 SDE 实现大于 $1/2$ 的强收敛阶。通过将 Milstein 格式成功推广到弯曲流形，TaSP-CM 方法为需要高数值精度和严格几何保真度的应用提供了实用工具，例如卫星控制、机器人位姿估计、目标跟踪以及流形上的随机控制。作者指出，这种方法为未来扩展到其他矩阵流形（如 Stiefel 流形、仿射群）以及黎曼扩散模型中的应用打开了大门。