Riemann–Liouville integral

Template:Short description {{#invoke:sidebar|collapsible | class = plainlist | titlestyle = padding-bottom:0.25em; | pretitle = Part of a series of articles about | title = Calculus | image = <math>\int_{a}^{b} f'(t) \, dt = f(b) - f(a)</math> | listtitlestyle = text-align:center; | liststyle = border-top:1px solid #aaa;padding-top:0.15em;border-bottom:1px solid #aaa; | expanded = Fractional calculus | abovestyle = padding:0.15em 0.25em 0.3em;font-weight:normal; | above =

| list2name = differential | list2titlestyle = display:block;margin-top:0.65em; | list2title = Template:Bigger | list2 ={{#invoke:sidebar|sidebar|child=yes

 |contentclass=hlist
 | heading1 = Definitions
 | content1 =

 | heading2 = Concepts
 | content2 =

 | heading3 = Rules and identities
 | content3 =

}}

| list3name = integral | list3title = Template:Bigger | list3 ={{#invoke:sidebar|sidebar|child=yes

 |contentclass=hlist
 | content1 =

| heading2 = Definitions

 | content2 =

 | heading3 = Integration by
 | content3 =

}}

| list4name = series | list4title = Template:Bigger | list4 ={{#invoke:sidebar|sidebar|child=yes

 |contentclass=hlist
 | content1 =

 | heading2 = Convergence tests
 | content2 =

}}

| list5name = vector | list5title = Template:Bigger | list5 ={{#invoke:sidebar|sidebar|child=yes

 |contentclass=hlist
 | content1 =

 | heading2 = Theorems
 | content2 =

}}

| list6name = multivariable | list6title = Template:Bigger | list6 ={{#invoke:sidebar|sidebar|child=yes

 |contentclass=hlist
 | heading1 = Formalisms
 | content1 =

 | heading2 = Definitions
 | content2 =

}}

| list7name = advanced | list7title = Template:Bigger | list7 ={{#invoke:sidebar|sidebar|child=yes

 |contentclass=hlist
 | content1 =

}}

| list8name = specialized | list8title = Template:Bigger | list8 =

| list9name = miscellanea | list9title = Template:Bigger | list9 =

}}

In mathematics, the Riemann–Liouville integral associates with a real function <math>f: \mathbb{R} \rightarrow \mathbb{R}</math> another function Template:Math of the same kind for each value of the parameter Template:Math. The integral is a manner of generalization of the repeated antiderivative of Template:Mvar in the sense that for positive integer values of Template:Mvar, Template:Math is an iterated antiderivative of Template:Mvar of order Template:Mvar. The Riemann–Liouville integral is named for Bernhard Riemann and Joseph Liouville, the latter of whom was the first to consider the possibility of fractional calculus in 1832.<ref name="Lizorkin">Template:Harvnb</ref><ref>Template:Citation.</ref><ref>Template:Citation.</ref><ref>Template:Citation.</ref> The operator agrees with the Euler transform, after Leonhard Euler, when applied to analytic functions.<ref>Template:Harvnb</ref> It was generalized to arbitrary dimensions by Marcel Riesz, who introduced the Riesz potential.

Motivation

The Riemann-Liouville integral is motivated from Cauchy formula for repeated integration. For a function Template:Mvar continuous on the interval [[[:Template:Mvar]],Template:Mvar], the Cauchy formula for Template:Mvar-fold repeated integration states that

<math display="block">I^n f(x) = f^{(-n)}(x) = \frac{1}{(n-1)!} \int_a^x\left(x-t\right)^{n-1} f(t)\,\mathrm{d}t.</math>

Now, this formula can be generalized to any positive real number by replacing positive integer Template:Mvar with Template:Mvar, Therefore we obtain the definition of Riemann-Liouville fractional Integral by

