Intégrale de Riemann-Liouville

L'intégrale de Riemann-Liouville est définie par

${\displaystyle \mathrm {I} ^{\alpha }f(x)={\frac {1}{\Gamma (\alpha )}}\int _{a}^{x}f(t)(x-t)^{\alpha -1}\,\mathrm {d} t}$

où Γ est la fonction gamma et a est un point de base arbitraire mais fixe. L'intégrale est bien définie pourvu que f soit une fonction localement intégrable, et α soit un nombre complexe dans le demi-plan Re(α) > 0 . La dépendance au point de base a est souvent supprimée, et représente une liberté en constante d'intégration. Il est clair que I¹f est une primitive de f (du premier ordre), et pour des valeurs entières positives de α, I^αf est une primitive d'ordre α par la formule de Cauchy pour l'intégration successive. Une autre notation, qui met l'accent sur le point de base, est [3]:

${\displaystyle {}_{a}\mathrm {D} _{x}^{-\alpha }f(x)={\frac {1}{\Gamma (\alpha )}}\int _{a}^{x}f(t)(x-t)^{\alpha -1}\,\mathrm {d} t.}$

Cela a également un sens si a = −∞, avec des restrictions appropriées sur f .

Les relations fondamentales sont vérifiées

${\displaystyle {\frac {\mathrm {d} }{\mathrm {d} x}}\mathrm {I} ^{\alpha +1}f(x)=\mathrm {I} ^{\alpha }f(x),\quad \mathrm {I} ^{\alpha }(\mathrm {I} ^{\beta }f)=\mathrm {I} ^{\alpha +\beta }f,}$

ce dernier est une propriété de semi-groupe[1]. Ces propriétés permettent non seulement de définir l'intégration fractionnaire, mais aussi la dérivation fractionnaire, en prenant suffisamment de dérivées de I^α f.

On fixe un intervalle borné (a , b). L'opérateur I^α associe à chaque fonction intégrable f sur (a , b) la fonction I^α f sur (a , b) qui est également intégrable par le théorème de Fubini. Ainsi I^α définit un opérateur linéaire sur L¹(a,b) :

${\displaystyle \mathrm {I} ^{\alpha }:L^{1}(a,b)\to L^{1}(a,b).}$

Le théorème de Fubini montre également que cet opérateur est continu par rapport à la structure de l'espace de Banach sur L¹, et que l'inégalité suivante est vérifiée :

${\displaystyle \left\|\mathrm {I} ^{\alpha }f\right\|_{1}\leq {\frac {|b-a|^{\Re (\alpha )}}{\Re (\alpha )|\Gamma (\alpha )|}}\|f\|_{1}.}$

Ici ‖ · ‖₁ désigne la norme sur L¹(a,b) .

Plus généralement, d'après l'inégalité de Hölder, il s'ensuit que si f ∈ L^p(a, b), alors I^α f ∈ L^p(a, b) également, et l'inégalité analogue est vraie :

${\displaystyle \left\|\mathrm {I} ^{\alpha }f\right\|_{p}\leq {\frac {|b-a|^{\Re (\alpha )/p}}{\Re (\alpha )|\Gamma (\alpha )|}}\|f\|_{p}}$

où ‖ · ‖_p est la norme L ^p sur l'intervalle (a, b ) . On a donc un opérateur linéaire borné I^α : L^p(a, b) → L^p(a, b). De plus, I^α f → f au sens L^p comme α → 0 le long de l'axe réel. C'est-à-dire

${\displaystyle \lim _{\alpha \to 0^{+}}\|\mathrm {I} ^{\alpha }f-f\|_{p}=0}$

pour tout p ≥ 1 . De plus, en estimant la fonction maximale (en) de I, on peut montrer que la limite I^α f → f est vraie ponctuellement presque partout.

L'opérateur I^α est bien défini sur l'ensemble des fonctions localement intégrables sur toute la droite réelle ${\displaystyle \mathbb {R} }$. Il définit une transformation bornée sur l'un des espaces de Banach des fonctions de type exponentielle ${\displaystyle X_{\sigma }=L^{1}(\mathrm {e} ^{-\sigma |t|}\mathrm {d} t),}$ constitué de fonctions localement intégrables dont la norme

