A numerical study of the Legendre-Galerkin method for the evaluation of the prolate spheroidal wave functions

作者：Schmutzhard Sebastian^*; Hrycak Tomasz; Feichtinger Hans G

来源：Numerical Algorithms, 2015, 68(4): 691-710.

DOI：10.1007/s11075-014-9867-3

• Prolate spheroidal wave functions
• Legendre-Galerkin method
• Discretization
• Condition number
• Accuracy
• Complexity

We conduct a numerical study of the Legendre-Galerkin method for the evaluation of the prolate spheroidal wave functions (PSWFs) viewed as the eigenfunctions of the prolate differential operator with boundary conditions of continuity. Our experiments indicate that the minimal dimension N of the Legendre-Galerkin matrix for the evaluation of the nth prolate with precision epsilon is , as , where c > 0 is the bandwidth parameter. The behavior of N, when either c or n is held constant, is also examined. As a consequence, we obtain an upper bound on the complexity of the evaluation of the prolates. We also study the condition number of the approximate Legendre coefficients, computed as an eigenvector of the Legendre-Galerkin matrix. We observe experimentally that for fixed precision epsilon, an error estimate based on this condition number is as . We conclude that the Legendre-Galerkin method is accurate for fairly large values of n and c.

• 出版日期2015-4