Table of Contents Author Guidelines Submit a Manuscript
Mathematical Problems in Engineering
Volume 2016 (2016), Article ID 1683849, 6 pages
http://dx.doi.org/10.1155/2016/1683849
Research Article

Numerical Solution of Time-Fractional Order Telegraph Equation by Bernstein Polynomials Operational Matrices

1Department of Mathematics, Islamic Azad University, Karaj Branch, Karaj, Iran
2Department of Mathematics, Islamic Azad University, Science and Research Branch, Tehran, Iran

Received 27 May 2016; Revised 5 September 2016; Accepted 18 October 2016

Academic Editor: Zhen-Lai Han

Copyright © 2016 M. Asgari 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

We present a new method to solve time-fractional order telegraph equation (TFOTE) by using Bernstein polynomials. By implementation of Bernstein polynomials operational matrices of fractional differential on TFOTE, we reduce the original problem to a linear system of algebraic equations. Also, we prove the convergence analysis. In order to show the efficiency of the proposed method, we present two numerical examples.

1. Introduction

Telegraph equations are hyperbolic partial differential equations that are applicable in modeling the reaction diffusion processes in various branches of engineering sciences and biological sciences. Those equations frequently arise in the study of wave propagation of electrical signal in a cable of transmission line and wave phenomena [14].

Many authors have used various numerical and analytical methods to solve the TFOTE. Chen and coworkers derived the analytical solution of the nonhomogeneous TFOTE by method of separation of variables [5]. Huang presented a new analytical solution for three basic problems of time-fractional telegraph equation. He solved Cauchy and signaling problems by Laplace and Fourier transforms and the boundary problem by spatial Sine transform [6]. Dehghan and Shokri developed a numerical method to solve the one-dimensional hyperbolic telegraph equation using the collocation points and approximated the solution by using thin plate spline radial basic functions [7]. Saadatmandi and Dehghan developed a numerical solution based on Chebyshev tau method [4]. Yousefi in [8] used Legendre multiwavelet Galerkin method for solving the hyperbolic telegraph equation. In [9], Das and Gupta used homotopy analysis method for solving fractional hyperbolic partial differential equation. In [10], Mollahasani et al. applied hybrid functions of Legendre polynomials and block pulse functions to obtain the solution of telegraph equation of fractional order.

In this paper, our study focuses on the time-fractional telegraph equation of order :with the initial and boundary conditions: The right-hand-side function is given, , and also , , and are arbitrary positive constants. If , we have one-dimensional hyperbolic telegraph equation.

The rest of this paper is organized as follows: First, we present some preliminaries in fractional calculus. In Section 3, we briefly review some general concepts concerning Bernstein polynomials and the Bernstein polynomials operational matrix for fractional derivative. In Section 4, the method is applied to solve linear TFOTE. Section 5 exhibits an error estimation for the presented method. Section 6 illustrates two numerical examples to show the convergence and accuracy of the proposed method.

2. Basic Definitions

In this section, we present some basic definitions and properties of the fractional calculus which are going to be used in this paper.

Definition 1 (see [11]). The Riemann-Liouville fractional integral operator of order on the Lebesgue space is given by where denotes Euler Gamma function.

Definition 2 (see [11]). The Caputo fractional derivative of order is defined by where and is the smallest integer greater than  . For the Caputo derivative, we have We use the ceiling function to denote the smallest integer greater than or equal to and . If , the Caputo differential operator coincides with the usual differential operator of an integer order.

For , , we have the following properties:

3. Bernstein Polynomials and Their Properties

3.1. The Definition of Bernstein Polynomials Basis

The Bernstein polynomials (BPs) of degree on the interval are defined by These polynomials have the following properties on [12]:(1), , ,(2),(3), ,(4), .

Theorem 3 (see [13]). Suppose that is a Hilbert space with the inner product and is a finite dimensional and closed subspace; therefore, is a complete subspace of . So, if is an arbitrary element in , it has a unique best approximation out of such as ; that is, where and . So, there exist unique coefficients such that where and

Lemma 4 (see [14]). If is a complete basis, then , where is an upper triangular matrix with for and .

3.2. Function Approximation

A function can be expressed in terms of the Bernstein polynomials basis as where with

3.3. Operational Matrix for Fractional Derivative

Theorem 5. Let be vector defined in (10); then, where , , and are matrices that is defined in (11); is a diagonal matrix with is called the Bernstein polynomials operational matrix of fractional derivative [14].

4. Description of the Method

Clearly, can be approximated by using Bernstein polynomials as where and are vectors defined in (10) and the unknown is matrix that can be shown as Now, we can write Substituting (19) and (21) into (1), we have Now, we collocate (22) in Newton-Cotes nodes as So, we have equations as Applying (17) and (20) in the initial and boundary conditions (2), we getBy collocating (25) in Newton-Cotes nodes and for and , we get equations. These equations together with (24) give equations, which can be solved for , . So, the unknown function can be approximated.

5. Convergence Analysis

