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 which 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 2N-dimensional phase space (q,p), where p is the wavevector space dual to q. Here, 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 1 small-angle MTs. The algorithm complexity scales as O(K N3 Np), where Np is the number of grid points. We present a numerical implementation of this algorithm and discuss how to mitigate the associated numerical instabilities.

0

Turn this paper into a lesson

ArcXiv compiles a structured reading guide from this paper's metadata: plain-English importance, contributions, prerequisite concepts, which sections to read first, flashcards, and a quiz. Grounded in the abstract, never invented.

Discussion (0)

Sign in to join the discussion.

Loading comments…