Pseudo-differential representation of the metaplectic transform and its application to fast algorithms
Abstract
The metaplectic transform (MT), also known as the linear canonical transform, is a unitary integral mapping that is widely used in signal processing and can be viewed as a generalization of the Fourier transform. For a given function $ψ$ on an $N$-dimensional continuous space q, the MT of $ψ$ is parameterized by a rotation (or more generally, a linear symplectic transformation) of the 2$N$-dimensional phase space (q,p), where p is the wavevector space dual to q. In this study, we derive a pseudo-differential form of the MT. For small-angle rotations, or near-identity transformations of the phase space, it readily yields asymptotic differential representations of the MT, which are easy to compute numerically. Rotations by larger angles are implemented as successive applications of $$K \gg 1$$ small-angle MTs. The algorithm complexity scales as $$O(KN^3N_p)$$ where $$N_p$$ is the number of grid points. Here, we present a numerical implementation of this algorithm and discuss how to mitigate the associated numerical instabilities.
- Authors:
-
- Princeton Univ., NJ (United States). Dept. of Astrophysical Sciences
- Princeton Univ., NJ (United States). Dept. of Astrophysical Sciences; Princeton Plasma Physics Lab. (PPPL), Princeton, NJ (United States)
- Publication Date:
- Research Org.:
- Princeton Plasma Physics Lab. (PPPL), Princeton, NJ (United States)
- Sponsoring Org.:
- USDOE
- OSTI Identifier:
- 1573502
- Alternate Identifier(s):
- OSTI ID: 1570879
- Grant/Contract Number:
- AC02-09CH11466
- Resource Type:
- Accepted Manuscript
- Journal Name:
- Journal of the Optical Society of America. A, Optics, Image Science, and Vision
- Additional Journal Information:
- Journal Volume: 36; Journal Issue: 11; Journal ID: ISSN 1084-7529
- Publisher:
- Optical Society of America (OSA)
- Country of Publication:
- United States
- Language:
- English
- Subject:
- 79 ASTRONOMY AND ASTROPHYSICS
Citation Formats
Lopez, N. A., and Dodin, I. Y. Pseudo-differential representation of the metaplectic transform and its application to fast algorithms. United States: N. p., 2019.
Web. doi:10.1364/JOSAA.36.001846.
Lopez, N. A., & Dodin, I. Y. Pseudo-differential representation of the metaplectic transform and its application to fast algorithms. United States. doi:10.1364/JOSAA.36.001846.
Lopez, N. A., and Dodin, I. Y. Tue .
"Pseudo-differential representation of the metaplectic transform and its application to fast algorithms". United States. doi:10.1364/JOSAA.36.001846. https://www.osti.gov/servlets/purl/1573502.
@article{osti_1573502,
title = {Pseudo-differential representation of the metaplectic transform and its application to fast algorithms},
author = {Lopez, N. A. and Dodin, I. Y.},
abstractNote = {The metaplectic transform (MT), also known as the linear canonical transform, is a unitary integral mapping that is widely used in signal processing and can be viewed as a generalization of the Fourier transform. For a given function $ψ$ on an $N$-dimensional continuous space q, the MT of $ψ$ is parameterized by a rotation (or more generally, a linear symplectic transformation) of the 2$N$-dimensional phase space (q,p), where p is the wavevector space dual to q. In this study, we derive a pseudo-differential form of the MT. For small-angle rotations, or near-identity transformations of the phase space, it readily yields asymptotic differential representations of the MT, which are easy to compute numerically. Rotations by larger angles are implemented as successive applications of $K \gg 1$ small-angle MTs. The algorithm complexity scales as $O(KN^3N_p)$ where $N_p$ is the number of grid points. Here, we present a numerical implementation of this algorithm and discuss how to mitigate the associated numerical instabilities.},
doi = {10.1364/JOSAA.36.001846},
journal = {Journal of the Optical Society of America. A, Optics, Image Science, and Vision},
number = 11,
volume = 36,
place = {United States},
year = {2019},
month = {10}
}
Figures / Tables:

Works referenced in this record:
Fast and accurate computation of two-dimensional non-separable quadratic-phase integrals
journal, January 2010
- Koç, Aykut; Ozaktas, Haldun M.; Hesselink, Lambertus
- Journal of the Optical Society of America A, Vol. 27, Issue 6
Lectures on nonlinear orbit dynamics
conference, January 1982
- Dragt, Alex J.
- AIP Conference Proceedings
Segmented fast linear canonical transform
journal, January 2018
- Sun, Yan-Nan; Li, Bing-Zhao
- Journal of the Optical Society of America A, Vol. 35, Issue 8
Quasioptical modeling of wave beams with and without mode conversion. I. Basic theory
journal, July 2019
- Dodin, I. Y.; Ruiz, D. E.; Yanagihara, K.
- Physics of Plasmas, Vol. 26, Issue 7
Fast numerical algorithm for the linear canonical transform
journal, January 2005
- Hennelly, Bryan M.; Sheridan, John T.
- Journal of the Optical Society of America A, Vol. 22, Issue 5
Synchrotron radiation representation in phase space
journal, May 2012
- Bazarov, Ivan V.
- Physical Review Special Topics - Accelerators and Beams, Vol. 15, Issue 5
Simulating first order optical systems—algorithms for and composition of discrete linear canonical transforms
journal, December 2017
- Healy, John J.
- Journal of Optics, Vol. 20, Issue 1
A non-negative Wigner-type distribution
journal, January 1976
- Cartwright, N. D.
- Physica A: Statistical Mechanics and its Applications, Vol. 83, Issue 1
Digital computation of the fractional Fourier transform
journal, January 1996
- Ozaktas, H. M.; Arikan, O.; Kutay, M. A.
- IEEE Transactions on Signal Processing, Vol. 44, Issue 9
The Symplectic Group and Classical Mechanics
journal, June 2005
- Dragt, Alex J.
- Annals of the New York Academy of Sciences, Vol. 1045, Issue 1
Noninterferometric phase retrieval using a fractional Fourier system
journal, December 2007
- Gopinathan, Unnikrishnan; Situ, Guohai; Naughton, Thomas J.
- Journal of the Optical Society of America A, Vol. 25, Issue 1
Linear Canonical Transformations and Their Unitary Representations
journal, August 1971
- Moshinsky, M.; Quesne, C.
- Journal of Mathematical Physics, Vol. 12, Issue 8
On the Quantum Correction For Thermodynamic Equilibrium
journal, June 1932
- Wigner, E.
- Physical Review, Vol. 40, Issue 5
Improved implementation algorithms of the two-dimensional nonseparable linear canonical transform
journal, January 2012
- Ding, Jian-Jiun; Pei, Soo-Chang; Liu, Chun-Lin
- Journal of the Optical Society of America A, Vol. 29, Issue 8
Fast linear canonical transforms
journal, December 2009
- Healy, John J.; Sheridan, John T.
- Journal of the Optical Society of America A, Vol. 27, Issue 1
Chirp slope keying for underwater communications
conference, May 2005
- Kaminsky, Edit J.; Simanjuntak, Lastri
- Defense and Security, SPIE Proceedings
Phase-space tomography with a programmable Radon–Wigner display
journal, January 2011
- Cámara, Alejandro; Alieva, Tatiana; Rodrigo, José A.
- Optics Letters, Vol. 36, Issue 13
Operator methods in physical optics
journal, January 1981
- Stoler, David
- Journal of the Optical Society of America, Vol. 71, Issue 3
Image rotation, Wigner rotation, and the fractional Fourier transform
journal, January 1993
- Lohmann, Adolf W.
- Journal of the Optical Society of America A, Vol. 10, Issue 10
Metaplectic formulation of linear mode conversion
journal, September 1993
- Tracy, E. R.; Kaufman, Allan N.
- Physical Review E, Vol. 48, Issue 3
Some properties of a non-negative quantum-mechanical distribution function
journal, September 1981
- O'Connell, R. F.; Wigner, E. P.
- Physics Letters A, Vol. 85, Issue 3
Lens-System Diffraction Integral Written in Terms of Matrix Optics*
journal, January 1970
- Collins, Stuart A.
- Journal of the Optical Society of America, Vol. 60, Issue 9
The semiclassical evolution of wave packets
journal, May 1986
- Littlejohn, R.
- Physics Reports, Vol. 138, Issue 4-5
Two-dimensional nonseparable discrete linear canonical transform based on CM-CC-CM-CC decomposition
journal, January 2016
- Pei, Soo-Chang; Huang, Shih-Gu
- Journal of the Optical Society of America A, Vol. 33, Issue 2
Local fields for asymptotic matching in multidimensional mode conversion
journal, August 2007
- Tracy, E. R.; Kaufman, A. N.; Jaun, A.
- Physics of Plasmas, Vol. 14, Issue 8
Symplectically Invariant WKB Wave Functions
journal, April 1985
- Littlejohn, Robert G.
- Physical Review Letters, Vol. 54, Issue 16
Precision measurement beyond the shot-noise limit
journal, July 1987
- Xiao, Min; Wu, Ling-An; Kimble, H. J.
- Physical Review Letters, Vol. 59, Issue 3
Smoothing and Differentiation of Data by Simplified Least Squares Procedures.
journal, July 1964
- Savitzky, Abraham.; Golay, M. J. E.
- Analytical Chemistry, Vol. 36, Issue 8
Figures / Tables found in this record: