JAMJournal of Applied Mathematics1687-00421110-757XHindawi Publishing Corporation15418910.1155/2010/154189154189Research ArticleA Note on Fractional Sumudu TransformGuptaV. G.1ShramaBhavna2KiliçmanAdem3GeurtsBernard1Department of MathematicsUniversity of RajasthanJaipur 302055Indiauniraj.ac.in2Department of MathematicsShri Balaji Institute of Engineering and TechnologyJaipur 302045Indiashantiniketanss.com3Department of Mathematics and Institute for Mathematical ResearchUniversiti Putra Malaysia(UPM), Serdang, Selangor 43400Malaysiaupm.edu.my20102605201020101301201025032010210420102010Copyright © 2010This 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.

We propose a new definition of a fractional-order Sumudu transform for fractional differentiable functions. In the development of the definition we use fractional analysis based on the modified Riemann-Liouville derivative that we name the fractional Sumudu transform. We also established a relationship between fractional Laplace and Sumudu duality with complex inversion formula for fractional Sumudu transform and apply new definition to solve fractional differential equations.

1. Introduction

In the literature there are numerous integral transforms that are widely used in physics, astronomy, as well as engineering. In order to solve the differential equations, the integral transforms were extensively used and thus there are several works on the theory and application of integral transforms such as the Laplace, Fourier, Mellin, and Hankel, to name but a few. In the sequence of these transforms in early 90s, Watugala  introduced a new integral transforms named the Sumudu transform and further applied it to the solution of ordinary differential equation in control engineering problems. For further detail and properties about Sumudu transforms see  and many others. Recently Kiliçman et al. applied this transform to solve the system of differential equations; see . The Sumudu transform is defined over the set of the functionsA={f(t):M,τ1,τ2>0,|f(t)|<Met/τj,if  t(-1)j×[0,)} by the following formula:G(u)=S[f(t);u]=:0f(ut)e-tdt,u(-τ1,τ2). The existence and the uniqueness were discussed in ; for further details and properties of the Sumudu transform and its derivatives we refer to . In , some fundamental properties of the Sumudu transform were established. In , this new transform was applied to the one-dimensional neutron transport equation. In fact, one can easily show that there is a strong relationship between double Sumudu and double Laplace transforms; see . Further in , the Sumudu transform was extended to the distributions and some of their properties were also studied in .

The function f(t) so involved is usually continuous and continuously differentiable. Suppose that the function is continuous but its fractional derivative exists of order α,0<α<1, but no derivative, and then (1.2) fails to apply. Thus we have to introduce a new definition of Sumudu transform. For the convenience of the reader, firstly we will give a brief background on the definition of the fractional derivative and basic notations for more details see  and .

1.1. Fractional Derivative via Fractional DifferenceDefinition 1.1.

Let f:, tf(t) denote a continuous (but not necessarily differentiable) function, and let h>0 denote a constant discretization span. Define the forward operator FW(h) by the equality FW(h)f(t):=f(t+h). Then the fractional difference of order α, 0<α<1 of f(t) is defined by the expression Δαf(t):=(FW-1)α=k=0(-1)k(αk)f[t+(α-k)h], and its fractional derivative of order α is defined by the limit f(α)(t)=limh0Δαf(t)hα. See the details in .

1.2. Modified Fractional Riemann-Liouville Derivative

Jumarie proposed an alternative way to the Riemann-Liouville definition of the fractional derivative; see .

Definition 1.2.

Let f: be a continuous but not necessarily differentiable function. Further, consider the following.

Assume that f(t) is a constant K. Then its fractional derivative of order α is DtαK=KΓ-1(1-α)t-α,α0,=0,α>0.

