Analyse fractionnaire

L'analyse fractionnaire est une branche de l'analyse mathématique qui étudie la possibilité de définir des puissances non entières des opérateurs de dérivation et d'intégration.

Ces dérivées ou intégrations fractionnaires rentrent dans le cadre plus général des opérateurs pseudo-différentiels.

Par exemple, on peut se demander comment interpréter convenablement la racine carrée

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

de l'opérateur de dérivation, c'est-à-dire une expression d'un certain opérateur qui, lorsqu'elle est appliquée deux fois à une fonction, aura le même effet que la dérivation. Plus généralement, on peut examiner le problème de définir

$\mathrm {D} ^{\alpha }$

pour des valeurs réelles de α, de telle sorte que lorsque α prend une valeur entière n, on récupère la dérivation n-ième usuelle pour n > 0 ou l'intégration itérée |n| fois pour n < 0. Le terme « fractionnaire » est utilisé de façon impropre : α n'est pas nécessairement un nombre rationnel, et l'on devrait donc plutôt parler de dérivation non entière. Cependant, le terme « analyse fractionnaire » est devenu traditionnel.

Les dérivées fractionnaires sont utilisées par exemple dans certains domaines de la physique faisant intervenir des phénomènes de diffusion comme l'électromagnétisme, l'acoustique ou la thermique, en définissant des opérateurs pseudo-différentiels diffusifs, avec conditions de bord à « géométrie fractale ».

Dérivée fractionnaire

Les fondations de ce sujet ont été jetées par Liouville dans un article de 1832. La dérivée fractionnaire d'ordre α d'une fonction en un point x est désormais souvent définie à partir de la transformée de Fourier ou de la transformée de Laplace.

Un point important est que la dérivée fractionnaire d'une fonction en un point x est une propriété locale seulement lorsque l'ordre α est entier ; dans les autres cas, on ne peut plus dire que la dérivée fractionnaire d'une fonction f en x ne dépend que du graphe de f au voisinage de x, comme c'est le cas en ce qui concerne les ordres de dérivation entiers.

Pour illustrer ceci, introduisons l'opérateur « de translation » $T:f(x)\mapsto f(x+h)$ et l'opérateur identité $\mathrm {Id}$. La limite, lorsque h tend vers 0, de l'opérateur

$\Delta ={\frac {T-{\rm {\mathrm {I} d}}}{h}}$

correspond bien à l'opérateur de dérivation au premier ordre.

Grâce à la formule du binôme généralisée, on peut alors élever cet opérateur à une puissance non entière. On obtient ainsi une série infinie :

$\Delta ^{\alpha }={\frac {(T-\mathrm {\mathrm {I} d} )^{\alpha }}{h^{\alpha }}}={\frac {1}{h^{\alpha }}}\sum _{k=0}^{\infty }{\alpha \choose k}$$(-T)^{k}$^{[Information douteuse]},

où ${\alpha \choose k}$ désigne le coefficient binomial généralisé ${\alpha \choose k}={\frac {\alpha (\alpha -1)(\alpha -2)\dots (\alpha -k+1)}{k!}}$.

Une telle définition induit un caractère non local de l'opération de dérivation à un ordre non entier.

Approche naturelle

Une question naturelle qui se pose est : existe-t-il un opérateur linéaire H tel que

$\mathrm {H} ^{2}f=\mathrm {D} f={\frac {\mathrm {d} }{\mathrm {d} x}}f=f'$ ?

Il apparaît qu'il existe un tel opérateur et même, pour tout α > 0 rationnel, il existe un opérateur P tel que $\mathrm {P} ^{\alpha }=\mathrm {D}$ (plus précisément, si α = p/q, $\mathrm {P} ^{p}=\mathrm {D} ^{q}$) ou, pour le formuler autrement, que ${\tfrac {\mathrm {d} ^{n}f}{\mathrm {d} x^{n}}}$ est défini pour toutes valeurs réelles n > 0.

Un résultat similaire s'applique pour l'intégration. Considérant une fonction f définie pour x > 0, nous pouvons former son intégrale définie de 0 à x :

$(\mathrm {I} f)(x)=\int _{0}^{x}f(t)\;\mathrm {d} t$.

En répétant ce processus, on obtient

$(\mathrm {I} ^{2}f)(x)=\int _{0}^{x}(\mathrm {I} f)(t)\;\mathrm {d} t=\int _{0}^{x}\left(\int _{0}^{t}f(u)\;\mathrm {d} u\right)\mathrm {d} t$,

et ceci peut être répété arbitrairement.

La formule suivante, appelée formule de Cauchy pour l'intégration successive,

$(\mathrm {I} ^{n}f)(x)={\frac {1}{(n-1)!}}\int _{0}^{x}(x-t)^{n-1}f(t)\;\mathrm {d} t$

exprime par une seule intégrale une primitive n-ième d'une fonction f. Ceci mène tout droit à une généralisation pour tout réel α > 0 et même, pour tout nombre complexe de partie réelle strictement positive.

La fonction gamma Γ, qui étend la factorielle aux valeurs complexes, est définie de telle sorte que :

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

En utilisant la fonction gamma pour se libérer de la nature discrète de la factorielle, nous obtenons un candidat naturel pour les puissances non entières de l'opérateur d'intégration :

$(\mathrm {I} ^{\alpha }f)(x)={\frac {1}{\Gamma (\alpha )}}\int _{0}^{x}(x-t)^{\alpha -1}f(t)\;\mathrm {d} t\quad (\operatorname {Re} (\alpha )>0)$.

Exemple :

Si f est une fonction puissance, de la forme $f(x)=x^{k}\quad (k\in \mathbb {C} ,\;x>0)$, alors $(\mathrm {I} ^{\alpha }f)(x)={\frac {\Gamma (k+1)}{\Gamma (\alpha +k+1)}}x^{\alpha +k}$.

Démonstration

$\int _{0}^{x}(x-t)^{\alpha -1}t^{k}\;\mathrm {d} t=x^{\alpha +k}\int _{0}^{1}(1-u)^{\alpha -1}u^{k}\;\mathrm {d} u=x^{\alpha +k}\mathrm {B} (\alpha ,k+1)=x^{\alpha +k}{\frac {\Gamma (\alpha )\,\Gamma (k+1)}{\Gamma (\alpha +k+1)}}$, où B est la fonction bêta.

La famille d'opérateurs intégraux $(\mathrm {I} ^{\alpha })$ vérifie :

$\mathrm {I} ^{\alpha }\circ \mathrm {I} ^{\beta }=\mathrm {I} ^{\alpha +\beta }$ (donc $\mathrm {I} ^{\alpha }\circ \mathrm {I} ^{\beta }=\mathrm {I} ^{\beta }\circ \mathrm {I} ^{\alpha }$).

Démonstration

${\begin{aligned}\left(I^{\alpha }\right)\left(I^{\beta }f\right)(x)&={\frac {1}{\Gamma (\alpha )}}\int _{0}^{x}(x-t)^{\alpha -1}\left(I^{\beta }f\right)(t)\,\mathrm {d} t\\&={\frac {1}{\Gamma (\alpha )\Gamma (\beta )}}\int _{0}^{x}\left(\int _{0}^{t}\left(x-t\right)^{\alpha -1}\left(t-s\right)^{\beta -1}f(s)\,\mathrm {d} s\right)\,\mathrm {d} t\\&={\frac {1}{\Gamma (\alpha )\Gamma (\beta )}}\int _{0}^{x}f(s)\left(\int _{s}^{x}\left(x-t\right)^{\alpha -1}\left(t-s\right)^{\beta -1}\,\mathrm {d} t\right)\,\mathrm {d} s\end{aligned}}$

On définit une nouvelle variable r par t = s + (x − s)r et l'on obtient :

$\left(I^{\alpha }\right)\left(I^{\beta }f\right)(x)={\frac {1}{\Gamma (\alpha )\Gamma (\beta )}}\int _{0}^{x}\left(x-s\right)^{\alpha +\beta -1}f(s)\left(\int _{0}^{1}\left(1-r\right)^{\alpha -1}r^{\beta -1}\,\mathrm {d} r\right)\,\mathrm {d} s$.

L'intégrale à l'intérieur est la fonction bêta, qui vérifie :

$\int _{0}^{1}\left(1-r\right)^{\alpha -1}r^{\beta -1}\,\mathrm {d} r=\mathrm {B} (\alpha ,\beta )={\frac {\Gamma (\alpha )\,\Gamma (\beta )}{\Gamma (\alpha +\beta )}}$,

si bien que

$\left(I^{\alpha }\right)\left(I^{\beta }f\right)(x)={\frac {1}{\Gamma (\alpha +\beta )}}\int _{0}^{x}\left(x-s\right)^{\alpha +\beta -1}f(s)\,\mathrm {d} s=\left(I^{\alpha +\beta }f\right)(x)$.

Malheureusement, le processus analogue pour l'opérateur de dérivation D est considérablement plus compliqué. La famille $(\mathrm {D} ^{\alpha })$ n'est ni additive (on n'a pas en général $\mathrm {D} ^{\alpha }\circ \mathrm {D} ^{\beta }=\mathrm {D} ^{\alpha +\beta }$), ni même commutative ($\mathrm {D} ^{\alpha }\circ \mathrm {D} ^{\beta }\neq \mathrm {D} ^{\beta }\circ \mathrm {D} ^{\alpha }$).

Une définition élémentaire

L'idée la plus simple est de partir de formules « régulières » pour la dérivée n-ième et de remplacer n par le réel α ; on obtient ainsi pour l'exponentielle $\mathrm {\mathrm {D} } ^{\alpha }(\mathrm {e} ^{\lambda x})=\lambda ^{\alpha }\mathrm {e} ^{\lambda x}$ et pour la fonction sinus : $(\mathrm {\mathrm {D} } ^{\alpha }\sin )(x)=\sin \left(x+{\frac {\alpha \pi }{2}}\right).$ La même idée pour les fonctions puissances oblige comme précédemment à introduire la fonction gamma : puisque ${\frac {\mathrm {d} ^{n}}{\mathrm {d} x^{n}}}x^{k}={\frac {k!}{(k-n)!}}x^{k-n}$, on aura $\mathrm {\mathrm {D} } ^{\alpha }(x^{k})={\frac {\Gamma (k+1)}{\Gamma (k+1-\alpha )}}x^{k-\alpha }$. Pour ces formules, on a bien l'additivité ($\mathrm {D} ^{\alpha }\circ \mathrm {D} ^{\beta }=\mathrm {D} ^{\alpha +\beta }$), qui permet d'obtenir, par exemple, une « racine carrée » de la dérivation en prenant α = 1/2.

Mais cette approche élémentaire, non seulement n'est pas généralisable, mais contredit les définitions plus générales construites à partir d'opérateurs intégraux.

Définitions générales des dérivées fractionnaires

Dérivées fractionnaires d'une fonction gaussienne, interpolant entre la fonction et sa dérivée usuelle.

Dérivée de Riemann-Liouville

L'idée, pour définir une dérivée d'ordre $\alpha \in \mathbb {C}$ avec Re(α) ≥ 0, est ici de calculer la dérivée n-ième usuelle de l'intégrale fractionnaire d'ordre n − α, pour n = ⌊Re(α)⌋ + 1. Deux variantes de la définition existent :

$\mathrm {D} _{a+}^{\alpha }f(x)={\frac {\mathrm {d} ^{n}}{\mathrm {d} x^{n}}}\mathrm {D} _{a+}^{-(n-\alpha )}f(x)={\frac {\mathrm {d} ^{n}}{\mathrm {d} x^{n}}}\mathrm {I} _{a+}^{n-\alpha }f(x)$ (pour $x>a$) ;

$\mathrm {D} _{b-}^{\alpha }f(x)=\left(-{\frac {\mathrm {d} }{\mathrm {d} x}}\right)^{n}\mathrm {D} _{b-}^{-(n-\alpha )}f(x)=\left(-{\frac {\mathrm {d} }{\mathrm {d} x}}\right)^{n}\mathrm {I} _{b-}^{n-\alpha }f(x)$ (pour $x<b$).

Les familles d'opérateurs d'intégration $(\mathrm {I} _{a+}^{\alpha })$ et $(\mathrm {I} _{b-}^{\alpha })$ (avec Re(α) > 0) utilisées ici généralisent la famille $(\mathrm {I} ^{\alpha })=(\mathrm {I} _{0+}^{\alpha })$ définie ci-dessus pour x > 0 :

$(\mathrm {I} _{a+}^{\alpha })f(x)={\frac {1}{\Gamma (\alpha )}}\int _{a}^{x}(x-t)^{\alpha -1}f(t)\;\mathrm {d} t\quad (x>a)$ ;

$(\mathrm {I} _{b-}^{\alpha })f(x)={\frac {1}{\Gamma (\alpha )}}\int _{x}^{b}(t-x)^{\alpha -1}f(t)\;\mathrm {d} t\quad (x<b)$.

Exemple :

Reprenons l'exemple de la fonction puissance $f(x)=x^{k}$, dont nous avions calculé l'image par $\mathrm {I} ^{\alpha }$. Supposons Re(k) > −1 et Re(α) ≥ 0. Alors, pour x > 0,

$\mathrm {D} _{0+}^{\alpha }f(x)={\frac {\mathrm {d} ^{n}}{\mathrm {d} x^{n}}}{\frac {\Gamma (k+1)}{\Gamma (n-\alpha +k+1)}}x^{n-\alpha +k}={\frac {\Gamma (k+1)}{\cancel {\Gamma (n-\alpha +k+1)}}}{\frac {\cancel {\Gamma (n-\alpha +k+1)}}{\Gamma (-\alpha +k+1)}}x^{-\alpha +k}$,

donc la définition ad hoc ci-dessus de $\mathrm {\mathrm {D} } ^{\alpha }(x^{k})$ est en fait un cas particulier de dérivée de Riemann-Liouville.

Dérivée de Liouville-Weyl

Une variante naturelle de la dérivée de Riemann-Liouville, pour une fonction définie sur $\mathbb {R}$ tout entier, consiste à prendre $a=-\infty$ et $b=+\infty$, donc à poser, pour tout $x\in \mathbb {R}$ :

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

$(\mathrm {I} _{-}^{\alpha })f(x)={\frac {1}{\Gamma (\alpha )}}\int _{x}^{+\infty }(t-x)^{\alpha -1}f(t)\;\mathrm {d} t$

puis, comme précédemment,

$\mathrm {D} _{+}^{\alpha }={\frac {\mathrm {d} ^{n}}{\mathrm {d} x^{n}}}\circ \mathrm {I} _{+}^{n-\alpha }$ et $\mathrm {D} _{-}^{\alpha }=\left(-{\frac {\mathrm {d} }{\mathrm {d} x}}\right)^{n}\circ \mathrm {I} _{-}^{n-\alpha }$.

Exemple :

Soit $f(x)=\mathrm {e} ^{\lambda x}$ avec Re(λ) > 0, et soit Re(α) ≥ 0. Alors,

$(\mathrm {I} _{+}^{\alpha })f(x)=\lambda ^{-\alpha }\mathrm {e} ^{\lambda x}$ et

$(\mathrm {D} _{+}^{\alpha })f(x)=\lambda ^{\alpha }\mathrm {e} ^{\lambda x}$,

donc la définition ad hoc ci-dessus de $\mathrm {\mathrm {D} } ^{\alpha }(\mathrm {e} ^{\lambda x})$ est en fait un cas particulier de dérivée de Liouville-Weyl.

Dérivée de Caputo

En 1967, Michele Caputo introduisit une nouvelle définition ne nécessitant pas de conditions aux bornes. La définition de Caputo diffère de celle de Riemann-Liouville en ce qu'elle effectue la dérivation n fois avant l'intégrale fractionnaire d'ordre n − α :

${}^{C}\mathrm {D} ^{\alpha }f=\mathrm {I} ^{n-\alpha }\left(f^{(n)}\right)\qquad (n-1<\alpha <n)$.

Elle a l'avantage d'être nulle pour f constante et d'avoir une transformée de Laplace exprimée à l'aide de celle de f et des valeurs initiales de f^(k) pour 0 ≤ k < n.

Plus généralement, pour α complexe non entier tel que Re(α) ≥ 0 et pour f de classe Cⁿ sur [a, b] avec n = ⌈Re(α)⌉, on définit

${}^{C}\mathrm {D} _{a+}^{\alpha }f=\mathrm {D} _{a+}^{\alpha }\left[f(t)-\sum _{k=0}^{n-1}{\frac {y^{(k)}(a)}{k!}}(t-a)^{k}\right]\quad {\text{et}}\quad ^{C}\mathrm {D} _{b-}^{\alpha }f=\mathrm {D} _{b-}^{\alpha }\left[f(t)-\sum _{k=0}^{n-1}{\frac {y^{(k)}(b)}{k!}}(b-t)^{k}\right]$,

et l'on démontre que

${}^{C}\mathrm {D} _{a+}^{\alpha }f=\mathrm {I} _{a+}^{n-\alpha }\left(f^{(n)}\right)\quad {\text{et}}\quad ^{C}\mathrm {D} _{b-}^{\alpha }f=(-1)^{n}\mathrm {I} _{b-}^{n-\alpha }\left(f^{(n)}\right)$.

Dérivée de Grünwald-Letnikov

On généralise d'abord la différence finie arrière d'ordre entier en posant, pour α > 0 :

$\Delta _{h}^{\alpha }f(x)=\sum _{k=0}^{\infty }(-1)^{k}{\alpha \choose k}f(x+(\alpha -k)h)$,

puis

${}^{GL}\mathrm {D} _{+}^{\alpha }f(x)=\lim _{h\to 0^{+}}{\frac {\Delta _{h}^{\alpha }f(x)}{h^{\alpha }}}\quad {\text{et}}\quad ^{GL}\mathrm {D} _{-}^{\alpha }f(x)=\lim _{h\to 0^{+}}{\frac {\Delta _{-h}^{\alpha }f(x)}{h^{\alpha }}}$.

Exemples :

1. Soit $f(x)=\mathrm {e} ^{\lambda x}$ avec Re(λ) ≥ 0, et soit α > 0. Alors,$\forall x\in \mathbb {R} \quad ^{GL}\mathrm {D} _{+}^{\alpha }f(x)=\lambda ^{\alpha }\mathrm {e} ^{\lambda x}$,donc la définition ad hoc ci-dessus de $\mathrm {\mathrm {D} } ^{\alpha }(\mathrm {e} ^{\lambda x})$ est en fait un cas particulier de dérivée (non seulement de Liouville-Weyl, mais aussi) de Grünwald-Letnikov.
2. Soit $g(x)=\sin(\mu x)$ avec μ ∈ ℝ, et soit α > 0. Alors, d'après l'exemple 1,$\forall x\in \mathbb {R} \quad ^{GL}\mathrm {D} _{+}^{\alpha }g(x)=\mu ^{\alpha }\sin \left(\mu x+{\frac {\alpha \pi }{2}}\right)$,donc la définition ad hoc ci-dessus de $\mathrm {\mathrm {D} } ^{\alpha }{\bigl (}\sin(x){\bigr )}$ est en fait un cas particulier de dérivée de Grünwald-Letnikov (ru).

Applications

Ce champ des mathématiques trouve une application pratique en automatique via la commande CRONE (Commande Robuste d'Ordre Non Entier).

Notes et références

(en) Cet article est partiellement ou en totalité issu des articles intitulés en anglais « Differintegral » (voir la liste des auteurs) et « Fractional calculus » (voir la liste des auteurs).

Voir aussi

Articles connexes

• Initialisation de l'analyse fractionnaire (en)
• Opérateur différentiel
• Théorème de Hille-Yosida

Liens externes

• (en) Eric W. Weisstein, « Fractional calculus », sur MathWorld
• (en) Eric W. Weisstein, « Fractional derivative », sur MathWorld
• (en) Igor Podlubny, « Fractional calculus: Resources »

Ouvrages

• (en) A. Carpinteri et F. Mainardi, Fractals and Fractional Calculus in Continuum Mechanics, Springer, 1998, 348 p., relié (ISBN 978-3-211-82913-4, lire en ligne)
• (en) Richard Herrmann, Fractional Calculus - An Introduction for Physicists, World Scientific, 2018, 3^e éd., 610 p. (DOI 10.1142/11107, lire en ligne)
• (en) F. Mainardi, Fractional Calculus and Waves in Linear Viscoelasticity: An Introduction to Mathematical Models, Imperial College Press, 2010, 368 p. (lire en ligne)
• (en) Kenneth S. Miller et Bertram Ross, An Introduction to the Fractional Calculus and Fractional Differential Equations, John Wiley & Sons, 1993, 1^re éd., 384 p., relié (ISBN 978-0-471-58884-9)
• (en) Keith B. Oldham et Jerome Spanier, The Fractional Calculus; Theory and Applications of Differentiation and Integration to Arbitrary Order, Academic Press, coll. « Mathematics in Science and Engineering » (n^o 111), 1974, relié (ISBN 978-0-12-525550-9, lire en ligne)
• (en) Igor Podlubny, Fractional Differential Equations: An Introduction to Fractional Derivatives, Fractional Differential Equations, to Methods of Their Solution and Some of Their Applications, Academic Press, coll. « Mathematics in Science and Engineering » (n^o 198), 1998 (ISBN 978-0-08-053198-4, lire en ligne)
• (en) Vasily E. Tarasov, Fractional Dynamics: Applications of Fractional Calculus to Dynamics of Particles, Fields and Media, Springer, 2010, 450 p. (présentation en ligne, lire en ligne)
• (en) Vladimir V. Uchaikin, Fractional Derivatives for Physicists and Engineers, Higher Education Press (en)/Springer, 2012, 385 p. (présentation en ligne, lire en ligne)
• (en) Bruce J. West, Mauro Bologna et Paolo Grigolini, Physics of Fractal Operators, Springer, 2003, 368 p., relié (ISBN 978-0-387-95554-4, Bibcode 2003PhT....56l..65W, lire en ligne)

Journaux spécialisés

• Fractional Calculus and Applied Analysis (FCAA) (journal de l'Institut de mathématiques de l'Académie bulgare des sciences)
• Fractional Differential Calculus (FDC)
• Journal of Fractional Calculus and Applications (JFCA)
• Progress in Fractional Differentiation and Applications (PFDA)