[论文解读] A new and efficient method for the computation of Legendre coefficients

Enrico De Micheli, G. A. Viano|arXiv (Cornell University)|Jun 2, 2011

Digital Filter Design and Implementation参考文献 7被引用 4

一句话总结

本文提出了一种新颖的 O(N log N) 算法，通过将计算勒让德系数的问题转化为计算函数 f(x) 的阿贝尔型积分变换的傅里叶系数，实现了高效且高精度的计算。该方法利用勒让德多项式的狄利克雷–默菲积分表示，并应用一次快速傅里叶变换，显著优于标准的 O(N²) 数值积分方法。

ABSTRACT

An efficient procedure for the computation of the coefficients of Legendre expansions is here presented. We prove that the Legendre coefficients associated with a function f(x) can be represented as the Fourier coefficients of an Abel-type transform of f(x). The computation of N Legendre coefficients can then be performed in O(N log N) operations with a single Fast Fourier Transform of the Abel-type transform of f(x).

研究动机与目标

• 开发一种更快、更准确的勒让德系数计算算法，该算法在逼近理论、偏微分方程和谱方法中具有关键作用。
• 克服标准数值积分方法在高指标勒让德系数计算中因快速振荡被积函数导致的 O(N²) 计算复杂度。
• 通过阿贝尔型积分变换，揭示勒让德展开与傅里叶级数之间的联系，以实现高效计算。
• 在通过标准数值积分计算阿贝尔变换后，仅使用一次快速傅里叶变换，实现 O(N log N) 的时间复杂度。

提出的方法

• 利用勒让德多项式的狄利克雷–默菲积分表示，将勒让德系数与 f(x) 的阿贝尔型变换的傅里叶系数联系起来。
• 通过涉及 f(x) 和平方根奇点的积分定义一个变换函数 f̂(y)，并采用非线性变量变换来处理弱奇点。
• 系数 cₙ 表示为 cₙ = (2n+1) × aₙ，其中 aₙ 是 f̂(y) 在 [−π, π] 上的非负傅里叶系数。
• 通过一次快速傅里叶变换高效计算傅里叶系数 aₙ，从而将整体复杂度降低至 O(N log N)。
• 利用高斯-勒让德积分和非线性变量替换，高精度地计算阿贝尔型积分，以处理代数奇点。
• 使用 GNU 科学库（GSL）以双精度实现该算法，并通过高精度参考值进行验证。

实验结果

研究问题

• RQ1 能否使用谱方法实现比 O(N²) 更快的勒让德系数计算？
• RQ2 勒让德展开与傅里叶级数之间是否存在可被利用以实现更快计算的联系？
• RQ3 f(x) 的阿贝尔型积分变换能否被高精度且高效地计算，从而支持 O(N log N) 算法？
• RQ4 与标准数值积分和现有的 O(N(log N)²) 算法相比，该方法在精度和速度上表现如何？

主要发现

• 该方法实现了 O(N log N) 的计算复杂度，相较于标准的 O(N²) 数值积分方法有显著改进。
• 对于 f(x) = |x|³/²，计算得到的勒让德系数与参考值的误差在 10⁻¹⁰ 到 10⁻¹³ 之间，证实了其高精度。
• 该算法的性能提升比与 N/log N 成正比，当 N 较大时，加速效果愈发显著。
• 在准确计算阿贝尔变换后仅使用一次 FFT，实现了速度与精度的兼顾，且在多项测试函数（包括多项式、指数函数和有理函数）上均得到验证。
• 该方法具有通用性和可扩展性，未来可推广至多变量勒让德变换和盖根鲍尔多项式展开。