When f(t) is not a constant, then we will set f(t)=f(0)+(f(t)-f(0)), and its fractional derivative will be defined by the expression f(α)(t)=Dtαf(0)+Dtα(f(t)-f(0)), in which, for negative α, one has Dtα(f(t)-f(0)):=1Γ(-α)ot(t-ξ)-α-1f(ξ)dξ,α<0, whilst for positive α, we will set Dtα(f(t)-f(0))=Dtαf(t)=Dt(fα-1(t)). When nα<n+1, we will set f(α)(t):=(f(α-n)(t))(n),nα<n+1,n1.

We will refer to this fractional derivative as the modified Riemann-Liouville derivative, and it is in order to point out that this definition is strictly equivalent to Definition 1.1, via (1.4).

1.3. Integration with respect to <inline-formula><mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M30"><mml:mo stretchy="false">(</mml:mo><mml:mi>d</mml:mi><mml:mi>t</mml:mi><mml:msup><mml:mrow><mml:mo stretchy="false">)</mml:mo></mml:mrow><mml:mrow><mml:mi>α</mml:mi></mml:mrow></mml:msup></mml:math></inline-formula>

The integral with respect to (dt)α is defined as the solution of the fractional differential equationdy=f(x)(dx)α,x0,y(0)=0, which is provided by the following results.

Lemma 1.3.

Let f(x) denote a continuous function; then the solution y(x) with y(0)=0, of (1.12), is defined by the equality y=oxf(ξ)(dξ)α=αox(x-ξ)α-1f(ξ)dξ,0<α<1.

2. Sumudu Transform of Fractional OrderDefinition 2.1.

Let f(t) denote a function which vanishes for negative values of t. Its Sumudu's transform of order α (or its fractional Sumudu's transform) is defined by the following expression, when it is finite: Sα{f(t)}:=:Gα(u):=0Eα(-tα)f(ut)(dt)α  :=limM0MEα(-tα)f(ut)(dt)α, where u𝒞, and Eα(x) is the Mittag-Leffler function k=0(xk/αk!).

Recently Tchuenche and Mbare introduced the double Sumudu transform . Analogously, we define the fractional double Sumudu transform in following way.

Definition 2.2.

Let f(x,t) denote a function which vanishes for negative values of x and t. Its double Sumudu transform of fractional order (or its fractional double Sumudu transform) is defined as Sα2{f(t,x)}:=:Gα2(u,v)=0Eα[-(t+x)α]f(ut,vx)(dt)α(dx)α, where u,v𝒞, and Eα(x) is the Mittag-Leffler function.

2.1. The Laplace-Sumudu Duality of Fractional Order

The following definition was given in .

Definition 2.3.

Let f(t) denote a function which vanishes for negative values of t. Its Laplace's transform of order α (or its αth fractional Laplace's transform) is defined by the following expression: Lα{f(t)}:=:Fα(u):=0Eα(-(ut)α)f(t)(dt)α=limM0MEα(-(ut)α)f(t)(dt)α provided that integral exists.

Theorem 2.4.

If the Laplace transform of fractional order of a function f(t) is Lα{f(t)}=Fα(u) and the Sumudu transform of this function is Sα{f(t)}=Gα(u), then Gα(u)=1uαFα(1u),0<α<1.

Proof.

By definition of fractional Sumudu transformation, Gα(u)=Sα{f(t)}:=limM0MEα(-tα)f(ut)(dt)α=limMα0M(M-t)α-1Eα(-tα)f(ut)dt. By using the change of variable utt=1uαlimMα0Mu(Mu-t)α-1Eα(-(tu)α)f(t)dt'=1uα0Eα(-(tu)α)f(t)(dt)α=1uαFα(1u). Similarly, on using the definition of fractional Sumudu transform, the following operational formulae can easily be obtained: (i)Sα{f(at)}=Gα(au),(ii)Sα{f(t-b)}=Eα(-bα)Gα(u),(iii)Sα{Eα(-cαtα)f(t)}=1(1+cu)αGα(u1+cu),(iv)Sα{0tf(t)(dt)α}=uαΓ(1+α)Gα(u),(v)Sα{fα(t)}=Gα(u)-Γ(1+α)f(0)uα.

