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

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.


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: (0, ) = 0 ( ) , (1, ) = 1 ( ) , where ( , ) denotes the voltage at position and time and ( , ) is the external source term. Suppose = / , = / and 2 = 1/ ; 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][7][8][9][10][11][12][13][14]. The Scientific World Journal 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 [0, 1) interval, with continuous orthogonal Sinc functions defined on (−∞, +∞) interval. It is essential to be said that [17][18][19][20][21][22][23][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][23][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.

Haar Functions.
The orthogonal set of Haar functions is a number of square waves with magnitude of ±1 in some intervals and zeros elsewhere [25]. The Haar functions are defined on the interval = [0, ) 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, ℎ 0 ( ) is defined for = = 0 and is given by We can expand any function ( ) ∈ 2 (0, 1) in first terms of Haar functions as where = 2 +1 , = 0, 1, 2, . . . .
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 The Scientific World Journal 3

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 R by [22] The Sinc functions form an interpolatory set of functions; that is, If a function ( ) is defined on the real axis, then for ℎ > 0 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 > 0, Approximations can be constructed for infinite, semiinfinite, and finite intervals. To construct approximations on the interval (0, 1) 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 (0, 1) 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]:
Proof. The error at th level may be defined as The Scientific World Journal where ( ) ≤ ∀ ∈ (0, 1) and is positive constant.

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

Lemma 4. Let be Sinc collocation points, given in
, .
The residual , ( , ) for (4) can be written as

Illustrative Examples
In this section, we apply the proposed method for solving (4) and show the efficiency of the method with the numerical with the following initial and boundary conditions: The exact solution to this problem is ( , ) = −2 sinh ( ) .
(48) Table 1 shows the absolute error function | exact ( , ) − , ( , )| obtained by the present method with = 3 and different values of . Also, Figure 1 displays the convergence rate of our method with = 4, 8, 16, 32 for = 0.5. 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 = 3.
The maximum absolute errors for = 8 and different values of (Sinc collocation points) are shown graphically in Figure 2 for = 0.5. In Figure 2, we observe that the values of maximum absolute error decay exponentially as expected from Theorem 3.
with the following initial and boundary conditions: The Scientific World Journal  The exact solution is given by Table 2 shows the absolute error values using the proposed method with = 3 and = 8, 16, 32. Figure 3 displays the values of maximum absolute error with = 4, 8, 16, 32 for = 0.5, 0.75. 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 = 0.75 is lower than = 0.5. Figure 4 displays the convergence rate of our method for different values of for = 0.75 and = 8. 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.
with the following initial and boundary conditions: The exact solution is given by  Mathematica software, we are able to establish more accurate results.    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 = 3 and = 8, 16, 32. 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 .

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.