Enhanced Fast Fractional Fourier Transform (FRFT) Scheme Based on Closed Newton-Cotes Rules

The paper presents an enhanced numerical framework for computing the one-dimensional fast Fractional Fourier Transform (FRFT) by integrating closed-form Composite Newton-Cotes quadrature rules. We show that a FRFT of a QN-length weighted sequence can be decomposed analytically into two mathematically commutative compositions: one involving the composition of a FRFT of an N-length sequence and a FRFT of a Q-length weighted sequence, and the other in reverse order. The composite FRFT approach is applied to the inversion of Fourier and Laplace transforms, with a focus on estimating probability densities for distributions with complex-valued characteristic functions. Numerical experiments on the Variance-Gamma (VG) and Generalized Tempered Stable (GTS) models show that the proposed scheme significantly improves accuracy over standard (non-weighted) fast FRFT and classical Newton-Cotes quadrature, while preserving computational efficiency. The findings suggest that the composite FRFT framework offers a robust and mathematically sound tool for transform-based numerical approximations, particularly in applications involving oscillatory integrals and complex-valued characteristic functions.