arxiv: v1 [math.ca] 6 Nov 2015

Transcription

1 rxiv:5.0207v mh.a 6 Nov 205 uo derivives of frcionl vrile order: numericl roximions Din Tvres, dvres@ileiri. Ricrdo Almeid ricrdo.lmeid@u. Delfim F. M. Torres delfim@u. ESES, Polyechnic Insiue of Leiri, Leiri, Porugl ener for Reserch nd Develomen in Mhemics nd Alicions IDMA), Dermen of Mhemics, Universiy of Aveiro, Aveiro, Porugl Asrc We resen new numericl ool o solve ril differenil equions involving uo derivives of frcionl vrile order. Three uo-ye frcionl oerors re considered, nd for ech one of hem n roximion formul is oined in erms of sndrd inegerorder) derivives only. Esimions for he error of he roximions re lso rovided. We hen comre he numericl roximion of some es funcion wih is exc frcionl derivive. We end wih n exemlificion of how he resened mehods cn e used o solve ril frcionl differenil equions of vrile order. Keywords: frcionl clculus, frcionl vrile order, frcionl differenil equions, roximion mehods. MS 200: 33F05, 35R. PAS 200: G, Lj. Inroducion As is well nown, severl hysicl henomen re ofen eer descried y frcionl derivives, 8, 36. This is minly due o he fc h frcionl oerors e ino considerion he evoluion of he sysem, y ing he glol correlion, nd no only locl chrcerisics. Moreover, ineger-order clculus someimes conrdic he exerimenl resuls nd herefore derivives of frcionl order my e more suile 2. An ineresing recen generlizion of he heory of frcionl clculus consiss o llow he frcionl order of he derivive o e non-consn, deending on ime 5, 9, 20. Wih his roch, he non-locl roeries re more eviden nd numerous licions hve een found in hysics, conrol nd signl rocessing 7, 3, 22, 2, 26, 27, 34. One difficul issue, h usully rises when deling wih such frcionl oerors, is he exreme difficuly in solving nlyiclly such rolems 2, 37. Thus, in mos cses, we do nonow he exc soluion for he rolem nd one needs o see numericl roximion. Severl numericl mehods cn e found in he lierure, yiclly lying some discreizion over ime or relcing he frcionl oerors y roer decomosiion 2, 37. Recenly, new roximion formuls were given for frcionl consn order oerors, wih he dvnge h higher-order derivives re no required o oin good ccurcy of he Pr of firs uhor s Ph.D., which is crried ou he Universiy of Aveiro under he Docorl Progrmme Mhemics nd Alicions of Universiies of Aveiro nd Minho. Thisis rerinofer whose finl nd definie formisin ommunicions in Nonliner Science nd Numericl Simulion, ISSN: Per sumied 27/My/205; revised 06/Oc/205; cceed 30/Oc/205.

2 mehod, 23, 24. These decomosiions only deend on ineger-order derivives, nd y relcing he frcionl oerors h er in he rolem y hem, one leves he frcionl conex ending u in he resence of sndrd rolem, where numerous ools re ville o solve hem. Here we exend such decomosiions o uo frcionl rolems of vrile order. The er is orgnized s follows. To sr, in Secion 2 we formule he needed definiions. Nmely, we resen hree yes of uo derivives of vrile frcionl order. Firs, we consider one indeenden vrile only Secion 2.); hen we generlize for severl indeenden vriles Secion 2.2). Secion 3 is he min core of he er: we rove roximion formuls for he given frcionl oerors nd uer ound formuls for he errors. To es he efficiency of he roosed mehod, in Secion 4 we comre he exc frcionl derivive of some es funcion wih he numericl roximions oined from he decomosiion formuls given in Secion 3. To end, in Secion 5 we ly our mehod o roxime wo hysicl rolems involving uo frcionl oerors of vrile order ime-frcionl diffusion equion in Secion 5. nd frcionl Burgers ril differenil equion in fluid mechnics in Secion 5.2) y clssicl rolems h my e solved y well-nown sndrd echniques. 2 Frcionl clculus of vrile order In he lierure of frcionl clculus, severl differen definiions of derivives re found 28. One of hose, inroduced y uo in nd sudied indeendenly y oher uhors, lie Džršjn nd Nersesjn in nd Ronov in , hs found mny licions nd seems o e more suile o model hysic henomen 6, 8, 9, 5, 6, 3, 33, 35. Before generlizing he uo derivive for vrile order of differeniion, we recll wo yes of secil funcions: he Gmm nd Psi funcions. The Gmm funcion is n exension of he fcoril funcion o rel numers, nd is defined y Γ) = 0 τ ex τ)dτ, > 0. We menion h oher definiions exis for he Gmm funcion, nd i is ossile o define i for comlex numers, exce he non-osiive inegers. A sic u fundmenl roery h we will use ler is he following: Γ+) = Γ). The Psi funcion is he derivive of he logrihm of he Gmm funcion: Ψ) = d d lnγ)) = Γ ) Γ). Given α 0,), he lef nd righ uo frcionl derivives of order α of funcion x :, R re defined y D α D α x) x)) nd D α D α x) x)), resecively, where D αx) nd D α x) denoe he lef nd righ Riemnn Liouville frcionl derivive of order α, h is, nd D α d Γ α) d D α d Γ α) d τ) α xτ)dτ τ ) α xτ)dτ. 2

3 If x is differenile, hen, inegring y rs, one cn rove he following equivlen definiions: nd Dα Γ α) Dα Γ α) τ) α x τ)dτ τ ) α x τ)dτ. From hese definiions, i is cler h he uo frcionl derivive of consn is zero, which is flse when we consider he Riemnn Liouville frcionl derivive. Also, he oundry condiions h er in he Llce rnsform of he uo derivive deend on ineger-order derivives, nd so coincide wih he clssicl cse. 2. Vrile order uo derivives for funcions of one vrile Our gol is o consider frcionl derivives of vrile order, wih α deending on ime. In fc, some henomen in hysics re eer descried when he order of he frcionl oeror is no consn, for exmle, in he diffusion rocess in n inhomogeneous or heerogeneous medium, or rocesses where he chnges in he environmen modify he dynmic of he ricle 4, 30, 32. Moived y he ove considerions, we inroduce hree yes of uo frcionl derivives. The order of he derivive is considered s funcion α) ing vlues on he oen inervl 0, ). To sr, we define wo differeninds of Riemnn Liouville frcionl derivives. Definiion Riemnn Liouville frcionl derivives of order α) yes I nd II). Given funcion x :, R,. he ye I lef Riemnn Liouville frcionl derivive of order α) is defined y D α) d Γ α)) d τ) α) xτ)dτ; 2. he ye I righ Riemnn Liouville frcionl derivive of order α) is defined y D α) d Γ α)) d τ ) α) xτ)dτ; 3. he ye II lef Riemnn Liouville frcionl derivive of order α) is defined y D α) d d Γ α)) ) τ) α) xτ)dτ ; 4. he ye II righ Riemnn Liouville frcionl derivive of order α) is defined y D α) d ) τ ) α) xτ)dτ. d Γ α)) The uo derivives re given using he revious Riemnn Liouville frcionl derivives. Definiion 2 uo frcionl derivives of order α) yes I, II nd III). Given funcion x :, R,. he ye I lef uo derivive of order α) is defined y D α) D α) x) x)) = d Γ α)) d τ) α) xτ) x)dτ; 3

4 2. he ye I righ uo derivive of order α) is defined y Dα) D α) x) x)) = d Γ α)) d τ ) α) xτ) x)dτ; 3. he ye II lef uo derivive of order α) is defined y D α) D α) x) x)) = d d Γ α)) τ) α) xτ) x)dτ 4. he ye II righ uo derivive of order α) is defined y D α) D α) x) x)) = d ) τ ) α) xτ) x)dτ ; d Γ α)) 5. he ye III lef uo derivive of order α) is defined y D α) Γ α)) τ) α) x τ)dτ; 6. he ye III righ uo derivive of order α) is defined y Dα) Γ α)) τ ) α) x τ)dτ. In conrs wih he cse when α is consn, definiions of differen yes do no coincide. Theorem 3. The following relions eween he lef frcionl oerors hold: D α) D α) α ) x) + τ) α) x τ) Γ2 α)) α) ln τ) dτ ) nd Dα) Dα) Proof. Inegring y rs, one ges Dα) x) α )Ψ α)) Γ α)) d Γ α)) d d = Γ α)) d Differeniing he inegrl, i follows h D α) α ) Γ α)) α)) 2 + α) ) ; τ) α) xτ) x)dτ. 2) τ) α) xτ) x)dτ α) τ) α) x τ) = α ) Dα) x) + Γ2 α)) The second formul follows from direc clculions. τ) α) x τ)dτ τ) α) x τ)dτ. α )ln τ)+ α) dτ τ τ) α) x τ) α) ln τ) dτ. 4

5 Therefore, when he order α) c is consn, or for consn funcions x), we hve Dα) Similrly, we oin he nex resul. Dα) Dα) x). Theorem 4. The following relions eween he righ frcionl oerors hold: nd Dα) Dα) x) + Dα) Dα) α ) Γ2 α)) x)+ α )Ψ α)) Γ α)) Theorem 5. Le x,,r). A = = D α) Dα) D α) Dα) τ ) α) x τ) lnτ ) dτ α) τ ) α) xτ) x)dτ. D α) 0; Dα) 0. Proof. We sr roving he hird equliy he iniil ime =. We simly noe h D α) x x) Γ α)) τ) α) dτ = x Γ2 α)) ) α), which is zero =. For he firs equliy =, using equion ), nd he wo nex relions τ) α) x τ) α) dτ x α))2 α)) )2 α) nd τ) α) x τ)ln τ)dτ x 2 α) )2 α) ln ) 2 α), his ler inequliy oined from inegrion y rs, we rove h Dα) 0 =. Finlly, we rove he second equliy = y considering equion 2): erforming n inegrion y rs, we ge τ) α) xτ) x)dτ x α))2 α)) )2 α) nd so D α) 0 =. The roof h he righ frcionl oerors lso vnish he end oin = follows y similr rgumens. Wih some comuions, relionshi eween he Riemnn Liouville nd he uo frcionl derivives is esily deduced: D α) D α) x) + = Dα) x) + x) d Γ α)) d x) Γ α)) ) α) + x)α ) Γ2 α)) ) α) τ) α) dτ α) ln ) 5

6 nd D α) D α) x)+x) d d For he righ frcionl oerors, we hve Γ α)) = x) Dα) x) + Γ α)) ) α) D α) D α) x) + ) τ) α) dτ + x)α ) Γ2 α)) ) α) Ψ2 α)) ln ). x) Γ α)) ) α) x)α ) Γ2 α)) ) α) α) ln ) nd D α) D α) x) x) + Γ α)) ) α) x)α ) Γ2 α)) ) α) Ψ2 α)) ln ). Thus, i is immedie o conclude h if x) = 0, hen nd if x) = 0, hen D α) D α) D α) x) nd D α) Dα) x) nd D α) D α) x) Dα) x). Nex we oin formuls for he uo frcionl derivives of ower funcion. Lemm 6. Le ) γ wih γ > 0. Then, D α) Γγ +) Γγ α)+) )γ α) α Γγ +) ) Γγ α)+2) )γ α)+ ln ) Ψγ α)+2)+ψ α)), Γγ +) Dα) Γγ α)+) )γ α) α Γγ +) ) Γγ α)+2) )γ α)+ ln ) Ψγ α)+2), Γγ +) Dα) Γγ α)+) )γ α). Proof. The formul for D α) x) follows immediely from 29. For he second equliy, one hs D α) d ) τ) α) τ ) γ dτ d Γ α)) = d ) α) τ ) α) τ ) dτ) γ. d Γ α)) 6

7 Wih he chnge of vriles τ = s ), nd wih he hel of he Be funcion B, ), we rove h Dα) d ) α) ) s) α) s γ ) γ+ ds d Γ α)) 0 = d ) ) γ α)+ Bγ +, α)) d Γ α)) = d ) Γγ +) d Γγ α)+2) )γ α)+. We oin he desired formul y differeniing his ler exression. The ls equliy follows in similr wy. Anlogous relions o hose of Lemm 6, for he righ uo frcionl derivives of vrile order, re esily oined. Lemm 7. Le ) γ wih γ > 0. Then, Dα) Dα) D α) Γγ +) Γγ α)+) )γ α) +α Γγ +) ) Γγ α)+2) )γ α)+ ln ) Ψγ α)+2)+ψ α)), Γγ +) Γγ α)+) )γ α) +α Γγ +) ) Γγ α)+2) )γ α)+ ln ) Ψγ α)+2), Γγ +) Γγ α)+) )γ α). Wih Lemm 6 in mind, we immediely see h Dα) x) Dα) x) Dα) x). Also, lesforheowerfuncion, isuggessh Dα) x) my e more suile inverse oerion of he frcionl inegrl when he order is vrile. For exmle, consider funcions 2 nd y) = ) 2, nd he frcionl order α) = 5+ 0, 0,. Then, 0. α) 0.6 for ll. Nex we comre he frcionl derivives of x nd y of order α) wih he frcionl derivives of consn order α = 0. nd α = 0.6. By Lemm 6, we now h he lef uo frcionl derivives of order α) of x re given y 0 Dα) 0 D α) 0 D α) 2 Γ3 α)) 2 α) 3 α) Γ4 α)) ln) Ψ4 α))+ψ α)), 2 Γ3 α)) 2 α) 3 α) 2 Γ3 α)) 2 α), Γ4 α)) ln) Ψ4 α)), while y Lemm 7, he righ uo frcionl derivives of order α) of y re given y D α) y) = D α) y) = Dα) y) = 2 Γ3 α)) )2 α) + )3 α) Γ4 α)) ln ) Ψ4 α))+ψ α)), 2 Γ3 α)) )2 α) + )3 α) 2 Γ3 α)) )2 α). Γ4 α)) ln ) Ψ4 α)), 7

