Appl. Comput. Harmon. Anal. 48 (2020) 539–569
Contents lists available at ScienceDirect
Applied and Computational Harmonic Analysis
www.elsevier.com/locate/acha
Ball prolate spheroidal wave functions in arbitrary dimensions
Jing Zhang a,1, Huiyuan Li b,∗,2, Li-Lian Wang c,3, Zhimin Zhang d,e,4 a School of Mathematics and Statistics & Hubei Key Laboratory of Mathematical Sciences, Central China Normal University, Wuhan 430079, China b State Key Laboratory of Computer Science/Laboratory of Parallel Computing, Institute of Software, Chinese Academy of Sciences, Beijing 100190, China c Division of Mathematical Sciences, School of Physical and Mathematical Sciences, Nanyang Technological University, 637371, Singapore d Beijing Computational Science Research Center, Beijing 100193, China e Department of Mathematics, Wayne State University, MI 48202, USA a r t i c l e i n f o a b s t r a c t
Article history: In this paper, we introduce the prolate spheroidal wave functions (PSWFs) of real Received 24 February 2018 order α > −1on the unit ball in arbitrary dimension, termed as ball PSWFs. Received in revised form 1 August They are eigenfunctions of both an integral operator, and a Sturm–Liouville 2018 differential operator. Different from existing works on multi-dimensional PSWFs, Accepted 3 August 2018 the ball PSWFs are defined as a generalization of orthogonal bal l polynomials in Available online 7 August 2018 Communicated by Gregory Beylkin primitive variables with a tuning parameter c > 0, through a “perturbation” of the Sturm–Liouville equation of the ball polynomials. From this perspective, we MSC: can explore some interesting intrinsic connections between the ball PSWFs and the 42B37 finite Fourier and Hankel transforms. We provide an efficient and accurate algorithm 33E30 for computing the ball PSWFs and the associated eigenvalues, and present various 33C47 numerical results to illustrate the efficiency of the method. Under this uniform 42C05 framework, we can recover the existing PSWFs by suitable variable substitutions. 65D20 © 2018 Elsevier Inc. All rights reserved. 41A10
Keywords: Generalized prolate spheroidal wave functions Arbitrary unit ball Sturm–Liouville differential equation Finite Fourier transform Bouwkamp spectral-algorithm
* Corresponding author. E-mail address: [email protected] (H. Li). 1 The work of this author is partially supported by the National Natural Science Foundation of China (NSFC 11671166). 2 The research of this author is partially supported by the National Natural Science Foundation of China (NSFC 11471312 and NSFC 91430216). 3 The research of this author is partially supported by Singapore MOE AcRF Tier 1 Grant (RG 27/15) and Singapore MOE AcRF Tier 2 Grant (MOE2017-T2-2-144). 4 The research of this author is supported in part by the National Natural Science Foundation of China (NSFC 11471031, NSFC 91430216 and NSAF U1530401) and the U.S. National Science Foundation (DMS-1419040). https://doi.org/10.1016/j.acha.2018.08.001 1063-5203/© 2018 Elsevier Inc. All rights reserved. 540 J. Zhang et al. / Appl. Comput. Harmon. Anal. 48 (2020) 539–569
1. Introduction
The PSWFs are a family of orthogonal bandlimited functions, originated from the investigation of time- frequency concentration problem in the 1960s (cf. [26,27,39,38]). In the study of time-frequency concentration ∞ problem, Slepian was the first to note that the PSWFs, denoted by ψn(x; c) n=0, are the eigenfunctions of an integral operator related to the finite Fourier transform: