Journal of Applied Mathematics

Volume 2014 (2014), Article ID 172327, 7 pages

http://dx.doi.org/10.1155/2014/172327

Research Article

The Approximate Solution of Fredholm Integral Equations with Oscillatory Trigonometric Kernels

Qinghua Wu

School of Mathematics and Statistics, Central South University, Changsha, Hunan 410083, China

Received 9 March 2014; Accepted 4 April 2014; Published 17 April 2014

Academic Editor: Turgut Öziş

Copyright © 2014 Qinghua Wu. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

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 [1–4]. 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 [7–10]. 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.

(a)

(b)

(c)

Figure 1: Example 1: the absolute error versus the -coordinate for (a), (b), and (c) with negative data ignored.

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.

(a)

(b)

(c)

Figure 2: Example 2: the absolute error versus the -coordinate for (a), (b), and (c) with negative data ignored.

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.

Table 1: Approximations for with .

Table 2: Approximations for with .

Table 3: Approximations for with .

Table 4: Approximations for with .

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.

Acknowledgment

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

References

V. M. Aleksandrov and E. V. Kovalenko, “Mathematical method in the displacement problem,” Inzhenernyj Zhurnal. Mekhanika Tverdogo Tela, vol. 2, pp. 77–89, 1984.

H. Brunner, On the Numerical Solution of First-Kind Volterra Integral Equations with Highly Oscillatory Kernels, Highly Oscillatory Problems: From Theory to Applications, Isaac Newton Institute, HOP 13-17, 2010.

R. Estrada and R. P. Kanwal, Singular Integral Equations, Birkhäauser, Boston, Mass, USA, 2000.

A. D. Polyanin and A. V. Manzhirov, Handbook of Integral Equations, CRC Press, 1998.

H. Brunner, Collocation Methods for Volterra Integral and Related Functional Equations, Cambridge University Press, Cambridge, Cambridge, UK, 2004.

I. G. Graham, “Galerkin methods for second kind integral equations with singularities,” Mathematics of Computation, vol. 39, no. 160, pp. 519–533, 1982.

S. Xiang and H. Brunner, “Efficient methods for Volterra integral equations with highly oscillatory Bessel kernels,” BIT. Numerical Mathematics, vol. 53, no. 1, pp. 241–263, 2013.

S. Xiang and K. He, “On the implementation of discontinuous Galerkin methods for Volterra integral equations with highly oscillatory Bessel kernels,” Applied Mathematics and Computation, vol. 219, no. 9, pp. 4884–4891, 2013.

S. Xiang and Q. Wu, “Numerical solutions to Volterra integral equations of the second kind with oscillatory trigonometric kernels,” Applied Mathematics and Computation, vol. 223, pp. 34–44, 2013.

Q. Wu, “On graded meshes for weakly singular Volterra integral equations with oscillatory trigonometric kernels,” Journal of Computational and Applied Mathematics, vol. 263, pp. 370–376, 2014.

A. Iserles and S. P. Nørsett, “On quadrature methods for highly oscillatory integrals and their implementation,” BIT. Numerical Mathematics, vol. 44, no. 4, pp. 755–772, 2004.

A. Iserles and S. P. Nørsett, “Efficient quadrature of highly oscillatory integrals using derivatives,” Proceedings of The Royal Society of London A, vol. 461, no. 2057, pp. 1383–1399, 2005.

S. Xiang, “Efficient Filon-type methods for ,” Numerische Mathematik, vol. 105, no. 4, pp. 633–658, 2007.

D. Levin, “Fast integration of rapidly oscillatory functions,” Journal of Computational and Applied Mathematics, vol. 67, no. 1, pp. 95–101, 1996.

S. Xiang, Y. J. Cho, H. Wang, and H. Brunner, “Clenshaw-Curtis-Filon-type methods for highly oscillatory Bessel transforms and applications,” IMA Journal of Numerical Analysis, vol. 31, no. 4, pp. 1281–1314, 2011.

G. A. Evans and K. C. Chung, “Some theoretical aspects of generalised quadrature methods,” Journal of Complexity, vol. 19, no. 3, pp. 272–285, 2003.

R. Piessens and F. Poleunis, “A numerical method for the integration of oscillatory functions,” vol. 11, pp. 317–327, 1971.

A. Alaylioglu, G. A. Evans, and J. Hyslop, “The use of Chebyshev series for the evaluation of oscillatory integrals,” The Computer Journal, vol. 19, no. 3, pp. 258–267, 1976.

M. Abramowitz and I. A. Stegun, Handbook of Mathematical Functions, National Bureau of Standards, Washington, DC, USA, 1964.

I. S. Gradshteyn and I. M. Ryzhik, Table of Integrals, Series, and Products, New York Academic Press, 1965.