${\displaystyle \|f\|=\int _{-\infty }^{\infty }|f(t)|\mathrm {e} ^{-\sigma |t|}\,\mathrm {d} t}$

est fini. Pour f ∈ X_σ, la transformée de Laplace de I^α f prend la forme particulièrement simple

${\displaystyle ({\mathcal {L}}\mathrm {I} ^{\alpha }f)(s)=s^{-\alpha }F(s)}$

pour Re(s) > σ . Ici F(s) désigne la transformée de Laplace de f, et cette propriété exprime que I^α est un multiplicateur de Fourier.

On peut également définir les dérivées d'ordre fractionnaire de f par

${\displaystyle {\frac {\mathrm {d} ^{\alpha }}{\mathrm {d} x^{\alpha }}}f\ {\overset {\text{def}}{=}}\ {\frac {\mathrm {d} ^{\lceil \alpha \rceil }}{\mathrm {d} x^{\lceil \alpha \rceil }}}\mathrm {I} ^{\lceil \alpha \rceil -\alpha }f}$

où ⌈ · ⌉ désigne la fonction partie entière par excès. On obtient également une différence intégrale interpolant entre différenciation et intégration en définissant

${\displaystyle \mathrm {D} _{x}^{\alpha }f(x)={\begin{cases}{\frac {\mathrm {d} ^{\lceil \alpha \rceil }}{\mathrm {d} x^{\lceil \alpha \rceil }}}\mathrm {I} ^{\lceil \alpha \rceil -\alpha }f(x)&\ {\textrm {pour}}\ \alpha >0\\f(x)&\ {\textrm {pour}}\ \alpha =0\\\mathrm {I} ^{-\alpha }f(x)&\ {\textrm {pour}}\ \alpha <0.\end{cases}}}$

Une dérivée fractionnaire alternative a été introduite par Caputo 1967, et produit une dérivée qui a des propriétés différentes : elle produit zéro à partir de fonctions constantes et, plus important encore, les termes de valeur initiale de la transformée de Laplace sont exprimés au moyen des valeurs de cette fonction et de sa dérivée d'ordre entier plutôt que les dérivées d'ordre fractionnaire comme dans la dérivée de Riemann-Liouville (Loverro 2004). La dérivée fractionnaire de Caputo de point de base x, est alors :