<math>I^\alpha f(x) = \frac{1}{\Gamma(\alpha)} \int_a^x\left(x-t\right)^{\alpha-1} f(t)\,\mathrm{d}t</math>

Definition

The Riemann–Liouville integral is defined by

<math>I^\alpha f(x) = \frac{1}{\Gamma(\alpha)}\int_a^xf(t)(x-t)^{\alpha-1}\,dt</math>

where Template:Math is the gamma function and Template:Mvar is an arbitrary but fixed base point. The integral is well-defined provided Template:Mvar is a locally integrable function, and Template:Mvar is a complex number in the half-plane Template:Math. The dependence on the base-point Template:Mvar is often suppressed, and represents a freedom in constant of integration. Clearly Template:Math is an antiderivative of Template:Mvar (of first order), and for positive integer values of Template:Mvar, Template:Math is an antiderivative of order Template:Mvar by Cauchy formula for repeated integration. Another notation, which emphasizes the base point, is<ref>Template:Harvnb</ref>

<math>{}_aD_x^{-\alpha}f(x) = \frac{1}{\Gamma(\alpha)}\int_a^x f(t)(x-t)^{\alpha-1}\,dt.</math>

This also makes sense if Template:Math, with suitable restrictions on Template:Mvar.

The fundamental relations hold

<math>\frac{d}{dx}I^{\alpha+1} f(x) = I^\alpha f(x),\quad I^\alpha(I^\beta f) = I^{\alpha+\beta}f,</math>

the latter of which is a semigroup property.<ref name="Lizorkin" /> These properties make possible not only the definition of fractional integration, but also of fractional differentiation, by taking enough derivatives of Template:Math.

Properties

Fix a bounded interval Template:Open-open. The operator Template:Math associates to each integrable function Template:Mvar on Template:Open-open the function Template:Math on Template:Open-open which is also integrable by Fubini's theorem. Thus Template:Math defines a linear operator on [[Lp space|Template:Math]]:

<math>I^\alpha : L^1(a,b) \to L^1(a,b).</math>

Fubini's theorem also shows that this operator is continuous with respect to the Banach space structure on Template:Mvar¹, and that the following inequality holds:

<math>\left \|I^\alpha f \right \|_1 \le \frac{|b-a|^{\Re(\alpha)}}{\Re(\alpha)|\Gamma(\alpha)|}\|f\|_1.</math>

Here Template:Math denotes the norm on Template:Math.

More generally, by Hölder's inequality, it follows that if Template:Math, then Template:Math as well, and the analogous inequality holds

<math>\left \|I^\alpha f \right \|_p \le \frac{|b-a|^{\Re(\alpha)/p}}{\Re(\alpha)|\Gamma(\alpha)|}\|f\|_p</math>

where Template:Math is the [[Lp space|Template:Mvar^{Template:Mvar} norm]] on the interval Template:Open-open. Thus we have a bounded linear operator Template:Math. Furthermore, Template:Math in the Template:Math sense as Template:Math along the real axis. That is

<math>\lim_{\alpha\to 0^+} \|I^\alpha f - f\|_p = 0</math>

for all Template:Math. Moreover, by estimating the maximal function of Template:Mvar, one can show that the limit Template:Math holds pointwise almost everywhere.

The operator Template:Math is well-defined on the set of locally integrable function on the whole real line <math>\mathbb{R}</math>. It defines a bounded transformation on any of the Banach spaces of functions of exponential type <math>X_{\sigma} = L^1(e^{-\sigma|t|}dt),</math> consisting of locally integrable functions for which the norm

<math>\|f\| = \int_{-\infty}^\infty |f(t)|e^{-\sigma|t|}\,dt</math>

is finite. For Template:Math, the Laplace transform of Template:Math takes the particularly simple form

<math>(\mathcal{L}I^\alpha f)(s) = s^{-\alpha}F(s)</math>

for Template:Math. Here Template:Math denotes the Laplace transform of Template:Mvar, and this property expresses that Template:Math is a Fourier multiplier.

