Calcul fractionnaire - Fractional calculus

Le calcul fractionnaire est une branche de l'analyse mathématique qui étudie les différentes possibilités de définir des puissances de nombres réels ou des puissances de nombres complexes de l' opérateur de différenciation D

et de l' opérateur d' intégration J

et développer un calcul pour de tels opérateurs généralisant le classique.

Dans ce contexte, le terme puissance fait référence à l'application itérative d'un opérateur linéaire D à une fonction f , c'est-à-dire à la composition répétée de D avec elle-même, comme dans .

Par exemple, on peut demander une interprétation significative de

comme un analogue de la racine carrée fonctionnelle de l'opérateur de différenciation, c'est-à-dire une expression pour un opérateur linéaire qui, appliqué deux fois à n'importe quelle fonction, aura le même effet que la différenciation . Plus généralement, on peut se poser la question de la définition d'un opérateur linéaire

pour chaque nombre réel de telle sorte que, quand un prend un nombre entier valeur nZ , elle coïncide avec l'habituel n différenciation -fois D si n > 0 , et avec le (- n ) pouvoir -ième de J lorsque n < 0 .

L'une des motivations derrière l'introduction et l'étude de ces sortes d'extensions de l'opérateur de différenciation D est que les ensembles d'opérateurs puissances { D a | aR } ainsi définis sont des semi-groupes continus de paramètre a , dont le semi-groupe discret original de { D n | n ∈ Z } pour entier n est un sous-groupe dénombrable : puisque les semi-groupes continus ont une théorie mathématique bien développée, ils peuvent être appliqués à d'autres branches des mathématiques.

Les équations différentielles fractionnaires , également connues sous le nom d'équations différentielles extraordinaires, sont une généralisation des équations différentielles grâce à l'application du calcul fractionnaire.

Notes historiques

En mathématiques appliquées et en analyse mathématique, une dérivée fractionnaire est une dérivée de tout ordre arbitraire, réel ou complexe. Sa première apparition est dans une lettre écrite à Guillaume de l'Hôpital par Gottfried Wilhelm Leibniz en 1695. Vers la même époque, Leibniz écrit à l'un des frères Bernoulli décrivant la similitude entre le théorème du binôme et la règle de Leibniz pour la dérivée fractionnaire de un produit de deux fonctions. Le calcul fractionnaire a été introduit dans l'un des premiers articles de Niels Henrik Abel où tous les éléments peuvent être trouvés : l'idée d'intégration et de différenciation d'ordre fractionnaire, la relation mutuellement inverse entre elles, la compréhension que la différenciation et l'intégration d'ordre fractionnaire peuvent être considérée comme la même opération généralisée, et même la notation unifiée pour la différenciation et l'intégration de l'ordre réel arbitraire. Indépendamment, les bases du sujet ont été posées par Liouville dans un article de 1832. L' autodidacte Oliver Heaviside a introduit l'utilisation pratique des opérateurs différentiels fractionnaires dans l'analyse des lignes de transmission électrique vers 1890. La théorie et les applications du calcul fractionnaire se sont considérablement développées au cours du 19e et 20e siècles, et de nombreux contributeurs ont donné des définitions pour les dérivées fractionnaires et les intégrales.

Nature de la dérivée fractionnaire

L' un ième dérivée d'une fonction f  ( x ) en un point x est une propriété locale uniquement lorsque a est un nombre entier; ce n'est pas le cas pour les dérivés de puissance non entiers. En d'autres termes, une dérivée fractionnaire non entière d'une fonction f  ( x ) à x = a dépend de toutes les valeurs de f , même celles qui sont éloignées de a . Par conséquent, on s'attend à ce que l'opération de dérivée fractionnaire implique une sorte de conditions aux limites , impliquant des informations sur la fonction plus loin.

La dérivée fractionnaire d'une fonction d'ordre a est souvent maintenant définie au moyen des transformées intégrales de Fourier ou de Mellin .

Heuristique

Une question assez naturelle à se poser est de savoir s'il existe un opérateur linéaire H , ou semi-dérivé, tel que

Il s'avère qu'il existe un tel opérateur, et en effet pour tout a > 0 , il existe un opérateur P tel que

ou pour le dire autrement, la définition de d n oui/dx npeut être étendu à toutes les valeurs réelles de n .

Soit f  ( x ) une fonction définie pour x > 0 . Former l'intégrale définie de 0 à x . Appelez ceci

La répétition de ce processus donne

et cela peut être étendu arbitrairement.

La formule de Cauchy pour l'intégration répétée , à savoir

conduit de manière directe à une généralisation pour le réel n .

