A flexible numerical calculation method of angular spectrum based on matrix product

Fast Fourier transform (FFT) is the most commonly used mathematical method in numerical calculation, and the FFT-based angular spectrum method (ASM) is also used widely in diffraction calculation. However, the frequency and spatial sampling rules in FFT limit the effective propagation distance and the observation window range of ASM. A novel method for calculating the angular spectrum based on the matrix product is proposed in this Letter. This method realizes the fast calculation of discrete Fourier transform (DFT) based on the matrix product, in which the sampling matrix is orthogonally decomposed into two vectors. Instead of FFT, angular spectrum diffraction calculation is carried out based on the matrix product, which is named the matrix product ASM. The method in this Letter uses a simple mathematical transformation to achieve maximum compression of the sampling interval in the frequency domain, which significantly increases the effective propagation distance of the angular spectrum. Additionally, the size of the observation window can be enlarged to obtain a wider calculation range by changing the spatial sampling of the output plane.