A New Algorithm for the Numerical Solution of Telegraph Equations by Using Fibonacci Polynomials

Abstract: In this study, we present a numerical scheme to solve the telegraph equation by using Fibonacci polynomials. This method is based on the Fibonacci collocation method which transforms the equation into a matrix equation, and the unknown of this equation is a Fibonacci coefficients matrix. Some numerical examples with comparisons are included to demonstrate the validity and applicability of the proposed method. The results show the efficiency and accuracy of this paper.


Introduction
Linear partial differential equations play an important role in the modelling of engineering problems and in many areas of science.Especially, hyperbolic type partial differential equations (PDEs) is a significant class of the PDEs because of a wide range of applications in many fields.
The hyperbolic type PDEs model atomic physics, aerospace, industry, biology and engineering problems such as vibrations of structures, beams and buildings [1,2].
In this paper, we show a new algorithm for the numerical solutions of the Telegraph equation which is a typical example of the hyperbolic type PDEs.This equation is commonly used in signal analysis for transmission and propagation of electrical signals and also has applications in other fields.Telegraph equation models mix diffusion and wave propagation by introducing a term that accounts for effects of finite velocity to standard heat or mass transport equation [2][3][4][5][6].
We consider the second-order linear hyperbolic partial differential Telegraph equation.
In this study, we introduce Fibonacci collocation method based on matrix relations which has been used to find the approximate solutions of differential, integral and integro-differential equations, differential-difference equations, Fredholm integro differential-difference equations, Pantograph-type functional differential equations and linear Volterra integro differential equations [11][12][13][14][15].Because Fibonacci polynomials are orthogonal polynomials and the collocation method is easily programmable by symbolic computation programs, this method is preferable [15].Also, recently, several researchers have used other collocation methods to solve different types of the partial differential equations by using various special polynomials without Fibonacci polynomials [16][17][18][19][20]. On the other hand, these matrix methods were used for the numerical analysis of the longitudinal vibration of rods [21].In addition to this, Shaban et al. [22] use the homotopy analysis method based on the operational matrix with Chebyshev polynomials.
In this paper, by developing the Fibonacci collocation method studied in [11][12][13][14][15], we will obtain the approximate solution of Equation ( 1) in the truncated Fibonacci series form.
where a mn ; m, n " 1, . . ., N are the unknown Fibonacci coefficients and F m pxq , n " 1, . . ., N are the Fibonacci functions defined by,

Fundamental Relations
To obtain the numerical solution of hyperbolic partial differential Telegraph equation with Fibonacci Collocation Method, we first evaluate the Fibonacci coefficients of the unknown function.For convenience, the solution function Equation (4) can be written in the matrix form where F pxq "

¨¨¨.
. .In there, the matrix form F pxq can be written as So that  C is the characteristic matrix of the matrix relations [11]. where

Method of Solution
We can compute the matrix forms of the terms u px, tq , u x px, tq , u t px, tq , u xx px, tq , u tt px, tq. of Equation (1) now.From the relations Equations ( 5)- (7) we obtain the matrix forms.u px, tq " X pxq C T X ptq C T A (8) where where all 0 matrix is N ˆN.
Shortly, Equation ( 13) can be written as where W px, tq " " w ij ‰ , i, j " 1, 2, . . ., N 2 .By using in Equation ( 14) the collocation points `xi , t j ˘defined by So, the system of the matrix equations can be written as Or, in summary, the fundamental matrix equation is By means of the relation Equations ( 8) and (10) for the conditions Equation ( 2), we can obtain the matrix forms.
By substituting the collocation points Equation (15) into Equations ( 17) and ( 18), we have the fundamental matrix equations of the initial conditions.
Similarly, by using Equation ( 8) for the conditions Equation (3), we can have the matrix forms as follow u p0, tq " On the other hand, by substituting the collocation points Equation ( 16) into Equations ( 19) and ( 20), we obtain the fundamental matrix equations of the Dirichlet boundary conditions To obtain the approximate solution of Equation ( 1) under conditions Equations ( 2) and (3), the augmented matrix is formed as follows So, the unknown Fibonacci coefficients are obtained as ı is generated by using the Gauss elimination method and then removing zero rows of the augmented matrix r | W ; q G s.

Error Estimation Algorithm
In this section, we introduce an error estimation algorithm for obtaining the Fibonacci polynomial solution giving last sections and use this technique to obtain the corrected solution of the problem.This algorithm is based on residual function [23][24][25], so we define residual function as follows: where From Equations ( 21) and ( 22), with the initial conditions u N px, 0q " k 1 pxq , 0 ă x ă l pu N q t px, 0q " k 2 pxq , 0 ă x ă l (24) and Dirichlet boundary conditions Also, the error function e N px, tq is defined as follows: e N px, tq " u px, tq ´uN px, tq Because of the fact that L is linear operator, error problem can be defined by using error function as with homogenous conditions e N px, 0q " 0, 0 ă x ă l pe N q t px, 0q " 0, 0 ă x ă l e N p0, tq " 0, 0 ă t ď τ e N pl, tq " 0, 0 ă t ď τ (28) The problem (27) under conditions (28) was solved by Fibonacci collocation method given by Section 3. When we solve the Equation (28) for M ą N which is a truncated number, we get the approximation e N,M px, tq to e N px, tq even if the exact solution of the problem is unknown.
Consequently, by means of the error estimation function e N,M px, tq, we obtain the corrected Fibonacci polynomial solution.

