arxiv: v1 [math.gm] 13 Nov 2015

Transcription

1 rxiv: v1 [mth.gm] 13 Nov 215 Modelling some rel phenomen by frctionl differentil equtions Ricrdo Almeid 1, Nuno R. O. Bstos 1,2 nd M. Teres T. Monteiro 3 1 enter for Reserch nd Development in Mthemtics nd Applictions (IDMA) Deprtment of Mthemtics, University of Aveiro, Aveiro, Portugl 2 entre for the Study of Eduction, Technologies nd Helth (I&DETS) Deprtment of Mthemtics, School of Technology nd Mngement of Viseu Polytechnic Institute of Viseu, Viseu, Portugl 3 Algoritmi R&D enter, Deprtment of Production nd Systems University of Minho, mpus de Gultr, Brg, Portugl Abstrct This pper dels with frctionl differentil equtions, with dependence on puto frctionl derivtive of rel order. The gol is to show, bsed on concrete exmples nd experimentl dt from severl experiments, tht frctionl differentil equtions my model more efficiently certin problems thn ordinry differentil equtions. A numericl optimiztion pproch bsed on lest squres pproximtion is used to determine the order of the frctionl opertor tht better describes rel dt, s well s other relted prmeters. keywords: Frctionl clculus, frctionl differentil eqution, numericl optimiztion 1 Introduction 1.1 Frctionl clculus Westrtwithreviewonfrctionlclculus, spresentedine.g. [4,9,1]. Frctionlclculusisn extension of ordinry clculus, in wy tht derivtives nd integrls re defined for rbitrry rel order. In some phenomen, frctionl opertors llow to model better thn ordinry derivtives nd ordinry integrls, nd cn represent more efficiently systems with high-order dynmics nd complex nonliner phenomen. This is due to two min resons; first, we hve the freedom to choose ny order for the derivtive nd integrl opertors, nd not be restricted to integer-order only. Secondly, frctionl order derivtives depend not only on locl conditions but lso on the pst, useful when the system hs long-term memory. In fmous letter dted 1695, L Hopitl s sked Leibniz wht would be the derivtive of order α = 1/2, nd the response of Leibniz An pprent prdox, from which one dy useful consequences will be drwn becme the birth of frctionl clculus. Along the centuries, severl ttempts were mde to define frctionl opertors. For exmple, in 173, using the formul d n x m dx n = m(m 1)...(m n+1)xm n = Γ(m+1) Γ(m n+1) xm n, Euler obtined the following expression d 1/2 x 4x dx = 1/2 π. 1

2 Here, Γ denotes the Gmm function, Γ(t) = s t 1 e s ds, t R\Z, which is n extension of the fctoril function to rel numbers. Two of the bsic properties of the Gmm function re Γ(m) = (m 1)!, for ll m N, nd Γ(t+1) = tγ(t), for ll t R\Z. However, the most fmous nd importnt definition is due to Riemnn. Strting with uchy s formul s1 sn 1 1 t ds 1 ds 2... y(s n )ds n = (t s) n 1 y(s)ds, (n 1)! Riemnn defined frctionl integrl type of order α > of function y : [,b] R s I α t y(t) = 1 Γ(α) (t s) α 1 y(s)ds. Frctionl derivtives re defined using the frctionl integrl ide. The Riemnn Liouville frctionl derivtive of order α > is given by ( ) n ( ) n d Dt α y(t) = It n α 1 d t y(t) = (t s) n α 1 y(s)ds, dt Γ(n α) dt where n N is such tht α (n 1,n). These definitions re probbly the most importnt with respect to frctionl opertors, nd until the 2th century the subject ws relevnt in pure mthemtics only. Nowdys, this is n importnt field not only in mthemtics but lso in other sciences, engineering, economics, etc. In fct, using these more generl concepts, we cn describe better certin rel world phenomen which re not possible using integer-order derivtives. For exmple, we cn find pplictions of frctionl clculus in mechnics [1], engineering [6], viscoelsticity [7], dynmicl systems [12], etc. Although the Riemnn Liouville frctionl derivtive is of gret importnce, from historil resons, it my not be suitble to model some rel world phenomen. So, different opertors re considered in the literture, for exmple, the puto frctionl derivtive. It hs proven its pplicbility due to two resons: the frctionl derivtive of constnt is zero nd the initil vlue problems depend on integer-order derivtives only. The puto frctionl derivtive of function y : [,b] R of order α > is defined by Dt α y(t) = It n α y (n) 1 (t) = Γ(n α) where n = [α]+1. In prticulr, when α (,1), we obtin Dα t y(t) = 1 Γ(1 α) (t s) n α 1 y (n) (s)ds, (t s) α y (s)ds. The frctionl integrl nd frctionl derivtive re inverse opertions, in the sense tht (see [4]): Lemm 1. Let α > nd n N be such tht α (n 1,n). If y A n [,b], or y n [,b], then n 1 It α Dα t y(t) = y(t) y (k) () (t ) k. k! k= 2

