# Generalized Differential Transform Method to Space-Time Fractional Telegraph Equation

Authors: Mridula Garg; Pratibha Manohar; Shyam L. Kalla
Journal: International Journal of Differential Equations
 (2011)

---

## Abstract

We use generalized differential transform method (GDTM) to derive the solution of space-time fractional telegraph equation in closed form. The space and time fractional derivatives are considered in Caputo sense and the solution is obtained in terms of Mittag-Leffler functions.


---

## Body

## 1. Introduction

Differential equations of fractional order have been successfully employed for modeling the so called anomalous phenomena during last two decades.  As a consequence, there has been an intensive development of the theory of fractional differential equations [1–4]. Recently, various analytical and numerical methods have been employed to solve linear and nonlinear fractional differential equations. A few to mention are Adomian decomposition method [5–7], homotopy perturbation method [8, 9], homotopy analysis method [10], variational iteration method [11, 12], matrix method [13], and differential transform method [14–16]. The differential transform method was proposed by Zhou [17] to solve linear and nonlinear initial value problems in electric circuit analysis. This method constructs an analytical solution in the form of a polynomial. It is different from the traditional higher order Taylor series method, which requires symbolic computation of the necessary derivatives of the data functions and takes long time in computation, whereas the differential transform is an iterative procedure for obtaining analytic Taylor series solution. The method is further developed by Momani, Odibat, and Erturk in their papers [14–16] for solving two-dimensional linear and nonlinear partial differential equations of fractional order. Recently, Biazar and Eslami [18] applied differential transform method to solve systems of Volterra integral equations of the first kind, El-Said et al. [19] developed extended Weierstrass transformation method for nonlinear evolution equations, and Keskin and Oturanc [20] developed the reduced differential transform method to solve fractional partial differential equations.In the present paper, we apply the method of generalized differential transform to solve space-time fractional telegraph equation. The classical telegraph equation is a partial differential equation with constant coefficients given by [21](1.1)utt-c2uxx+aut+bu=0,
where a, b and c are constants. This equation is used in modeling reaction diffusion and signal analysis for propagation of electrical signals in a cable of transmission line [21, 22]. Both current Iand voltage Vsatisfy an equation of the form (1.1). This equation also arises in the propagation of pressure waves in the study of pulsatile blood flow in arteries and in one-dimensional random motion of bugs along a hedge. Compared to the heat equation, the telegraph equation is found to be a superior model for describing certain fluid flow problems involving suspensions [23]. This equation is used in modeling reaction diffusion and signal analysis for transmission and propagation of electrical signals.The classical telegraph equation and space or time fractional telegraph equations have been studied by a number of researchers namely Biazar et al. [24], Cascaval et al. [25], Kaya [26], Momani [5], Odibat and Momani [27], Sevimlican [12], and Yıldırim [9]. Orsingher and Zhao [28] have shown that the law of the iterated Brownian motion and the telegraph processes with Brownian time are governed by time-fractional telegraph equations. Orsingher and Beghin [29] presented that the transition function of a symmetric process with discontinuous trajectories satisfies the space-fractional telegraph equation. Several techniques such as transform method, Adomian decomposition method, juxtaposition of transforms, generalized differential transform method, variational iteration method, and homotopy perturbation method have been used to solve space or time fractional telegraph equation.In the present paper, we make an attempt to solve homogeneous and nonhomogeneous space-time fractional telegraph equation by means of generalized differential transform method.


## 2. Preliminaries

