[论文解读] Galerkin--Chebyshev approximation of Gaussian random fields on compact Riemannian manifolds

Annika Lang, Mike Pereira|arXiv (Cornell University)|Jul 6, 2021

Soil Geostatistics and Mapping参考文献 57被引用 4

一句话总结

该论文提出了一种Galerkin–Chebyshev逼近方法，通过拉普拉斯–贝尔特拉米算子的谱分解和切比雪夫级数的多项式逼近，高效地在紧致黎曼流形上采样高斯随机场。该方法在网格尺寸上实现线性计算和存储成本，并在Galerkin和Galerkin–Chebyshev逼近中均达到最优收敛速率，数值验证在球面和双曲面上均有效。

ABSTRACT

A new numerical approximation method for a class of Gaussian random fields on compact connected oriented Riemannian manifolds is introduced. This class of random fields is characterized by the Laplace--Beltrami operator on the manifold. A Galerkin approximation is combined with a polynomial approximation using Chebyshev series. This so-called Galerkin--Chebyshev approximation scheme yields efficient and generic sampling algorithms for Gaussian random fields on manifolds. Strong and weak orders of convergence for the Galerkin approximation and strong convergence orders for the Galerkin--Chebyshev approximation are shown and confirmed through numerical experiments.

研究动机与目标

• 开发一种通用、高效且可扩展的方法，用于在紧致黎曼流形（不仅限于球面）上采样高斯随机场。
• 将Galerkin逼近与切比雪夫多项式展开相结合，以实现有限维表示中相关随机系数的快速计算。
• 建立Galerkin逼近的强收敛和弱收敛速率，以及Galerkin–Chebyshev方案的强收敛速率。
• 通过蒙特卡洛协方差估计，在球面和双曲面上数值验证理论收敛速率。
• 在二维流形上，实现对Matérn型场ε-精度采样的计算复杂度𝒪(ε⁻²⁄ρ)。

提出的方法

• 该方法通过拉普拉斯–贝尔特拉米算子−Δℳ的函数展开，在紧致黎曼流形ℳ上构造高斯随机场𝒁。
• 通过Galerkin方法将𝒁投影到有限维空间Vₙ（如线性有限元空间）上，形成Galerkin逼近ℤₙ。
• 推导出ℤₙ在基表示下系数向量的协方差矩阵的闭式表达，从而实现相关系数的显式采样。
• 对谱函数应用切比雪夫多项式逼近，通过矩阵-向量乘积加速场实现的计算。
• 该方法利用质量矩阵对角化和多项式求根技术，在保持收敛速率的同时维持计算效率。
• 理论收敛性通过均方误差和协方差误差进行分析，数值实验验证了预测的收敛速率。

实验结果

研究问题

• RQ1 能否为一般紧致黎曼流形（不仅限于球面）上的高斯随机场开发一种通用且高效的采样算法？
• RQ2 在均方误差和协方差误差方面，此类场的Galerkin逼近收敛速率如何？
• RQ3 对谱函数的切比雪夫逼近如何影响场采样的精度和计算成本？
• RQ4 所提出的Galerkin–Chebyshev方案能否在自由度数量上实现线性扩展，同时保持最优收敛性？
• RQ5 理论收敛速率在非球面流形（如双曲面）上是否在实践中成立？

主要发现

• 在标准假设下，Galerkin逼近ℤₙ以𝒪(n⁻¹)阶的均方误差收敛于真实场𝒁。
• ℤₙ与𝒁之间的协方差误差收敛速率为2αβ − 1 = ν，与Whittle–Matérn场的理论预测一致。
• 在球面上，Galerkin–Chebyshev逼近对具有光滑参数ν的Matérn场实现了𝒪(n⁻ν)阶的均方根误差收敛速率。
• 对于二维流形，以精度ε > 0采样场的计算成本为𝒪(ε⁻²⁄ρ)，其中ρ = ν对应Matérn场。
• 在球面和双曲面上的数值实验验证了理论收敛速率：当ν = 1时，协方差误差以𝒪(n⁻¹)衰减。
• 为达到|cₖ/cₘₐₓ| < 10⁻¹²的精度，所需多项式逼近阶数K保持较低，即使在高K下，每次采样的计算时间也低于0.5秒。