Equivalence of the Initialized Riemann-Liouville Derivatives and the Initialized Caputo Derivatives arxiv: v1 [math.

Transcription

1 Equivalence of the Initialized Riemann-Liouville Derivatives and the Initialized Caputo Derivatives arxiv: v1 [math.gm] 14 Nov 218 Jian Yuan 1, College of Mathematic and Information Science, Shandong Technology and Business University, Yantai 2645, P.R.China Guozhong Xiu, Bao Shi Institute of System Science and Mathematics, Naval Aeronautical University, Yantai 2641, P.R.China Abstract Initialization of fractional differential equations remains an ongoing problem. In recent years, the initialization function approach and the infinite state approach provide two effective ways to deal with this problem. The purpose of this paper is to prove the equivalence of the initialized Riemann- Liouville derivatives and the initialized Caputo derivatives with arbitrary orders. By synthesizing the above two initialization theories, the diffusive representations of the two initialized derivatives with arbitrary orders are derived. Laplace transforms of the two initialized derivatives are shown to be equal. As a result, the two most commonly used derivatives are proved to be equivalent when initial conditions are properly imposed. Keywords: Fractional calculus, initialization problems, diffusive representation, equivalence of fractional derivatives 1. Introduction Fractional calculus provides a powerful tool to model real-world phenomena exhibiting memory and hereditary properties [1]. Initial conditions of fractional differential equations characterize the past histories of dynamical systems from the starting time to the initial time. As a result, initial Corresponding author. yuanjianscar@gmail.com Preprint submitted to Elsevier November 29, 218

2 conditions of fractional differential equations should be imposed in a different way from that of integer-order differential equations. However, proper initialization of fractional-order systems remains an ongoing problem [2-4]. This issue can date back as far as Riemann s complementary function theory, in which many mathematicians were made confused, including Liouville, Peacoch, Cayley, and Riemann himself [3]. In recent years, the initialization function approach and the infinite state approach provide two effective ways to impose physically coherent initial conditions to fractional systems. In the initialization function theory [3, 5], initialization functions are required to initialize fractional differential equations. The initialization function is a time-varying function and a generalization of the constant of integration required for the order-one integral. The Riemann-Liouville fractional derivatives and the Caputo derivatives can be initialized using the initialization function approach. In the infinite state theory, the Riemann-Liouville fractional integral with order between and 1 is viewed a linear system which is characterized by the impulse function and excited by the integrand function. The linear system is termed as the fractional integrator, which is also an exactly equivalent frequency distributed differential system. As a result, initial conditions of fractional systems can be represented by the distributed initial conditions of the fractional integrator [6-8]. Equivalence and compatibility of the above initialization theories have been proved in [9, 1]. The Riemann-Liouville derivative and the Caputo derivative are the two most commonly used derivatives in the real-world modeling of factional systems. However, which derivative to choose is a trial and error process. In many applications, the two derivatives are expected to be equivalent as long as the initial conditions are properly taken into account [11, 12]. In [3] and [13], Laplace transforms of the Riemann-Liouville derivative and the Caputo derivative have been derived and shown to be equal in a special case where the history function is a constant and the fractional order lies between and 1. In this paper, we present the equivalence and compatibility of the initialized Riemann-Liouville derivatives and the initialized Caputo derivatives with arbitrary fractional orders. By synthesizing the initialization function approach and the infinite state approach, the diffusive representations of the initialized Riemann-Liouville derivatives and the initialized Caputo derivatives are derived. Laplace transforms of the two initialized derivatives are shown to be equal. Consequently, the Riemann-Liouville derivatives and the Caputo derivatives are proved to equivalent when initial conditions are 2

3 properly imposed. The rest of this paper is organized as follows. Section 2 revisits the diffusive representation for the initialized Riemann-Liouville fractional integral. Section 3 presents the equivalence of the initialized Riemann-Liouville derivatives and the initialized Caputo derivatives with order between and 1. Section 4 shows the equivalence of the two initialized derivatives with order between 1 and 2. Section 5 proves the equivalence of the two initialized derivatives with arbitrary orders. Finally, the paper is concluded in Section Preliminaries 2.1. Diffusive model of the fractional integrator [7] The Riemann-Liouville fractional integral of a function f (t)with order < α < 1 is defined as t I α t f (t) = 1 Γ(α) t t (t τ) α 1 f (τ)dτ (1) Eq.(1)can also be viewed as a convolution between f (t) and the impulse response h α (t) = tα 1 Γ(α), i.e., t I α t f(t) = h α(t) f (t). From this viewpoint, the fractional integral can be obtained as the output of a linear system, which is characterized by the impulse response h α (t) and excited by f (t): x(t) = h α (t) f (t). This linear system is called the fractional integrator. The frequency distributed state z(ω, t)of the fractional integrator verifies the following ordinary differential equation: z(ω, t) t and the output of the fractional integrator is x(t) = = ωz(ω,t)+f(t) (2) + µ α (ω)z(ω,t), (3) 3

4 where the elementary frequency ω is ranging from to, and µ α (ω) = sinαπ ω α. π The relations (2) and (3) are called the frequency distributed model or diffusive model of the fractional integrator The initialized fractional integrals [14, 15] In the initialization theory for fractional calculus, the integrand of Riemann- Liouville fractional integral is represented as a composite function, t < a v(t) = f in (t), a t f(t), t > where t = a is the starting time of integral, t = is the initial time, f i (t) is the function describing the behavior during the initialization period [ a, ], f(t) is the function of primary interest after the initial time t =. The initialized Riemann-Liouville fractional integral of order α is defined as Dt α f (t) = 1 t (t τ) α 1 f (τ)dτ +ψ(t), t > Γ(α) or D α t f (t) = d α t f (t)+ψ(t), t >. where D α t f (t) is called the initialized fractional integral, d α t f (t) is called the uninitialized fractional integral, ψ(t) is called the initialization function describing the hereditary effect of the past Diffusive representation of the initialized Riemann-Liouville fractional integrals [1] Lemma 1. TheuninitializedRiemann-Liouvillefractionalintegral d α t f(t) with order < α < 1 is equivalent to the diffusive model with zero initial condition: z(ω, t) = ωz(ω,t)+f(t) t z(ω,t ) = and d α t f(t) = + 4 µ α (ω)z(ω,t).

5 Lemma 2. TheinitializedRiemann-Liouvillefractionalintegral Dt α f(t) with order < α < 1 is equivalent to the diffusive model with distributed initial condition: and z(ω, t) t z(ω,) = D α t f(t) = = ωz(ω,t)+f(t), a + e ωτ f in (τ)dτ µ α (ω)z(ω,t). 3. Equivalence of the two Derivatives with order between and Diffusive model for the initialized Riemann-Liouville derivatives with order between and 1 TheinitializedRiemann-Liouvillefractionalderivative D α t f(t)withorder < α < 1 is defined in [14] as Dt α f (t) = d dt Dt α 1 f (t). By virtue of Lemma 2, the diffusive representation of the initialized Riemann- Liouville fractional integral Dt α 1 f(t) is Dt α 1 f (t) = µ 1-α (ω)z (ω,t) where z (ω,t) satisfies { z (ω,t) = ωz t (ω,t)+f(t) z (ω,) = a e ωτ f in (τ)dτ As a result, the diffusive model for the initialized Riemann-Liouville derivatives is D α t f (t) = d dt (4) µ 1-α (ω)z (ω,t) (5) Taking Laplace transform of the first equation of Eq.(4), we have Z (ω,t) = z (ω,)+f (s) (6) 5

6 Taking Laplace transform of the initialized Riemann-Liouville derivative from Eq.(5), we have Dt α f (t) } { d } = L µ 1-α (ω)z (ω,t) dt = s µ 1-α (ω)z (ω,t) Substituting Eq. (6) into Eq.(7), we obtain Dt α f (t) } = sf (s) Because Eq.(8) becomes µ 1-α (ω) +s µ 1-α (ω)z (ω,) Dt α f (t)} = s α F (s)+s µ α (ω) = s α, µ 1-α (ω)z (ω,) µ 1-α (ω)z (ω,) (7) µ 1-α (ω)z (ω,) µ 1-α (ω)z (ω,) ( ) s = s α F (s)+ 1 µ 1-α (ω)z (ω,) = s α F (s) ωµ 1-α (ω)z (ω,) 3.2. Diffusive model for initialized Caputo derivatives with order between and 1 The initialized Caputofractional derivative Dt α f(t)withorder < α < 1 is defined in [14] as C Dα t f (t) = Dt α 1 f (t). By virtue of Lemma 2, the diffusive representation of the initialized Caputo fractional derivative is C Dα t f (t) = (8) (9) µ 1-α (ω)z C (ω,t) (1) 6

7 where z C (ω,t) satisfies { zc (ω,t) = ωz t C (ω,t)+f(t) z C (ω,) = a e ωτ f in (τ)dτ Taking Laplace transform of Eq.(11), we have Z C (ω,t) = sf (s) f ()+z C (ω,) Taking Laplace transform of Eq.(1), we have L { { } C Dt α f (t)} = L µ 1-α (ω)z C (ω,t) = Substituting Eq.(12) into Eq.(13), we obtain L { C D α t f (t) } = sf (s) f () = s α F (s)+ µ 1-α (ω)z C (ω,t) µ 1-α (ω) + µ 1-α (ω) µ 1-α (ω)z C (ω,) µ 1-α (ω)z C (ω,) s α 1 f () In terms of Eq.(11) and integrating by parts, we have z C (ω,) = a In terms of Eq.(4), we have Then µ 1-α (ω)z C (ω,) = f () e ωτ f in (τ)dτ = e ωτ f in (τ) τ= τ= a ω a e ωτ f in (τ)dτ (11) (12) (13) (14) z C (ω,) = f () ωz Rl (ω,) (15) = s α 1 f () µ 1 α (ω) 7 µ 1-α (ω)z (ω,) ωµ 1-α (ω)z (ω,) (16)

8 Substituting Eq.(16) into Eq.(14), we get L { C D α t f (t)} = s α F (s) From Eq.(9) and Eq.(17), it is clear that ωµ 1-α (ω)z (ω,) (17) D α t f (t) } = L { C D α t f (t) }. (18) Thus Dα t f (t)=c Dα t f (t). (19) Eq.(19) shows the equivalence of the initialized Riemann-Liouville derivatives and Caputo derivatives with order between and Equivalence of the two Derivatives with order between 1 and Diffusive model for the initialized Riemann-Liouville derivatives with order between 1 and 2 TheinitializedRiemann-Liouvillefractionalderivative D α t f(t)withorder 1 < α < 2 is defined as Dt α d2 f (t) = D α 2 dt 2 t f (t). By virtue of Lemma 2, the diffusive representation of the initialized Riemann- Liouville fractional integral D α 2 f(t) is D α 2 t f (t) = t µ 2 α (ω)z (ω,t) where z (ω,t) satisfies Eq.(4). As a result, the diffusive model for the initialized Riemann-Liouville derivatives is Dt α f (t) = d2 dt 2 Taking Laplace transform of Eq.(2), we have L { { d Dt α 2 } f (t)} = L µ dt 2 1-α (ω)z (ω,t) = s 2 µ 2-α (ω)z (ω,t) s µ 2-α (ω) z (ω,t) t t= 8 µ 1-α (ω)z (ω,t) (2) µ 2-α (ω)z (ω,) (21)

