The Scientific World Journal
Volume 2014 (2014), Article ID 340752, 9 pages
http://dx.doi.org/10.1155/2014/340752

Research Article

Numerical Algorithm Based on Haar-Sinc Collocation Method for Solving the Hyperbolic PDEs

A. Pirkhedri,¹ H. H. S. Javadi,² and H. R. Navidi²

¹Department of Computer Engineering, Islamic Azad University, Science and Research Branch, Tehran, Iran
²Department of Applied Mathematics, Faculty of Mathematics and Computer Science, Shahed University, Tehran, Iran

Received 13 April 2014; Revised 20 October 2014; Accepted 27 October 2014; Published 16 November 2014

Academic Editor: Zacharias Anastassi

Copyright © 2014 A. Pirkhedri et al. 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.

The present study investigates the Haar-Sinc collocation method for the solution of the hyperbolic partial telegraph equations. The advantages of this technique are that not only is the convergence rate of Sinc approximation exponential but the computational speed also is high due to the use of the Haar operational matrices. This technique is used to convert the problem to the solution of linear algebraic equations via expanding the required approximation based on the elements of Sinc functions in space and Haar functions in time with unknown coefficients. To analyze the efficiency, precision, and performance of the proposed method, we presented four examples through which our claim was confirmed.

1. Introduction

In recent years, the field of hyperbolic partial differential equations has attracted the attention of scientists in several areas and has also been used to solve many actual problems modeled in mathematical physics, such as the vibrations of structures (e.g., buildings, beams, and machines), fluid mechanics, and atomic physics [1].

The telegraph equation has typically been used for transmission and propagation of electrical signals [2], wave propagation model [3], random walk theory [4], and so forth. Let us consider an infinitesimal piece of a telegraph wire which consists of resistor of resistance , a coil of inductance , a resistor of conductance , and a capacitor of capacitance . The telegraph equation is concerned with the changes in voltage between the ends of the piece of the wire which can be formulated as follows [5]:

with initial and boundary conditions:

where denotes the voltage at position and time and is the external source term. Suppose , and ; thus, we have

Because of the lack of appropriate mathematical methods, most of the analytical solutions for linear and nonlinear partial differential equations are challenging to acquire. Therefore, approximation and numerical techniques such as Adomian decomposition method, local radial basis collocation method, variational iteration, homotopy perturbation, and Laplace transform method have been applied [6–14]. The numerical methods for solving the second-order telegraph equation are well developed. Lakestani and Saray [1] solved this problem through expanding the obtained approximate solution as the elements of interpolating scaling functions. Furthermore, the numerical approximation based on differential transform method (DTM) was considered to solve telegraph equation [5]. Using DTM, it is possible to find the exact solution or a closed approximate solution for an equation. The fully discrete local discontinuous Galerkin finite element method based on a finite difference scheme in time was introduced by [15]. Chen et al. [16] used the method of separation of variables for the analytical solution of the nonhomogeneous telegraph equation under three types of nonhomogeneous boundary conditions.

We are inspired to have an algorithm which not only is appropriate for long-time calculations but also reflects the global behavior of exact solutions. The novelty of the present paper is that we investigate the behaviour of the combination of two different groups of orthogonal functions from two different intervals. The combination of piecewise orthogonal Haar functions defined on interval, with continuous orthogonal Sinc functions defined on interval. It is essential to be said that [17–24] have previously used the Haar and Sinc functions separately for solving optimal control problems and some nonlinear ordinary differential equations. For example, Karimi and Lohmann [17] applied the Haar functions for modeling and robust control of bounce and pitch vibration for the enginebody vibration structure. Also, the Haar wavelet method has been investigated for optimal control of time-varying state-delayed [18], linear singularly perturbed systems [19], and second-order linear systems [20]. In addition, based on the properties of orthogonal Sinc functions, it is apparent that the convergence rate of approximation is exponential [21]. By using this property, the authors of [22–24] studied the Sinc collocation method for solving nonlinear singular equations like Thomas-Fermi, Lane-Emden, and Blasius equations.

Another advantage of proposed method is that it transforms the problem into a system of algebraic equations, so that the computation becomes simple and computer oriented. In the new proposed algorithm, we extend the solution of the problem to the sum of basis functions and take good advantage of the orthogonality of Haar and Sinc functions to build a set of equations for the coefficients of the solution. Although the test model is a very simple one, the proposed method is also applicable to many other problems such as fractional and two-dimensional nonlinear PDEs.

The organization of the rest of the paper is as follows. In Section 2, we present a brief introduction to the essential definitions of the Haar and Sinc functions from which are derived some tools for developing our method. In Section 3, the convergence rate analysis of the Haar and Sinc functions is given. In Section 4, we apply the method of Haar-Sinc collocation for solving the model equation. In Section 5, the proposed method is used in some types of telegraph equations and it is compared with the current analytic solutions revealed in different published works within the literature. The conclusion is presented in the final section.