Fractional derivatives

One can define fractional-order derivatives of Template:Mvar as well by

<math>\frac{d^\alpha}{dx^\alpha} f \, \overset{\text{def}}{=} \frac{d^{\lceil\alpha\rceil}}{dx^{\lceil\alpha\rceil}} I^{\lceil \alpha \rceil-\alpha}f</math>

where Template:Math denotes the ceiling function. One also obtains a differintegral interpolating between differentiation and integration by defining

<math>D^\alpha_x f(x) = \begin{cases} \frac{d^{\lceil\alpha\rceil}}{dx^{\lceil\alpha\rceil}} I^{\lceil\alpha\rceil-\alpha}f(x)& \alpha>0\\ f(x) & \alpha=0\\ I^{-\alpha}f(x) & \alpha<0. \end{cases}</math>

An alternative fractional derivative was introduced by Caputo in 1967,<ref>Template:Harvnb</ref> and produces a derivative that has different properties: it produces zero from constant functions and, more importantly, the initial value terms of the Laplace Transform are expressed by means of the values of that function and of its derivative of integer order rather than the derivatives of fractional order as in the Riemann–Liouville derivative.<ref>Template:Harvnb</ref> The Caputo fractional derivative with base point Template:Mvar, is then:

<math>D_x^{\alpha}f(y)=\frac{1}{\Gamma(1-\alpha)}\int_x^y f'(y-u)(u-x)^{-\alpha}du.</math>

Another representation is:

<math>{}_a\tilde{D}^\alpha_x f(x)=I^{\lceil \alpha\rceil-\alpha}\left(\frac{d^{\lceil \alpha\rceil}f}{dx^{\lceil \alpha\rceil}}\right).</math>

Fractional derivative of a basic power function

The half derivative (purple curve) of the function Template:Math (blue curve) together with the first derivative (red curve).

Let us assume that Template:Math is a monomial of the form

The first derivative is as usual

Repeating this gives the more general result that

<math>\frac{d^a}{dx^a}x^k = \dfrac{k!}{(k-a)!}x^{k-a}\,,</math>

which, after replacing the factorials with the gamma function, leads to

<math>\frac{d^a}{dx^a}x^k = \dfrac{\Gamma(k+1)}{\Gamma(k-a+1)}x^{k-a}, \quad\ k > 0.</math>

For Template:Math and Template:Math, we obtain the half-derivative of the function <math>x \mapsto x</math> as

<math>\frac{d^\frac12}{dx^\frac12}x=\frac{\Gamma(1+1)}{\Gamma\left(1-\frac12 + 1\right)} x^{1-\frac12}=\frac{\Gamma(2)}{\Gamma\left(\frac{3}{2}\right)}x^\frac12 = \frac{1}{\frac{\sqrt{\pi}}{2}}x^\frac12.</math>

To demonstrate that this is, in fact, the "half derivative" (where Template:Math), we repeat the process to get:

<math>\dfrac{d^\frac12}{dx^\frac12} \dfrac{2x^{\frac12}}{\sqrt{\pi}}

=\frac{2}{\sqrt{\pi}}\dfrac{\Gamma(1+\frac12)}{\Gamma(\frac12-\frac12+1)}x^{\frac12-\frac12} =\frac{2}{\sqrt{\pi}} \frac{\Gamma\left(\frac{3}{2}\right)}{\Gamma(1)} x^0 =\frac{2 \frac{\sqrt{\pi}}{2} x^0}{\sqrt{\pi}}=1\,,</math>

(because <math display="inline">\Gamma\!\left(\frac{3}{2}\right) = \frac{\sqrt{\pi}}{2}</math> and Template:Math) which is indeed the expected result of

<math>\left(\frac{d^\frac12}{dx^\frac12} \frac{d^\frac12}{dx^\frac12}\right)\!x = \frac{d}{dx} x = 1\,.</math>

