Pseudo-differential representation of the metaplectic transform and its application to fast algorithms
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.
- Research Organization:
- Princeton Plasma Physics Laboratory (PPPL), Princeton, NJ (United States)
- Sponsoring Organization:
- USDOE
- Grant/Contract Number:
- AC02-09CH11466
- OSTI ID:
- 1573502
- Alternate ID(s):
- OSTI ID: 1570879
- Journal Information:
- Journal of the Optical Society of America. A, Optics, Image Science, and Vision, Vol. 36, Issue 11; ISSN 1084-7529
- Publisher:
- Optical Society of America (OSA)Copyright Statement
- Country of Publication:
- United States
- Language:
- English
Web of Science
Similar Records
Restoring geometrical optics near caustics using sequenced metaplectic transforms
Extension of the Conley-Zehnder index, a product formula, and an application to the Weyl representation of metaplectic operators