Proof.

It can easily be proved by using Definition 2.1.

Proof.

We start with the following equality by using (2.1): 0MEα(-tα)f(u(t-b))(dt)α=0M(M-t)α-1Eα(-tα)f(u(t-b))dt=0M(M-t-b)α-1Eα(-(b+t)α)f(ut)dt on using the change of variable t-bt. Then it follows that Eα(x+y)α=Eα(xα)Eα(yα).

Proof.

We start from equality (1.13): 0MEα(-tα)Eα(-cαuαtα)f(ut)(dt)α=0M(M-t)α-1Eα(-(1+cu)αtα)f(ut)dt using the change of variable (1+cu)t=1(1+cu)α0M(1+cu)(M(1+cu)-t)α-1Eα(-tα)f(u1+cut)dt.

Proof.

Using fractional Laplace-Sumudu duality and using the result of Jumarie (see ), we can easily obtain these results.

Now we will obtain very similar properties for the fractional double Sumudu transform. Since proofs of these properties are straight, due to this reason, we will give only statements of these properties:(vi)Sα2{f(at)g(bx)}=Gα(au)Hα(bv),(vii)Sα2{f(at,bx)}=Gα2(au,bv),(viii)Sα2{f(t-a,x-b)}=Eα(-(a+b)α)Gα2(au,bv),(ix)Sα2{tαf(t,x)}=Gα2(u,v)-Γ(1+α)f(0,x)uα, where tα is the fractional partial derivative of order α(0<α<1) (see ).

3. The Convolution Theorem and Complex Inversion FormulaProposition 3.1.

If one defines the convolution of order of the two functions f(t) and g(t) by the expression (f(x)*g(x))α:=0xf(x-v)g(v)(dv)α, then Sα{(f(t)*g(t))α}=uαGα(u)Hα(u), where Gα(u)=Sα{f(t)} and Hα(u)=Sα{g(t)}.

Proof.

First recall that the Laplace transform of fractional order of (f*g)α is given by Lα{(f(t)*g(t))α}=Lα{f(t)}Lα{g(t)}. Now, by the fractional Laplace-Sumudu duality relation, Sα{(f(t)*g(t))α}=1uαLα{(f(t)*g(t))α}=1uαLα{f(t)}Lα{g(t)}=uαLα{f(t)}uαLα{g(t)}uα=uαGα(u)Hα(u).

Proposition 3.2.

Given Sumudus transforms that one recalls here for convenience: Gα(u)=0Eα(-xα)f(ux)dx,0<α<1, one has the inversion formula f(x)=1(Mα)α-iiEα((xu)α)uαG((1u)α)(du)α, where Mα is the period of the Mittag-Leffler function.

Proof.

By using complex inversion formula of fractional Laplace transform, see , if Fα(u)=0Eα(-(ux)α)f(x)dx,0<α<1, then inversion formula is given as f(x)=1(Mα)α-iiEα((xu)α)Fα(u)(du)α According to fractional Sumudu-Laplace duality, we can easily yield the desired result.

4. An Application of Fractional Sumudu TransformExample 4.1.

Solution of the equation y(α)+y=f(x),y(0)=0,  0<α<1, is given by f(x)=1(Mα)α-iiEα((xu)α)uαGα((1u)α)(du)α.

Proof.

Taking Sumudu transform of (4.1) both side, we can easily get yα(u)=uα1+uαGα(u) on using y(0)=0. Then by applying the complex inversion formula of fractional Sumudu transforms we get the following result: y(x)=1(Mα)α-iiEα((xu)α)uα(1+uα)Gα((1u)α)(du)α. Now we apply the fractional double Sumudu transform to solve fractional partial differential equation.

Example 4.2.

Consider the linear fractional partial differential equation (see ) tαz(x,t)=cxβz(x,t),x,t+, with the boundary condition z(0,t)=f(t),z(x,0)=g(x), where c is a positive coefficient, and 0<α, β<1.