8 For consn order α, we hve 0 Dα 2 Γ3 α) 2 α nd Dα y) = 2 Γ3 α) )2 α. The resuls cn e seen in Figure. 2.2 Vrile order uo derivives for funcions of severl vriles Pril frcionl derivives re nurl exension nd re defined in similr wy. Le m N, m {,...,m}, nd consider funcion x : i, i R wih m vriles. For simliciy, we define he vecors nd i= τ =,...,,τ,+,..., m ) R m ) =,..., m ) R m. Definiion 8 Pril uo frcionl derivives of vrile order yes I, II nd III). Given funcion x : m i= i, i R nd frcionl orders α :, 0,), {,...,m},. he ye I ril lef uo derivive of order α ) is defined y D α ) Γ α )) τ) α ) xτ x )dτ; 2. he ye I ril righ uo derivive of order α ) is defined y D α ) Γ α )) τ ) α ) xτ x )dτ; 3. he ye II ril lef uo derivive of order α ) is defined y D α ) Γ α )) τ) α ) xτ x )dτ 4. he ye II ril righ uo derivive of order α ) is defined y D α ) ) τ ) α ) xτ x )dτ ; Γ α )) 5. he ye III ril lef uo derivive of order α ) is defined y D α ) τ) α ) x τ dτ; Γ α )) 6. he ye III ril righ uo derivive of order α ) is defined y D α ) τ ) α ) x τ dτ. Γ α )) Similrly s done efore, relions eween hese definiions cn e roven. ) ; 8

9 ) 0 Dα) x) ) 0 Dα) x) c) 0 Dα) x) d) Dα) y) e) Dα) y) f) Dα) y) Figure : omrison eween vrile order nd consn order frcionl derivives. Theorem 9. The following four formuls hold: D α ) D α ) x) + α ) Γ2 α )) τ) α ) x τ α ) ln τ) dτ, 3) 9

10 D α ) D α ) x) α )Ψ α )) τ) α ) xτ x dτ, 4) Γ α )) nd D α ) D α ) x) + α ) Γ2 α )) τ ) α ) x τ D α ) D α ) x)+ α )Ψ α )) Γ α )) α ) lnτ ) dτ τ ) α ) xτ x dτ. 3 Aroximion of vrile order uo derivives Le N. We define A = + Γ+ α )) l=n + Γα ) n+) B = Γ α ))Γα )) n)!, n x V ) = τ ) τ dτ, L ) = mx x τ. τ, Γα ) n+l), Γα ) )l n+)! Theorem 0. Le x n+ m i= i, i,r) wih n N. Then, for ll {,...,m} nd for ll N N such h N n, we hve D α ) n = A ) α ) x The roximion error E) is ounded y + B ) n α ) V )+E). exn α )) 2 +n α )) E) L n+ ) Γn+ α ))N n α ) n α )) ) n+ α ). Proof. By definiion, D α ) =n τ) α ) x τ dτ Γ α )) nd, inegring y rs wih u τ) = τ) α ) nd vτ) = x τ, we deduce h D α ) ) α ) x + Γ2 α )) Γ2 α )) τ) α ) 2 x 2 Inegring gin y rs, ing u τ) = τ) α ) nd vτ) = 2 xτ 2, we ge D α ) ) α ) x + ) 2 α) Γ2 α )) Γ3 α )) + Γ3 α )) 2 x 2 τ) 2 α ) 3 x 3 τ dτ. τ dτ. 0

11 Reeing he sme rocedure n 2 more imes, we ge he exnsion formul D α ) Using he equliies wih we rrive D α ) n = ) α ) x Γ+ α )) + Γn+ α )) τ) n α ) = ) n α ) τ = n = = ) n α ) E) = =N+ ) α ) x Γ+ α )) + ) n α ) Γn+ α )) n = N =0 n α ) τ) n α ) n+ x n+ ) n α ) n α ) ) ) τ ) ), =0 ) α ) x Γ+ α )) + ) n α ) Γn+ α )) =0 ) ) τ ) ) +E) τ dτ. n α ) ) ) τ ) n+ x ) n+ τ dτ +E) ) n α ) ) ) τ ) n+ x n+ τ dτ +E) wih E) = ) n α ) E) n+ x Γn+ α )) n+ τ dτ. Now, we sli he ls sum ino = 0 nd he remining erms =,...,N nd inegre y rs wih uτ) = τ ) nd v τ) = n+ xτ. Oserving h we oin: ) n α ) Γn+ α )) =0 ) n α ) n+ ) n α ) ) = Γα ) n+), Γα ) n)! ) ) = ) n α ) n x Γn+ α )) n n x n ) n x n τ ) n x n τ ) n+ x n+ + ) n α ) Γn+ α )) τ dτ τ dτ = Γα ) n+) Γα ) n)! )