Numerical Examples
In this section, we will give several examples to illustrate the applicability of the Fibonacci matrix method and all of them are performed on the computer.Maximum absolute error L 8 is defined as with the initial conditions u px, 0q " x 2 , 0 ă x ă 1 u t px, 0q " 1, 0 ă x ă 1 and Dirichlet boundary conditions u p0, tq " t, 0 ă t ď 1 u p1, tq " 1 `t, 0 ă t ď 1.
Here l " 1, τ " 1, α " 1, β " 1, k 1 pxq " x 2 , k 2 pxq " 1, g px, tq " x 2 `t ´1.By using the solution algorithm in given Section 3, we find the Fibonacci polynomial solution of the problem is u 5 px, tq " x 2 `t, which is the exact solution of the problem.Further, we can obtain the exact solution of the problem for the any value of N.
The exact solution is given by u px, tq " e p´tq sin pxq .
g px, tq " p2 ´α `βq e p´tq sin pxq , h 1 ptq " h 2 ptq " 0 In Table 1 shows a comparison of the absolute errors values of Fibonacci polynomial solution and corrected Fibonacci polynomial solution with Taylor Matrix Method and Legendre Multiwavelet Method.Also, Figure 1 clearly shows when the N values increase, and when the absolute error function values rapidly decrease.
g px, tq " p2 ´2t `t2 q ´x ´x2 ¯ep´tq 2t   Exact solution and numerical solution for  = 11 are plotted in Figure 2. On the other hand, taking  = 7, 11, 14, we compared the absolute error functions in Figure 3 and Table 2. Exact solution and numerical solution for N " 11 are plotted in Figure 2. On the other hand, taking N " 7, 11, 14, we compared the absolute error functions in Figure 3 and Table 2.

Conclusions
In this article, we have presented a new approach using the Fibonacci polynomials to solve the Telegraph equation which is a typical example of the hyperbolic type PDEs.The solution of this equation is expressed as a truncated series of Fibonacci polynomials, which are then transformed from algebraic form to matrix form.The new method has been applied to three illustrative examples.Illustrated examples are contained to demonstrate the applicability and the validity of the method and the obtained results are compared with the different methods in literature.The results show the efficiently and the accuracy of the present work.Moreover, it is well seen from these examples that the method yields either the exact solution or a high accuracy approximate solution for the telegraph equations.Also, as N increases, the errors decrease more rapidly.One of the considerable advantages of this technique is that the approximate solutions are found very easily by using the computer programs.As seen from this study, the new approach with Fibonacci polynomial is a good approximation for solving the Telegraph equations.

Figure 2 .
Figure 2. Exact solution and Fibonacci polynomial solution of Example 3 for N " 11.(a) Exact solution; (b) Present method for N " 11.

Figure 3 .
Figure 3. Graphs of the absolute errors functions of Example 3 for  = 7, 11, 14.(a) The absolute error function for  = 7; The absolute error function for  = 11; The absolute error function for  = 14.

Figure 3 .
Figure 3. Graphs of the absolute errors functions of Example 3 for N " 7, 11, 14.(a) The absolute error function for N = 7; (b) The absolute error function for N = 11; (c) The absolute error function for N = 14.

Table 1 .
Comparison of the absolute errors values of Fibonacci polynomial solution and corrected Fibonacci polynomial solution with Taylor Matrix Method and Legendre Multiwavelet Method for Example 1 and α " 4, β " 2.

Table 3
shows a comparison of the  ∞ values of Fibonacci polynomial solution and corrected Fibonacci polynomial solution with B-spline Collocation Method and Tau Method for Example 3.

Table 3 .
Comparison of the  ∞ values of Fibonacci polynomial solution and corrected Fibonacci polynomial solution with B-spline Collocation Method and Tau Method for Example 3.

Table 2 .
Comparison of the absolute error values E N px r , 1q of Example 3 for N " 7, 11, 14.

Table 3
shows a comparison of the L 8 values of Fibonacci polynomial solution and corrected Fibonacci polynomial solution with B-spline Collocation Method and Tau Method for Example 3.

Table 3 .
Comparison of the L 8 values of Fibonacci polynomial solution and corrected Fibonacci polynomial solution with B-spline Collocation Method and Tau Method for Example 3.