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.
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.
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.
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 .
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.