9 Substituting Eq.(6) into Eq.(21), we have Dt α f (t)} = sf (s) Because s µ 2-α (ω) µ 2-α (ω)z (ω,) +s2 µ 2-α (ω)z (ω,) µ 2 α (ω) z (ω,t) t t= µ 2 α (ω) = sα 2, By substituting the first equation of Eq.(4) into Eq.(22), we obtain L { ( ) s Dα t f (t)} = s α 2 F (s)+ µ 2-α (ω)z (ω,) f () = s α F (s)+ f () µ 2 α (ω) ω 2 µ 2-α (ω)z (ω,) µ 2-α (ω) (22) (23) 4.2. Diffusive model for the initialized Caputo derivatives with order between 1 and 2 The initialized Caputofractional derivative Dt α f(t)withorder 1 < α < 2 is defined as C Dα t f (t) = Dt α 2 f (t) By virtue of Lemma 2, the diffusive representation of the initialized Caputo fractional derivative is C Dα t f (t) = where z C (ω,t) satisfies { zc (ω,t) = ωz t C (ω,t)+f (t) z C (ω,) = a e ωτ f in (τ)dτ µ 2-α (ω)z C (ω,t) (24) 9 (25)

10 Taking Laplace transform of the first equation in Eq.(25), we have Z C (ω,t) = s2 F (s) sf () f ()+z C (ω,). (26) In terms of the initial conditions of Eq.(25), z C (ω,) = a e ωτ f in (τ)dτ = e ωτ f in (τ) τ= ω e ωτ f τ= a in (τ)dτ a = f () ω e ωτ f in (τ)dτ a = f () ωf ()+ω 2 z R (ω,) (27) Substituting Eq.(27) into Eq.(26), we have Z C (ω,t) = s2 F (s) ()f ()+ω 2 z (ω,). (28) Taking Laplace transform of Eq.(24)and substituting Eq.(28) into it, we obtain L { C Dt α f (t) } { = L = = s 2 F (s) f () } µ 2-α (ω)z C (ω,t) µ 2-α (ω)z C (ω,t) = s α F (s)+ f () µ 2 α (ω) + µ 2-α (ω) ω 2 µ 2-α (ω)z (ω,) µ 2-α (ω) ω 2 µ 2-α (ω)z (ω,) (29) 1

11 From Eq.(23) and Eq.(29), it is clear that Thus, D α t f (t)} = L { C D α t f (t)}. (3) Dα t f (t)=c Dα t f (t). (31) Eq.(31) shows the equivalence of the initialized Riemann-Liouville derivatives and Caputo derivatives with order between 1 and Equivalence of the two Derivatives with arbitrary orders 5.1. Diffusive model for initialized Riemann-Liouville derivatives with order between n and n 1 TheinitializedRiemann-Liouvillefractionalderivative D α t f(t)withorder n < α < n 1 is defined as Dt α f (t) = dn D α n dt n t f (t). By virtue of Lemma 2, the diffusive representation of the initialized Riemann- Liouville fractional integral Dt α n f (t) is D α n t f (t) = µ n α (ω)z (ω,t) where z (ω,t) satisfies Eq.(4). As a result, the diffusive model for the initialized Riemann-Liouville derivative is dn Dα t f (t) = µ dt n n-α (ω)z (ω,t) (32) Taking Laplace transform of Eq.(32), we have 11

