Complex generalized Gauss-Radau quadrature rules for Hankel transforms of integer order

Haiyong Wang,Menghan Wu
2024-03-28
Abstract:Complex Gaussian quadrature rules for oscillatory integral transforms have the advantage that they can achieve optimal asymptotic order. However, their existence for Hankel transform can only be guaranteed when the order of the transform belongs to $[0,1/2]$. In this paper we consider the construction of generalized Gauss-Radau quadrature rules for Hankel transform. We show that, if adding certain value and derivative information at the left endpoint, then complex generalized Gauss-Radau quadrature rules for Hankel transform of integer order can be constructed with theoretical guarantees. Orthogonal polynomials that are closely related to such quadrature rules are investigated and their existence for even degrees is proved. Numerical experiments are presented to confirm our findings.
Numerical Analysis,Classical Analysis and ODEs
What problem does this paper attempt to address?
The problem that this paper attempts to solve is to construct complex generalized Gauss - Radau quadrature rules for the integer - order Hankel transform and ensure the theoretical feasibility of these rules. Specifically, the paper has solved the following problems: 1. **Existence problem**: The traditional complex Gauss quadrature rules can only guarantee existence for the Hankel transform when the transform order ν belongs to [0, 1/2]. In this paper, by introducing specific values and derivative information, this range is extended to integer orders, and the existence of such rules is proved. 2. **Optimal asymptotic order**: The paper shows that the proposed complex generalized Gauss - Radau quadrature rules can not only be constructed successfully but also reach the optimal asymptotic order, that is, their errors decay at the fastest algebraic rate as the frequency ω increases. 3. **Orthogonal polynomial research**: The paper explores the properties of orthogonal polynomials closely related to these quadrature rules and proves that for even - order cases, these polynomials do exist. 4. **Numerical experiment verification**: The correctness and effectiveness of the above - mentioned theoretical results are verified through numerical experiments. ### Main contributions - **Construction method**: A method of adding specific values and derivative information at the left - end point is proposed, so that the complex generalized Gauss - Radau quadrature rules can be applied to the integer - order Hankel transform. - **Theoretical guarantee**: The theoretical feasibility of the proposed method is proved, including the existence of quadrature nodes and weights and the realization of the optimal asymptotic order. - **Application expansion**: The discussion is extended to the fractional - order Hankel transform and the Fourier sine transform, and an application example of the oscillatory Hilbert transform is given. ### Key formulas - Definition of Hankel transform: \[ (H_\nu f)(\omega) := \int_0^\infty f(x) J_\nu(\omega x) dx \] where \( J_\nu(x) \) is the Bessel function of the first kind, and \( \omega \) is the oscillation frequency. - Complex generalized Gauss - Radau quadrature rules: \[ (QHI_{2n,\mu} f)(\omega) = \frac{1}{\omega} \left( \sum_{k = 0}^{\mu - 1} \hat{w}_0^k \frac{f^{(k)}(\hat{x}_0)}{\omega^k} + 2 \sum_{j = 1}^{2n} \hat{w}_j f\left(\frac{\hat{x}_j}{\omega}\right) \right) \] These contributions provide a new and efficient numerical method for dealing with the Hankel transform and lay a theoretical foundation for further research in related fields.