The Bessel Expansion of Fourier Integral on Finite Interval

: In this paper, we further extend the Filon-type method to the Bessel function expansion for calculating Fourier integral. By means of complex analysis, this expansion is effective for all the oscillation frequencies. Namely, the errors of the expansion not only decrease as the order of the derivative increases, but also decrease rapidly as the frequency increases. Some numerical experiments are also presented to verify the effectiveness of the method.


Introduction
Highly oscillatory problems widely appear in many fields of mathematics and engineering, such as asymptotic analysis and electromagnetics [1][2][3][4][5][6][7]. Especially, the Fourier integral 1 −1 f (x)e iωx dx is often encountered. While the parameter |ω| 1, the integrand is highly oscillatory, and this kind of integral is often called highly oscillatory integral. This kind of integration often appears in high frequency acoustic scattering, which is why we study this kind of problem. However, when using the traditional method for calculating the highly oscillatory integrals, such as Gaussian quadrature rule, usually very small discrete steps have to be taken before they can receive a satisfactory numerical results, which makes the computational cost quite high, thus the traditional numerical integration method usually cannot be used for highly oscillatory integral calculation [3][4][5][6][7]. In 1928, Filon first proposed the Filon method [8], which laid the foundation for the numerical method in the future. In recent years, many numerical analysts have developed various Filon-type methods [3][4][5]. An excellent reference on this topic is [3]. For a simple but classical example, the Filon-type method can be defined as However, the expression and error analysis of Filon-type method are not clear enough. Therefore, we extend the Filon-type method to the Bessel function expansion in this paper. In addition, the upper bounds for reminders of the expansions are given by complex analysis.
The structure of the paper is as follows. Section 2 gives some basic theoretical results. The Bessel expansion of normal Fourier integral can be found in Section 3 and the Bessel expansion of Fourier integral with symmetric branch points are shown in Section 4. The numerical experiments are provided in Section 5 to show the effectiveness of the new method.

Some Basic Theoretical Results
We first give some theoretical results in this section. Firstly, we introduce Equations (1) and (2) (also see [9], Equation (8) and [10]).According to e iωx = ∑ ∞ j=0 (iωx) j j! and integrating item by item, Thus, according to the series expression of Bessel function of the first kind [11] We only focus on the case of ω > 0 for no segmentation. If ω > 0, we have Similarly, Now, we derive some results about Cauchy-type integral kernel, function approximation and its remainder estimation.
Assuming that U is a region, f is analytic in the region U and continuous in the closed regionŪ, i.e., f ∈ O(U) C(Ū), V ⊂ U, ∂U = Γ, and P is the polynomial. We rewrite Cauchy integral formula where we define 1 P(ζ)−P(z) as Cauchy-type integral kernel. If |P(z)| < |P(ζ)| for z ∈ V ⊂ U, ζ ∈ Γ, then the geometric series 1 is absolutely and consistently convergent such that we can exchange the order of integration and summation, where coefficients c k (z) where L is the length of Γ, the distance d = inf{|z − ζ||z ∈ V, ζ ∈ Γ}, and Proof. By the integral value theorem, the conclusion is obviously true according to the formula dζ.

The Bessel Expansion of Normal Fourier Integral
In this section, we focus on the normal Fourier integral of the form

Hermite Series of Analytic Functions
Let P(z) = 1 − z 2 , then according to Equation (3), where Γ is a Jordan curve [12] with the points ±1 in its interior, then we derive the series of f , and we define it as Hermite series, in which the coefficients are and they can be computed by the residue theorem. Let by Taylor expansion, then we obtain the coefficients Similarly, let According to Theorem 1, we estimate the truncation errors of Equation (5) as follows. where

The Bessel Expansion of Fourier Integral With Analytic Function
where the remainder satisfies Proof. The path γ is shown in Figure 1, and can be represented as according to Cauchy theorem, the error estimation According to the inequality in Equation (6), In the integral on the right side of Equation (9), since |e iωz | = e −ωa(1−t 2 ) , dz dt ≤ 1 + 4a 2 1/2 on the path γ, it is true that According to Equations (1) and (2), the conclusion is true. Letting a = 0, the right side of the inequality in Equation (8) behaves as then we have the following.

Corollary 2.
Under the assuming condition in Theorem 2, it is true that

Remark 2.
We further enlarge the inequality in Equation (8): Firstly, if n → ∞ with a = 0 and ω fixed, the right side of Equation (8) behaves as const 1 √ n(R 2 −1) n → 0. Secondly, if ω → +∞ with a = 1 and n fixed, In brief, the errors of Equation (7) not only decrease as the order n of the derivatives increases, but also decrease rapidly as the frequency ω increases.

Hermite Series with Symmetric Branch Points
The function ϕ can be expressed by the multiple Cauchy integral formula where Γ is the skeleton of U, and, if then the above integral kernel can be written as p and noting that f (z) = ϕ(z 1 , z 2 ), we have Letting the coefficients we can rewrite Equation (11) as Hermite series of f , where the remainder is The coefficients can be computed by the multi-dimension residue theorem. Letting by multi-Taylor expansion, then we obtain the coefficients Similarly, letting (12), if f (z) = ϕ(z 1 , z 2 ), ϕ ∈ O(U) C(Ū) with z 1 = z, z 2 = (1 − z 2 ) 1 p , p ∈ Z + and r 1 > 1, r 2 > (1 + r 2 1 ) 1/p , then we have the error estimation of the remainder

The Bessel Expansion of Fourier Integral with Symmetric Branch Points
Similar to Theorem 2, we have the Bessel expansion of Fourier integral with symmetric branch points, and the error estimation is as follows.

The Numerical Experiments
In this section, for the two cases of analytic functions and branch functions, we present two examples that verify the validity of Bessel expansion to calculate Fourier integrals. These numerical experiments were done in Maple 16.

The Conclusions and Acknowledgement
In this article, we make the expression of the Filon-type method more concise by the Bessel expansions for calculating Fourier integral. Moreover, we also derive the convergence rate of the Filon-type method by means of complex analysis. In fact, both error analysis and numerical tests have shown that the Filon-type method is efficient. On the road of frequency analysis and calculation, we will continue to study.
Funding: This research was funded by project of enhancing school with innovation of Guangdong Ocean University grant number GDOU2016050228 and the APC was funded by project of enhancing school with innovation of Guangdong Ocean University(Q18306).