3 We cn obtin ordinry derivtives by computing the limit when α n N. In fct, we hve the following property: lim α n D α t y(t) = y (n) (t) nd lim α n + D α t y(t) = y (n) (t) y (n) (). For exmple, let y (t) = t, with t. Then, for α (,1), we hve D α t y (t) = while for α > 1, we hve D α t y (t) = (see Fig. 1). 1 Γ(2 α) t1 α, D α t y (t).5 α=1 α=.9 α=.95 α= t Figure 1: Derivtive of y nd its frctionl derivtives. However, we hve the following result. Theorem 1. Let f : [,b] R R be continuous function nd α > rel. If yα is solution of the system { Dt α y(t) = f(t,y), α (,1) y() = y, nd if the limit lim α 1 ±y α(t) := y (t) exists for ll t [,b], then y is solution of the uchy problem { y (t) = f(t,y) y() = y. Proof. Strting with eqution Dα t y α (t) = f(t,y α ), pplying the Iα opertor to both sides nd using Lemm 1, we obtin tht yα stisfies the Volterr eqution Tking the limit α 1 ±, we get y α (t) = y + 1 Γ(α) y (t) = y + tht is, y stisfies the uchy problem. (t s) α 1 f(s,y α (s))ds. f(s,y (s))ds, 3

4 Remrk 1. The ide of Theorem 1 is the following. We first consider frctionl differentil eqution of order α (,1), nd s α 1, the frctionl system converts into n ordinry uchy problem. Then, in order to obtin higher ccurcy for the method, we consider the solution of the frctionl system y, depending on the prmeter α >, nd compute the limit α 1 ±. 1.2 omputtionl issues The MATLAB numericl experiments to find the model prmeters s well s the non-integer derivtives tht fits better in the model is crried out. The numericl experiments were done in MATLAB [8] using the routine lsqcurvefit tht solves nonliner dt-fitting problems in lest-squres sense, tht is, given dt (y i ) i=1...m, nd model (M i ) i=1...m, the error is: E = m (y i M i ) 2. i=1 This routine is bsed on n itertive method with locl convergence, i.e., depending on the initil pproximtion to the prmeters to estimte. The initil pproximtions were found, bsed on the dt nlysis of ech model. In the computtionl tests the trust-region-reflective lgorithm ws selected. In ll the studied cses, we will see tht the frctionl pproch is more close to the dt thn the clssicl one. To test the efficiency of the model, we compre the error of the clssicl model E clssicl with the error given by the frctionl model E frctionl : 1.3 Outline E clssicl E frctionl E clssicl This rticle is orgnized s follows. In ech section we present rel problem, described by n ordinry or system of ordinry differentil equtions. We then consider the sme problem, but modeled by frctionl or system of frctionl differentil equtions, nd compre which of the two models re more suitble to describe the process, bsed on rel experimentl dt. Given n ordinry differentil eqution y (t) = f(t,y), we replce it by the frctionl differentil eqution Dα t y(t) = f(t,y), with α (,1). When we consider the limit α 1, we obtin the initil one. Ifweconsiderα (1,2)swell,ndthentkethelimitα 1 +, wewouldgety (t) y () = f(t,y), nd this is the reson why we will neglect this cse firstly. After computing the solution of the frctionl differentil eqution, which depends on α, we consider the function with frctionl order α on the intervl (,2) in order to increse the ccurcy of the method. Three distinct problems re studied in this pper. Section 2 dels with the exponentil lw of growth, nd we replce it by the Mittg Leffler function, which is generliztion of the exponentil function. In Section 3 we del with n ppliction relted to Blood Alcohol Level using integer-order derivtives. In the finl Section 4 we study model of tpe counter redings t n instnt time. In ech cse, we compute the vlues of the prmeters tht better fit with the given dt, nd lso the error in ech of both pproches. 2 World Popultion Growth 2.1 Stndrd pproch There exist severl ttempts to describe the World Popultion Growth [11]. The simplest model is the following, known s the Mlthusin lw of popultion growth, which is used to predict popultions under idel conditions. Let N(t) be the number of individuls in popultion t time t, B nd M the birth nd mortlity rtes, respectively, so tht the net growth rte is given by. N (t) = (B M)N(t) = PN(t), (1) 4

5 where P := B M is the production rte. Here, we ssume tht B nd M re constnt, nd thus P is lso constnt. The solution of this differentil eqution is the function N(t) = N e Pt, t, (2) where N is the popultion t t =. Becuse of the solution (2), this model is lso known s the exponentil growth model. 2.2 Frctionl pproch onsidered now tht the World Popultion Growth model is ruled by the frctionl differentil eqution Dt α N(t) = PN(t), t, α (,1). (3) Observe tht, tking the limit α 1, Eq (3) converts into Eq. 1, but if we consider α (1,2) nd tke the limit α 1 +, we obtin N (t) N () = PN(t). Using Theorem 7.2 in [2], the solution of this frctionl differentil eqution is the function where E α is the Mittg Leffler function E α (t) = N(t) = N E α (Pt α ), (4) k= t k Γ(αk +1), t R. onsider now function (4), with α (,2). Then, s α 1 ±, we recover the solution for the clssicl problem (2). 2.3 Numericl experiments For our numericl tretment of the problem, we find severl dtbses with the world popultion through the centuries. Here, we use the one provided by the United Ntions [13], from yer 191 until 21, consisting in 11 vlues, where the initil vlue is N = 175. For the clssicl pproch, the production rte is P nd the error from the dt with respect to the nlytic solution (2) is given by E clssicl If we tke into considertion the frctionl model (4), for α (,2) we obtin tht the best vlues re α = nd P with error E frctionl The gin of the efficiency in this procedures is The grph of Figure 2 illustrtes the World Popultion Growth model with the dt, the clssicl model nd the frctionl model. 5

6 8 7 Experiment lssicl Model, Error= x1 5 Frctionl Model, Error= x1 5 World Popultion (in Millions) Yer Figure 2: World Popultion Growth model. 3 Blood Alcohol Level problem 3.1 Stndrd pproch In [5] we find simple model to determine the level of Blood Alcohol, described by system of two differentil equtions. Let A represent the concentrtion of lcohol in the stomch nd B the concentrtion of lcohol in the blood. The problem is described by the following uchy system: A (t) = k 1 A(t) B (t) = k 1 A(t) k 2 B(t) A() = A B() =, where A is the initil lcohol ingested by the subject nd k 1,k 2 some rel constnts. The solution of this system is given by the two functions A(t) = A e k1t nd B(t) = A k 1 k 2 k 1 (e k1t e k2t ). Also, in [5] some experimentl dt is obtined in order to determine the rbitrry constnts k 1 nd k 2 (Tble 1). The time is in minutes nd the Blood Alcohol Level (BAL) in mg/l. Time BAL Tble 1: Experimentl dt for Blood Alcohol Level. Using the dt from the tble, the vlues tht minimize the Men Absolute Error when the vlues from B re fitted with the experimentl dt re: A = , k 1 =.19456, k 2 =.17727, 6

7 nd the error in this pproximtion is E clssicl = , nd for comprison, nd get the following results in Tble 2 Time Error BAL (Experiment) BAL (lssicl) Frctionl pproch Tble 2: Experimentl dt vs theoreticl model. Now we show tht, if we consider the problem modeled by system of frctionl differentil equtions, we obtin curve tht better fit with the experimentl results. Let α,β (,1) nd consider the system of frctionl differentil equtions Dα t A(t) = k 1A(t) Dβ t B(t) = k 1A(t) k 2 B(t) A() = A B() =. Using Theorem 7.2 in [2], the solution with respect to A is A(t) = A E α ( k 1 t α ). (5) To determine B, it cn be found s the solution of the frctionl differentil liner eqution Dβ t B(t) = k 2B(t)+k 1 A E α ( k 1 t α ). By Theorem 7.2 in [2], we obtin the solution with respect to B: B(t) = β k 1 A E α ( k 1 (t s) α )s β 1 E β ( k 2s β )ds ( k 1 ) m (t s) mα = βk 1 A s β 1 (n+1)( k 2 ) n s nβ ds Γ(αm+1) Γ(nβ +β +1) m= n= ( k 1 ) m ( k 2 ) n (n+1) t = βk 1 A (t s) mα s nβ+β 1 ds. Γ(αm+1)Γ(nβ +β +1) m=n= We remrk tht, s the expression inside the series is continuous nd uniformly convergent, the previous clcultions re vlid. Let us re-write function B( ). Using the Bet function, the following property B(,b) = 1 t 1 (1 t) b 1 dt,,b >, B(,b) = Γ()Γ(b) Γ(+b), nd doing the chnge of vribles u = s/t, we obtin nd so 1 (t s) mα s nβ+β 1 ds = t mα+nβ+β (1 u) mα u nβ+β 1 du = t mα+nβ+β B(nβ +β,mα+1) = t mα+nβ+βγ(nβ +β)γ(mα+1) Γ(nβ +β +mα+1), B(t) = k 1 A m=n= ( k 1 ) m ( k 2 ) n Γ(nβ +β +mα+1) tmα+nβ+β. (6) 7

8 3.3 Numericl experiments For computtionl purposes we consider the upper bounds n = m = 45 in Eq. (6), nd we intend to determine the frctionl orders α nd β on the intervl (,2), the initil vlue A, nd the prmeters k 1 nd k 2 tht best fit with the dt. To tht purpose, we obtin the vlues with frctionl orders nd the error in this pproximtion is By this method, the efficiency hs incresed A , k 1.643, k 2.88, α , β 1.52, E frctionl In Tble 3 we summrize the results, compring the stndrd pproch, using first-order derivtives, with ours using frctionl-order derivtives. Time Error BAL (Experiment) BAL (lssicl) BAL (Frctionl) Tble 3: Results for the Blood Alcohol Level problem. The grph of Figure 3 illustrtes the Blood Alcohol Level with experimentl dt, clssicl model nd frctionl model Experiment lssicl Model, Error= Frctionl Model, Error= Blood Alcohol Level (mg/l) Time (minutes) Figure 3: The blood lcohol level problem. 8

9 4 Video Tpe problem 4.1 Stndrd pproch For the video tpe problem, we intend to determine the tpe counter redings n(t) t n instnt time t >. For tht, the problem is modeled in the following wy. Let c represents the thickness of the tpe nd v the velocity of the tpe, which is constnt in time. Also, R(t) indictes the rdius of the wheel s is filled by the tpe, nd A(t) mesures the ngle t moment t. The number of revolutions is proportionl to the ngle, tht is, n(t) = ka(t), for some positive constnt k. In [3], n ordinry differentil eqution is obtined to describe the behvior of the model: A v (t) = R() cv, whereb = bt+1 πr 2 (). (7) Using the initil condition A() =, the solution is obtined: A(t) = 2v br() ( bt+1 1). In conclusion, the tpe counter redings is given by the expression 4.2 Frctionl pproch n(t) = ( bt+1 1), where = 2kv br(). (8) In this section we consider the video tpe problem modeled by frctionl differentil eqution of order α (,1): D α t A(t) = A() = v R() bt+1, whereb = cv πr 2 () Applying the frctionl integrl opertor I α t to both sides of the frctionl differentil eqution, nd using Lemm 1, tht reds s we obtin tht the solution is given by n(t) = p Γ(α) I α t D α t A(t) = A(t) A(), ifα (,1), 4.3 Numericl experiments (t s) α 1 bs+1 ds, wherep = kv R(), b = cv πr 2 (). (9) Using the experimentl dt vilble t tht counts the number of revolutions every 5 minutes, from t = until t = 24, we first determine the constnts nd b for which the model given in Eq. (8) best fit the dt. nd the error is given by = nd b =.219, E clssicl When we consider the frctionl pproch, with frctionl order α (,2) for Function (9), we obtin the vlues α =.9917, p = nd b =.199 with error E frctionl

10 For this lst cse, we hve The grph of Figure 4 illustrtes Video Tpe model with experimentl dt, clssicl model nd frctionl model. 15 Experiment lssicl Model, Error= Frctionl Model, Error= Figure 4: The Video Tpe problem. 5 onclusions Frctionl differentil equtions cn describe better certin rel world phenomen. This is understndble since systems re not usully perfect, nd cn be perturbed (like friction, mnipultion, externl forces, etc) nd becuse of it integer-order derivtives my not be dequte to understnd the trjectories of the stte vribles. By considering frctionl derivtives, we hve n infinite choices of derivtive orders tht we cn consider, nd with it determine wht is the frctionl differentil eqution tht better describes the dynmics of the model. We hve seen this in our experimentl dt, where non-integer derivtives llow the solution curve to model more efficiently the problems. Acknowledgements This work ws supported by Portuguese funds through the IDMA - enter for Reserch nd Development in Mthemtics nd Applictions, nd the Portuguese Foundtion for Science nd Technology(FT-Fundção pr iênci e Tecnologi), within project UID/MAT/416/213; third uthor by the ALGORITMI R&D enter nd project PEst-UID/E/319/213. The uthors re very grteful to two referees for mny constructive comments nd remrks. References [1] T.M. Atnckovic, S. Pilipovic, B. Stnkovic nd D. Zoric, Frctionl lculus with Applictions in Mechnics: Vibrtions nd Diffusion Processes. Wiley-ISTE,

11 [2] K. Diethelm, The Anlysis of Frctionl Differentil Equtions. An Appliction-Oriented Exposition Using Differentil Opertors of puto Type, LNM Springer, Volume 24, 21. [3]. Glphin, J. Glembocki nd J. Tompkins, Video Tpe ounters: The Exponentil Mesure of Time (vilble online t [4] A.A. Kilbs, H.M. Srivstv nd J.J. Trujillo, Theory nd Applictions of Frctionl Differentil Equtions. North-Hollnd Mthemtics Studies, 24. Elsevier Science B.V., Amsterdm, 26. [5]. Ludwin, Blood Alcohol ontent, Undergrd. J. Mth. Model.: One + Two 3 (2) Article 1 (211). [6] J.A.T. Mchdo, M.F. Silv, R.S. Brbos, I.S. Jesus,.M. Reis, M.G. Mrcos nd A.F. Glhno. Some Applictions of Frctionl lculus in Engineering, Mthemticl Problems in Engineering, vol. 21, Article ID 63981, 34 pges, 21. [7] F.. Merl, T.J. Royston nd R. Mgin, Frctionl clculus in viscoelsticity: An experimentl study, ommun. Nonliner Sci. Numer. Simul. 15 (4), (21). [8]. Moler, J. Little, nd S. Bngert, Mtlb User s Guide - The Lnguge of Technicl omputing, The MthWorks, Sherborn, Mss. (21). [9] I. Podlubny, Frctionl differentil equtions, Mthemtics in Science nd Engineering, 198. Acdemic Press, Inc., Sn Diego, A, [1] S.G. Smko, A.A. Kilbs nd O.I. Mrichev, Frctionl integrls nd derivtives, trnslted from the 1987 Russin originl, Gordon nd Brech, Yverdon, [11] D.A. Smith, Humn popultion growth: Stbility or explosion? Mth Mg 5 (4), (1977). [12] V.E. Trsov, Applictions of Frctionl lculus to Dynmics of Prticles, Fields nd Medi, Nonliner Physicl Science, Springer, 21. [13] United Ntions, The World t Six Billion Off Site, Tble 1, World Popultion From Yer to Stbiliztion, 5,