2. Basic Definitions

2.1. Haar Functions

The orthogonal set of Haar functions is a number of square waves with magnitude of in some intervals and zeros elsewhere [25]. The Haar functions are defined on the interval by

where,

The value of is defined by two parameters and as

The integer indicates the level of the wavelet and the maximal level of resolution is the integer . Also, is defined for and is given by

We can expand any function in first terms of Haar functions as

where

The Haar functions coefficient vector and Haar functions vector are defined as

The matrix can be expressed as

Furthermore, the integration of the defined in (12) is given by

where is the operational matrix for integration and is given in [26] as

where , .

2.2. Sinc Functions

The Sinc function is defined on the whole real line by

For each integer and the mesh size , the Sinc basis functions are defined on by [22]

The Sinc functions form an interpolatory set of functions; that is,

If a function is defined on the real axis, then for the series,

is called the Whittaker cardinal expansion of whenever this series converges. The properties of the Whittaker cardinal expansion have been extensively studied in [27]. These properties are derived in the infinite strip of the complex -plane, where for ,

Approximations can be constructed for infinite, semiinfinite, and finite intervals. To construct approximations on the interval which is used in this paper, the eye-shaped domain in the z-plane,

is mapped conformally onto the infinite strip via

The basis functions on are taken to be the composite translated Sinc functions:

where is defined by . The inverse map of is

Thus we may define the inverse images of the real line and of the evenly spaced nodes as

Also, the nth derivative of the function at some points can be approximated [24]:

3. Convergence Rate Analysis

3.1. Haar Functions

Theorem 1. Assume that with the bounded first derivative on ; then, the error norm at th level satisfies the following inequality:

where are some real constants [28, 29].

Proof. The error at th level may be defined as

where . But , where and . Then

where and is positive constant.

3.2. Sinc Functions

The following theorem for which the proof can be found in [21] shows that the convergence rate of Sinc approximation is exponential.

Definition 2. Let be the class of functions which are analytic in (the eye-shaped domain defined in (21)) satisfy

where and the function satisfies the following equation on the boundary of :

Theorem 3. Assume that ; then, for all ,

Moreover, if for some positive constants , if the selection , then

where depends only on , , and .

4. Haar-Sinc Collocation Method

A discrete approximation to the can be expanded into Sinc functions and Haar functions as

Lemma 4. Let be Sinc collocation points, given in (26). Then the following relations hold:

where

Proof. Employing (2), (14), (27), and (37) we have

Also, using (40) we get

In addition, using (28), (29), and (41) we have

which completes our proof.

The residual for (4) can be written as

The equations for obtaining the coefficients arise from equalizing to zero at Sinc points and Haar points are defined by

By substitution collocation points in and equalizing to zero we have

Equation (45) gives linear algebraic equations which can be solved for the unknown coefficients by using the Newton’s method. Consequently, given in (4) can be calculated.

5. Illustrative Examples

In this section, we apply the proposed method for solving (4) and show the efficiency of the method with the numerical results of some examples. In all examples we choose .

Example 1. Consider the linear telegraph equation [1]:

with the following initial and boundary conditions:

The exact solution to this problem is

Table 1 shows the absolute error function obtained by the present method with and different values of .
Also, Figure 1 displays the convergence rate of our method with for . It is seen from the Figure 1 that for each fixed point the absolute errors get smaller and smaller as increases. Furthermore, we can see that the presented method provides accurate results even by using .
The maximum absolute errors for and different values of (Sinc collocation points) are shown graphically in Figure 2 for . In Figure 2, we observe that the values of maximum absolute error decay exponentially as expected from Theorem 3.

Table 1: The values of absolute error with for Example 1.

0.3
0.25

0.50

0.75

0.6
0.25

0.50

0.75

0.9
0.25

0.50

0.75

Figure 1: Plot of the absolute errors for different values of with and for Example 1.

Figure 2: Plot of the maximum absolute errors for different values of with and for Example 1.

Example 2. Consider the linear telegraph equation [5]:

with the following initial and boundary conditions:

The exact solution is given by

Table 2 shows the absolute error values using the proposed method with and .
Figure 3 displays the values of maximum absolute error with for . This figure demonstrates the validity and applicability of the present technique for this problem. Also, we observe that the convergence rate of our method for is lower than .
Figure 4 displays the convergence rate of our method for different values of for and . We see from Figure 4 that our method is in good agreement with the actual rate of convergence and the values of maximum absolute error decay exponentially.

Table 2: The values of absolute error with and for Example 2.

0.1

0.2