12 = ) n α ) Γn+ α )) + ) n α ) Γn α )) Thus, we ge n D α ) n x n + ) n α ) Γn+ α )) = = n x n + = Γα ) n+) Γα ) n)! Γα ) n+) Γα )+ n) )! ) τ ) n x n τ dτ. ) α ) Γ+ α )) x + ) n α ) n x Γn+ α )) n + + ) n α ) Γn α )) τ ) n x n = = Γα ) n+) Γα ) n)! Γα ) n+) Γα )+ n) )! ) τ dτ +E). Reeing he rocess n more imes wih resec o he ls sum, h is, sliing he firs erm of he sum nd inegring y rs he oined resul, we rrive o n D α ) ) α ) x Γα ) n+l) Γ+ α )) + Γα ) )l n+)! = + =n l=n + Γα ) n+) Γ α ))Γα )) n)! ) n α ) τ ) n x τ dτ +E). We now see he uer ound formul for E). Using he wo relions τ ), if τ, nd n α ) exn α )) 2 +n α )), n+ α ) we ge Then, E) =N+ N exn α )) 2 +n α )) n+ α ) exn α )) 2 +n α )) n+ α ) d = exn α )) 2 +n α )). N n α ) n α )) exn α )) 2 +n α )) E) L n+ ) Γn+ α ))N n α ) n α )) ) n+ α ). This concludes he roof. Remr. In Theorem 0 we hve for ll m i= i, i nd n N. lim E) = 0 N 2

13 Theorem 2. Le x n+ m i= i, i,r) wih n N. Then, for ll {,...,m} nd for ll N N such h N n, we hve D α ) n A ) α ) x + B ) n α ) V ) = + α ) ) α ) Γ2 α )) =n α ) ln ) + =0 The roximion error E) is ounded y ) N =0 ) α ) ) ) V n+) α ) ) ) N r r= ) +rv n++r) exn α )) 2 +n α )) E) L n+ ) Γn+ α ))N n α ) n α )) ) n+ α ) + α ex α )) 2 + α )) ) L ) Γ2 α ))N α ) α )) α ) ln ) + ) 2 α ). N Proof. Ting ino ccoun relion 3) nd Theorem 0, we only need o exnd he erm +E). α ) τ) α ) x τ Γ2 α )) α ) ln τ) dτ. 5) Sliing he inegrl, nd using he exnsion formuls τ) α ) = ) α ) τ = ) α ) N =0 α ) ) α ) ) ) τ ) ) +E ) wih nd wih E ) = =N+ α ) ) ) τ ) ) ln τ) = ln )+ln τ ) τ ) r = ln ) r ) r E 2) E 2 ) = r=n+ r= τ ) r r ) r, 3

14 we conclude h 5) is equivlen o ) α ) Γ2 α )) α ) ln ) τ) α ) x τ dτ τ) α ) ln τ ) x τ dτ α = ) ) Γ2 α )) α ) ln ) ) α ) α ) ) ) τ ) x =0 ) τ dτ + ) α ) α ) ) ) τ ) τ ) r x =0 ) r r= ) r τ dτ α + ) ) Γ2 α )) α ) ln ) ) α ) E ) x τ dτ + ) α ) E )E 2 ) x τ dτ = α ) ) α ) ) Γ2 α )) α ) ln N ) ) α ) ) =0 ) V n+) + α ) ) ) N r =0 r= ) +rv n++r) + α ) ) α ) Γ2 α )) ) α ) ln ). E ) x τ dτ + For he error nlysis, we now from Theorem 0 h E ) ex α )) 2 + α )). N α ) α )) E )E 2 ) x τ dτ Then, ) α ) E ) x τ dτ L ) ex α )) 2 + α )) N α ) 2 α ). ) α )) 6) On he oher hnd, we hve ) α ) E )E 2 ) x τ dτ L ) ex α )) 2 + α )) N α ) α ) ) α )) r ) r τ ) r dτ = L ) ex α )) 2 + α )) N α ) α )) ) α ) L ) ex α )) 2 + α )) N 2 α ) 2 α ). ) α )) r=n+ r=n+ We ge he desired resul y comining inequliies 6) nd 7). rr +) 7) 4

15 Theorem 3. Le x n+ m i= i, i,r) wih n N. Then, for ll {,...,m} nd for ll N N such h N n, we hve D α ) n A ) α ) x + B ) n α ) V ) = + α ) ) α ) Γ2 α )) + =0 α ) ) ) N r= The roximion error E) is ounded y =n Ψ2 α )) ln )) =0 r ) +rv n++r) +E). ) α ) ) ) V n+) exn α )) 2 +n α )) E) L n+ ) Γn+ α ))N n α ) n α )) ) n+ α ) + α ex α )) 2 + α )) ) L ) Γ2 α ))N α ) α )) Ψ2 α )) ln ) + ) 2 α ). N Proof. Sring wih relion 4), nd inegring y rs he inegrl, we oin h D α ) D α ) x)+ α )Ψ α )) Γ2 α )) The res of he roof is similr o he one of Theorem 2. τ) α ) x τ dτ. Remr 4. As riculr cses of Theorems 0, 2 nd 3, we oin exnsion formuls for x), D α) x) nd D α) x). D α) Wih resec o he hree righ frcionl oerors of Definiion 8, we se, for N, ) Γα ) n+l) = +, Γ+ α )) Γα ) )l n+)! l=n + Γα ) n+) D = Γ α ))Γα )) n)!, n x W ) = τ) τ dτ, M ) = mx x τ. τ, The exnsion formuls re given in Theorems 5, 6 nd 7. We omi he roofs since hey re similr o he corresonding lef ones. Theorem 5. Le x n+ m i= i, i,r) wih n N. Then, for ll {,...,m} nd for ll N N such h N n, we hve D α ) n = ) α ) x The roximion error E) is ounded y + D ) n α ) W )+E). exn α )) 2 +n α )) E) M n+ ) Γn+ α ))N n α ) n α )) ) n+ α ). =n 5

