Table of Contents Author Guidelines Submit a Manuscript
Mathematical Problems in Engineering
Volume 2009, Article ID 786368, 10 pages
http://dx.doi.org/10.1155/2009/786368
Research Article

Fourier Approximation for Integral Equations on the Real Line

Department of Mathematics, Faculty of Science, K. N. Toosi University of Technology, P.O. Box 16765-165, Tehran 19697, Iran

Received 8 February 2009; Accepted 22 May 2009

Academic Editor: Francesco Pellicano

Copyright © 2009 S. M. Hashemiparast and H. Fallahgoul. 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

Based on Guass quadradure method a class of integral equations having unknown periodic solution on the real line is investigated, by using Fourier series expansion for the solution of the integral equation and applying a process for changing the interval to the finite interval (; 1), the Chebychev weights become appropriate and examples indicate the high accuracy and very good approximation to the solution of the integral.

1. Introduction

Many numerical methods for approximating the solution of integral equations are based on Galerkin and Collocation methods (see [1, 2]). Some modifications in these methods such as iterated Galerkin and iterated Collocation methods give more accurate approximations to the solution [3], some others have extended the technics to Petrov-Galerkin and Hammersteion equations degenerate kernel methods (see [47]), in some of recent methods further improvement in some steps is not possible [3, 6], the computation of sufficiently accurate solution may require the use of much partition of the domain and may thus involve the solution of a correspondingly large linear system. On the other hand, it is well known that the coefficient matrix of an integral equation is dense, and it requires heavy work to generate the coefficient matrix if the partition is much finer (see [8, 9]). To improve the difficulties in this methods, improved iterated Galerkin methods are recommended (see [10]), other improvement with less additional computational efforts for each step of iteration are introduced and shows that every step of iteration the coefficient matrix of the linear system to be inverted remains the same as that in the original approximation method (see [11]), the additional order of convergence for every step of iteration is which is less than in Galerkin, Collocation and degenerate kernel methods, where is the norm of the partition and is order of the piecewise polynomials used in the approximation. Less works consider the integral equations on the real line (see [12, 13]), integral equations with singularity (see [14, 15]), integral equations with periodic solutions (see [16]).

In this paper we consider the approximate solutions of Fredholm integral equations on the real line having periodic solutions.

We consider the following Fredholm integral equations of second kind with periodic solution on the real line the main part of the first section considers the several admissibility results and existence principle on real line, however before proceeding to other sections we recall some facts about periodic functions, we discuss the general approach that will be taken to obtain the various existence principles [16].

2. Fredholm Integral Equation on the Real Line

Let and suppose is an interval of that contains at least one compact subinterval of length , which we denote by . We define to be the subspace of (bounded continuous functions on with values in ) consisting of all -periodic mappings, that is, if , then f is continuous on and , for all t such that The norm on is the same as the norm on and it is clear that for we have We consider Fredholm integral equations over the real line that is separable where and is a polynomial of degree at most , which is nonnegative on the half line and is an arbitrary function.

Problem (2.2) can be simplified by first obtaining the partial fraction decomposition of (2.4) in the form where the sum is over all pairs of conjugate complex poles of with corresponding multiplicities Here,

Thus, without loss of generality we can consider only weights of the form In (2.2) the function is -periodic:

Let

3. Reduction of the Integral to a Finite Interval

Milovanovic and et al. [17] show how to reduce the integral (2.9) to an integral on the finite interval. For this purpose we need the sum of the following series: Since (cf. [9]) in the simplest, but the most important case for and we obtain We bring a lemma and two theorems in the following which have been proved in [17].

Lemma 3.1. Let be given by (2.7), and Then

The proof of this result can be done by an integration of the function over the rectangular contour with vertices at the points where is such that the poles of the function are inside of Then, taking the corresponding integral over tends to zero, because when Then by Cauchy’s residue theorem, we get For (3.4) reduce to (3.3). When where We can suspect the following form of our sum: where is an algebraic polynomial. Indeed, we can prove the following result.

Theorem 3.2. Let and Then where is a nonnegative polynomial on of degree These polynomials satisfy the recurrence relation where

Proof. We start with (3.3) written in the form , where Thus, the formula (3.9) is true for .
Suppose that (3.9) holds for some . Then, differentiating with respect to , we get that is, Thus, the result is proved.

We are ready now to give a transformation of the integral (2.9) to one on a finite interval. Putting and using the periodicity of the function , we have

because of the uniform convergence of the series (3.1). Thus, where is defined by (3.1) and given by Theorem (3.2). We see that that is, is an even weight function.

Because of the last property of the weight function, we have

Changing the variables and putting we get the following result.

Theorem 3.3. The integral (2.9) can be transformed to the form where is a polynomial determined by the recurrence relation Theorem 3.2, and is defined by (3.20).

4. Gaussian Type Formula for Chebychev Weights

In order to evaluate the integral (3.20) it would seem more natural and simpler to apply the Gauss-Chebychev quadrature formula, that is, taking as in integrating function with respect to the Chebychev weight (Chw)

In this case, when for some the function satisfied the condition the error of the -point Gauss-Chebychev quadrature can be estimated as follow: where is a constant independent on and . Hence, when is very close to 1, even if the integrand is bounded, it gives a very large bound.

Thus, for evaluating the integral (3.20) it is more convenient to construct the Gaussian quadratures for the measure where the function includes the algebraic polynomial that is, Here, denotes the set of all polynomials of degree at most and is the root of Chebychev polynomials.

5. Solving Integral Equation

Let this integral equation which is -periodic, BE separable, that this integral equation is simplified according to Sections 2 and 3: where and is obtained from Theorem 3.2.