0.3

0.4

0.5

0.6

0.7

0.8

0.9

Figure 3: Plot of the maximum absolute errors for different values of with and for Example 2.

Figure 4: Plot of the maximum absolute errors for different values of with and for Example 2.

Example 3. Consider the nonlinear telegraph equation [5]:

with the following initial and boundary conditions:

The exact solution is given by

Figures 5, 6, 7, and 8 show the absolute error function obtained by the present method with and and . We can see clearly that better accuracy can be achieved by increasing the number of Haar collocation points and using the arbitrary precision ability of Mathematica software, we are able to establish more accurate results.

Figure 5: Plot of the absolute error with for Example 3.

Figure 6: Plot of the absolute error with for Example 3.

Figure 7: Plot of the absolute error with for Example 3.

Figure 8: Plot of the absolute error with for Example 3.

Example 4. Consider the numerical computation of the nonlinear telegraph equations [31]:

The exact solution is given by

We extract the initial and boundary conditions from the exact solution. Table 3 shows the comparison of the error between approximations obtained by the radial basis functions (IMQ, TSP) [30] and the Haar-Sinc collocation method with and . The results obtained in the table show that the error between the numerical and the exact solution can be reduced by increasing the values of .

Table 3: Results for the value of error obtained by Haar-Sinc collocation method with and for Example 4.

error [30]
error (present method)

TPS
IMQ

0.5

1.0

1.5

2.0

3.0

6. Conclusion

A numerical method for solving the telegraph equations based on the combination of two orthogonal Haar and Sinc functions was proposed. Also, [32, 33] have previously applied the combination of the Sinc functions with the other different groups of orthogonal Legendre and Chebyshev polynomials for solving similar PDEs. It is worth mentioning that since the convergence rate of the Sinc approximation is exponential, we expand the problem in space with the Sinc basis functions. Furthermore, according to the initial conditions of the problem and the Haar operational matrix for integration, the approximation in time was expanded via the elements of Haar functions. This method can be applied to solve similar problems in physics and provides a powerful alternative for physicians to investigate such types of nonlinear PDEs. The effectiveness of the method was examined via comparing the obtained results with the exact solutions. Based on the numerical results, it is obvious that it can more entirely simulate the global property of the exact solution and can provide more information about the structures of the problem. Also, the absolute error may be decreased if we take more collocation points.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

References

M. Lakestani and B. N. Saray, “Numerical solution of telegraph equation using interpolating scaling functions,” Computers & Mathematics with Applications, vol. 60, no. 7, pp. 1964–1972, 2010.

M. S. El-Azab and M. El-Gamel, “A numerical algorithm for the solution of telegraph equations,” Applied Mathematics and Computation, vol. 190, no. 1, pp. 757–764, 2007.

A. C. Metaxas and R. J. Meredith, Industrial Microwave Heating, Peter Peregrinus, London, UK, 1993.

G. Roussy and J. A. Pearcy, Foundations and Industrial Applications of Microwaves and Radio Frequency Fields, John Wiley & Sons, New York, NY, USA, 1995.

J. Biazar and M. Eslami, “Analytic solution for Telegraph equation by differential transform method,” Physics Letters A, vol. 374, no. 29, pp. 2904–2906, 2010.

S. Momani, “Analytic and approximate solutions of the space- and time-fractional telegraph equations,” Applied Mathematics and Computation, vol. 170, no. 2, pp. 1126–1134, 2005.

Siraj-ul-Islam, B. Sarler, and R. Vertnik, “Local radial basis function collocation method along with explicit time stepping for hyperbolic partial differential equations,” Applied Numerical Mathematics, vol. 67, pp. 136–151, 2013.

J. Biazar, H. Ebrahimi, and Z. Ayati, “An approximation to the solution of telegraph equation by variational iteration method,” Numerical Methods for Partial Differential Equations, vol. 25, no. 4, pp. 797–801, 2009.

A. Yıldırım, “He's homotopy perturbation method for solving the space- and time-fractional telegraph equations,” International Journal of Computer Mathematics, vol. 87, no. 13, pp. 2998–3006, 2010.

S. Kumar, “A new analytical modelling for fractional telegraph equation via Laplace transform,” Applied Mathematical Modelling, vol. 38, no. 13, pp. 3154–3163, 2014.

S. Kumar, “An analytical algorithm for nonlinear fractional Fornberg-Whitham equation arising in wave breaking based on a new iterative method,” Alexandria Engineering Journal, vol. 53, no. 1, pp. 225–231, 2014.

S. Kumar, A. Yildirim, Y. Khan, and L. Wei, “A fractional model of the diffusion equation and its analytical solution using Laplace transform,” Scientia Iranica, vol. 19, no. 4, pp. 1117–1123, 2012.

