Related papers: Computation of 2D Fourier transforms and diffracti…
We propose a method based on sinc series approximations for computing the Rayleigh-Sommerfeld and Fresnel diffraction integrals of optics. The diffraction integrals are given in terms of a convolution, and our proposed numerical approach is…
We present a method of Gaussian-beam-mode decomposition to calculate classical diffraction of optical beams by apertures. This method offers a entirely different approach to examine the classic problem. Although our method is based on a…
Accurately estimating the point spread function (PSF) of an optical system requires solving free-space wave propagation, which entails evaluating a diffraction integral. This integral is traditionally computed numerically using Fast Fourier…
Based on diffraction theory and the propagation of the light, Fourier optics is a powerful tool allowing the estimation of a visible-range imaging system to transfer the spatial frequency components of an object. The analyses of the imaging…
Computing accurate estimates of the Fourier transform of analog signals from discrete data points is important in many fields of science and engineering. The conventional approach of performing the discrete Fourier transform of the data…
We present a fast algorithm for computing the diffracted field from arbitrary binary (sharp-edged) planar apertures and occulters in the scalar Fresnel approximation, for up to moderately high Fresnel numbers ($\lesssim 10^3$). It uses a…
We discuss efficient algorithms for the accurate forward and reverse evaluation of the discrete Fourier-Bessel transform (dFBT) as numerical tools to assist in the 2D polar convolution of two radially symmetric functions, relevant, e.g., to…
We introduce two efficient algorithms for computing the partial Fourier transforms in one and two dimensions. Our study is motivated by the wave extrapolation procedure in reflection seismology. In both algorithms, the main idea is to…
In this article, we derive a theoretical formalism that unifies the rigorous coupled wave analysis and the dynamical diffraction theory. Based on this formalism, we design a computational approach for the diffraction calculation for the…
We derive and showcase a novel approach to approximating Fourier transforms in higher dimensions, focusing specifically on the case of 2D radially concentrated ('ring-like') functions. We first reduce the problem to that of evaluating the…
We introduce a versatile numerical method for modeling light diffraction in periodically patterned photonic structures containing quadratically nonlinear non-centrosymmetric optical materials. Our approach extends the generalized source…
A two-dimensional (2D) mathematical model of quadratically distorted (QD) grating is established with the principles of Fraunhofer diffraction and Fourier optics. Discrete sampling and bisection algorithm are applied for finding numerical…
The paper develops a method for discrete computational Fourier analysis of functions defined on quasicrystals and other almost periodic sets. A key point is to build the analysis around the emerging theory of quasicrystals and diffraction…
Within the expansive domain of optical sciences, achieving the precise characterization of light beams stands as a fundamental pursuit, pivotal for various applications, including telecommunications and imaging technologies. This study…
Image computation is a fundamental tool for performance assessment of astronomical instrumentation, usually implemented by Fourier transform techniques. We review the numerical implementation, evaluating a direct implementation of the…
The design of starshades, i.e. external occulters for stellar coronography, relies on the fast and precise computation of their associated diffraction patterns of incoming plane waves in the telescope aperture plane. We present here a…
It was demonstrated in earlier work that, by approximating its range kernel using shiftable functions, the non-linear bilateral filter can be computed using a series of fast convolutions. Previous approaches based on shiftable approximation…
Building on recent advances in studying the co-homological properties of Feynman integrals, we apply intersection theory to the computation of Fourier integrals. We discuss applications pertinent to gravitational bremsstrahlung and deep…
Recasting the $N$-point one loop scalar integral from Feynman to Schwinger parameters gives an integrand with a Gaussian form. By application of a Fourier transform, it is easy to derive explicit expressions for the two, three and…
In this paper, we propose a numerical method of computing an integral whose integrand is a slowly decaying oscillatory function. In the proposed method, we consider a complex analytic function in the upper-half complex plane, which is…