For negative integer power Template:Math, 1/<math display="inline">\Gamma</math> is 0, so it is convenient to use the following relation:<ref>Template:Citation</ref>

<math>\frac{d^a}{dx^a}x^{-k} = \left(-1\right)^a\dfrac{\Gamma(k+a)}{\Gamma(k)}x^{-(k+a)} \quad\text{ for } k \ge 0.</math>

This extension of the above differential operator need not be constrained only to real powers; it also applies for complex powers. For example, the Template:Math-th derivative of the Template:Math-th derivative yields the second derivative. Also setting negative values for Template:Mvar yields integrals.

For a general function Template:Math and Template:Math, the complete fractional derivative is

<math>D^\alpha f(x)=\frac{1}{\Gamma(1-\alpha)}\frac{d}{dx} \int_0^x \frac{f(t)}{\left(x-t\right)^\alpha} \, dt.</math>

For arbitrary Template:Mvar, since the gamma function is infinite for negative (real) integers, it is necessary to apply the fractional derivative after the integer derivative has been performed. For example,

Laplace transform

{{ safesubst:#invoke:Unsubst||date=__DATE__ |$B= {{ safesubst:#invoke:Unsubst||date=__DATE__ |$B= Template:Ambox }} }} We can also come at the question via the Laplace transform. Knowing that

<math>\mathcal L \left\{Jf\right\}(s) = \mathcal L \left\{\int_0^t f(\tau)\,d\tau\right\}(s) = \frac1 s \bigl(\mathcal L\left\{f\right\}\bigr)(s)</math>

and

<math>\mathcal L \left\{J^2f\right\}=\frac1s\bigl(\mathcal L \left\{Jf\right\} \bigr)(s)=\frac1{s^2}\bigl(\mathcal L\left\{f\right\}\bigr)(s)</math>

and so on, we assert

<math>J^\alpha f=\mathcal L^{-1}\left\{s^{-\alpha}\bigl(\mathcal L\{f\}\bigr)(s)\right\}</math>.

For example,

<math>J^\alpha(t^k) = \mathcal L^{-1}\left\{\frac{\Gamma(k+1)}{s^{\alpha+k+1}}\right\} = \frac{\Gamma(k+1)}{\Gamma(\alpha+k+1)} t^{\alpha+k} </math>

as expected. Indeed, given the convolution rule

<math>\mathcal L\{f*g\}=\bigl(\mathcal L\{f\}\bigr)\bigl(\mathcal L\{g\}\bigr)</math>

and shorthanding Template:Math for clarity, we find that

<math>\begin{align}

\left(J^\alpha f\right)(t) &= \frac{1}{\Gamma(\alpha)}\mathcal L^{-1}\left\{\bigl(\mathcal L\{p\}\bigr)\bigl(\mathcal L\{f\}\bigr)\right\}\\ &=\frac{1}{\Gamma(\alpha)}(p*f)\\ &=\frac{1}{\Gamma(\alpha)}\int_0^t p(t-\tau)f(\tau)\,d\tau\\ &=\frac{1}{\Gamma(\alpha)}\int_0^t\left(t-\tau\right)^{\alpha-1}f(\tau)\,d\tau\\ \end{align}</math>

which is what Cauchy gave us above.

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

Notes

References

Template:Sfn whitelist

External links

{{#invoke:citation/CS1|citation

|CitationClass=web }}

{{#invoke:citation/CS1|citation

|CitationClass=web }}

Template:Bernhard Riemann

Riemann–Liouville integral

Contents

Motivation

Definition

Properties

Fractional derivatives

Fractional derivative of a basic power function

Laplace transform

See also

Notes

References

External links

Navigation menu

Riemann–Liouville integral

Motivation

Definition

Properties

Fractional derivatives

Fractional derivative of a basic power function

Laplace transform

See also

Notes

References

External links

Navigation menu

Search