16 Theorem 6. Le x n+ m i= i, i,r) wih n N. Then, for ll {,...,m} nd for ll N N such h N n, we hve D α ) n ) α ) x + D ) n α ) W ) = + α ) ) α ) Γ2 α )) + =0 =n α ) ln ) α ) ) ) N r r= ) +rw n++r) The roximion error E) is ounded y ) N =0 +E). ) α ) ) ) W n+) exn α )) 2 +n α )) E) M n+ ) Γn+ α ))N n α ) n α )) ) n+ α ) + α ex α )) 2 + α )) ) M ) Γ2 α ))N α ) α )) α ) ln ) + ) 2 α ). N Theorem 7. Le x n+ m i= i, i,r) wih n N. Then, for ll {,...,m} nd for ll N N such h N n, we hve D α ) n ) α ) x + D ) n α ) W ) = + α ) ) α ) Γ2 α )) + =0 α ) ) ) N r= The roximion error E) is ounded y =n Ψ2 α )) ln )) =0 r ) +rw n++r) +E). ) α ) ) ) W n+) exn α )) 2 +n α )) E) M n+ ) Γn+ α ))N n α ) n α )) ) n+ α ) + α ex α )) 2 + α )) ) M ) Γ2 α ))N α ) α )) Ψ2 α )) ln ) + ) 2 α ). N 4 An exmle To es he ccurcy of he roosed mehod, we comre he frcionl derivive of concree given funcion wih some numericl roximions of i. For 0,, le 2 e he es funcion. For he order of he frcionl derivives we consider wo cses: α) = nd β) = We consider he roximions given in Theorems 0, 2 nd 3, wih fixed n = nd N {2,4,6}. The error of roximing f) y f) is mesured y f) f). See Figures

17 ) 0 Dα) x) ) Error Figure 2: Tye III lef uo derivive of order α) for he exmle of Secion 4 nlyic versus numericl roximions oined from Theorem 0. ) 0 Dα) x) ) Error Figure 3: Tye I lef uo derivive of order α) for he exmle of Secion 4 nlyic versus numericl roximions oined from Theorem 2. 5 Alicions In his secion we ly he roosed echnique o some concree frcionl differenil equions of hysicl relevnce. 5. A ime-frcionl diffusion equion We exend he one-dimensionl ime-frcionl diffusion equion 4 o he vrile order cse. onsider u = ux,) wih domin 0, 2. The ril frcionl differenil equion of order α) is he following: 0 D α) ux,) 2 u x2x,) = fx,) for x 0,, 0,, 8) 7

18 ) 0 Dα) x) ) Error Figure 4: Tye II lef uo derivive of order α) for he exmle of Secion 4 nlyic versus numericl roximions oined from Theorem 3. ) 0 Dβ) x) ) Error Figure 5: Tye III lef uo derivive of order β) for he exmle of Secion 4 nlyic versus numericl roximions oined from Theorem 0. sujec o he oundry condiions nd ux,0) = gx), for x 0,), 9) u0,) = u,) = 0, for 0,. 0) We menion h when α), one oins he clssicl diffusion equion, nd when α) 0 one ges he clssicl Helmholz elliic equion. Using Lemm 6, i is esy o chec h ux,) = 2 sin2πx) is soluion o 8) 0) wih ) 2 fx,) = Γ3 α)) 2 α) +4π 2 2 sin2πx) 8

19 ) 0 Dβ) x) ) Error Figure 6: Tye I lef uo derivive of order β) for he exmle of Secion 4 nlyic versus numericl roximions oined from Theorem 2. ) 0 Dβ) x) ) Error Figure 7: Tye II lef uo derivive of order β) for he exmle of Secion 4 nlyic versus numericl roximions oined from Theorem 3. nd gx) = 0 comre wih Exmle in 4). The numericl rocedure is he following: relce 0 Dα) u wih he roximion given in Theorem 0, ing n = nd n rirry N, h is, 0 Dα) ux,) A α) u N x,)+ B α) V x,) = 9