L'utilisation de la fonction gamma pour supprimer la nature discrète de la fonction factorielle nous donne un candidat naturel pour les applications fractionnaires de l'opérateur intégral.

Il s'agit en fait d'un opérateur bien défini.

Il est simple de montrer que l' opérateur J satisfait

Preuve  —

où, à la dernière étape, nous avons échangé l'ordre d'intégration et extrait le facteur f  ( s ) de l' intégration t . Changer les variables en r défini par t = s + ( xs ) r ,

L'intégrale interne est la fonction bêta qui vérifie la propriété suivante :

Remplacement dans l'équation

Interchangeant α et ß montre que l'ordre dans lequel la J opérateur est appliqué est sans importance et achève la démonstration.

Cette relation est appelée la propriété de semi-groupe des opérateurs différentiels fractionnaires . Malheureusement, le processus comparable pour l'opérateur dérivé D est significativement plus complexe, mais on peut montrer que D n'est ni commutatif ni additif en général.

Dérivée fractionnaire d'une fonction puissance de base

La demi-dérivée (courbe violette) de la fonction f  ( x ) = x (courbe bleue) avec la dérivée première (courbe rouge).
L'animation montre l'opérateur dérivé oscillant entre la primitive ( α = −1 : y =1/2x 2 ) et la dérivée ( α = +1 : y = 1 ) de la fonction puissance simple y = x en continu.

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

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

La répétition de ceci donne le résultat plus général que

qui, après avoir remplacé les factorielles par la fonction gamma , conduit à

Pour k = 1 et a =1/2, on obtient la demi-dérivée de la fonction x comme

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

(car et (1) = 1 ) qui est bien le résultat attendu de

Pour une puissance entière négative k , la fonction gamma est indéfinie et nous devons utiliser la relation suivante :

Cette extension de l'opérateur différentiel ci-dessus n'a pas besoin d'être limitée aux seules puissances réelles ; elle s'applique également aux puissances complexes. Par exemple, la dérivée (1 + i ) -ième de la dérivée (1 − i ) -ième donne la dérivée seconde. La définition de valeurs négatives pour a donne également des intégrales.

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

Pour α arbitraire , puisque la fonction gamma est indéfinie pour les arguments dont la partie réelle est un entier négatif et dont la partie imaginaire est nulle, 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,

transformation de Laplace

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

et

et ainsi de suite, nous affirmons

.

Par exemple,

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

et en abrégé p ( x ) = x α − 1 pour plus de clarté, nous trouvons que

c'est ce que Cauchy nous a donné plus haut.

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

Intégrales fractionnaires

Intégrale fractionnaire de Riemann-Liouville

La forme classique du calcul fractionnaire est donnée par l' intégrale de Riemann-Liouville , qui est essentiellement ce qui a été décrit ci-dessus. La théorie des fonctions périodiques (incluant donc la "condition aux limites" de répétition après une période) est l' intégrale de Weyl . Il est défini sur les séries de Fourier , et nécessite que le coefficient de Fourier constant s'annule (ainsi, il s'applique aux fonctions sur le cercle unité dont les intégrales s'évaluent à zéro). L'intégrale de Riemann-Liouville existe sous deux formes, supérieure et inférieure. Compte tenu de l' intervalle [ a , b ] , les intégrales sont définies comme

Où le premier est valable pour t > a et le dernier est valable pour t < b .

En revanche, la dérivée de Grünwald-Letnikov commence par la dérivée au lieu de l'intégrale.

Intégrale fractionnaire d'Hadamard

L' intégrale fractionnaire d'Hadamard est introduite par Jacques Hadamard et est donnée par la formule suivante,

Intégrale fractionnaire Atangana-Baleanu

L'intégrale fractionnaire d'Atangana-Baleanu d'une fonction continue est définie comme :

Dérivés fractionnaires

Contrairement aux dérivées newtoniennes classiques, une dérivée fractionnaire est définie via une intégrale fractionnaire.

Dérivées fractionnaires d'une gaussienne, interpolant en continu entre la fonction et sa première dérivée.

Dérivée fractionnaire de Riemann-Liouville

La dérivée correspondante est calculée en utilisant la règle de Lagrange pour les opérateurs différentiels. Le calcul de n ième dérivée d'ordre au-dessus de l'intégrale d'ordre ( n - α ) , l' α dérivée d'ordre est obtenue. Il est important de remarquer que n est le plus petit entier supérieur à α ( c'est-à-dire n = ⌈ α ). Semblable aux définitions de l'intégrale de Riemann-Liouville, la dérivée a des variantes supérieure et inférieure.