Theorem 6. Suppose that is a continuous function and all partial derivatives of exist and are continuous. Let , . If is the best approximation for out of and also all th-order partial derivatives of are bounded in magnitude by , then the error bound is presented as follows: where and .

Proof. By applying the Taylor expansion in two variables for , we have Hence, from Taylor expansion, we have Since is the best approximation for out of , we conclude that Clearly, we obtain the following result:

Remark 7. Equation (26) shows that if , then

6. Numerical Examples

To demonstrate the validity and applicability of the numerical scheme, we apply the present method for the following illustrative examples.

Example 1. Consider the time-fractional telegraph equation of order : The exact solution of this example is . In Table 1, we compare the obtained numerical results with the method of interpolating scaling functions operational matrix [15] (ISF), for . Also, in this table, we present the maximum absolute errors in for and various values of .

Table 1: The maximum absolute errors in Example 1.

Example 2. Consider the time-fractional telegraph equation of order :Table 2 shows the absolute errors using the proposed method in the previous section for and different values of . The exact solution of this equation is .

Table 2: The maximum absolute errors for .

7. Conclusion

This paper proposed a numerical approach for solving TFOTE by using the operational matrices of Bernstein polynomials. The operational matrix of fractional derivative, , as well as collocation method was used to transform the TFOTE to a linear system of algebraic equations. The numerical results show that the proposed method in this paper can be a suitable method for solving these equations. To obtain the numerical results, we applied Mathematica 9 software.

Competing Interests

The authors declare that they have no competing interests.

References

  1. J. Banasiak and J. R. Mika, “Singularly perturbed telegraph equations with applications in the random walk theory,” Journal of Applied Mathematics and Stochastic Analysis, vol. 11, no. 1, pp. 9–28, 1998. View at Publisher · View at Google Scholar · View at MathSciNet · View at Scopus
  2. V. H. Weston and S. He, “Wave splitting of the telegraph equation in R3 and its application to inverse scattering,” Inverse Problems, vol. 9, no. 6, pp. 789–812, 1993. View at Publisher · View at Google Scholar
  3. P. M. Jordan and A. Puri, “Digital signal propagation in dispersive media,” Journal of Applied Physics, vol. 85, no. 3, pp. 1273–1282, 1999. View at Publisher · View at Google Scholar · View at Scopus
  4. A. Saadatmandi and M. Dehghan, “Numerical solution of hyperbolic telegraph equation using the Chebyshev tau method,” Numerical Methods for Partial Differential Equations, vol. 26, no. 1, pp. 239–252, 2010. View at Publisher · View at Google Scholar · View at MathSciNet · View at Scopus
  5. 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. View at Publisher · View at Google Scholar · View at MathSciNet · View at Scopus
  6. F. Huang, “Analytical solution for the time-fractional telegraph equation,” Journal of Applied Mathematics, vol. 2009, Article ID 890158, 9 pages, 2009. View at Publisher · View at Google Scholar · View at MathSciNet · View at Scopus
  7. M. Dehghan and A. Shokri, “A numerical method for solving the hyperbolic telegraph equation,” Numerical Methods for Partial Differential Equations, vol. 24, no. 4, pp. 1080–1093, 2008. View at Publisher · View at Google Scholar · View at MathSciNet · View at Scopus
  8. S. A. Yousefi, “Legendre multiwavelet Galerkin method for solving the hyperbolic telegraph equation,” Numerical Methods for Partial Differential Equations, vol. 26, no. 3, pp. 535–543, 2010. View at Publisher · View at Google Scholar · View at MathSciNet · View at Scopus
  9. S. Das and P. K. Gupta, “Homotopy analysis method for solving fractional hyperbolic partial differential equations,” International Journal of Computer Mathematics, vol. 88, no. 3, pp. 578–588, 2011. View at Publisher · View at Google Scholar · View at MathSciNet · View at Scopus
  10. N. Mollahasani, M. Mohseni Moghadam, and K. Afrooz, “A new treatment based on hybrid functions to the solution of telegraph equations of fractional order,” Applied Mathematical Modelling, vol. 40, no. 4, pp. 2804–2814, 2015. View at Publisher · View at Google Scholar
  11. I. Podlubeny, Fractional Differential Equations, Academic Press, San Diego, Calif, USA, 1990.
  12. H. Jafari and H. Tajadodi, “Fractional order optimal control problems via the operational matrices of Bernstein polynomials,” UPB Scientific Bulletin, Series A, vol. 76, no. 3, pp. 115–128, 2014. View at Google Scholar
  13. E. Kreyszig, Introductory Functional Analysis with Applications, John Wiley & Sons, New York, NY, USA, 1978. View at MathSciNet
  14. M. Alipour, D. Rostamy, and D. Baleanu, “Solving multidimensional fractional optimal control problems with inequality constraint by Bernstein polynomials operational matrices,” Journal of Vibration and Control, vol. 19, no. 16, pp. 2523–2540, 2012. View at Google Scholar
  15. 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. View at Publisher · View at Google Scholar · View at MathSciNet · View at Scopus