20 wih A = + Γ2 α)) l= Γα) +) B = Γ α))γα)) )!, V x,) = 0 τ u x,τ)dτ. Γα) +l), Γα) )l! Then, he iniil frcionl rolem 8) 0) is roximed y he following sysem of secondorder ril differenil equions: nd A α) u x,)+ N V = B α) V x,) 2 u x2x,) = fx,) x,) = ux,), =,...,N, for x 0, nd for 0,, sujec o he oundry condiions ux,0) = 0, for x 0,), nd u0,) = u,) = 0, for 0,, V x,0) = 0, for x 0,, =,...,N. Remr 8. As ws menioned in Theorem 0, s N increses, he error of he roximion decreses nd he given roximion formul converges o he frcionl derivive. Thus, in order o hve good ccurcy for he mehod, one should e higher vlues for N. Remr 9. We re no wre of similr mehods o our, concerning vrile frcionl clculus, in order o comre he erformnce of he roosed mehod o oher numericl roximion mehods. For his reson, we decided o comre wih he exc soluion. In he ville lierure, using discreizion rocess, FDEs re solved s finie differences. Our echnique is quie differen: we rewrie he FDE s sysem of ordinry differenil equions, nd hen we cn ly ny nown echnique o solve i. Noe h he reson why we soed here wih N = 6 ws o hve n roximion h is enough close o he exc soluion u sill visully disinguishly when we increse N more, he roximion nd he exc soluion er o e he sme in he los). In erms of erformnce of he mehod, i is roughly he sme o u N = 6 or igger. 5.2 A frcionl ril differenil equion in fluid mechnics We now ly our roximion echniques o he following one-dimensionl liner inhomogeneous frcionl Burgers equion of vrile order see 7, Exmle 5.2): 0 Dα) ux,)+ u x x,) 2 u 22 α) x2x,) = +2x 2, for x 0,, 0,, ) Γ3 α)) sujec o he oundry condiion ux,0) = x 2, for x 0,). 2) Here, Fx,) = 22 α) Γ3 α)) +2x 2 20

21 is he exernl force field. Burgers equion is used o model gs dynmics, rffic flow, urulence, fluid mechnics, ec. The exc soluion is ux,) = x The frcionl rolem ) 2) cn e roximed y A α) u N x,)+ B α) V x,)+ u x x,) 2 u 22 α) x2x,) = Γ3 α)) +2x 2 = wih A, B nd V, {,...,N}, s in Secion 5.. The roximion error cn e decresed s much s desired y incresing he vlue of N. Acnowledgmens This wor ws suored y Poruguese funds hrough he ener for Reserch nd Develomen in Mhemics nd Alicions IDMA) nd The Poruguese Foundion for Science nd Technology FT), wihin rojec UID/MAT/0406/203. Tvres ws lso suored y FT hrough he Ph.D. fellowshi SFRH/BD/42557/2007; Torres y rojec PTD/EEI-AUT/450/202, co-finnced y FEDER under POF-QREN wih OMPETE reference FOMP FEDER The uhors re very greful o hree nonymous referees, for severl commens nd imrovemen suggesions. References T. M. Ancovic, M. Jnev, S. Piliovic nd D. Zoric, An exnsion formul for frcionl derivives of vrile order, en. Eur. J. Phys. 203), no. 0, A. Angn nd A. H. loo, Siliy nd convergence of he sce frcionl vrile-order Schrödinger equion, Adv. Difference Equ. 203, 203:80, 0. 3 M. uo, Liner model of dissiion whose Q is lmos frequency indeenden II, Geohys. J. R. Asr. Soc ), A. V. hechin, R. Gorenflo nd I. M. Soolov, Frcionl diffusion in inhomogeneous medi, J. Phys. A ), no. 42, L679 L S. hen, F. Liu nd K. Burrge, Numericl simulion of new wo-dimensionl vrile-order frcionl ercolion equion in non-homogeneous orous medi, omu. Mh. Al ), no. 9, F. M. oimr, Mechnics wih vrile-order differenil oerors, Ann. Phys. 8) ), no. -2, F. M. oimr,. M. Soon nd M. H. Koyshi, The vrile viscoelsiciy oeror, Annlen der Physi ), M. Dlir nd M. Bshour, Alicions of frcionl clculus, Al. Mh. Sci. Ruse) 4 200), no. 2-24, K. Diehelm, The nlysis of frcionl differenil equions, Lecure Noes in Mhemics, 2004, Sringer, Berlin, M. M. Džršjn nd A. B. Nersesjn, Frcionl derivives nd he uchy rolem for differenil equions of frcionl order, Izv. Ad. Nu Armjn. SSR Ser. M ), no., R. Herrmnn, Folded oenils in cluser hysics comrison of Yuw nd oulom oenils wih Riesz frcionl inegrls, J. Phys. A ), no. 40, , 2. 2 R. Hilfer, Alicions of frcionl clculus in hysics, World Sci. Pulishing, River Edge, NJ, D. Ingmn nd J. Suzdlnisy, onrol of dming oscillions y frcionl differenil oeror wih ime-deenden order, omu. Mehods Al. Mech. Engrg ), no. 52, Y. Lin nd. Xu, Finie difference/secrl roximions for he ime-frcionl diffusion equion, J. omu. Phys ), no. 2,