12 L { { d Dt α n } f (t)} = L µ dt n n α (ω)z (ω,t) = s n µ n α (ω)z (ω,t) s n 1 s n 2 µ n α (ω)z (ω,) s n 3 µ n α (ω)z (ω,) µ n α (ω)z (n 1) (ω,) Substituting Eq.(6) into Eq.(33), andapplying Dα t f (t)} = s α F (s)+s n µ n α (ω)z (ω,) (33) µ 2 α (ω) = sα 2, we have µ n α (ω)z (ω,) s n 1 µ n α (ω)z (ω,) s n 2 µ n α (ω)z (ω,) s n 3 µ n α (ω)z (ω,) µ n α (ω)z (n 1) (ω,) (34) From Eq.(4), we can calculate high-order derivatives of z (ω,t): 3 z (ω,t) t 3 n 1 z (ω,t) t n 1 2 z (ω,t) t 2 = ω 2 z (ω,t) ωf (t)+f (t) = ω 3 z (ω,t)+ω 2 f (t) ωf (t)+f (t) = ( ω) n 1 z (ω,t)+( ω) n 2 f (t)+( ω) n 3 f (t) + +f (n 2) (t) 12

13 Substituting the above derivatives into Eq.(34), we obtain D α t f (t)} = s α F (s)+s n where µ n α (ω)z (ω,) s n 1 µ n α (ω)z (ω,) s n 2 µ n α (ω)[ ωz (ω,)+f()] s n 3 µ n α (ω) [ ω 2 z (ω,) ωf ()+f () ] µ n α (ω) [ ( ω) n 1 z (ω,)+( ω) n 2 f () +( ω) n 3 f ()+ +f (n 2) () ] [ s = s α n F (s)+ µ n α (ω)z (ω,) sn 1 +ωs n 2 ω 2 s n 3 + +( ω) n 2 s+( ω) n 1] + = s α F (s)+ µ n α (ω){ 1 } µ n α (ω)z (ω,) ( ω)n µ n α (ω){ 1 } 1 = [ s n 2 +( ω)s n 3 +( ω) 2 s n 4 + ( ω) n 2] f () (35) + [ s n 3 +( ω)s n 4 + ( ω) n 3] f ()+ +f (n 2) () (36) 5.2. Diffusive model for initialized Caputo derivatives with order between n and n 1 The initialized Caputo fractional derivative Dt α f(t) with order n 1 < α < n is defined as C Dα t f (t) = Dt α n f (n) (t). 13

