Saved in:
Bibliographic Details
Main Authors: Lopez, N. A., Dodin, I. Y.
Format: Preprint
Published: 2020
Subjects:
Online Access:https://arxiv.org/abs/2012.06028
Tags: Add Tag
No Tags, Be the first to tag this record!
_version_ 1866912317662298112
author Lopez, N. A.
Dodin, I. Y.
author_facet Lopez, N. A.
Dodin, I. Y.
contents The metaplectic transform (MT), a generalization of the Fourier transform sometimes called the linear canonical transform, is a tool used ubiquitously in modern optics, for example, when calculating the transformations of light beams in paraxial optical systems. The MT is also an essential ingredient of the geometrical-optics modeling of caustics that was recently proposed by the authors. In particular, this application relies on the near-identity MT (NIMT); however, the NIMT approximation used so far is not exactly unitary and leads to numerical instability. Here, we develop a discrete MT that is exactly unitary, and approximate it to obtain a discrete NIMT that is also unitary and can be computed in linear time. We prove that the discrete NIMT converges to the discrete MT when iterated, thereby allowing the NIMT to compute MTs that are not necessarily near-identity. We then demonstrate the new algorithms with a series of examples.
format Preprint
id arxiv_https___arxiv_org_abs_2012_06028
institution arXiv
publishDate 2020
record_format arxiv
spellingShingle Exactly unitary discrete representations of the metaplectic transform for linear-time algorithms
Lopez, N. A.
Dodin, I. Y.
Optics
Mathematical Physics
The metaplectic transform (MT), a generalization of the Fourier transform sometimes called the linear canonical transform, is a tool used ubiquitously in modern optics, for example, when calculating the transformations of light beams in paraxial optical systems. The MT is also an essential ingredient of the geometrical-optics modeling of caustics that was recently proposed by the authors. In particular, this application relies on the near-identity MT (NIMT); however, the NIMT approximation used so far is not exactly unitary and leads to numerical instability. Here, we develop a discrete MT that is exactly unitary, and approximate it to obtain a discrete NIMT that is also unitary and can be computed in linear time. We prove that the discrete NIMT converges to the discrete MT when iterated, thereby allowing the NIMT to compute MTs that are not necessarily near-identity. We then demonstrate the new algorithms with a series of examples.
title Exactly unitary discrete representations of the metaplectic transform for linear-time algorithms
topic Optics
Mathematical Physics
url https://arxiv.org/abs/2012.06028