Definition 2.1.
Caputo fractional derivative of orderα is defined as [30]:
(2.1)Daαf(x)=1Γ(m-α)∫axf(m)(ξ)(x-ξ)α-m+1dξ,(m-1<α≤m),m∈N.Definition 2.2.
The Mittag-Leffler function which is a generalization of exponential function is defined as [31]:
(2.2)Eα(z)=∑n=0∞znΓ(αn+1)(α∈C,R(α)>0).
A further generalization of (2.2) is given in the form [32]
(2.3)Eα,β(z)=∑n=0∞znΓ(αn+β);(α,β∈C,R(α)>0,R(β)>0).
For α=1,Eα(z) reduces to ez.Definition 2.3.
Generalized two-dimensional differential transform [14–16] is as given below: Consider a function of two variables u(x,y) and suppose that it can be represented as a product of two single-variable functions, that is, u(x,y)=f(x)g(y). If function u(x,y) is analytic and differentiated continuously with respect to x and y in the domain of interest, then the generalized two-dimensional differential transform of the function u(x,y) is given by
(2.4)Uα,β(k,h)=1Γ(αk+1)Γ(βh+1)[(Dx0α)k(Dy0β)hu(x,y)](x0,y0),
where 0<α, β≤1(Dx0α)k=Dx0α⋯Dx0α⋯Dx0α (k times), Dx0α is defined by (2.1) and Uα,β(k,h) is the transformed function.The generalized differential transform inverse ofUα,β(k,h) is given by(2.5)u(x,y)=∑k=0∞∑h=0∞Uα,β(k,h)(x-x0)kα(y-y0)hβ.Some basic properties of the generalized two-dimensional differential transform are as given below.LetUα,β(k,h), Vα,β(k,h) and Wα,β(k,h) be generalized two-dimensional differential transform of the functions u(x,y), v(x,y), and w(x,y), respectively, then(a)
ifu(x,y)=v(x,y)±w(x,y), then Uα,β(k,h)=Vα,β(k,h)±Wα,β(k,h),(b)
ifu(x,y)=av(x,y), a is constant, then Uα,β(k,h)=aVα,β(k,h),(c)
ifu(x,y)=Dx0γv(x,y) where m-1<γ≤m, m∈ℕ then Uα,β(k,h)=(Γ(αk+γ+1)/Γ(αk+1))Vα,β(k+γ/α,h),(d)
ifu(x,y)=Dy0μv(x,y) where n-1<μ≤n, n∈ℕ, then Uα,β(k,h)=(Γ(βh+μ+1)/Γ(βh+1))Vα,β(k,h+μ/β).


## 3. Solution of Space-Time Fractional Telegraph Equations by Generalized Two-Dimensional Differential Transform Method

