An Alternative Derivation of the Method of Stationary Phase for Multivariate Signals

The method of stationary phase (msp) is a commonly established technique to asymptotically determine integrals (such as Fourier integrals) with strongly oscillating integrands. Basically, the method states that an integral over a strongly oscillating integrand only delivers contributions at those (“stationary”) points or regions, where the integrand does not change, i.e., where the phase derivative of its oscillation vanishes. In the engineering literature, various approaches for the derivation exist for univariate (1-D) signals, which, however, find their limitations for multivariate (multidimensional) signals, which cannot be factorized. On the other hand, from the field of applied mathematics, there are derivations, formulating and solving the problem in a rigorous and sophisticated way. This paper tries to close what the author interprets as a gap between these two worlds, aiming to substitute intuitive insights by some basic Dirac impulse modeling concepts, which are familiar to signal processing engineers. We first start with the general multivariate (multidimensional) case, introducing vectorial time and frequency representations. Here, we replace 1-D Taylor-series-based quadratic phase signals by multivariate signals with phase histories, being given by more general quadratic forms. The exponential, containing the quadratic form, is then related to an n-variate Dirac impulse, which provides the sifting property needed to establish the basic msp property that the integration over the entire integration range yields only one value of the integrand, evaluated in the point of stationary phase. This sifting capacity in connection with (time limiting) window functions in the integrand migrates time limitations of the signals to frequency limitations of the spectra-a multiplicative time window is transformed in a multiplicative frequency window of the same form, rather than into a convolution of the Fourier transforms. In this respect, this paper provides a closed-form comprehensive treatment of the subject.

[1]  Otmar Loffeld,et al.  Models and useful relations for bistatic SAR processing , 2004, IEEE Transactions on Geoscience and Remote Sensing.

[2]  Morris Kline,et al.  Asymptotic expansion of multiple integrals and the method of stationary phase , 2011 .

[3]  A. Papoulis Systems and transforms with applications in optics , 1981 .

[4]  A Large Time-Bandwidth Product Pulse-Compression Technique , 1962, IRE Transactions on Military Electronics.

[5]  Riccardo Lanari,et al.  Synthetic Aperture Radar Processing , 1999 .

[6]  N. G. Van Kampen,et al.  The method of stationary phase and the method of Fresnel zones , 1958 .

[7]  Stewart A. Greenhalgh,et al.  Modeling and Migration of 2-D Georadar Data: A Stationary Phase Approach , 2006, IEEE Transactions on Geoscience and Remote Sensing.

[8]  Roderick Wong,et al.  Asymptotic approximations of integrals , 1989, Classics in applied mathematics.

[9]  Giorgio Franceschetti,et al.  A SAR processor based on two-dimensional FFT codes , 1990 .

[10]  R.N. Bracewell,et al.  Signal analysis , 1978, Proceedings of the IEEE.

[11]  Otmar Loffeld,et al.  Focusing of General Bistatic SAR Configuration Data With 2-D Inverse Scaled FFT , 2006, IEEE Transactions on Geoscience and Remote Sensing.

[12]  Graeme L. James,et al.  Geometrical Theory of Diffraction for Electromagnetic Waves , 1980 .

[13]  Joachim H. G. Ender,et al.  Processing the Azimuth-Variant Bistatic SAR Data by Using Monostatic Imaging Algorithms Based on Two-Dimensional Principle of Stationary Phase , 2011, IEEE Transactions on Geoscience and Remote Sensing.

[14]  Roderick Wong,et al.  On a method of asymptotic evaluation of multiple integrals , 1981 .

[15]  N. G. Van Kampen,et al.  An asymptotic treatment of diffraction problems. II , 1949 .

[16]  A. Sieber,et al.  Fast algorithm for a near-field synthetic aperture radar processor , 1994 .

[17]  T. Wolff The stationary phase method , 2003 .

[18]  Masoud Alimohammadi Quantum chains with GLq(2) symmetry , 1994 .

[19]  S. Lang,et al.  An Introduction to Fourier Analysis and Generalised Functions , 1959 .

[20]  Alberto Moreira,et al.  Spotlight SAR data processing using the frequency scaling algorithm , 1999, IEEE Trans. Geosci. Remote. Sens..

[21]  R. Mager,et al.  An examination of the limited aperture problem of physical optics inverse scattering , 1978 .

[22]  Norman Bleistein,et al.  Mathematical Methods for Wave Phenomena , 1984 .