[论文解读] How to compute the length of a geodesic on a Riemannian manifold with small error in arbitrary Sobolev norms

Joerg Kampen|ArXiv.org|Nov 11, 2008

Differential Equations and Boundary Problems参考文献 5被引用 4

一句话总结

本文提出了一种基于多项式插值的稳定算法，可在任意 Sobolev 范数下以任意小的误差计算黎曼流形上的测地线长度。通过使用全局解的规则多项式插值求解 eikonal 方程，该方法实现了长度泛函的高精度逼近，并支持抛物型 PDE 的 WKB 展开的精确计算，在数学物理和金融学中具有应用价值。

ABSTRACT

We compute the length of geodesics on a Riemannian manifold by regular polynomial interpolation of the global solution of the eikonal equation related to the line element $ds^2=g_{ij}dx^idx^j$ of the manifold. Our algorithm approximates the length functional in arbitrarily strong Sobolev norms. Error estimates are obtained where the geometric information is used. It is pointed out how the algorithm can be used to get accurate approximation of solutions of parabolic partial differential equations leading obvious applications to finance and physics.

研究动机与目标

• 开发一种稳定、高精度的数值方法，用于计算黎曼流形上的测地线长度。
• 在任意强 Sobolev 范数（$H^{s,p}$）下，以受控误差逼近 eikonal 方程的解。
• 通过 WKB 类型展开，实现对具有变系数的线性抛物型 PDE 的基本解的精确计算。
• 提供基于几何数据和 Sobolev 范数的逼近误差的严格估计。
• 支持在数学物理、金融学（如期权定价）和统计学（如最大似然估计）中的应用。

提出的方法

• 该方法使用 eikonal 方程 $ d^2 = \frac{1}{4} \sum_{ij} a_{ij}(x) \partial_{x_i} \partial_{x_j} $ 的全局解的规则多项式插值，其中 $ a_{ij} $ 是黎曼度量 $ g_{ij} $ 的逆矩阵。
• 该算法构建了对真实 $ d^2 $ 的 $ H^{s,p} $-逼近 $ q_{s,p}^2 $，确保对任意给定的 $ \epsilon > 0 $，有 $ \|d^2 - q_{s,p}^2\|_{s,p} \leq \epsilon $。
• 插值在有界区域 $ \Omega $ 上进行，误差界由 $ d^2 $ 及其导数的光滑性导出。
• 该方法利用了 $ d^2 $ 在紧致区域上为 $ C^\infty $ 的事实，从而可通过泰勒展开和插值实现高阶局部逼近。
• 误差估计使用多指标记号，并基于高阶导数的有界性，其中常数 $ C $ 依赖于 $ \Omega $ 上 $ \partial^\gamma d^2 $ 的最大值。
• 该方法被扩展用于计算抛物型 PDE 的 WKB 系数 $ c_k $，要求对 $ d^2 $ 实现 $ H^{s,p} $-逼近，且满足 $ s \geq 2k $。

实验结果

研究问题

• RQ1 如何在任意 Sobolev 范数下，以任意小的误差逼近黎曼流形上测地线的长度？
• RQ2 在有界区域上，将多项式插值应用于 eikonal 方程解时，其收敛行为和误差界如何？
• RQ3 对 $ d^2 $ 的逼近如何影响抛物型 PDE 基本解的 WKB 展开的精度？
• RQ4 为确保在 Sobolev 范数下实现高阶收敛，度量 $ g_{ij} $ 所需的最小正则性和几何条件是什么？
• RQ5 该算法能否在物理和金融学的实际应用中实现数值稳定和高效？

主要发现

• 该算法对任意 $ \epsilon > 0 $ 实现 $ \|d^2 - q_{s,p}^2\|_{s,p} \leq \epsilon $，其中 $ q_{s,p}^2 $ 是 eikonal 解在任意 Sobolev 范数下的多项式逼近。
• 第 $ m $ 阶插值的局部误差满足 $ |d^2(x,y) - d^2_{M(\beta_m)}(x,y)| \leq C P \frac{h^m}{m!} $，其中 $ P $ 为 $ m $ 阶多指标的数量。
• 导数误差满足 $ \partial^\beta d^2 - \partial^\beta d^2_{M(\beta_m)} \leq C \frac{h^{m-|\beta|}}{(m-|\beta|)!} $，表明所有导数均具有高阶收敛性。
• 常数 $ C $ 可用 $ \Omega $ 上 $ \partial^\gamma d^2 $ 的最大值有界表示，该值可基于 eikonal 方程和度量性质预先估计。
• 对于抛物型 PDE 的 WKB 近似至 $ k $ 阶，只需对 $ d^2 $ 实现满足 $ s \geq 2k $ 的 $ H^{s,p} $-逼近，即可确保基本解展开的高精度。
• 该方法可实现扩散过程中过渡密度的精确计算，直接应用于期权定价和随机过程中的统计推断。