S. Kumar, H. Kocak, and A. Yildirim, “A fractional model of gas dynamics equations and its analytical approximate solution using laplace transform,” Zeitschrift fur Naturforschung, vol. 67, no. 6-7, pp. 389–396, 2012.

S. Kumar, “Numerical computation of time-fractional equation arising in solid state physics and circuit theory,” Zeitschrift für Naturforschung, vol. 68a, pp. 1–8, 2013.

L. Wei, H. Dai, D. Zhang, and Z. Si, “Fully discrete local discontinuous Galerkin method for solving the fractional telegraph equation,” Calcolo, vol. 51, no. 1, pp. 175–192, 2014.

J. Chen, F. Liu, and V. Anh, “Analytical solution for the time-fractional telegraph equation by the method of separating variables,” Journal of Mathematical Analysis and Applications, vol. 338, no. 2, pp. 1364–1377, 2008.

H. R. Karimi and B. Lohmann, “Haar wavelet-based robust optimal control for vibration reduction of vehicle engine-body system,” Electrical Engineering, vol. 89, no. 6, pp. 469–478, 2007.

H. R. Karimi, “A computational method for optimal control problem of time-varying state-delayed systems by Haar wavelets,” International Journal of Computer Mathematics, vol. 83, no. 2, pp. 235–246, 2006.

H. R. Karimi, P. Jabedar Maralani, B. Moshiri, and B. Lohmann, “Numerically efficient approximations to the optimal control of linear singularly perturbed systems based on Haar wavelets,” International Journal of Computer Mathematics, vol. 82, no. 4, pp. 495–507, 2005.

H. R. Karimi, B. Moshiri, B. Lohmann, and P. J. Maralani, “Haar wavelet-based approach for optimal control of second-order linear systems in time domain,” Journal of Dynamical and Control Systems, vol. 11, no. 2, pp. 237–252, 2005.

F. Stenger, “Integration formulae based on the trapezoidal formula,” Journal of the Institute of Mathematics and its Applications, vol. 12, pp. 103–114, 1973.

K. Parand, M. Dehghan, and A. Pirkhedri, “The sinc-collocation method for solving the Thomas-Fermi equation,” Journal of Computational and Applied Mathematics, vol. 237, no. 1, pp. 244–252, 2013.

K. Parand and A. Pirkhedri, “Sinc-collocation method for solving astrophysics equations,” New Astronomy, vol. 15, no. 6, pp. 533–537, 2010.

K. Parand, M. Dehghan, and A. Pirkhedri, “Sinc-collocation method for solving the Blasius equation,” Physics Letters A, vol. 373, no. 44, pp. 4060–4065, 2009.

C. F. Chen and C. H. Hsiao, “Haar wavelet method for solving lumped and distributed-parameter systems,” IEE Proceedings Control Theory and Applications, vol. 14, no. 1, pp. 87–94, 1997.

M. Razzaghi and Y. Ordokhani, “Solution of differential equations via rationalized Haar functions,” Mathematics and Computers in Simulation, vol. 56, no. 3, pp. 235–246, 2001.

J. Lund and K. L. Bowers, Sinc Methods for Quadrature and Differential Equations, SIAM, Philadelphia, Pa, USA, 1992.

Siraj-ul-Islam, B. Šarler, I. Aziz, and Fazal-i-Haq, “Haar wavelet collocation method for the numerical solution of boundary layer fluid flow problems,” International Journal of Thermal Sciences, vol. 50, no. 5, pp. 686–697, 2011.

Siraj-ul-Islam, I. Aziz, and B. Šarler, “The numerical solution of second-order boundary-value problems by collocation method with the Haar wavelets,” Mathematical and Computer Modelling, vol. 52, no. 9-10, pp. 1577–1590, 2010.

L. de Su, Z. W. Jiang, and T. S. Jiang, “Numerical solution for a kind of nonlinear telegraph equations using radial basis functions,” Communications in Computer and Information Science, vol. 391, pp. 140–149, 2013.

Y. D. Shang, “Explicit and exact solutions for a class of nonlinear wave equations,” Acta Mathematicae Applicatae Sinica, vol. 23, no. 1, pp. 21–30, 2000.

A. Saadatmandi, M. Dehghan, and M.-R. Azizi, “The sinc-Legendre collocation method for a class of fractional convection-diffusion equations with variable coefficients,” Communications in Nonlinear Science and Numerical Simulation, vol. 17, no. 11, pp. 4125–4136, 2012.

Z. Mao, A. Xiao, Z. Yu, and L. Shi, “Sinc-chebyshev collocation method for a class of fractional diffusion-wave equations,” The Scientific World Journal, vol. 2014, Article ID 143983, 7 pages, 2014.