Fractional calculus

From Wikipedia, the free encyclopedia

In mathematics, fractional calculus is a branch of mathematical analysis that studies the possibility of taking real number powers of the differential operator

$D = \frac{d}{dx} \,$

and the integration operator J. (Usually not I, to avoid confusion with other I-like glyphs, or identities, but J must not be confused with Bessel functions, which often come up in study of differential equations.)

In this context powers refer to iterative application, in the same sense that f²(x) = f(f(x)).
For example, one may pose the question of interpreting meaningfully

$\sqrt{D} = D^{1/2} \,$

as a square root of the differentiation operator (an operator half iterate), i.e., an expression for some operator that when applied twice to a function will have the same effect as differentiation. More generally, one can look at the question of defining

$D^s \,$

for real number values of s in such a way that when s takes an integer value n, the usual power of n-fold differentiation is recovered for n > 0, and the −nth power of J when n < 0.

There are various reasons for looking at this question. One is that in this way the semigroup of powers Dⁿ in the discrete variable n is seen inside a continuous semigroup (one hopes) with parameter s which is a real number. Continuous semigroups are prevalent in mathematics, and have an interesting theory. Notice here that fraction is then a misnomer for the exponent, since it need not be rational, but the term fractional calculus has become traditional.

1 Fractional derivative
2 Heuristics
3 Half derivative of a simple function
4 Laplace transform
5 Riemann-Liouville differintegral
6 Functional calculus
7 References
8 See also
9 External links

[edit] Fractional derivative

As far as the existence of such a theory is concerned, the foundations of the subject were laid by Liouville in a paper from 1832. The fractional derivative of a function to order a is often now defined by means of the Fourier or Mellin integral transforms. An important point is that the fractional derivative at a point x is a local property only when a is an integer; in non-integral cases we cannot say that the fractional derivative at x of a function f depends only on the graph of f very near x, in the way that integer-power derivatives certainly do. Therefore it is expected that the theory involves some sort of boundary conditions, involving information on the function further out. To use a metaphor, the fractional derivative requires some peripheral vision.