${\displaystyle \mathrm {D} _{x}^{\alpha }f(y)={\frac {1}{\Gamma (1-\alpha )}}\int _{x}^{y}f'(y-u)(u-x)^{-\alpha }\mathrm {d} u.}$

Une autre représentation est :

${\displaystyle {}_{a}{\widetilde {\mathrm {D} }}_{x}^{\alpha }f(x)=\mathrm {I} ^{\lceil \alpha \rceil -\alpha }\left({\frac {\mathrm {d} ^{\lceil \alpha \rceil }f}{\mathrm {d} x^{\lceil \alpha \rceil }}}\right).}$

La demi-dérivée (courbe violette) de la fonction f(x) = x (courbe bleue) avec la première dérivée (courbe rouge).

L'animation montre l'opérateur dérivé oscillant entre la primitive ( α = −1 : y = 1/2x²

Supposons que f(x) soit un monôme de la forme

${\displaystyle f(x)=x^{k}\,.}$

La dérivée première est comme d'habitude

${\displaystyle f'(x)={\frac {\mathrm {d} }{\mathrm {d} x}}f(x)=kx^{k-1}\,.}$

Répéter l'opération donne le résultat plus général :

${\displaystyle {\frac {\mathrm {d} ^{a}}{\mathrm {d} x^{a}}}x^{k}={\dfrac {k!}{(k-a)!}}x^{k-a}\,,}$

ce qui, après remplacement des factorielles par la fonction gamma, conduit à

${\displaystyle {\frac {\mathrm {d} ^{a}}{\mathrm {d} x^{a}}}x^{k}={\dfrac {\Gamma (k+1)}{\Gamma (k-a+1)}}x^{k-a},\quad \ \forall k>0.}$

Pour k = 1 et a = 1/2, on obtient la demi-dérivée de la fonction ${\displaystyle x\mapsto x}$ par :

${\displaystyle {\frac {\mathrm {d} ^{\frac {1}{2}}}{\mathrm {d} x^{\frac {1}{2}}}}x={\frac {\Gamma (1+1)}{\Gamma \left(1-{\frac {1}{2}}+1\right)}}x^{1-{\frac {1}{2}}}={\frac {\Gamma (2)}{\Gamma \left({\frac {3}{2}}\right)}}x^{\frac {1}{2}}={\frac {1}{\frac {\sqrt {\pi }}{2}}}x^{\frac {1}{2}}.}$

Pour démontrer qu'il s'agit en fait de la "demi-dérivée" (où H²f(x) = Df(x) ), on répète le processus pour obtenir :

${\displaystyle {\dfrac {\mathrm {d} ^{\frac {1}{2}}}{\mathrm {d} x^{\frac {1}{2}}}}{\dfrac {2x^{\frac {1}{2}}}{\sqrt {\pi }}}={\frac {2}{\sqrt {\pi }}}{\dfrac {\Gamma (1+{\frac {1}{2}})}{\Gamma ({\frac {1}{2}}-{\frac {1}{2}}+1)}}x^{{\frac {1}{2}}-{\frac {1}{2}}}={\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\,,}$

(car ${\textstyle \Gamma \!\left({\frac {3}{2}}\right)={\frac {\sqrt {\pi }}{2}}}$ et Γ(1) = 1 ) qui est bien le résultat attendu de

${\displaystyle \left({\frac {\mathrm {d} ^{\frac {1}{2}}}{\mathrm {d} x^{\frac {1}{2}}}}{\frac {\mathrm {d} ^{\frac {1}{2}}}{\mathrm {d} x^{\frac {1}{2}}}}\right)\!x={\frac {\mathrm {d} }{\mathrm {d} x}}x=1\,.}$

Pour une puissance entière négative k, ${\textstyle 1/\Gamma }$ est nulle, il est donc intéressant d'utiliser la relation suivante :

${\displaystyle {\frac {\mathrm {d} ^{a}}{\mathrm {d} x^{a}}}x^{-k}=\left(-1\right)^{a}{\dfrac {\Gamma (k+a)}{\Gamma (k)}}x^{-(k+a)}\quad {\text{ pour }}k\geq 0.}$

Cette extension de l'opérateur différentiel ci-dessus n'a pas besoin d'être limitée uniquement aux puissances réelles, et peut s'appliquer également aux pouvoirs complexes. Par exemple, la (1 + i) -ième dérivée de la (1 − i) -ième dérivée donne la seconde dérivée. Définir également des valeurs négatives pour a donne des intégrales.

Pour une fonction générale f(x) et 0 < α < 1, la dérivée fractionnaire complète est

${\displaystyle \mathrm {D} ^{\alpha }f(x)={\frac {1}{\Gamma (1-\alpha )}}{\frac {\mathrm {d} }{\mathrm {d} x}}\int _{0}^{x}{\frac {f(t)}{\left(x-t\right)^{\alpha }}}\,\mathrm {d} t.}$

Pour α arbitraire, puisque la fonction gamma est infinie pour les entiers négatifs, il est nécessaire d'appliquer la dérivée fractionnaire après que la dérivée entière a été effectuée. Par exemple,

${\displaystyle \mathrm {D} ^{\frac {3}{2}}f(x)=\mathrm {D} ^{\frac {1}{2}}\mathrm {D} ^{1}f(x)=\mathrm {D} ^{\frac {1}{2}}{\frac {\mathrm {d} }{\mathrm {d} x}}f(x).}$

On peut aussi venir à la question via la transformée de Laplace. Sachant que

${\displaystyle {\mathcal {L}}\left\{\mathrm {I} f\right\}(s)={\mathcal {L}}\left\{\int _{0}^{t}f(\tau )\,\mathrm {d} \tau \right\}(s)={\frac {1}{s}}\left({\mathcal {L}}\left\{f\right\}\right)(s)}$

et

${\displaystyle {\mathcal {L}}\left\{\mathrm {I} ^{2}f\right\}={\frac {1}{s}}\left({\mathcal {L}}\left\{\mathrm {I} f\right\}\right)(s)={\frac {1}{s^{2}}}\left({\mathcal {L}}\left\{f\right\}\right)(s)}$

et ainsi de suite, on a

${\displaystyle \mathrm {I} ^{\alpha }f={\mathcal {L}}^{-1}\left\{s^{-\alpha }\left({\mathcal {L}}\{f\}\right)(s)\right\}}$ .

Par exemple,

${\displaystyle \mathrm {I} ^{\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}}$

comme prévu. En effet, étant donné la règle de convolution

${\displaystyle {\mathcal {L}}\{f*g\}={\bigl (}{\mathcal {L}}\{f\}{\bigr )}{\bigl (}{\mathcal {L}}\{g\}{\bigr )}}$

et en raccourcissant p(x) = x^{α − 1} pour plus de clarté, on trouve

${\displaystyle {\begin{aligned}\left(\mathrm {I} ^{\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 )\,\mathrm {d} \tau \\&={\frac {1}{\Gamma (\alpha )}}\int _{0}^{t}\left(t-\tau \right)^{\alpha -1}f(\tau )\,\mathrm {d} \tau \\\end{aligned}}}$

ce qui revient au résultat obtenu par la formule de Cauchy donnée plus haut.

Les transformations de Laplace "fonctionnent" sur relativement peu de fonctions, mais elles sont souvent utiles pour résoudre des équations différentielles fractionnaires.

• Fonctions puissances

L'intégrale de Riemann-Liouville de ${\textstyle f:x\mapsto x^{\beta }1\!\!1_{\mathbb {R} _{+}}(x)}$ est

${\displaystyle \mathrm {I} ^{\alpha }f(x)={\frac {\Gamma (\beta +1)}{\Gamma (\alpha +\beta +1)}}x^{\alpha +\beta }}$

• Fonctions exponentielles

L'intégrale de Riemann-Liouville de ${\textstyle f:x\mapsto \mathrm {e} ^{kx}1\!\!1_{\mathbb {R} _{+}}(x)}$ est

${\displaystyle \mathrm {I} ^{\alpha }f(x)={\frac {1}{k^{\alpha }}}\sum _{p=0}^{+\infty }{\frac {(kx)^{p+a}}{\Gamma (\alpha +k+1)}}={\frac {(kx)^{a}}{k^{\alpha }}}E_{1,\alpha +1}(kx)}$

où E désigne la fonction de Mittag-Leffler.

• (en) Yu. A. Brychkov et A.P. Prudnikov, « Intégrale de Riemann-Liouville », dans Michiel Hazewinkel, Encyclopædia of Mathematics, Springer, 2002 (ISBN 978-1556080104, lire en ligne).
• (en) P.I. Lizorkin, « Intégrale de Riemann-Liouville », dans Michiel Hazewinkel, Encyclopædia of Mathematics, Springer, 2002 (ISBN 978-1556080104, lire en ligne).
• (en) Michele Caputo, « Linear model of dissipation whose Q is almost frequency independent. II », Geophysical Journal International, vol. 13, n^o 5,‎ 1967, p. 529–539 (DOI 10.1111/j.1365-246x.1967.tb02303.x , Bibcode 1967GeoJ...13..529C).
• (en) Einar Hille et Ralph S. Phillips, Functional analysis and semi-groups, Providence, R.I., American Mathematical Society, 1974 (MR 0423094).
• (en) Adam Loverro, Fractional Calculus: History, Definitions and Applications for the Engineer, Notre Dame, IN, University of Notre Dame, 8 mai 2004 (lire en ligne [archive du 29 octobre 2005])
• (en) Kenneth S. Miller et Bertram Ross, An Introduction to the Fractional Calculus and Fractional Differential Equations, John Wiley & Sons, 1993 (ISBN 0-471-58884-9).
• Marcel Riesz, « L'intégrale de Riemann-Liouville et le problème de Cauchy », Acta Mathematica, vol. 81, n^o 1,‎ 1949, p. 1–223 (ISSN 0001-5962, DOI 10.1007/BF02395016 , MR 0030102).
• Marcel Riesz, « L’intégrale de Riemann-Liouville et le problème de Cauchy pour l’équation des ondes », Bulletin de la S. M. F., vol. 67, n^o 1,‎ 1939, p. 153-170 (ISSN 0001-5962, lire en ligne)

• Symbole d'un opérateur différentiel

• Alan Beardon, « Fractional calculus II », University of Cambridge, 2000
• Alan Beardon, « Fractional calculus III », University of Cambridge, 2000