For getting unknown function in relation (5.1) which, is -periodic. We have the following Fourier expansion of be -periodic: indeed we will obtain unknowns coefficients and by substituting in (5.2) and using Chebychev integration the relation can be solved for these unknowns, then by substituting the Fourier coefficients in (5.3) will be determined.

6. Numerical Examples

In this section we consider some numerical examples to illustrate the presented transformation method. All computations were done by Maple 11, and the runtime for this calculation is too short. It should be noted that the real line is the integration interval of the following examples, but the integrands are periodic functions and they have been evaluated in the period .

Example 6.1. Let us consider the following integral equation: where and are arbitrary knowns and is unknown -periodic function on which should be determined. Also and are equal 1 in (2.8).Applying Gaussian quadrature with the ChW for , we get the approximate Fourier expantion of with relative error given in Table 1. The exact solution of (6.1) is the approximated solution of according to the values of and is which shows the difference of the solutions for equals and other values of the errors are given in Table 1, which indicates a super convergence of the solution even for small values of .

tab1
Table 1

Example 6.2. Consider the integral equation where and are arbitrary and is -periodic function on also and are equal 1 in (2.8).

Applying Gaussian quadrature with the Chw, for .

The exact and approximated solutions are considered to illustrate in Table 2 with comparison is exact solution the for the error is almost zero.

tab2
Table 2

Example 6.3. Consider the following integral equation: where is arbitrary amount and is any function of . If and by getting (Chw) quadrature we will obtain the approximations function for . The exact solution is if then and for and So the approximation and exact solution almost equal, or The approximation solution is almost exact.

7. Conclusion

A Fourier series approximation for the integral equations having periodic solution are considered in the real line, the procedure produced almost an exact approximation for the solution of the integral even by taking few number of the elements in the Fourier series expansion of the solution, the numerical examples indicate the high accuracy of the approximation.

References

  1. K. E. Atkinson, The Numerical Solution of Integral Equations of the Second Kind, vol. 4 of Cambridge Monographs on Applied and Computational Mathematics, Cambridge University Press, Cambridge, UK, 1997. View at Zentralblatt MATH · View at MathSciNet
  2. R. Kress, Linear Integral Equations, vol. 82 of Applied Mathematical Sciences, Springer, Berlin, Germany, 1989. View at MathSciNet
  3. I. H. Sloan, “Four variants of the Galerkin method for integral equations of the second kind,” IMA Journal of Numerical Analysis, vol. 4, no. 1, pp. 9–17, 1984. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet
  4. Z. Chen and Y. Xu, “The Petrov-Galerkin and iterated Petrov-Galerkin methods for second-kind integral equations,” SIAM Journal on Numerical Analysis, vol. 35, no. 1, pp. 406–434, 1998. View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet
  5. Z. Chen, Y. Xu, and J. Zhao, “The discrete Petrov-Galerkin method for weakly singular integral equations,” Journal of Integral Equations and Applications, vol. 11, no. 1, pp. 1–35, 1999. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet
  6. I. G. Graham, S. Joe, and I. H. Sloan, “Iterated Galerkin versus iterated collocation for integral equations of the second kind,” IMA Journal of Numerical Analysis, vol. 5, no. 3, pp. 355–369, 1985. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet
  7. H. Kaneko and Y. Xu, “Superconvergence of the iterated Galerkin methods for Hammerstein equations,” SIAM Journal on Numerical Analysis, vol. 33, no. 3, pp. 1048–1064, 1996. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet
  8. Z. Chen, C. A. Micchelli, and Y. Xu, “A multilevel method for solving operator equations,” Journal of Mathematical Analysis and Applications, vol. 262, no. 2, pp. 688–699, 2001. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet
  9. Z. Chen, C. A. Micchelli, and Y. Xu, “Fast collocation methods for second kind integral equations,” SIAM Journal on Numerical Analysis, vol. 40, no. 1, pp. 344–375, 2002. View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet
  10. D. Porter and D. S. G. Stirling, “The reiterated Galerkin method,” IMA Journal of Numerical Analysis, vol. 13, no. 1, pp. 125–139, 1993. View at Publisher · View at Google Scholar · View at MathSciNet
  11. G. Long and G. Nelakanti, “Iteration methods for Fredholm integral equations of the second kind,” Computers & Mathematics with Applications, vol. 53, no. 6, pp. 886–894, 2007. View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet
  12. C. Frammartino, C. Laurita, and G. Mastroianni, “On the numerical solution of Fredholm integral equations on unbounded intervals,” Journal of Computational and Applied Mathematics, vol. 158, no. 2, pp. 355–378, 2003. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet
  13. M. C. De Bonis, C. Frammartino, and G. Mastroianni, “Numerical methods for some special Fredholm integral equations on the real line,” Journal of Computational and Applied Mathematics, vol. 164-165, pp. 225–243, 2004. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet
  14. D. O'Regan, R. P. Agarwal, and K. Perera, “Nonlinear integral equations singular in the dependent variable,” Applied Mathematics Letters, vol. 20, no. 11, pp. 1137–1141, 2007. View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet
  15. Y. Cao, M. Huang, L. Liu, and Y. Xu, “Hybrid collocation methods for Fredholm integral equations with weakly singular kernels,” Applied Numerical Mathematics, vol. 57, no. 5–7, pp. 549–561, 2007. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet
  16. D. O'Regan and M. Meehan, “Periodic and almost periodic solutions of integral equations,” Applied Mathematics and Computation, vol. 105, no. 2-3, pp. 121–136, 1999. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet
  17. G. Mastroianni and G. V. Milovanović, “Weighted integration of periodic functions on the real line,” Applied Mathematics and Computation, vol. 128, no. 2-3, pp. 365–378, 2002. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet