Prolate spheroidal wave functions

The prolate spheroidal wave functions are a set of functions derived by timelimiting and lowpassing, and a second timelimit operation. Let $Q\_T$ denote the time truncation operator, such that $x=Q\_T\; x$ iff x is timelimited within $[-T/2;T/2]$. Similarly, let $P\_W$ denote an ideal low-pass filtering operator, such that $x=P\_W\; x$ iff x is bandlimited within $[-W;W]$. The operator $Q\_T\; P\_W\; Q\_T$ turns out to be linear, bounded and self-adjoint. For $n=1,2,ldots$ we denote with $psi\_n$ the n-th eigenfunction, defined as

: $Q\_T\; P\_W\; Q\_T\; psi\_n=lambda\_npsi\_n,$

where $\{lambda\_n\}\_n$ are the associated eigenvalues. The timelimited functions $\{psi\_n\}\_n$ are the Prolate Spheroidal Wave Functions (PSWFs).

These functions are also encountered in a different context. When solving the Helmholtz equation,$Delta\; Phi\; +\; k^2\; Phi$, by the method of separation of variables in prolate spheroidal coordinates, $(xi,eta,phi)$, with:

:$x=f\; xi\; eta,$

:$y=f\; sqrt\{(xi^2-1)(1-eta^2)\}\; cos\; phi,$

:$z=f\; sqrt\{(xi^2-1)(1-eta^2)\}\; sin\; phi,$

:$|xi|>1$ and . the solution $Phi(xi,eta,phi)$ can be writtenas the product of a radial spheroidal wavefunction $R\_\{mn\}(c,xi)$ and an angular spheroidal wavefunction $S\_\{mn\}(c,eta)$ by $e^\{i\; m\; phi\}$ with $c=fk/2$.

The radial wavefunction $R\_\{mn\}(c,xi)$ satisfies the linear ordinary differential equation:

:$(xi^2\; -1)\; frac\{d^2\; R\_\{mn\}(c,xi)\}\{d^2\; xi\}\; +\; 2xi\; frac\{d\; R\_\{mn\}(c,xi)\}\{d\; xi\}\; -left(lambda\_\{mn\}(c)\; -c^2\; xi^2\; +frac\{m^2\}\{xi^2-1\}\; ight)\; \{R\_\{mn\}(c,xi)\}\; =\; 0$

The eigenvalue $lambda\_\{mn\}(c)$ of this Sturm-Liouville differential equation is fixed by the requirement that $\{R\_\{mn\}(c,xi)\}$ must be finite for $|xi|\; o\; 1\_+$.

The angular wavefunction satisfies the differential equation:

:$(eta^2\; -1)\; frac\{d^2\; S\_\{mn\}(c,eta)\}\{d^2\; eta\}\; +\; 2eta\; frac\{d\; S\_\{mn\}(c,eta)\}\{d\; eta\}\; -left(lambda\_\{mn\}(c)\; -c^2\; eta^2\; +frac\{m^2\}\{eta^2-1\}\; ight)\; \{S\_\{mn\}(c,eta)\}\; =\; 0$

It is the same differential equation as in the case of the radial wavefunction. However, the range of the variable is different (in the radial wavefunction, $|xi|>1$) in the angular wavefunction ).

For $c=0$ these two differential equations reduce to the equations satisfied by the Legendre functions. For $c\; e\; 0$, the angular spheroidal wavefunctions can be expanded as a series of Legendre functions.

Let us note that if one writes $S\_\{mn\}(c,eta)=(1-eta^2)^\{m/2\}\; Y\_\{mn\}(c,eta)$, the function $Y\_\{mn\}(c,eta)$ satisfies the following linear ordinary differential equation:

$(1-eta^2)\; frac\{d^2\; Y\_\{mn\}(c,eta)\}\{d^2\; eta\}\; -2\; (m+1)\; eta\; frac\{d\; Y\_\{mn\}(c,eta)\}\{d\; eta\}\; -\; +left(c^2\; eta^2\; +m(m+1)-lambda\_\{mn\}(c)\; ight)\; \{Y\_\{mn\}(c,eta)\}\; =\; 0,$

which is known as the spheroidal wave equation. This auxiliary equation is used for instance by Stratton in his 1935 article.

There are different normalization schemes for spheroidal functions. A table of the different schemes can be found in Abramowitz and Stegun p. 758. Abramowitz and Stegun (and the present article) follow the notation of Flammer.

In the case of oblate spheroidal coordinates the solution of the Helmholtz equation yields oblate spheroidal wavefunctions.

Originally, the spheroidal wave functions were introduced by C. Niven in 1880 when studying the conduction of heat in an ellipsoid of revolution, which lead to a Helmholtz equation in spheroidal coordinates.

References

* W. J. Thomson [http://www.ece.nus.edu.sg/stfpage/elelilw/Software/00764220.pdf Spheroidal Wave functions] Computing in Science & Engineering p. 84, May-June 1999.
* I. Daubechies, "Ten Lectures on Wavelets"
* C. Niven " [http://gallica.bnf.fr/notice?N=FRBNF37571969 On the Conduction of Heat in Ellipsoids of Revolution.] " Philosophical transactions of the Royal Society of London, 171 p. 117 (1880)
* J. A. Stratton " [http://www.pnas.org/cgi/reprint/21/1/51?maxtoshow=&HITS=10&hits=10&RESULTFORMAT=1&title=spheroidal&andorexacttitle=and&andorexacttitleabs=and&andorexactfulltext=and&searchid=1&FIRSTINDEX=0&sortspec=relevance&resourcetype=HWCIT Spheroidal functions] " Proceedings of the National Academy of Sciences (USA) 21, 51 (1935)
* C. Flammer "Spheroidal Wave Functions." Stanford, CA: Stanford University Press, 1957.
* J. Meixner and F. W. Schäfke, "Mathieusche Funktionen und Sphäroidfunktionen." Berlin: Springer-Verlag, 1954.
* J. A. Stratton , P. M. Morse, J. L. Chu, J. D. C. Little, and F. J. Corbató, "Spheroidal Wave Functions." New York: Wiley, 1956.
* M. Abramowitz and I. Stegun "Handbook of Mathematical Functions" [http://www.math.sfu.ca/~cbm/aands/page_751.htm pp. 751-759] (Dover, New York, 1972)
* F. Sleator " [http://hdl.handle.net/2027.42/21231 Studies in Radar Cross-Sections -- XLIX. Diffraction and scattering by regular bodies III: the prolate spheroid] " (1964)
* H. E. Hunter " [http://hdl.handle.net/2027.42/5662 Tables of prolate spheroidal functions for m=0: Volume I.] " (1965)
* H. E. Hunter " [http://hdl.handle.net/2027.42/5663 Tables of prolate spheroidal functions for m=0 : Volume II.] " (1965)
* P. E. Falloon (2002) [http://ftp.physics.uwa.edu.au/pub/Theses/2002/Falloon/Spheroidal/ Thesis on numerical computation of spheroidal functions] University of Western Australia

External links

* MathWorld [http://mathworld.wolfram.com/SpheroidalWaveFunction.html Spheroidal Wave functions]
* MathWorld [http://mathworld.wolfram.com/ProlateSpheroidalWaveFunction.html Prolate Spheroidal Wave Function]
* MathWorld [http://mathworld.wolfram.com/OblateSpheroidalWaveFunction.html Oblate Spheroidal Wave function]