A method for approximating the solution of weakly singular Fredholm integral equation of the second kind with highly oscillatory trigonometric kernel is presented. The unknown function is approximated by expansion of Chebychev polynomial and the coefficients are determinated by classical collocation method. Due to the highly oscillatory kernels of integral equation, the discretised collocation equation will give rise to the computation of oscillatory integrals. These integrals are calculated by using recursion formula derived from the fundamental recurrence relation of Chebyshev polynomial. The effectiveness and accuracy of the proposed method are tested by numerical examples.

1. Introduction

The Fredholm integral equations of second kind, where is a continuous function and is a given function, have many applications in mathematical physics and engineering, such as heat conduction problem, potential problems, quantum mechanics, and seismology image processing [14]. Particularly, when and , (1) is called weakly singular.

In most of the cases, the integral equation cannot be done analytically and one has to resort to numerical methods. Many numerical methods, such as collocation method and Galerkin method, have been developed to solve (1); for details see [2, 5, 6]. These methods are well-established numerical algorithms; however, standard version of these classical methods may suffer from difficulty for computation of (1), containing highly oscillatory kernels since the computation of the highly oscillatory integrals by standard quadrature methods is exceedingly difficult and the cost steeply increases with the frequency. Furthermore, Galerkin method requires many double integrals when approximating the solution of integral equation. Specially, when the kernel is highly oscillatory, it requires the evaluation of many highly oscillatory double integrals, which can become computationally expensive.

Recently, for weakly singular Volterra integral equations of the second kind with highly oscillatory Bessel kernels, it was found that the collocation methods are much more easily implemented and can get higher accuracy than discontinuous Galerkin methods under the same piecewise polynomials space; for details see [710]. In addition, collocation method only involves single integrals which are a little easier to evaluate. Motivated by this fact, here we investigate the application of collocation method for the solution of Fredholm integral equation Here is the unknown function and is a given function and assume that is not the eigenvalue of integral equation, since for any finite interval it can be transformed to interval by linear transformation. In this paper, we only consider the case of .

On the other hand, when the unknown function is approximated by Chebychev polynomial, the discretization of the integral equation (2) by collocation method will give rise to highly oscillatory integral which can be formulated as where denotes the Chebychev polynomials of the first kind.

In last few years many efficient methods have been devised for the evaluation of oscillatory integral, such as asymptotic method [11], Filon-type method [12, 13], Levin’s collocation method [14], modified Clenshaw-Curtis method, Clenshaw-Curtis-Filon-type method [15], and generalized quadrature rule [16], although some of these methods do not involve Chebyshev polynomial. Piessens and Poleunis [17] consider a simpler case of and use Chebyshev polynomials to evaluate the integral by somewhat indirect method involving a truncated infinite series of Bessel functions. The evaluation of this infinite summation suggests a defect of this method. An alternative procedure which avoids this infinite series is the original Bakhvalo and Vasil'eva-Legendre work involved in the evaluation of the integral by recurrence relation, where denotes Legendre polynomial. Sadly, this relation proved to be unstable in the forwards direction for small . In [18], Alaylioglu et al proposed a simple alternative approach analogous to Newton-Cotes based formulae. This method avoids the instability to a certain degree compared with these two methods mentioned above. We used the same idea and generalized it to the case of weakly singular and calculated oscillatory integral (3).

This paper is organized as follows: in Section 2 we derive some basic formulae and introduce some mathematical preliminaries of the proposed method. In Section 3 we discuss the evaluation of the integrals occurring in collocation equation. In Section 4 numerical experiments are conducted to illustrate the performance of the proposed method.

2. Fundamental Relations

Firstly, by separation of real and imaginary part of and , we transform the integral equation (2) into equivalent systems of two linear integral equations of Fredholm in the forms where , , and .

Assume that and ; then, systems of linear integral equations (5) can be written in the matrix form where Set , , where , .

Hence, unknown functions can be expressed by where

Then, the aim is to find Chebyshev coefficients, that is, the matrix . We first substitute the Chebyshev collocation points, which are defined by ,  , into (6) and then rearrange a new matrix form to determine : in which is the integral part of (6) and

By substituting (8) into (10), the unknown coefficients can be easily computed from this linear algebraic equations and therefore we find the solution of integral equation (2).

3. Evaluation of the Integral

The discretization of integral equation will lead to the calculation of integral . In this section, we will derive the recursion formula to compute it efficiently from the fundamental recurrence relation of Chebyshev polynomial.

The Chebyshev polynomials are of the form and the coefficients of in can be easily calculated by the means of the recurrence relation By expanding the Chebyshev polynomial in terms of powers of , the integral would be transformed into the form of where .

By separation of real and imaginary part of , (14) is transformed into whose main difficulty now is turning to the evaluation of the following basic integrals:

With substitution for definite integral and binomial theorem, (16) is equivalent to where is number of combination and can be formulated analogy.

Efficient evaluation of (18) is based on accurate calculation of integrals in the formula which can be computed explicitly by the incomplete Gamma function [19]; that is, in which ; for details see ([20], pp 215).

Once the integral is obtained, by substituting into (10), the coefficient matrix is derived; then, we can compute the solution of integral by solving these linear algebraic equations.

4. Numerical Examples

In this section, we give some numerical examples to illustrate the performance of proposed method. In all the following examples, is the number of mesh points, denotes the approximate solution, where is the number of terms of the Chebyshev series, and denotes the exact solution, respectively. All the computations have been performed by using Matlab R2012a on a 2.5 GHz PC with 2 GB of RAM.

Example 1. We consider and set ; then, the solution of which is We plot the absolute error by Matlab 2012a internal function ezplot. Specially, for , with and , solve linear algebraic equations (10); it can be found that the approximate solutions are

It is easy to see from Figure 1 that the proposed method is converging and we only need so we could achieve high accuracy.

Example 2. We consider and set ; then, the solution of which is

In this case, for , with and , the approximate solutions are

It is obvious to see from Figure 2 that the proposed method is efficient and could achieve high accuracy with little number of collocation points.

Example 3. For a general case, we consider and set . Although the solution of this equation is unknown, we can verify the calculation precision of the method to a certain degree by comparing the approximate evaluation at mesh points.

In this case, for , with and , the approximate solutions are

It is also easy to see from Tables 1, 2, 3, and 4 that the presented method is efficient and accurate, although the exact solution is unknown.

5. Conclusion

In this paper, we explore quadrature methods for weakly singular Fredholm integral equation of the second kind with oscillatory trigonometric kernels and present collocation methods with Chebyshev series for calculation of the solution. For integral equation with highly oscillatory kernels, the standard collocation methods with classical quadrature methods are not suitable for the numerical approximation of the solution of integral equation, since the computation of the highly oscillatory integrals by standard quadrature methods is exceedingly difficult and the cost steeply increases with the frequency. Based on the recursion formula derived in Section 3, we compute the highly oscillatory integrals occurring in collocation equation, directly and efficiently. Numerical examples demonstrate the performance of algorithm.

Conflict of Interests

The author declares that there is no conflict of interests regarding the publication of this paper.


This paper was supported by NSF of China (no. 11371376).