Generalizing, optimizing, and inventing numerical algorithms for the fractional Fourier, Fresnel, and linear canonical transforms.

By use of matrix-based techniques it is shown how the space-bandwidth product (SBP) of a signal, as indicated by the location of the signal energy in the Wigner distribution function, can be tracked through any quadratic-phase optical system whose operation is described by the linear canonical transform. Then, applying the regular uniform sampling criteria imposed by the SBP and linking the criteria explicitly to a decomposition of the optical matrix of the system, it is shown how numerical algorithms (employing interpolation and decimation), which exhibit both invertibility and additivity, can be implemented. Algorithms appearing in the literature for a variety of transforms (Fresnel, fractional Fourier) are shown to be special cases of our general approach. The method is shown to allow the existing algorithms to be optimized and is also shown to permit the invention of many new algorithms.

[1]  J. Mayer,et al.  On the Quantum Correction for Thermodynamic Equilibrium , 1947 .

[2]  J. Tukey,et al.  An algorithm for the machine calculation of complex Fourier series , 1965 .

[3]  J. Goodman Introduction to Fourier optics , 1969 .

[4]  A. Papoulis Ambiguity function in Fourier optics , 1974 .

[5]  Mj Martin Bastiaans Wigner distribution function and its application to first-order optics , 1979 .

[6]  L.R. Rabiner,et al.  Interpolation and decimation of digital signals—A tutorial review , 1981, Proceedings of the IEEE.

[7]  Joseph Shamir,et al.  First-order optics—a canonical operator representation: lossless systems , 1982 .

[8]  A. Lohmann Image rotation, Wigner rotation, and the fractional Fourier transform , 1993 .

[9]  John T. Sheridan,et al.  Optical operations on wave functions as the Abelian subgroups of the special affine Fourier transformation. , 1994, Optics letters.

[10]  J. Sheridan,et al.  Generalization of the fractional Fourier transformation to an arbitrary linear lossless transformation an operator approach , 1994 .

[11]  Maciej Sypek,et al.  Light propagation in the Fresnel region. New numerical approach , 1995 .

[12]  R. Dorsch,et al.  Fractional-Fourier-transform calculation through the fast-Fourier-transform algorithm. , 1996, Applied optics.

[13]  Zeev Zalevsky,et al.  Space–bandwidth product of optical signals and systems , 1996 .

[14]  Gozde Bozdagi Akar,et al.  Digital computation of the fractional Fourier transform , 1996, IEEE Trans. Signal Process..

[15]  Zeev Zalevsky,et al.  Space–bandwidth product adaptation and its application to superresolution: examples , 1997 .

[16]  Mj Martin Bastiaans Application of the Wigner distribution function in optics , 1997 .

[17]  Zeev Zalevsky,et al.  Computation considerations and fast algorithms for calculating the diffraction integral , 1997 .

[18]  David Mendlovic,et al.  Space–bandwidth product adaptation and its application to superresolution: fundamentals , 1997 .

[19]  Francisco J. Marinho,et al.  Numerical calculation of fractional Fourier transforms with a single fast-Fourier-transform algorithm , 1998 .

[20]  Tomaso Erseghe,et al.  Unified fractional Fourier transform and sampling theorem , 1999, IEEE Trans. Signal Process..

[21]  Carlos Ferreira,et al.  Fast algorithms for free-space diffraction patterns calculation , 1999 .

[22]  J Gan,et al.  Fast algorithm for chirp transforms with zooming-in ability and its applications. , 2000, Journal of the Optical Society of America. A, Optics, image science, and vision.

[23]  D Mendlovic,et al.  Understanding superresolution in Wigner space. , 2000, Journal of the Optical Society of America. A, Optics, image science, and vision.

[24]  William T. Rhodes Numerical simulation of Fresnel-regime wave propagation: the light-tube model , 2001, SPIE Optics + Photonics.

[25]  Z. Zalevsky,et al.  The Fractional Fourier Transform: with Applications in Optics and Signal Processing , 2001 .

[26]  David Mas,et al.  Fast numerical calculation of Fresnel patterns in convergent systems , 2003 .

[27]  Bryan M Hennelly,et al.  Fast numerical algorithm for the linear canonical transform. , 2005, Journal of the Optical Society of America. A, Optics, image science, and vision.