In this section, we consider space-time fractional telegraph equations in the following form:(3.1)c2Dx2αu(x,t)=Dtpβu(x,t)+aDtrβu(x,t)+bu(x,t)+f(x,t),0<x<1,t>0,
where β=1/q, p,q,r∈ℕ, 1<2α≤2, 1<pβ≤2, 0<rβ≤1, Dx2α≡DxαDxα, Dtpβ≡DtβDtβ⋯Dtβ (p times), Dtrβ≡DtβDtβ⋯Dtβ (r times), Dxα, Dtβ are Caputo fractional derivatives defined by (2.1), a, b, and c are constants, f(x,t) is given function.Particularly forα=1, q=1, p=2, r=1, f=0, space-time fractional telegraph (3.1) reduces to classical telegraph (1.1).To give a clear overview of the methodology, we have selected three illustrative examples, the first is a homogeneous space-time fractional telegraph equation with conditions involving ordinary derivative with respect to space, the second is a homogeneous space-time fractional telegraph equation with conditions involving fractional derivative with respect to space, and the third is a nonhomogeneous space-time fractional telegraph equation with conditions involving fractional derivative with respect to space.Example 3.1.
Consider the following homogeneous space-time fractional telegraph equation:(3.2)Dx3/2u(x,t)=Dtpβu(x,t)+Dtrβu(x,t)+u(x,t),0<x<1,t>0,
where β=1/q, p,q,r∈ℕ, 1<pβ≤2, 0<rβ≤1, Dx3/2≡(Dx1/2)3, Dtpβ≡DtβDtβ⋯Dtβ (p times), Dtrβ≡DtβDtβ⋯Dtβ (r times), Dxα, Dtβ are Caputo fractional derivatives defined by (2.1), p+r is odd and
(3.3)u(0,t)=Eβ(-tβ),ux(0,t)=Eβ(-tβ).
Applying generalized two-dimensional differential transform (2.4) with x0=0=y0, α=1/2, to both sides of (3.2) and (3.3) and using properties (c) and (d), we obtain
(3.4)U1/2,β(k+3,h)=Γ((k/2)+1)Γ(((k+3)/2)+1)[Γ(β(h+p)+1)Γ(βh+1)U1/2,β(k,h+p)+Γ(β(h+r)+1)Γ(βh+1)U1/2,β(k,h+r)+U1/2,β(k,h)],(3.5)U1/2,β(0,h)=(-1)hΓ(βh+1),U1/2,β(1,h)=0,U1/2,β(2,h)=(-1)hΓ(βh+1),h=0,1,2,….
Utilizing recurrence relation (3.4), the transformed conditions (3.5) and the condition p+r is odd, we can easily obtain, for l,h=0,1,2,…(3.6)U1/2,β(3l,h)=(-1)hΓ((3/2)l+1)Γ(βh+1),(3.7)U1/2,β(3l+1,h)=0,(3.8)U1/2,β(3l+2,h)=(-1)hΓ((3/2)l+2)Γ(βh+1).
Now, from (2.5), we have
(3.9)u(x,t)=∑k=0∞∑h=0∞U1/2,β(k,h)xk/2tβh.
Using the values of U1/2,β(k,h) from (3.5)–(3.8) in (3.9), the exact solution of space-time fractional telegraph (3.2) is obtained as
(3.10)u(x,t)=[E3/2(x3/2)+xE3/2,2(x3/2)]Eβ(-tβ),
which is same as obtained by Garg and Sharma [33] using Adomian decomposition method.Settingp=2, q=r=1, the space-time fractional telegraph (3.2) reduces to space fractional telegraph equation and the solution is same as obtained by Momani [5], Odibat and Momani [27], and Yıldırim [9] using Adomian decomposition method, generalized differential transform method, and homotopy perturbation method, respectively.Further, settingα=1, it reduces to classical telegraph equation and the solution is same as obtained by Kaya [26] using Adomian decomposition method.Example 3.2.
Consider the following homogeneous space-time fractional telegraph equation:(3.11)Dx2αu(x,t)=Dtpβu(x,t)+Dtrβu(x,t)+u(x,t),0<x<1,t>0,
where β=1/q, p,q,r∈ℕ, 1<2α≤2, 1<pβ≤2, 0<rβ≤1, Dx2α≡DxαDxα, Dtpβ≡DtβDtβ⋯Dtβ (p times), Dtrβ≡DtβDtβ⋯Dtβ (r times), Dxα, Dtβ are Caputo fractional derivatives defined by (2.1), p+r is odd and
(3.12)u(0,t)=Eβ(-tβ),Dxαu(x,t)∣x=0=Eβ(-tβ).
Applying generalized two-dimensional differential transform (2.4) with x0=0=y0 to both sides of (3.11), (3.12) and using properties (c) and (d) we obtain
(3.13)Uα,β(k+2,h)=Γ(αk+1)Γ(α(k+2)+1)[Γ(β(h+p)+1)Γ(βh+1)Uα,β(k,h+p)+Γ(β(h+r)+1)Γ(βh+1)Uα,β(k,h+r)+Uα,β(k,h)],(3.14)Uα,β(0,h)=(-1)hΓ(βh+1),Uα,β(1,h)=(-1)hΓ(α+1)Γ(βh+1),h=0,1,2,….
Utilizing the recurrence relation (3.13), the transformed conditions (3.14) and the condition p+r is odd, we obtain
(3.15)Uα,β(k,h)=(-1)hΓ(kα+1)Γ(hβ+1),fork,h=0,1,2,….
Now, from (2.5), we have
(3.16)u(x,t)=∑k=0∞∑h=0∞Uα,β(k,h)xαktβh.
Using the values of Uα,β(k,h) from (3.15) in (3.16), the exact solution of homogeneous space-time fractional telegraph (3.11) is obtained as
(3.17)u(x,t)=Eα(xα)Eβ(-tβ).Remark 3.3.
(1) Settingq=1, p=2, r=1, (3.11) reduces to space fractional telegraph equation
(3.18)Dx2αu(x,t)=Dt2u(x,t)+Dtu(x,t)+u(x,t),0<x<1,t>0,
with solution
(3.19)u(x,t)=Eα(xα)e-t.
(2) Settingα=1, (3.11) reduces to time fractional telegraph equation:
(3.20)Dx2u(x,t)=Dtpβu(x,t)+Dtrβu(x,t)+u(x,t),0<x<1,t>0,
with solution
(3.21)u(x,t)=exEβ(-tβ).
(3) Settingα=1, q=1, p=2, r=1, (3.11) reduces to classical telegraph equation:
(3.22)Dx2u(x,t)=Dt2u(x,t)+Dtu(x,t)+u(x,t),0<x<1,t>0,
with solution
(3.23)u(x,t)=ex-t,
which is same as obtained by Kaya [26] using Adomian decomposition method.Example 3.4.
Consider the following non-homogeneous space-time fractional telegraph equation:(3.24)Dx2αu(x,t)=Dtpβu(x,t)+Dtrβu(x,t)+u(x,t)-2Eα(xα)Eβ(-tβ),0<x<1,t>0,
where β=1/q, p,q,r∈ℕ, 1<2α≤2, 1<pβ≤2, 0<rβ≤1, Dx2α≡DxαDxα, Dtpβ≡DtβDtβ⋯Dtβ (p times), Dtrβ≡DtβDtβ⋯Dtβ (r times), Dxα, Dtβ are Caputo fractional derivatives defined by (2.1), p and r are even and
(3.25)u(0,t)=Eβ(-tβ),Dxαu(x,t)∣x=0=Eβ(-tβ).
Applying generalized two-dimensional differential transform (2.4) with x0=0=y0 to both sides of (3.24), (3.25), and using properties (c) and (d) we obtain
(3.26)Uα,β(k+2,h)=Γ(αk+1)Γ(α(k+2)+1)[Γ(β(h+p)+1)Γ(βh+1)Uα,β(k,h+p)+Γ(β(h+r)+1)Γ(βh+1)Uα,β(k,h+r)+Uα,β(k,h)-2(-1)hΓ(αk+1)Γ(βh+1)],(3.27)Uα,β(0,h)=(-1)hΓ(βh+1),Uα,β(1,h)=(-1)hΓ(α+1)Γ(βh+1),h=0,1,2,….
Utilizing the recurrence relation (3.26) and the transformed conditions (3.27), we obtain
(3.28)Uα,β(k,h)=(-1)hΓ(kα+1)Γ(hβ+1),k,h=0,1,2,….
Now from (2.5), we have
(3.29)u(x,t)=∑k=0∞∑h=0∞Uα,β(k,h)xαktβh.
Using the values of Uα,β(k,h) from (3.28) in (3.29), the exact solution of non-homogeneous space-time fractional telegraph (3.24) is obtained as
(3.30)u(x,t)=Eα(xα)Eβ(-tβ).Remark 3.5.
(1) Settingq=2, p=4, r=2, (3.24) reduces to non-homogeneous space fractional telegraph equation:
(3.31)Dx2αu(x,t)=Dt2u(x,t)+Dtu(x,t)+u(x,t)-2Eα(xα)e-t,0<x<1,t>0,
with solution
(3.32)u(x,t)=Eα(xα)e-t.
(2) Settingα=1, (3.24) reduces to non-homogeneous time fractional telegraph equation:
(3.33)Dx2u(x,t)=Dtpβu(x,t)+Dtrβu(x,t)+u(x,t)-2exEβ(-tβ),0<x<1,t>0,
with solution
(3.34)u(x,t)=exEβ(-tβ).
(3) Settingα=1, q=2, p=4, r=2, (3.24) reduces to non-homogeneous telegraph equation:
(3.35)Dx2u(x,t)=Dt2u(x,t)+Dtu(x,t)+u(x,t)-2exE1/2(-t1/2),0<x<1,t>0,
with solution
(3.36)u(x,t)=exE1/2(-t1/2).



---

*Source: 548982-2011-10-11.xml*