Dérivée fractionnaire Caputo

Une autre option pour calculer les dérivées fractionnaires est la dérivée fractionnaire Caputo. Il a été introduit par Michele Caputo dans son article de 1967. Contrairement à la dérivée fractionnaire de Riemann-Liouville, lors de la résolution d'équations différentielles en utilisant la définition de Caputo, il n'est pas nécessaire de définir les conditions initiales d'ordre fractionnaire. La définition de Caputo est illustrée comme suit, où encore n = ⌈ α :

Il existe la dérivée fractionnaire de Caputo définie comme :

qui a l'avantage d'être nul lorsque f  ( t ) est constante et sa transformée de Laplace s'exprime à l'aide des valeurs initiales de la fonction et de sa dérivée. De plus, il existe la dérivée fractionnaire de Caputo d'ordre distribué définie comme
φ ( ν ) est une fonction de pondération et qui est utilisée pour représenter mathématiquement la présence de plusieurs formalismes de mémoire.

Dérivé fractionnaire de Caputo-Fabrizio

Dans un article de 2015, M. Caputo et M. Fabrizio ont présenté une définition de dérivée fractionnaire à noyau non singulier, pour une fonction de donnée par :

Dérivé fractionnaire Atangana-Baleanu

En 2016, Atangana et Baleanu ont proposé des opérateurs différentiels basés sur la fonction généralisée de Mittag-Leffler . L'objectif était d'introduire des opérateurs différentiels fractionnaires à noyau non local non singulier. Leurs opérateurs différentiels fractionnaires sont donnés ci-dessous respectivement au sens de Riemann-Liouville et au sens de Caputo. Pour une fonction de donnée par

Si la fonction est continue, la dérivée d'Atangana-Baleanu au sens de Riemann-Liouville est donnée par :

Le noyau utilisé dans la dérivée fractionnaire d'Atangana-Baleanu a certaines propriétés d'une fonction de distribution cumulative. Par exemple, pour tous , la fonction est de plus en plus sur la ligne réelle, converge vers dans et . Par conséquent, nous avons que la fonction est la fonction de distribution cumulative d'une mesure de probabilité sur les nombres réels positifs. La distribution est donc définie, et n'importe lequel de ses multiples, est appelé une distribution d'ordre Mittag-Leffler . Il est également bien connu que toutes ces distributions de probabilités sont absolument continues . En particulier, la fonction Mittag-Leffler a un cas particulier , qui est la fonction exponentielle, la distribution d'ordre Mittag-Leffler est donc une distribution exponentielle . Cependant, pour , les distributions de Mittag-Leffler sont à queue lourde . Leur transformée de Laplace est donnée par :

Cela implique directement que, pour , l'espérance est infinie. De plus, ces distributions sont des distributions géométriques stables .

Dérivé de Riesz

Le dérivé de Riesz est défini comme

où désigne la transformée de Fourier .

Autres types

Les dérivés fractionnaires classiques comprennent :

Les nouveaux dérivés fractionnaires comprennent :

  • Dérivé de Coimbra
  • Dérivé de Katugampola
  • Dérivé de Hilfer
  • Dérivé de Davidson
  • Dérivé de Chen
  • Dérivé de Caputo Fabrizio
  • Dérivé Atangana-Baeanu

Généralisations

Opérateur Erdélyi-Kober

L' opérateur d'Erdélyi-Kober est un opérateur intégral introduit par Arthur Erdélyi (1940). et Hermann Kober (1940) et est donnée par

qui généralise l' intégrale fractionnaire de Riemann-Liouville et l'intégrale de Weyl.

Calcul fonctionnel

Dans le cadre de l'analyse fonctionnelle , les fonctions f  ( D ) plus générales que les puissances sont étudiées dans le calcul fonctionnel de la théorie spectrale . La théorie des opérateurs pseudo-différentiels permet aussi de considérer des puissances de D . Les opérateurs apparaissant sont des exemples d' opérateurs intégraux singuliers ; et la généralisation de la théorie classique à des dimensions supérieures s'appelle la théorie des potentiels de Riesz . Il existe donc un certain nombre de théories contemporaines disponibles, au sein desquelles le calcul fractionnaire peut être discuté. Voir aussi l' opérateur d'Erdélyi-Kober , important dans la théorie des fonctions spéciales ( Kober 1940 ), ( Erdélyi 1950–51 ).

Applications

Conservation fractionnée de la masse

