[论文解读] Detailed analysis of prolate quadratures and interpolation formulas

Andrei Osipov, Vladimir Rokhlin|arXiv (Cornell University)|Aug 23, 2012

Electromagnetic Scattering and Analysis参考文献 20被引用 6

一句话总结

本文提出了计算抛物球面波函数（PSWFs）、特征值及带限函数求积规则的高效数值算法。其计算成本从与带宽 $c$ 的平方或立方关系降低至线性 $O(c)$，从而实现了在 $[-1,1]$ 上利用基于 PSWF 的求积规则对带限函数进行快速且高精度的积分与插值。

ABSTRACT

As demonstrated by Slepian et. al. in a sequence of classical papers, prolate spheroidal wave functions (PSWFs) provide a natural and efficient tool for computing with bandlimited functions defined on an interval. As a result, PSWFs are becoming increasing popular in various areas in which such function occur - this includes physics (e.g. wave phenomena, fluid dynamics), engineering (e.g. signal processing, filter design), etc. To use PSWFs as a computational tool, one needs fast and accurate numerical algorithms for the evaluation of PSWFs and related quantities, as well as for the construction of quadratures, interpolation formulas, etc. Even though, for the last half a century, substantial progress has been made in design of such algorithms, the complexity of many of the existing algorithms, however, is at least quadratic in the band limit $c$. For example, the evaluation of the $n$th eigenvalue of the prolate integral operator requires at least $O(c^2)$ operations. Therefore, while the existing algorithms are quite satisfactory for moderate values of $c$ (e.g. $c \leq 10^3$), they tend to be relatively slow when $c$ is large (e.g. $c \geq 10^4$). In this paper, we describe several numerical algorithms for the evaluation of PSWFs and related quantities, and design a class of PSWF-based quadratures for the integration of bandlimited functions. Also, we perform detailed analysis of the related properties of PSWFs. While the analysis is somewhat involved, the resulting numerical algorithms are quite simple and efficient in practice. For example, the evaluation of the $n$th eigenvalue of the prolate integral operator requires $O(n+c)$ operations; also, the construction of related accurate quadrature rules requires $O(c)$ operations. Our results are illustrated via several numerical experiments.

研究动机与目标

• 开发用于计算抛物球面波函数（PSWFs）及相关量（如特征值和求积权重）的快速且精确的算法。
• 设计一类基于 PSWF 的求积方法，使得对带宽为 $c$ 的前 $n$ 个 PSWF 进行积分的误差 $\leq \varepsilon$。
• 将现有算法的计算复杂度从 $O(c^2)$ 或 $O(c^3)$ 降低至求积构造的 $O(c)$ 和特征值计算的 $O(n + c)$。
• 对 PSWF 的性质进行详细的理论与数值分析，包括节点结构、特征值衰减规律以及勒让德展开系数。
• 通过大量数值实验验证算法的性能，结果表明其具有机器精度级别的准确性，并且在 $c$ 上呈现线性可扩展性。

提出的方法

• 本文通过渐近展开与递推关系，推导出一种新算法，可在 $O(n + c)$ 次操作内计算第 $n$ 个 PSWF 特征值。
• 将求积节点构造为第 $n$ 个 PSWF 的根，并通过一种新公式计算权重，该公式涉及导数 $\psi_n'(t_j)$ 和第二类勒让德函数。
• 利用包含勒让德系数 $\alpha_k^{(n)}$ 的级数展开计算权重，截断至 $O(n)$ 项以提高效率。
• 该方法利用了 PSWF 的正交性及其渐近行为，并结合了对特殊函数（如 $Q_k(t)$）的快速求值。
• 采用双精度算术实现稳健的数值方案，并利用对称性减少 $j > (n+1)/2$ 时的计算量。
• 理论结果通过数值实验得到支持，实验表明算法在 $c$ 上具有线性可扩展性，且精度达到机器精度。

实验结果

研究问题

• RQ1 能否将计算 PSWF 特征值的计算成本降低至 $O(c^2)$ 以下？若能，应采用何种算法方法？
• RQ2 为使带宽为 $c$ 的前 $n$ 个 PSWF 的积分误差 $\leq \varepsilon$，所需的最少求积节点数 $n$ 是多少？
• RQ3 如何仅利用 PSWF 的根与导数，高效且准确地计算求积权重？
• RQ4 勒让德函数展开在推导稳定且快速的权重计算公式中起什么作用？
• RQ5 对于大 $c$（例如 $c \geq 10^4$）的情形，所提出的算法在保持精度与效率方面能达到何种程度？

主要发现

• 第 $n$ 个 PSWF 特征值的计算可在 $O(n + c)$ 次操作内完成，相比以往方法的 $O(c^2)$ 复杂度有显著提升。
• 构造带宽为 $c$ 的带限函数的高精度求积规则仅需 $O(c)$ 次操作，相比以往的 $O(c^3)$ 成本大幅降低。
• 数值实验表明，通过新公式计算的求积权重与理论预期一致，残差量级为 $10^{-15}$，达到机器精度。
• 求积权重严格为正，并向中心单调递增，与定理 73 和注记 13 的理论预期一致。
• 计算权重与解析公式 (673) 之间的差异在机器精度范围内为零，验证了权重计算算法的正确性。
• 当 $c = 40$ 且 $n = 41$ 时，算法可准确计算所有权重，且备选权重公式 (674) 的误差量级为 $|\lambda_n| \approx 7 \times 10^{-9}$，与注记 14 的预测一致。