14 By virtue of Lemma 2, the diffusive representation of the initialized Caputo fractional derivative is C Dα t f (t) = µ n α (ω)z C (ω,t) (37) where z C (ω,t) satisfies { zc (ω,t) = ωz t C (ω,t)+f (n) (t) z C (ω,) = a e ωτ f (n) in (τ)dτ (38) Taking Laplace transform of Eq.(38), we have Z C (ω,s) = 1 [ zc (ω,)+s n F (s) s n 1 f () s n 1 f () sf (n 2) () f (n 1) () ] (39) In terms of the initial condition of Eq.(38) and integrating by parts, we have z C (ω,) = f (n 1) ()+( ω)f (n 2) ()+( ω) 2 f (n 3) ()+( ω) 3 f (n 4) () + +( ω) n 2 f ()+( ω) n 1 f ()+( ω) n z (ω,) Substituting Eq.(4) into Eq.(39), we have Z C (ω,s) = 1 { s n F (s)+ [ ( ω) n 1 s n 1] f () + [ ( ω) n 2 s n 2] f ()+ (4) +[( ω) s]f (n 2) ()+( ω) n Z (ω,) } (41) Taking Laplace transform of Eq.(37) and substituting Eq.(41) into it, we obtain L { { } C Dt α f (t)} = L µ n α (ω)z C (ω,t) = = s α F (s)+ µ n α (ω)z C (ω,t) ( ω) n µ n α(ω)z (ω,) µ n α (ω){ 2 } 14 (42)

15 where [ ] [ ] s n 1 ( ω) n 1 s n 2 ( ω) n 2 2 = f ()+ f ()+ +f (n 2) (). (43) From Eqs. (36)and(43), it is clear that 1 = 2, and Thus, D α t f (t)} = L { C D α t f (t)}. (44) Dt α f (t) = C Dα t f (t). (45) Eq.(45) shows the equivalence of the initialized Riemann-Liouville derivatives and Caputo derivatives with arbitrary orders α. 6. Conclusions This paper has presented the equivalence of the initialized Riemann- Liouville derivatives and the initialized Caputo derivatives with arbitrary orders. By synthesizing the initialization function theory and the infinite state theory, the diffusive representations of the two initialized derivatives have been obtained. Laplace transforms of the two initialized derivatives with arbitrary orders have been shown to be equal from Eqs.(18) (3) and (44). As a result, the two most commonly used derivatives have been shown equivalent from Eqs.(19) (31) and (45). Acknowledgements All the authors acknowledge the valuable suggestions from the peer reviewers. This work was supported by the National Natural Science Foundation of China (Grant No ). References [1] Podlubny, I. Fractional Differential Equations. in Mathematics in Science and Engineering [2] Fukunaga, M. and N.J.N.D. Shimizu, Role of Prehistories in the Initial Value Problems of Fractional Viscoelastic Equations. Nonlinear Dynamics, (1-4): p [3] Lorenzo, C.F. and T.T. Hartley, Initialization of Fractional-Order Operators and Fractional Differential Equations. Journal of Computational and Nonlinear Dynamics, 28. 3(2): p

16 [4] Trigeassou, J.C. and N. Maamri, Initial conditions and initialization of linear fractional differential equations. Signal Processing, (3): p [5] Lorenzo, C.F. and T.T. Hartley, Dynamics and Control of Initialized Fractional-Order Systems (1-4): p [6] Trigeassou, J.C., N. Maamri, and A. Oustaloup, The infinite state approach: Origin and necessity. Computers & Mathematics with Applications, (5): p [7] Trigeassou, J. C., Maamri, N., Sabatier, J., Oustaloup, A., State variables and transients of fractional order differential systems. Computers & Mathematics with Applications, (1): p [8] Trigeassou, J. C., Maamri, N., Sabatier, J., Oustaloup, A., A lyapunov approach to the stability of fractional differential equations. Signal Processing, (3): p [9] T.T. Hartley, Lorenzo, C.F., Trigeassou, J. C., Maamri, N., Equivalence of History-Function Based and Infinite-Dimensional-State Initializations for Fractional-Order Operators. Journal of Computational and Nonlinear Dynamics, (4): p [1] Yuan, J., Zhang, Y., Liu, J., Shi, B., Equivalence of initialized fractional integrals and the diffusive model. Journal of Computational and Nonlinear Dynamics, (3): p [11] Paola, M. D., Fiore, V., Pinnola, F. P., Valenza, A., On the influence of the initial ramp for a correct definition of the parameters of fractional viscoelastic materials. Mechanics of Materials, (1): p [12] Bagley, R., On the equivalence of the Riemann-Liouville and the Caputo fractional order derivatives in modeling of linear viscoelastic materials. Fractional Calculus and Applied Analysis, 27. 1(2): p [13] Lorenzo, C.F. and T.T. Hartley, Time-varying initialization and Laplace transform of the Caputo derivative : with order between zero and one. ASME 211 International Design Engineering Technical Conferences and Computers and Information in Engineering Conference

17 [14] Lorenzo, C.F. and T.T. Hartley, Initialized fractional calculus. 2, NASA/TP [15] Lorenzo, C.F. and T.T. Hartley, Initialization, conceptualization, and application in the generalized (fractional) calculus. Critical Reviews in Biomedical Engineering, 27, 35(6): p