Comme décrit par Wheatcraft et Meerschaert (2008), une équation de conservation fractionnaire de la masse est nécessaire pour modéliser l'écoulement de fluide lorsque le volume de contrôle n'est pas assez grand par rapport à l'échelle d' hétérogénéité et lorsque le flux dans le volume de contrôle est non linéaire. Dans l'article référencé, l'équation de conservation fractionnaire de la masse pour l'écoulement de fluide est :

Analyse électrochimique

Lors de l'étude du comportement redox d'un substrat en solution, une tension est appliquée à une surface d'électrode pour forcer le transfert d'électrons entre l'électrode et le substrat. Le transfert d'électrons qui en résulte est mesuré sous forme de courant. Le courant dépend de la concentration du substrat à la surface de l'électrode. Au fur et à mesure que le substrat est consommé, le substrat frais diffuse vers l'électrode comme décrit par la loi de Fick . Prendre la transformée de Laplace de la deuxième loi de Fick donne une équation différentielle ordinaire du second ordre (ici sous forme sans dimension):

dont la solution C(x,s) contient une demi-puissance dépendante de s. Prendre la dérivée de C(x,s) puis la transformée de Laplace inverse donne la relation suivante :

qui relie la concentration de substrat à la surface de l'électrode au courant. Cette relation est appliquée en cinétique électrochimique pour élucider le comportement mécaniste. Par exemple, il a été utilisé pour étudier le taux de dimérisation de substrats lors d'une réduction électrochimique.

Problème d'écoulement des eaux souterraines

En 2013-2014, Atangana et al. décrit quelques problèmes d'écoulement des eaux souterraines en utilisant le concept de dérivée avec ordre fractionnaire. Dans ces travaux, la loi de Darcy classique est généralisée en considérant le débit d'eau en fonction d'une dérivée d'ordre non entier de la hauteur piézométrique. Cette loi généralisée et la loi de conservation de la masse sont ensuite utilisées pour dériver une nouvelle équation pour l'écoulement des eaux souterraines.

Équation de dispersion d'advection fractionnaire

Cette équation s'est avérée utile pour modéliser le flux de contaminants dans des milieux poreux hétérogènes.

Atangana et Kilicman ont étendu l'équation de dispersion d'advection fractionnaire à une équation d'ordre variable. Dans leur travail, l'équation de dispersion hydrodynamique a été généralisée en utilisant le concept de dérivée d'ordre variationnel . L'équation modifiée a été résolue numériquement via la méthode Crank-Nicolson . La stabilité et la convergence dans les simulations numériques ont montré que l'équation modifiée est plus fiable pour prédire le mouvement de la pollution dans les aquifères déformables que les équations à dérivées fractionnaires et entières constantes

Modèles d'équations de diffusion fractionnaire dans l'espace-temps

Les processus de diffusion anormaux dans les milieux complexes peuvent être bien caractérisés en utilisant des modèles d'équation de diffusion d'ordre fractionnaire. Le terme dérivé dans le temps correspond à la décroissance de queue lourde à long terme et à la dérivée spatiale pour la non-localité de diffusion. L'équation régissant la diffusion fractionnaire dans l'espace-temps peut s'écrire sous la forme

Une simple extension de la dérivée fractionnaire est la variable dérivée d'ordre fractionnaire, α et β sont modifiés en α ( x , t ) et β ( x , t ) . Ses applications dans la modélisation de la diffusion anormale peuvent être trouvées dans la référence.

Modèles d'amortissement structurel

Les dérivés fractionnaires sont utilisés pour modéliser l' amortissement viscoélastique dans certains types de matériaux comme les polymères.

Contrôleurs PID

Généraliser les contrôleurs PID pour utiliser des ordres fractionnaires peut augmenter leur degré de liberté. La nouvelle équation reliant la variable de contrôle u ( t ) en termes d' une valeur d' erreur mesurée e ( t ) peut être écrite sous la forme

α et β sont des ordres fractionnaires positifs et K p , K i et K d , non-négatif, désignent les coefficients pour les proportionnels , intégraux et dérivés termes, respectivement (parfois notée P , I et D ).

Équations d'ondes acoustiques pour milieux complexes

La propagation des ondes acoustiques dans des milieux complexes, comme dans les tissus biologiques, implique couramment une atténuation obéissant à une loi de puissance en fréquence. Ce type de phénomène peut être décrit à l'aide d'une équation d'onde causale qui intègre des dérivées fractionnaires en temps :

Voir également Holm & Näsholm (2011) et les références qui s'y trouvent. De tels modèles sont liés à l'hypothèse communément admise que de multiples phénomènes de relaxation donnent lieu à l'atténuation mesurée dans les milieux complexes. Ce lien est décrit plus en détail dans Näsholm & Holm (2011b) et dans le document d'enquête, ainsi que dans l' article sur l' atténuation acoustique . Voir Holm & Nasholm (2013) pour un article qui compare les équations d'ondes fractionnaires qui modélisent l'atténuation de la loi de puissance. Ce livre sur l'atténuation de la loi de puissance couvre également le sujet plus en détail.

Pandey et Holm ont donné un sens physique aux équations différentielles fractionnaires en les dérivant de principes physiques et en interprétant l'ordre fractionnaire en termes de paramètres du milieu acoustique, par exemple dans les sédiments marins non consolidés granulaires saturés de fluide. Fait intéressant, Pandey et Holm ont dérivé la loi de Lomnitz en sismologie et la loi de Nutting en rhéologie non newtonienne en utilisant le cadre du calcul fractionnaire. La loi de Nutting a été utilisée pour modéliser la propagation des ondes dans les sédiments marins à l'aide de dérivés fractionnaires.

Équation fractionnaire de Schrödinger en théorie quantique

L' équation fractionnaire de Schrödinger , une équation fondamentale de la mécanique quantique fractionnaire , a la forme suivante :

où la solution de l'équation est la fonction d' onde ψ ( r , t ) - la mécanique quantique amplitude de probabilité pour que la particule d'avoir une donnée vecteur de position r à un moment donné t , et ħ est la constante de Planck réduite . La fonction d' énergie potentielle V ( r , t ) dépend du système.

De plus, =2/r 2est l' opérateur de Laplace , et D α est une constante d'échelle de dimension physique [ D α ] = J 1 − α ·m α ·s α = kg 1 − α ·m 2 − α ·s α − 2 , (à α = 2 , D 2 =1/2 mètrespour une particule de masse m ), et l'opérateur (− ħ 2 Δ) α /2 est la dérivée fractionnaire de Riesz tridimensionnelle définie par

L'indice α dans l'équation fractionnaire de Schrödinger est l'indice de Lévy, 1 < α ≤ 2 .

Équation de Schrödinger fractionnaire d'ordre variable

En tant que généralisation naturelle de l' équation fractionnaire de Schrödinger , l'équation fractionnaire de Schrödinger d'ordre variable a été exploitée pour étudier les phénomènes quantiques fractionnaires :

=2/r 2est l' opérateur de Laplace et l'opérateur (− ħ 2 Δ) β ( t )/2 est la dérivée fractionnaire de Riesz quantique d'ordre variable.

Voir également

Autres théories fractionnaires

Remarques

Les références

Sources

  • Kilbas, Anatolii Alexandrovitch ; Srivastava, Hari Mohan ; Trujillo, Juan J. (2006). Théorie et applications des équations différentielles fractionnaires . Amsterdam, Pays-Bas : Elsevier. ISBN 978-0-444-51832-3.

Lectures complémentaires

Articles concernant l'histoire du calcul fractionnaire

  • Ross, B. (1975). « Une brève histoire et une exposition de la théorie fondamentale du calcul fractionnaire ». Le calcul fractionnaire et ses applications . Calcul fractionnaire et ses applications. Notes de cours en mathématiques . Notes de cours en mathématiques. 457 . p. 1–36. doi : 10.1007/BFb0067096 . ISBN 978-3-540-07161-7.
  • Debnath, L. (2004). « Une brève introduction historique au calcul fractionnaire ». Journal international de l'enseignement des mathématiques en science et technologie . 35 (4) : 487–501. doi : 10.1080/00207390410001686571 . S2CID  122198977 .
  • Tenreiro Machado, J. ; Kiryakova, V.; Mainardi, F. (2011). « Histoire récente du calcul fractionnaire ». Communications en sciences non linéaires et simulation numérique . 16 (3) : 1140-1153. Bibcode : 2011CNSNS..16.1140M . doi : 10.1016/j.cnsns.2010.05.027 . hdl : 10400.22/4149 .
  • Tenreiro Machado, juge d'instruction ; Galhano, AM; Trujillo, JJ (2013). « La métrique scientifique sur le développement du calcul fractionnaire depuis 1966 ». Calcul fractionnaire et analyse appliquée . 16 (2) : 479-500. doi : 10.2478/s13540-013-0030-y . hdl : 10400.22/3773 . S2CID  122487513 .
  • Tenreiro Machado, juge d'instruction ; Galhano, AMSF; Trujillo, JJ (2014). « Sur le développement du calcul fractionnaire au cours des cinquante dernières années ». Scientométrie . 98 (1) : 577-582. doi : 10.1007/s11192-013-1032-6 . hdl : 10400.22/3769 . S2CID  16501850 .

Livres

Liens externes