Two-dimensional wavelets for numerical solution of integral equations

Purpose

In this paper, we shall investigate the numerical solution of two-dimensional Fredholm integral equations (2D-FIEs).

Methods

In this work, we apply two-dimensional Haar wavelets, to solve linear two dimensional Fredholm integral equations (2D-FIEs). Using 2D Haar wavelets and their properties, 2D-FIEs of the second kind reduce to a system of algebraic equations.

Results

The numerical examples illustrate the efficiency and accuracy of the method.

Conclusions

In comparison with other bases (for example, polynomial bases), one of the advantages of this method is, although the involved matrices have a large dimension, they contain a large percentage of zero entries, which keeps computational effort within reasonable limits.

The integral equations provide an important tool for modeling a numerous phenomena and processes, and for solving boundary value problems for both ordinary and partial differential equations. Their historical development is closely related to the solution of boundary value problems in potential theory. In the last decades, there has been much interest in numerical solutions of integral equations. The Nystrom and collocation methods are probably the two most important approaches for the numerical solution of these integral equations [1,2]. While several numerical methods are known for one-dimensional integral equations, fewer methods are known for two-dimensional integral equations [3-6].

Recently, many different basic functions have been used to estimate the solution of integral equations, such as orthogonal functions and wavelets. Haar wavelets are the simplest orthogonal wavelet with compact support, and they have been used in different numerical approximation problems.

In this work, we apply two-dimensional Haar wavelets, constructed on , to solve linear two-dimensional Fredholm integral equations (2D-FIEs) of the form:

where is an unknown function to be found and the functions and are given continuous functions defined on D and D², respectively. The existence and uniqueness results for Equation 1 can be found in the classical theory of Fredholm integral equations.

Two-dimensional Haar wavelets

We usually call the Haar wavelets containing one variable as one-dimensional, and those containing two variables as two-dimensional. One-dimensional Haar wavelets have been widely used for solving different problems [6-8]. Complete details for one-dimensional Haar wavelets is found in [9,10]. These discussions can also be extended to the two-dimensional one.

Definitions and properties

Definition 2.1

The orthogonal basis of one-dimensional Haar wavelets for the Hilbert space consists of

where

The integer 2^j indicates the level of the wavelet and k is the translation parameter.

Simple calculations show that

Also, it can be shown that any function can be expressed as , where [11].

Definition 2.2

Let be the one-dimensional Haar wavelets on . We call the two-dimensional Haar wavelets on as:

where , with and .

The family is orthogonal on and forms a basis for :

Theorem 2.3

The basis is orthonormal on .

Proof

Let or

Theorem 2.4

Proof

The expansion of a function

A function defined over may be expanded by the two-dimensional Haar wavelets as

where the wavelet coefficients, , are obtained as

If the infinite series in Equation 6 is truncated up to their k terms, then it can be written as

where , and is a nonnegative integer. Here, and are the Haar wavelet coefficients and Haar wavelet functions vectors, respectively, and defined as:

Similarly, a function of four variables, , on may be approximated with respect to Haar wavelets such as:

where and are two-dimensional Haar wavelets vectors of dimension k², and K is the two-dimensional Haar coefficient matrix.

Solution of 2D-FIEs of the second kind

Now, consider the second kind Fredholm integral equation of the form in Equation 1. Our goal is to reduce this equation to a linear system of algebraic equations by the method presented in this paper.

In order to approximate the solution of integral equation (Equation 1), we approximate functions and with respect to 2D-Haar wavelets by the way mentioned in Two-dimensional Haar wavelets section as

where is as defined in Equation 10, the vectors and matrix are Haar wavelets coefficients of and , respectively.

By substituting the approximations (Equation 12) into Equation 1, we obtain

which gives

However, the orthonormality property of the sequence implies that

By substituting Equation 15 shown in Equation 14, we get the Equation below:

By considering the inner product of the both sides of Equation 16 with and using the orthonormality property of the sequence , we obtain

which is a linear system of algebraic equations that can be easily solved by direct or iterative methods.

Numerical examples

In this section, we applied the method presented in this paper for solving integral equation (Equation 1) and solved some examples. The computations associated with the examples were performed in a personal computer using Mathematica 7.

Example 1. Consider the following two-dimensional Fredholm integral equation of the second kind [12]

where

and the exact solution is . Table 1 shows the absolute values of error for using the present method in selected grid points. Better approximation is expected by choosing the optimal value .

Example 2. As the second example, consider the following linear two-dimensional integral equation

where

and the exact solution Numerical results are shown in Table 2. Better approximation is expected by choosing the optimal value .

Finding exact solutions for two-dimensional integral equations is often difficult, so approximating these solutions is very important. In this work, a computational method has been presented for numerical solution of 2D-FIEs based on Haar wavelet series. In comparison with other bases (for example, polynomial bases), one of the advantages of this method is, although the involved matrices have a large dimension, they contain a large percentage of zero entries, which keeps computational effort within reasonable limits. We can modify this method for the numerical solution of linear and nonlinear two-dimensional Volterra and Fredholm integral equations in the future.

We can modify this method for the numerical solution of linear and nonlinear two- Dimensional Volterra and Fredholm integral equations in the future.

The authors declare that they have no competing interests.

HD carried out the two-dimensional wavelets studies, participated in the sequence alignment and drafted the manuscript. SS carried out the necessary programing. SS participated in the sequence alignment. HH and SS participated in the design of the study and performed the error analysis. AA conceived of the study, and participated in its design and programing. All authors read and approved the final manuscript.

The authors would like to thank both the referees for the valuable comments and Islamic Azad University-Karaj Branch for supporting this work.

1. Atkinson, KE: The Numerical Solution of Integral Equations of the Second Kind, Cambridge University Press, UK (1997)

2. Baker, CT: The Numerical Treatment of Integral Equations, Clarendon Press, Oxford (1977)

3. Carutasu, V: Numerical solution of two-dimensional nonlinear Fredholm integral equations of the second kind by spline functions. Gen. Mathematics. 9, 31 (2001)

4. Derili, H, Sohrabi, S: Numerical solution of singular integral equations using orthogonal functions. Mathematical Sci. (QJMS). 2(3), 261 (2008)

5. Hadizadeh, M, Moatamedi, N: A new differential transformation approach for two-dimensional Volterra integral equations. Int. J. Comp. Math. 84, 515 (2007). Publisher Full Text

6. Babolian, E, Shahsavaran, A: Numerical solution of nonlinear Fredholm integral equations of the second kind using Haar wavelets. J. Comp. Appl. Math. 225, 87 (2009). Publisher Full Text

7. Lepik, U, Tamme, E: Application of the Haar wavelets for solution of linear integral equations, In: Dynamical Systems and Applications, pp. 395407, 510 July, Antalya, Turkey (2005)

8. Hsiao, CH, Chen, CF: Haar wavelet method for solving lumped and distributed-parameter systems. IEE Proc. Control Theory Appl. 144, 87 (1997). Publisher Full Text

9. Wojtaszczyk, P: A Mathematical Introduction to Wavelets, Cambridge Univ Press, UK (1997)

10. Chui, CK: An Introduction to Wavelets, Wavelet Analysis and Its Applications, Academic press, Massachusetts (1992)

11. Rudin, W: Principles of Mathematical Analysis, McGraw-Hill, Singapore (1976)

12. Tari, A, Shahmorad, S: A computational method for solving two-dimensional linear Fredholm integral equations of the second kind. ANZIAM J. 49(4), 543 (2008). Publisher Full Text