Efficient method for the computation of oscillatory Bessel transform and Bessel Hilbert transform

详细信息    查看全文

• 作者：Zhenhua Xu^a ; ^{xuzhenhua19860536@163.com" class="auth_mail" title="E-mail the corresponding author} ; Gradimir V. Milovanović^b ; ^{gvm@mi.sanu.ac.rs" class="auth_mail" title="E-mail the corresponding author}
• 关键词：65D32 ; 65D30
• 刊名：Journal of Computational and Applied Mathematics
• 出版年：2016
• 出版时间：15 December 2016
• 年：2016
• 卷：308
• 期：Complete
• 页码：117-137
• 全文大小：740 K

文摘

In this paper, we study the numerical methods for the evaluation of two kinds of highly oscillatory Bessel transforms. Firstly, we rewrite both integrals as the sum of two integrals. By rewriting the Bessel function as a linear combination of Whittaker W function, we then transform one of integrals into the Fourier type, which can be transformed into the integrals on [0,+∞)$[0,+\infty )$, and can be computed by some proper Gaussian quadrature, which take into account the asymptotic property of Whittaker W function as x→0$x\to 0$. The other can be efficiently computed based on the evaluation of special functions. In addition, error analysis for the presented methods is given. Moreover, we also give an explicit formula for the integral ${⨍}_0^{+\infty }\frac{J_{\nu }\left(\omega x\right)}{x-\tau }dx$ in terms of Meijer G$G$-function, and then apply the method for the oscillatory Bessel transforms to the computation of highly oscillatory Bessel Hilbert transforms. Theoretical results and numerical examples demonstrate the efficiency and accuracy of the proposed methods.