Abstract

The beam propagation method (BPM) can be viewed as a chain of alternating convolutions and multiplications, as filtering operations alternately in the space and frequency domains or as multiplication operations sandwiched between linear canonical or fractional Fourier transforms. These structures provide alternative models of inhomogeneous media and potentially allow mathematical tools and algorithms associated with these transforms to be applied to the BPM. As an example, in the case where quadratic approximation is possible, it is shown that the BPM can be represented as a single LCT system, leading to significantly faster computation of the output field.

© 2022 Optica Publishing Group

Equivalence of linear canonical transform domains to fractional Fourier domains and the bicanonical width product: a generalization of the space–bandwidth product

Figen S. Oktem and Haldun M. Ozaktas
J. Opt. Soc. Am. A 27(8) 1885-1895 (2010)

Fast numerical algorithm for the linear canonical transform

Bryan M. Hennelly and John T. Sheridan
J. Opt. Soc. Am. A 22(5) 928-937 (2005)

Fast linear canonical transforms

John J. Healy and John T. Sheridan
J. Opt. Soc. Am. A 27(1) 21-30 (2010)

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

Bryan M. Hennelly and John T. Sheridan
J. Opt. Soc. Am. A 22(5) 917-927 (2005)

Eigenfunctions of the offset Fourier, fractional Fourier, and linear canonical transforms

Soo-Chang Pei and Jian-Jiun Ding
J. Opt. Soc. Am. A 20(3) 522-532 (2003)