Rafael G. Campos,Jared Figueroa
Abstract:In recent years there has been a renewed interest in finding fast algorithms to compute accurately the linear canonical transform (LCT) of a given function. This is driven by the large number of applications of the LCT in optics and signal processing. The well-known integral transforms: Fourier, fractional Fourier, bilateral Laplace and Fresnel transforms are special cases of the LCT. In this paper we obtain an O(N*Log N) algorithm to compute the LCT by using a chirp-FFT-chirp transformation yielded by a convergent quadrature formula for the fractional Fourier transform. This formula gives a unitary discrete LCT in closed form. In the case of the fractional Fourier transform the algorithm computes this transform for arbitrary complex values inside the unitary circle and not only at the boundary. In the case of the ordinary Fourier transform the algorithm improves the output of the FFT.
What problem does this paper attempt to address?
The problem that this paper attempts to solve is to develop a fast algorithm to calculate the Linear Canonical Transform (LCT) efficiently and accurately. The LCT is a three - parameter integral transform, which is related to the classical transform in quantum mechanics and has wide applications in the fields of optics and signal processing. Common integral transforms such as the Fourier transform, the fractional - order Fourier transform, the bilateral Laplace transform and the Fresnel transform are all special cases of the LCT.
Specifically, the goals of the paper include:
1. **Improve computational efficiency**: By proposing a fast algorithm with a time complexity of \(O(N \log N)\) to calculate the LCT, thereby significantly reducing the computation time and resource consumption.
2. **Expand the scope of application**: This algorithm is not only applicable to the standard Fourier transform, but also can handle any complex values in the fractional - order Fourier transform and perform calculations within the unit circle (not just on the boundary).
3. **Improve existing methods**: For the ordinary Fourier transform, this algorithm can improve the results of the Fast Fourier Transform (FFT).
To achieve these goals, the author uses a convergent quadrature formula derived from the asymptotic properties of Hermite polynomials to decompose the LCT into a scaling - modulation - DFT - modulation - scaling transform, and finally obtains a discrete form of the LCT, namely the Extended Fourier Transform (XFT). This discretization method makes the calculation of the LCT more efficient and accurate.
### Formula summary
The linear canonical transform is defined as:
\[
L_{a,b,c,d}\{f(x), y\} =
\begin{cases}
\frac{1}{\sqrt{2\pi i b}} \int_{-\infty}^{\infty} e^{i\left(\frac{a}{2b}x^2 - \frac{x y}{b} + \frac{d}{2b}y^2\right)} f(x) \, dx & \text{if } b \neq 0 \\
\sqrt{d} e^{i \frac{c}{2d} y^2} f(dy) & \text{if } b = 0
\end{cases}
\]
where \( a, b, c, d \) satisfy \( ad - bc = 1 \), so only three parameters are free.
Through the above methods, the author successfully simplifies the calculation of the LCT into an efficient matrix operation problem, thereby achieving fast calculation.