22 5 J. A. T. Mchdo, M. F. Silv, R. S. Bros, I. S. Jesus,. M. Reis, M. G. Mrcos nd A. F. Glhno, Some licions of frcionl clculus in engineering, Mh. Prol. Eng ), Ar. ID 63980, D. A. Murio nd. E. Mejí, Generlized ime frcionl IHP wih uo frcionl derivives, J. Phys. onf. Ser ), Ar. ID 02074, 8. 7 Z. Odi nd S. Momni, The vriionl ierion mehod: n efficien scheme for hndling frcionl ril differenil equions in fluid mechnics, omu. Mh. Al ), no. -2, T. Odzijewicz, A. B. Mlinows nd D. F. M. Torres, Frcionl clculus of vriions in erms of generlized frcionl inegrl wih licions o hysics, Asr. Al. Anl. 202, Ar. ID 8792, 24. rxiv: T. Odzijewicz, A. B. Mlinows nd D. F. M. Torres, Vrile order frcionl vriionl clculus for doule inegrls, Proceedings of he 5s IEEE onference on Decision nd onrol, Decemer 0 3, 202, Mui, Hwii, Ar. no ), rxiv: T. Odzijewicz, A. B. Mlinows nd D. F. M. Torres, Frcionl vriionl clculus of vrile order, in Advnces in hrmonic nlysis nd oeror heory, 29 30, Oer. Theory Adv. Al., 229, Birhäuser/Sringer Bsel AG, Bsel, 203. rxiv: T. Odzijewicz, A. B. Mlinows nd D. F. M. Torres, Noeher s heorem for frcionl vriionl rolems of vrile order, en. Eur. J. Phys. 203), no. 6, rxiv: P. W. Oslczy, P. Duch, D. W. Brzezińsi nd D. Snowsi, Order Funcions Selecion in he Vrile-, Frcionl-Order PID onroller, Advnces in Modelling nd onrol of Non-ineger-Order Sysems, Lecure Noes in Elecricl Engineering ), S. Pooseh, R. Almeid nd D. F. M. Torres, Aroximion of frcionl inegrls y mens of derivives, omu. Mh. Al ), no. 0, rxiv: S. Pooseh, R. Almeid nd D. F. M. Torres, Numericl roximions of frcionl derivives wih licions, Asin J. onrol 5 203), no. 3, rxiv: Yu. N. Ronov, ree rolems in srucurl memers, Norh-Hollnd Series in Alied Mhemics nd Mechnics, Amserdm/London, L. E. S. Rmirez nd. F. M. oimr, On he vrile order dynmics of he nonliner we cused y sedimening ricle, Phys. D ), no. 3, M. R. Rić nd A. Pisno, Vrile-order frcionl oerors for dive order nd rmeer esimion, IEEE Trns. Auom. onrol ), no. 3, S. G. Smo, A. A. Kils nd O. I. Mrichev, Frcionl inegrls nd derivives, rnsled from he 987 Russin originl, Gordon nd Brech, Yverdon, S. G. Smo nd B. Ross, Inegrion nd differeniion o vrile frcionl order, Inegrl Trnsform. Sec. Func. 993), no. 4, F. Snmri, S. Wils, E. de Schuer nd G. J. Augusine, Anomlous diffusion in Purinje cell dendries cused y sines, Neuron ), K. Singh, R. Sxen nd S. Kumr, uo-sed frcionl derivive in frcionl Fourier rnsform domin, IEEE Journl on Emerging nd Seleced Toics in ircuis nd Sysems 3 203), H. G. Sun, W. hen nd Y. Q. hen, Vrile order frcionl differenil oerors in nomlous diffusion modeling, Physic A ) N. H. Sweilm nd H. M. AL-Mrwm, On he Numericl Soluions of he Vrile Order Frcionl He Equion, Sudies in Nonliner Sciences 2 20), D. Vlério nd J. Sá d os, Vrile order frcionl conrollers, Asin J. onrol 5 203), no. 3, T. Yjim nd K. Ymsi, Geomery of surfces wih uo frcionl derivives nd licions o incomressile wo-dimensionl flows, J. Phys. A ), no. 6, 06520, B. Zheng, G /G)-exnsion mehod for solving frcionl ril differenil equions in he heory of mhemicl hysics, ommun. Theor. Phys. Beijing) ), no. 5, P. Zhung, F. Liu, V. Anh nd I. Turner, Numericl mehods for he vrile-order frcionl dvecion-diffusion equion wih nonliner source erm, SIAM J. Numer. Anl ), no. 3,