首页资源分类电子电路模拟及混合电路 > introduction to fractional derivative

introduction to fractional derivative

已有 456417个资源

下载专区

上传者其他资源

文档信息举报收藏

标    签: derivative

分    享:

文档简介

introduction to fractional derivative pdf

文档预览

Introductory Notes on Fractional Calculus xuru.org xuru@xuru.org July 31, 2006 Abstract Our understanding of Nature relies on calculus, which in turn relies on the intuitive concept of the derivative. It’s descriptive power comes from the fact that it analyses the behavior at scales small enough that its properties change linearly, so avoiding complexities that arise at larger ones. Fractional Calculus generalizes this concept from integer to non-integer order. Despite it seems not to have significant applications in fundamental physics, research on this core concept could be valuable in understanding Nature. These notes comprise an introduction to the field. 1 Introduction Fractional Calculus is the branch of calculus that generalizes the derivative of a function to non-integer order, allowing calculations such as deriving a function to 1/2 order. Despite “generalized” would be a better option, the name “fractional” is used for denoting this kind of derivative. !Note: These notes are not a summary of standard methods or an academic text, so some conventional developments are omitted while some other are introduced. The same happens with notation. The derivative of a function f is defined as D1 f (x) = lim f (x + h) − f (x) h→0 h (1.1) Iterating this operation yields an expression for the n-st derivative of a function. As can be easily seen -and proved by induction- for any natural number n, Dn f (x) = lim h−n n (−1)m n f (x + (n − m)h) h→0 m m=0 (1.2) 2 1 Introduction where n n! = m m!(n − m)! (1.3) or equivalently, Dn f (x) = lim h−n n (−1)m n f (x − mh) h→0 m m=0 (1.4) The case of n = 0 can be included as well. Such an expression could be valuable for instance in a simple program for plotting the n-st derivative of a function. Viewing this expression one asks immediately if it can be generalized to any non-integer, real or complex number n. There are some reasons that can make us think so, 1. The fact that for any natural number n the calculation of the n-st derivative is given by an explicit formula (1.2) or (1.4). 2. That the generalization of the factorial by the gamma function allows n n! Γ(n + 1) = = m m!(n − m)! Γ(m + 1)Γ(n − m + 1) (1.5) which also is valid for non-integer values. 3. The likeness of (1.2) to the binomial formula n (a + b)n = n an−m bm m m=0 (1.6) which can be generalized to any complex number α by ∞ (a + b)α = Γ(α + 1) aα−n bn n! Γ(α − n + 1) n=0 which is convergent if (1.7) |b| < a (1.8) There are some desirable properties that could be required to the fractional derivative, 1. Existance and continuity for m times derivable functions, for any n which modulus is equal or less than m. 2 Exponentials 3 2. For n = 0 the result should be the function itself; for n > 0 integer values it should be equal to the ordinary derivative and for n < 0 integer values it should be equal to ordinary integration -regardless the integration constant. 3. Iterating should not give problems, Dα+β f (x) = Dα Dβ f (x) (1.9) 4. Linearity, Dα(af (x) + bg(x)) = a Dα f (x) + b Dα g(x) (1.10) 5. Allowing Taylor’s expansion in some other way. 6. Its characteristic property should be preserved for the exponential function, Dα ex = ex (1.11) 2 Exponentials The case of the exponential function is specially simple and gives some clues about the generalization of the derivatives. Following (1.2), Dα eax = lim h−α α (−1)n α ea(x+(α−n)h) = h→0 n n=0 = eax lim h−α α (−1)n α (eah)α−n = h→0 n n=0 = eax lim h−α (eah − 1)α = h→0 = aα eax (2.1) the above limit exists for any complex number α. However, it should be noted that in the substitution of the binomial formula a natural number has been considered. We shall deal with this problem later to get our first generalization of the derivative. Applying this to the imaginary unit, Dα cos(x) + i Dα sin(x) = Dα eix = iαeix = = απi e2 eix = ei (x + απ 2 ) = απ απ = cos x + + i sin x + 2 2 (2.2) 4 2 Exponentials and Dα cos(x) − i Dα sin(x) = Dα e−ix = (−i)αe−ix = = −απ i e2 e−ix = e−i(x + απ 2 ) = απ απ = cos x + − i sin x + 2 2 (2.3) solving this system we have the next definition for the sine and cosine derivatives, Dα sin(x) = sin απ x+ 2 (2.4) and Dα cos(x) = cos απ x+ 2 (2.5) We could expect these relations for the sine and cosine derivatives to be maintained in the generalization of the derivative. Applying the above method we also can calculate the following, Dα cos(ax) + i Dα sin(ax) = Dα eiax = (ai)αeiax = = aα e απ 2 i eiax = aα ei (a x+ απ 2 ) = = aα cos απ ax + + iaα sin απ ax + 2 2 (2.6) and Dα cos(ax) − i Dα sin(ax) = Dα e−iax = (−ai)αe−iax = = aα e− απ 2 i e−iax = aα e−i (a x+ απ 2 ) = = aα cos απ ax + − iaα sin απ ax + 2 2 (2.7) Thus, Dα sin(ax) = aα sin απ ax + 2 (2.8) and Dα cos(ax) = aα cos απ ax + 2 (2.9) 3 Powers 5 Indeed, the above result of the exponential can be applied to any function that can be expanded in exponentials ∞ f (x) = an einx =⇒ n=−∞ ∞ ∞ =⇒ Dα f (x) = an Dα einx = an (ni)α einx = n=−∞ n=−∞ ∞ ∞ = an nα e απ 2 i einx = an nα ei (n x + απ 2 ) n=−∞ n=−∞ (2.10) Expanding the function in Fourier series, ∞ f (x) = π einx f (t)e−int dt =⇒ 2π n=−∞ −π ∞ =⇒ Dα f (x) = nα ei (n x+ απ 2 ) π f (t)e−int dt 2π n=−∞ −π (2.11) This method can be useful for calculating fractional derivatives of trigonometric functions. 3 Powers The case of powers of x also has some simplicity that allows its generalization. The case of integer order derivatives n−1 D1 xa = a xa−1 =⇒ Dn xa = xa−n (a − m) = a! xa−n (a − n)! m=0 (3.1) can be easily generalized to non-integer order derivatives Dα xa = Γ(a + 1) xa−α Γ(a − α + 1) (3.2) which can be applied to any function that can be expanded in powers of x ∞ f (x) = an xn =⇒ n=−∞ =⇒ Dα f (x) = ∞ an Dα xn = ∞ an Γ(n Γ(n − + α 1) + 1) xn−α n=−∞ n=−∞ (3.3) 6 4 Binomial Formula Expanding the function in Taylor series, f (x + a) = ∞ Dn f (a) xn =⇒ n! n=0 ∞ =⇒ Dα f (x + a) = Dn f (a) xn−α Γ(n − α + 1) n=0 (3.4) or expanding it in Laurent series, 1 f (x + a) = ∞ xn 2πi f (t + a) tn+1 dt =⇒ n=−∞ C =⇒ Dα f (x + a) = 1 ∞ Γ(n + 1) xn−α 2πi Γ(n − α + 1) f (t + a) tn+1 dt n=−∞ C (3.5) This can be an useful tool for calculating fractional derivatives. However, we should compare these results of powers with the previous results of exponentials to see if they agree. With the result of exponentials (2.1), Dα eax = aα eax = aα ∞ (a x)n = ∞ an+α xn n! n! n=0 n=0 (3.6) but with the result of powers (3.4), Dα eax = Dα ∞ (a x)n = ∞ an xn−α n! Γ(n − α + 1) n=0 n=0 (3.7) If we compare both results, we see that they only agree for integer values of α. We shall see later where these discrepancies come from, and how they can be avoided. 4 Binomial Formula In the expression (1.2) the exponentian function allows the substitution of the binomial formula as done in (2.1), but this is not possible for any given function. For applying this substitution we require the following displacement operator, whose iteration yields dh f (x) = f (x + h) dhα f (x) = f (x + αh) (4.1) (4.2) 5 Functions of the Derivative 7 what allows the application of the binomial formula (1.6) for natural numbers and the generalized binomial formula (1.7) for complex numbers, Dn f (x) = lim h−n n (−1)m n f (x + (n − m)h) = h→0 m m=0 = lim h−n h→0 n (−1)m n m dhn−m f (x) = m=0 = lim dh − 1 n f (x) h→0 h (4.3) so that for any complex number α it can be generalized Dα f (x) = lim dh − 1 α f (x) h→0 h (4.4) This sheds more light on the derivative and its generalization. Using the expression of the generalized binomial formula (1.7) for non-integer numbers, Dα f (x) = lim dh − 1 α f (x) = h→0 h = lim h−α h→0 ∞ n! Γ(α + 1) Γ(α − n + 1) (−1)n dhα−n f (x) = n=0 ∞ = lim h−α (−1)n Γ(α + 1) f (x + (α − n)h) h→0 n! Γ(α − n + 1) n=0 (4.5) In the case of integer values the summatory only extends α terms and it is equal to the ordinary derivative. Finally, it is obvious that as h goes to 0 the last equation is equivalent to the following ∞ Dα f (x) = lim h−α (−1)n Γ(α + 1) f (x − nh) h→0 n! Γ(α − n + 1) n=0 (4.6) 5 Functions of the Derivative It is worth to mention that being (4.4) the expression of the α-st derivative of a function the derivative itself is D = lim dh − 1 h→0 h (5.1) And being the α-st derivative or iterating the differentiation α times powering it to α, applying other functions to the derivative could be also 8 5 Functions of the Derivative considered. If the function applied to the derivative can be expanded in powers of x, ∞ ∞ g(x) = an xn =⇒ g(D)f (x) = an Dn f (x) n=−∞ n=−∞ (5.2) The result is a weighted sum of different order derivatives. These functions of the derivative are usually known as “formal differential operators”. As an example, the exponential of the derivative applied to the exponential would give the following result that could be valuable for calculating functions of the derivative when both f and g can be expanded in exponentials eax = ∞ an xn =⇒ n! n=0 =⇒ eaD ebx = ∞ an Dn ebx = ∞ an bn ebx = eab ebx = eb(x+a) n! n! n=0 n=0 (5.3) If both functions f and g can be expanded in positive powers of x, ∞ ∞ g(x) = an xn , f (x) = bn xn =⇒ n=0 n=0 ∞ ∞ ∞ ∞ =⇒ g(D)f (x) = an Dn bm xm = an bm Dn xm = n=0 m=0 n=0 m=0 = ∞ an ∞ bm Γ(m Γ(m − + n 1) + 1) xm−n = n=0 m=0 = ∞ an ∞ bl+n Γ(l + Γ(l n + + 1) 1) xl = n=0 l=−n = ∞ an ∞ bl+n (l + n)! l! xl = n=0 l=0 ∞ xl ∞ = l! (l + n)! an bl+n l=0 n=0 (5.4) Thus, ∞ g(D)f (x) = cn xn n=0 1∞ cn = n! (n + m)! am bn+m m=0 (5.5) At a first glance this seems not very interesting, but interesting properties could be hidden under this apparent mess. We shall return later to this point with the help of other tools. 6 Grunwald-Letnikov Derivative 9 6 Grunwald-Letnikov Derivative Grunwald-Letnikov derivative or also named Grunwald-Letnikov differinte- gral, is a generalization of the derivative analogous to our generalization by the binomial formula (4.6), but it is based on the direct generalization of the equation (1.4). The idea behind is that h should approach 0 as n approaches infinity, Dn f (x) = lim h−n n (−1)m n f (x − mh) = h→0 m m=0 n = lim h−n (−1)m n! f (x − mh) = h→0 m! Γ(n − m + 1) m=0 ∞ = lim h−n (−1)m n! f (x − mh) = h→0 m! Γ(n − m + 1) m=0 x−a h = lim h−n (−1)m n! f (x − mh) h→0 m! Γ(n − m + 1) m=0 (6.1) being a negative infinity. However, the generalization is done so that any wished a less than x can be chosen. The above equation can be generalized now to get a formula equivalent to (4.6) x−a h Dα f (x) = lim h−α (−1)m Γ(α + 1) f (x − mh) h→0 m=0 m! Γ(α − m + 1) or equivalently, (6.2) Dα f (x) = lim n→∞ n αn (−1)m Γ(α + 1) x−a f x−m x−a m! Γ(α − m + 1) n m=0 (6.3) These generalizations have the inconvenient that if α is a negative integer they are not valid for the gamma function is infinity in negative integers and zero. However, for negative α values the following can be used −n n+m−1 n+m−1 l −n = l=−n−m+1 = m m! (−l) (−1)m l l=n = m! l=n = m! (−1)m (n + m − 1)! (−1)m Γ(n + m) = = m!(n − 1)! m! Γ(n) (6.4) so that for negative values equations (6.2) and (6.3) can be changed to x−a D−α f (x) = lim h−α h Γ(α + m) f (x − mh) h→0 m! Γ(α) m=0 (6.5) 10 7 Riemann-Liouville Derivative and D−α f (x) = lim n α n Γ(α + m) x−a f x−m n→∞ x − a m! Γ(α) n m=0 (6.6) The reason why x − a is considered in the above equations instead of a single positive number can be seen considering the case of α = −1, which must give the same result that integrating the function, x f (t)dt = D−1 f (x) = lim h n Γ(m + 1) f (x − mh) = a n→∞ h→0 m! Γ(1) m=0 n nh = lim hf (x − mh) = lim f (x − t)dt = n→∞ n→∞ h→0 m=0 h→0 0 x = lim f (t)dt =⇒ n→∞ h→0 x−nh x−a x−a =⇒ a = x − nh =⇒ n = , h= h n (6.7) Yes, in fractional derivatives limits must be considered as they must be in integrals. Limits only vanish with integer order derivatives. This also means that fractional derivatives are nonlocal, which may be the reason that makes this kind of derivatives less useful in describing Nature. We shall see later more about the integration -or better say differentiation- limits. In the case of the derivative generalized by the binomial formula (4.6), since n goes to infinity regardless of h, x − a must be infinity, so that the derivative defined in (4.4) and (4.6) is equivalent to the Grunwald-Letnikov derivative with a lower limit of negative infinity. 7 Riemann-Liouville Derivative Riemann-Liouville derivative is the most used generalization of the derivative. It is based on Cauchy’s formula for calculating iterated integrals. If the first integral of a function, which must equal to deriving it to −1, is as follows x D−1 f (x) = f (t)dt 0 (7.1) 7 Riemann-Liouville Derivative 11 the calculation of the second can be simplified by interchanging the integration order x t2 xx x x D−2 f (x) = f (t1)dt1 dt2 = f (t1)dt2 dt1 = f (t1) dt2 dt1 = 00 0 t1 0 t1 x = f (t)(x − t)dt (7.2) 0 This method can be applied repeatedly, resulting in the following formula for calculating iterated integrals, x D−n f (x) = 1 f (t)(x − t)n−1 dt (n − 1)! 0 (7.3) Now this can be easily generalized to non-integer values, in what is the Riemann-Liouville derivative, x Dα f (x) = 1 Γ(−α) f (t) (x − t)α+1 dt 0 (7.4) Note however that in the above formulas the election of 0 as the lower limit of integration has been arbitrary, and any other number could be chosen. Generally, the election of the integration limits in this and other generalizations of the derivative is indicated with subscripts. The RiemannLiouville derivative with the lower integration limit a would be x aDxα f (x) = 1 Γ(−α) f (t) (x − t)α+1 dt a (7.5) The problem with this generalization is that if the real part of α is positive or zero the integral diverges. So it only can be used to calculate generalized integrals. However, this can be solved easily by deriving first by ordinary derivative more than the amount necessary, thus making the remaining necessary differentiation negative and then applying the generalized derivative for completing the rest in which will be a negative differentiation, aDxα f (x) = Dn aDxα−n f (x) n = (α) + 1 (7.6) Once seen how integration limits must be specified in fractional derivatives, we can return to the derivative of the exponential (2.1) and the derivative of powers (3.2) to see how their disagreement (3.7) was caused by having 12 8 Domain Transforms different limits. With the exponential, supposing b a positive number, x b−1 ebx = aDx−1 ebx = ebx dx = ebx − eba =⇒ eba =0 =⇒ b b a =⇒ a = −∞ (7.7) while with the powers x xb+1 b+1 = aDx−1 xb = xb dx = xb+1 − ab+1 =⇒ ab+1 = 0 =⇒ b+1 b+1 a =⇒ a = 0 (7.8) So both results were different in the case of the exponential function because the limits of the derivatives were different. And the case of this discrepancy is not casual. Differences in the results of derivatives with different differentiation limits are indeed important. Finally, it can be proved that the Grunwald-Letnikov derivative with any given integration limits is equal to the Riemann-Liouville derivative with the same limits for any complex number α with a negative real part. This important result means that our derivative generalized by the binomial formula in (4.4) and (4.6) is equivalent to the Riemann-Liouville derivative with a lower limit of negative infinity provided that the real part of α is negative, −∞Dxα f (x) = lim h→0 dh − 1 α f (x) = 1 h Γ(−α) x f (t) (x − t)α+1 dt −∞ (7.9) This kind of Riemann-Liouville derivative with a lower limit of negative infinity is known as the Weyl derivative. 8 Domain Transforms The Laplace and Fourier transforms that serve to transform to the frequency domain can be used to get generalizations of the derivative valid for functions that allow such transformations. The Laplace transform is defined by ∞ L f (x) = e−txf (x)dx (8.1) 0 while its inverse transform is a+i∞ L−1 f (x) = 1 2πi etxf (x)dx a−i∞ (8.2) 8 Domain Transforms 13 where a is chosen so that it is greather than the real part of any of the singularities of f (x). An important property of the Laplace transform is related to the transform of the n-st derivative of a function, n−1 L Dn f (x) = tn L f (x) − tm (Dn−m−1 f )(0) m=0 (8.3) In the case that the terms in the summatory are zero the relation is particularly simple, and for those kind of functions the derivative can be generalized so that this property holds true for non-integer values of α (Dl 0Dxα−m f )(0) = 0 l = 0 . . . m − 1 m = (α) + 1 =⇒ =⇒ L Dα f (x) = tα L f (x) (8.4) for which the generalized derivative can be defined as Dα f (x) = L−1 tα L f (x) (8.5) Keeping in mind the result of the generalized derivative of the exponential (2.1), the following development provides an easy understanding of the reasons involved in the above generalization, Dα f (x) = Dα L−1 L f (x) =   a+i∞ ∞ = Dα  1 2πi etx e−txf (x)dxdt =  a−i∞ 0 a+i∞ ∞ 1 = Dα etx e−tx f (x) dxdt = 2πi a−i∞ 0 a+i∞ ∞ 1 = tα etx e−tx f (x) dxdt = 2πi a−i∞ 0 = L−1 tα L f (x) (8.6) On the other hand, the Fourier transform is defined by ∞ F f (x) = e−itxf (x)dx −∞ while its inverse transform is ∞ F −1 f (x) = 1 eitxf (x)dx 2π −∞ (8.7) (8.8) 14 9 Convolution This transform also has an analogous property related to the transform of the n-st derivative of a function, F Dn f (x) = (it)n F f (x) (8.9) and the derivative can be generalized so that this property holds true for non-integer values of α F Dα f (x) = (it)α F f (x) (8.10) yielding the following definition of the generalized derivative Dα f (x) = F −1 (it)α F f (x) (8.11) The application of (2.1) can provide again an easy understanding of the reasons involved in the generalization, Dα f (x) = Dα F −1 F f (x) =   ∞ ∞ = Dα  1 eitx e−itx f (x) dxdt = 2π  −∞ −∞ ∞ ∞ 1 = 2π Dα eitx e−itx f (x) dx dt = −∞ −∞ ∞ ∞ 1 = (it)α eitx e−itx f (x) dx dt = 2π −∞ −∞ = F −1 (it)α F f (x) (8.12) In these two generalizations the implicit limits of differentiation should be determined. In the case of Laplace transform, the generalized derivative is a Riemann-Liouville derivative with a lower limit of 0, whereas in the case of Fourier transform it is a Weyl derivative. Indeed, if we check for instance the derivative generalized by the Fourier transform in the cases of the sine and cosine functions -calculated in (2.4), (2.5), (2.8) and (2.9) with the generalized derivative of the exponential that we have seen in (7.7) that also is a Weyl derivative- we will find that they match perfectly. 9 Convolution The generalizations of the derivative as expressed in (4.6) and (7.4) suggest that they can be formulated in terms of the convolution, which would be important for the convolution is a very simple operation in the frequency 9 Convolution 15 spaces achieved by Laplace and Fourier transforms. The following development shows how this is the case, and how after all the derivative of a function is its convolution with certain function, xα−1 Φα(x) = Γ(α) (9.1) which Laplace convolution with f (x) yields the Riemann-Liouville derivative of order −α x x (x − t)α−1 Φα(x) ∗ f (x) = Φα(x − t)f (t)dt = f (t)dt = Γ(−α) 0 0 = D−α f (x) (9.2) And since the transform of the function expressed in (9.1) gives the following simple result L xα Γ(α + 1) = tα+1 =⇒ L {Φα} = L xα−1 Γ(α) = t−α (9.3) and the transform of the convolution is the multiplication of the transforms, in the case that the function f fulfills all the requirements given previously for the simplicity of the Laplace transform of the derivative, we again get the equation (8.4), which was the key for the generalization by Laplace transform L D−α f (x) = L Φα(x) ∗ f (x) = L Φα(x) L f (x) = = t−α L f (x) (9.4) These last results showing that the generalized derivative is a convolution with certain function, open the interesting question of what other kind of operators would be defined if functions other than (9.1) had been chosen. The answer is that the new operators defined would be functions of the derivative. To see this, we can exploit the linearity of the convolution, 16 9 Convolution supposing that the function g can be expanded in powers of x ∞ g(x) = an xn =⇒ n=−∞ ∞ ∞ =⇒ g(D)f (x) = an Dn f (x) = an Φ−n(x) ∗ f (x) = n=−∞ n=−∞   ∞ = an x−n−1 ∗ f (x) = Γ(−n) n=−∞   ∞ = a−n−1 xn ∗ f (x) = Γ(n + 1) n=−∞ = h(x) ∗ f (x) (9.5) where ∞ h(x) = bn xn n=−∞ bn = a−n−1 Γ(n + 1) (9.6) This proves the equivalence between functions of the derivative and convolutions. For any convolution there is an equivalent function of the derivative if the function implied can be expanded in powers of x and vice versa. Now, the Laplace transform shows what was expectable, L g(D)f (x) = L h(x) ∗ f (x) = L h(x) L f (x) =   ∞  =L a−n−1 xn  L f (x) = n=−∞ Γ(n + 1)  ∞ = a−n−1 L {xn} L f (x) = Γ(n + 1) n=−∞ ∞ = a−n−1 tn+1 L f (x) = n=−∞ ∞ = an tn L f (x) = n=−∞ = g(t)L f (x) (9.7) since the linearity of the Laplace transform combined with the result of (8.5) 10 Cauchy Integral Formula 17 shows clearly that ∞ ∞ g(D)f (x) = an Dn f (x) = an L−1 tn L f (x) = n=−∞ n=−∞   ∞   = L−1 an tn L f (x) = n=−∞  = L−1 g(t)L f (x) (9.8) as well as (8.11) shows for the Fourier transform that ∞ ∞ g(D)f (x) = an Dn f (x) = an F −1 (it)n F n=−∞ n=−∞   ∞   = F −1 an (it)n F f (x) = n=−∞  f (x) = F −1 g(it)F f (x) = (9.9) These are useful tools for the calculation of functions of the derivative. As an example, the following case is considered with the help of Fourier transforms cos(D) sin(x) = F −1 cos(it)F sin(x) = = F −1 cos(it)πi δ(t + 1) − δ(t − 1) = = cosh(1) sin(x) (9.10) and directly as in (5.2) cos(D) sin(x) = ∞ (−1)n D2n sin(x) = ∞ (−1)n (−1)n sin(x) = (2n)! (2n)! n=0 n=0 ∞1 = sin(x) = cosh(1) sin(x) (2n)! n=0 (9.11) yielding both methods the same result. The case of (5.3) also matches when calculated with Fourier transforms. 10 Cauchy Integral Formula Another way of generalizing the derivative to non-integer order is given by the Cauchy integral formula that plays a key role in complex analysis, Dn f (z) = n! 2πi f (t) (t − z)n+1 dt (10.1) 18 11 Properties Despite its generalization to any complex number α seems straightforward, it must be taken into account that while being n integer there is an isolated singularity at t = z, being it non-integer there is a branch point, what means that the integration contour has to be chosen carefully. Otherwise, the generalization only involves changing the factorial to the gamma function, so defining what is known as the Cauchy-type fractional derivative z0 Dzα f (z) = Γ(α + 1) 2πi f (t) (t − z)α+1 dt C (z0 , z+) (10.2) where supposing that the branch line starts at t = z and passes through z0, the contour C starts at t = z0, encircles t = z once in the positive sense and returns to t = z0 where now the integrand has a different value. It can be proved that this generalization of the derivative is equivalent to the Riemann-Liouville derivative with a lower limit of z0 for the appropriate values of α in which both derivatives are defined. 11 Properties Fractional derivatives satisfy quite well all the properties that one could expect from them, despite some of them are only characteristic of integer order differentiation and some other have restrictions. For instance, the property of linearity (1.10) is fulfilled, while that of the iteration (1.9) has some restrictions in the cases that positive differentiation orders are present. Some properties that include summatories can be generalized changing the summatories into integrals. One such property is the expansion in Taylor series f (z) = ∞ aDzt Γ(t f (z0) + 1) (z − z0)t dt −∞ (11.1) and other is the Leibniz rule ∞ aDzα f (z)g(z) = Γ(t + Γ(α + 1) Γ(α 1) − t + 1) aDzα−t f (z) aDzt g(z) dt −∞ (11.2) These and other generalized properties can be applied to the study of special functions, which often can be expressed in terms of simple formulas involving fractional derivatives. For instance, Gauss’s hypergeometric function can be expressed as 2 F1(a, b; c; x) = Γ(a) x1−c Γ(b) 0Dxb−c xb−1 (1 − x)a (11.3) 12 Local Operators 19 12 Local Operators Despite of some interesting results, we have left unanswered important questions about the functions of the derivative. For instance, it would be very interesting to know if integer order derivatives are or not the only kind of functions of the derivative that are local, and if are or not other kind of local operators that are not functions of the derivative. Looking at (5.2) it seems that in fact the only local functions of the derivative are integer order derivatives and their finite sums, being nonlocality the result of adding infinite terms with displacement functions steadily increasing the distance to the point in which the function is evaluated. However, this impediment can be overcome by defining the derivative in the following way Dα f (x) = lim a→x aDxα f (x) (12.1) so that the limit assures locality. In order to avoid another feature of fractional derivatives that sometimes is not wanted -the fact that the derivative of a constant is not zero- it is usual to define the local fractional derivative as Dα f (x) = lim y→x xDyα f (y) − f (x) (12.2) or similar forms. So there can be operators being local and taking an infinite number of evaluations of the function, provided that the points in which these evaluations are carried out -or the value of their displacement functionremains being infinitesimal. The question now is if there are other local operators besides the local fractional derivative and functions of it. It could be argued that if a function admits its expansion in Taylor series, all the information about the function is indeed in its derivatives, so that any operator giving information about the function will be unnecessary and reducible to a sum of integer order derivatives. However, on the one hand there are functions that are not differentiable but are local fractional differentiable, and on the other hand, even if local operators could be reduced to ordinary derivatives for differentiable functions, that does not assure that ordinary derivatives are the most appropriate local operators for all the tasks. Characterizing all local operators is important, for locality makes of them tools that could be useful in understanding Nature.

Top_arrow
回到顶部
EEWORLD下载中心所有资源均来自网友分享,如有侵权,请发送举报邮件到客服邮箱bbs_service@eeworld.com.cn 或通过站内短信息或QQ:273568022联系管理员 高进,我们会尽快处理。