For the history of the subject, see the thesis (in French): Stéphane Dugowson, Les différentielles métaphysiques (histoire et philosophie de la généralisation de l'ordre de dérivation), Thèse, Université Paris Nord (1994)

[edit] Heuristics

A fairly natural question to ask is, does there exist an operator $H$ , or half-derivative, such that

$H^2 f(x) = D f(x) = \frac{d}{dx} f(x) = f'(x)$ ?

It turns out that there is such an operator, and indeed for any $a > 0$ , there exists an operator $P$ such that

$(P ^ a f)(x) = f'(x) \,$ ,

or to put it another way, $\frac{d^ny}{dx^n}$ is well-defined for all real values of n > 0. A similar result applies to integration.

To delve into a little detail, start with the Gamma function $\Gamma \,$ , which is defined such that

$n! = \Gamma(n+1) \,$ .

Assuming a function $f (x)$ that is well defined where $x > 0$ , we can form the definite integral from 0 to x. Let's call this

$( J f ) ( x ) = \int_0^x f(t) \; dt$ .

Repeating this process gives

$( J^2 f ) ( x ) = \int_0^x ( J f ) ( t ) dt = \int_0^x \left( \int_0^t f(s) \; ds \right) \; dt$ ,

and this can be extended arbitrarily.

The Cauchy formula for repeated integration, namely

$(J^n f) ( x ) = { 1 \over (n-1) ! } \int_0^x (x-t)^{n-1} f(t) \; dt,$

leads to a straightforward way to a generalization for real n.

Simply using the Gamma function to remove the discrete nature of the factorial function gives us a natural candidate for fractional applications of the integral operator.

$(J^\alpha f) ( x ) = { 1 \over \Gamma ( \alpha ) } \int_0^x (x-t)^{\alpha-1} f(t) \; dt$

This is in fact a well-defined operator.

It can be shown that the J operator is both commutative and additive. That is,

$(J^\alpha) (J^\beta) f = (J^\beta) (J^\alpha) f = (J^{\alpha+\beta} ) f = { 1 \over \Gamma ( \alpha + \beta) } \int_0^x (x-t)^{\alpha+\beta-1} f(t) \; dt$

This property is called the Semi-Group property of fractional differintegral operators. Unfortunately the comparable process for the derivative operator D is significantly more complex, but it can be shown that D is neither commutative, nor additive in general.

[edit] Half derivative of a simple function

The half derivative (maroon curve) of the function y=x.

Let us assume that $f (x)$ is a monomial of the form

$f(x) = x^k\;.$

The first derivative is as usual

$f'(x) = {d \over dx } f(x) = k x^{k-1}\;.$

Repeating this gives the more general result that

${d^a \over dx^a } x^k = { k! \over (k - a) ! } x^{k-a}\;,$

Which, after replacing the factorials with the Gamma function, leads us to

${d^a \over dx^a } x^k = { \Gamma(k+1) \over \Gamma(k - a + 1) } x^{k-a}\;.$

So, for example, the half-derivative of $x$ is

${ d^{1 \over 2} \over dx^{1 \over 2} } x = { \Gamma(1 + 1) \over \Gamma ( 1 - {1 \over 2} + 1 ) } x^{1-{1 \over 2}} = { \Gamma( 2 ) \over \Gamma ( { 3 \over 2 } ) } x^{1 \over 2} = {2 \pi^{-{1 \over 2}}} x^{1 \over 2}\;.$

Repeating this process gives

${ d^{1 \over 2} \over dx^{1 \over 2} } {2 \pi^{-{1 \over 2}}} x^{1 \over 2} = {2 \pi^{-{1 \over 2}}} { \Gamma ( 1 + {1 \over 2} ) \over \Gamma ( {1 \over 2} - { 1 \over 2 } + 1 ) } x^{{1 \over 2} - {1 \over 2}} = {2 \pi^{-{1 \over 2}}} { \Gamma( { 3 \over 2 } ) \over \Gamma ( 1 ) } x^0 = { 1 \over \Gamma (1) } = 1\;,$

which is indeed the expected result of

$\left( \frac{d^{1/2}}{dx^{1/2}} \frac{d^{1/2}}{dx^{1/2}} \right) x = { d \over dx } x = 1 \,$

[edit] Laplace transform

One can also come at the question via the Laplace transform. Noting that

$\mathcal L\left(t\mapsto\int_0^t f(\tau)\,d\tau\right)=\mathcal LJf=s\mapsto\frac1s(\mathcal Lf)(s)$

and

$\mathcal LJ^2f=s\mapsto\frac1s(\mathcal LJf)(s)=s\mapsto\frac1{s^2}(\mathcal Lf)(s)$

etc., we assert

$J^\alpha f=\mathcal L^{-1}\left(s\mapsto s^{-\alpha}(\mathcal Lf)(s)\right)$ .

For example

$J^\alpha\left(t\mapsto t^k\right)$ $=\mathcal L^{-1}\left(s\mapsto{\Gamma(k+1)\over s^{\alpha+k+1}}\right)$ $=t\mapsto{\Gamma(k+1)\over\Gamma(\alpha+k+1)}t^{\alpha+k}$

as expected. Indeed, given the convolution rule $\mathcal L(f*g)=(\mathcal Lf)(\mathcal Lg)$ (and shorthanding $p (x) = x α - 1$ for clarity) we find that

$J^\alpha f=\frac1{\Gamma(\alpha)}\mathcal L^{-1}\left(\left(\mathcal Lp\right)(\mathcal Lf)\right)$ $=\frac1{\Gamma(\alpha)}(p*f)$ $=x\mapsto\frac1{\Gamma(\alpha)}\int_0^xp(x-t)f(t)\,dt$ $=x\mapsto\frac1{\Gamma(\alpha)}\int_0^x(x-t)^{\alpha-1}f(t)\,dt$

which is what Cauchy gave us above.

Laplace transforms "work" on relatively few functions, but they are often useful for solving fractional differential equations.

[edit] Riemann-Liouville differintegral

The classical form of fractional calculus is given by the Riemann-Liouville differintegral, essentially what has been described above. The theory for periodic functions, therefore including the 'boundary condition' of repeating after a period, is the Weyl differintegral. It is defined on Fourier series, and requires the constant Fourier coefficient to vanish (so, applies to functions on the unit circle integrating to 0).

By contrast the Grunwald-Letnikov differintegral starts with the derivative.

[edit] Functional calculus

In the context of functional analysis, functions f(D) more general than powers are studied in the functional calculus of spectral theory. The theory of pseudo-differential operators also allows one to consider powers of D. The operators arising are examples of singular integral operators; and the generalisation of the classical theory to higher dimensions is called the theory of Riesz potentials. So there are a number of contemporary theories available, within which fractional calculus can be discussed. See also Erdelyi-Kober operator, important in special function theory.

For possible geometric and physical interpretation of fractional-order integration and fractional-order differentiation, see:

Podlubny, I., Geometric and physical interpretation of fractional integration and fractional differentiation. Fractional Calculus and Applied Analysis, vol. 5, no. 4, 2002, 367–386. (available as original article, or preprint at Arxiv.org)

[edit] References

Theory and Applications of Fractional Differential Equations, by Kilbas, A. A.; Srivastava, H. M.; and Trujillo, J. J. Amsterdam, Netherlands, Elsevier, Febrary 2006. ISBN 0-444-51832-0 (http://www.elsevier.com/wps/find/bookdescription.cws_home/707212/description#description)
An Introduction to the Fractional Calculus and Fractional Differential Equations, by Kenneth S. Miller, Bertram Ross (Editor). Hardcover: 384 pages. Publisher: John Wiley & Sons; 1 edition (May 19, 1993). ISBN 0-471-58884-9
The Fractional Calculus; Theory and Applications of Differentiation and Integration to Arbitrary Order (Mathematics in Science and Engineering, V), by Keith B. Oldham, Jerome Spanier. Hardcover. Publisher: Academic Press; (November 1974). ISBN 0-12-525550-0
Fractional Differential Equations. An Introduction to Fractional Derivatives, Fractional Differential Equations, Some Methods of Their Solution and Some of Their Applications., (Mathematics in Science and Engineering, vol. 198), by Igor Podlubny. Hardcover. Publisher: Academic Press; (October 1998) ISBN 0-12-558840-2
Fractals and Fractional Calculus in Continuum Mechanics, by A. Carpinteri (Editor), F. Mainardi (Editor). Paperback: 348 pages. Publisher: Springer-Verlag Telos; (January 1998). ISBN 3-211-82913-X
Physics of Fractal Operators, by Bruce J. West, Mauro Bologna, Paolo Grigolini. Hardcover: 368 pages. Publisher: Springer Verlag; (January 14, 2003). ISBN 0-387-95554-2
Fractional Calculus and the Taylor-Riemann Series, Rose-Hulman Undergrad. J. Math. Vol.6(1) (2005).