Proof.

Taking fractional double Sumudu transform of (4.5) both side, we can easily get (1uα-1vβ)Gα2(u,v)=Γ(1+α)uαf(t)-Γ(1+β)vαg(x), which gives Gα2(u,v)=Γ(1+α)(vβvβ-uα)f(t)-Γ(1+β)(uαvβ-uα)g(x).

Acknowledgments

The third author gratefully acknowledges that this research was partially supported by the University Putra Malaysia under the Research University Grant Scheme 05-01-09-0720RU and Fundamental Research Grant Scheme 01-11-09-723FR. The authors also thank the referee(s) for the very constructive comments and suggestions.

WatugalaG. K.Sumudu transform: a new integral transform to solve differential equations and control engineering problemsInternational Journal of Mathematical Education in Science and Technology19932413543MR120684710.1080/0020739930240105ZBL0768.44003AsiruM. A.Sumudu transform and the solution of integral equations of convolution typeInternational Journal of Mathematical Education in Science and Technology2001326906910MR187254310.1080/002073901317147870ZBL1008.45003AşiruM. A.Further properties of the Sumudu transform and its applicationsInternational Journal of Mathematical Education in Science and Technology2002333441449MR191244610.1080/002073902760047940ZBL1013.44001AsiruM. A.Classroom note: application of the Sumudu transform to discrete dynamic systemsInternational Journal of Mathematical Education in Science and Technology2003346944949BelgacemF. B. M.KaraballiA. A.KallaS. L.Analytical investigations of the Sumudu transform and applications to integral production equationsMathematical Problems in Engineering20033-410311810.1155/S1024123X03207018MR2032184EltayebH.KiliçmanA.aKılıçman@putra.upm.edu.myFisherB.A new integral transform and associated distributionsIntegral Transforms and Special Functions20102153673792-s2.0-7034949091910.1080/10652460903335061GuptaV. G.SharmaB.Application of Sumudu transform in reaction-diffusion systems and nonlinear wavesApplied Mathematical Sciences201049–12435446MR2580638KılıçmanA.EltayebH.AgarwalP. R.On Sumudu transform and system of differential equationsAbstract and Applied Analysis201020101159870210.1155/2010/598702KılıçmanA.EltayebH.A note on integral transforms and partial differential equationsApplied Mathematical Sciences201041–4109118MR2576040KademA.Solving the one-dimensional neutron transport equation using Chebyshev polynomials and the Sumudu transformAnalele Universitatii din Oradea. Fascicola Matematica200512153171MR2210623ZBL1164.82331KılıçmanA.GadainH. E.On the applications of Laplace and Sumudu transformsJournal of the Franklin Institute. In press10.1016/j.jfranklin.2010.03.008JumarieG.Fractional partial differential equations and modified Riemann-Liouville derivative new methods for solutionJournal of Applied Mathematics and Computing2007241-23148MR231194710.1007/BF02832299ZBL1145.26302JumarieG.Table of some basic fractional calculus formulae derived from a modified Riemann-Liouville derivative for non-differentiable functionsApplied Mathematics Letters200922337838510.1016/j.aml.2008.06.003MR2483503ZBL1171.26305JumarieG.Laplace's transform of fractional order via the Mittag-Leffler function and modified Riemann-Liouville derivativeApplied Mathematics Letters200922111659166410.1016/j.aml.2009.05.011MR2569059ZBL1181.44001AlmeidaR.ricardo.almeida@ua.ptMalinowskaA. B.abmalinowska@ua.ptTorresD. F. M.delfim@ua.ptA fractional calculus of variations for multiple integrals with application to vibrating stringJournal of Mathematical Physics20105131120335032-s2.0-5304908843610.1063/1.3319559TchuencheJ. M.MbareN. S.An application of the double Sumudu transformApplied Mathematical Sciences200711–43139MR2300973ZBL1154.44001