Transformation de Laplace
En mathématiques, la transformation de Laplace est une transformation intégrale, c'est-à-dire une opération associant à une fonction ƒ (définie sur les réels positifs et à valeurs réelles) une nouvelle fonction dite transformée de Laplace de ƒ (notée traditionnellement F et définie et à valeurs complexes), via une intégrale.
Remarque : on note traditionnellement t le paramètre générique de ƒ (formant ainsi ƒ(t)), tandis que l'on note plutôt p celui de sa transformée F (on écrit donc F(p)).
La transformation de Laplace est injective et par calcul (ou par usage de tables) il est possible d'inverser la transformation. Le grand avantage de la transformation de Laplace est que la plupart des opérations courantes sur la fonction originale ƒ(t), telle que la dérivation, ou une translation sur la variable t, ont une traduction (plus) simple sur la transformée F(p). Ainsi :
- la transformée de Laplace de la dérivée ƒ '(t) est simplement pF(p) – ƒ(0–) ;
- la transformée de la fonction ƒ(t – τ) (translation) est simplement e–pτ F(p).
Cette transformation fut introduite pour la première fois sous une forme proche de celle utilisée par Laplace en 1774, dans le cadre de la théorie des probabilités.
La transformation de Laplace généralise la transformation de Fourier qui est également utilisée pour résoudre les équations différentielles : contrairement à cette dernière, elle tient compte des conditions initiales et peut ainsi être utilisée en théorie des vibrations mécaniques ou en électricité dans l'étude des régimes forcés sans négliger le régime transitoire. Elle converge pour toutes les fonctions qui, pondérées par une exponentielle, admettent une transformée de Fourier ; par conséquent les fonctions admettant une transformée de Fourier admettent toutes une transformée de Laplace, mais la réciproque n'est pas vraie. De manière générale, ses propriétés vis-à-vis de la dérivation permettent un traitement plus simple de certaines équations différentielles, et elle est de ce fait très utilisée en automatique.
Dans ce type d'analyse, la transformation de Laplace est souvent interprétée comme un passage du domaine temps, dans lequel les entrées et sorties sont des fonctions du temps, dans le domaine des fréquences, dans lequel les mêmes entrées et sorties sont des fonctions de la « fréquence » (complexe) p. Ainsi; il est possible d'analyser simplement l'effet du système sur l'entrée pour donner la sortie en matière d'opérations algébriques simples (cf. théorie des fonctions de transfert en électronique ou en mécanique).
Définition
En mathématiques et en particulier en analyse fonctionnelle, la transformée de Laplace monolatérale d'une fonction ƒ (éventuellement généralisée, telle que la « fonction de Dirac ») d'une variable réelle t, à support positif, est la fonction F de la variable complexe p, définie par :
Plus précisément, cette formule est valide lorsque :
- Re(p) > α, où α est l'abscisse de convergence (définie plus bas), –∞ ≤ α ≤ +∞ ;
- et ƒ est une fonction localement intégrable à support positif, c'est-à-dire nulle en dehors de l'intervalle I = [0, +∞[, ou plus généralement un « germe » de distributions définies dans un voisinage ouvert (et borné inférieurement) de l'intervalle I = [0, +∞[ dont la restriction au complémentaire de I dans ce voisinage est une fonction indéfiniment dérivable (voir l'article Transformation bilatérale de Laplace)[1].
C'est un tel germe appelé ici, par abus de langage, une fonction généralisée à support positif, et la transformation de Laplace est injective appliquée à ces fonctions généralisées.
L'abscisse de convergence α se définit comme suit :
- soit, pour un réel β, . Alors α est la borne inférieure dans de l'ensemble B des β pour lesquels ƒβ est une distribution tempérée (donc α = +∞ si B est vide).
La « fonction de Dirac » est de cette nature. Sa transformée de Laplace vaut 1 avec une abscisse de convergence de –∞.
Les propriétés de cette transformation lui confèrent une grande utilité dans l'analyse des systèmes dynamiques linéaires. La plus intéressante de ces propriétés est que l'intégration et la dérivation sont transformées en division et multiplication par p, de la même manière que le logarithme transforme la multiplication en addition. Elle permet ainsi de ramener la résolution des équations différentielles linéaires à coefficients constants à la résolution d'équations affines (dont les solutions sont des fonctions rationnelles de p).
La transformation de Laplace est très utilisée par les ingénieurs pour résoudre des équations différentielles et déterminer la fonction de transfert d'un système linéaire. Par exemple, en électronique, contrairement à la décomposition de Fourier qui est utilisée pour la détermination du spectre d'un signal périodique ou même quelconque, elle tient compte de l'existence d'un régime transitoire précédant le régime permanent (exemple : la prise en compte de l'allure du signal avant et après la mise en marche d'un générateur de fréquence).
Il suffit en effet de transposer l'équation différentielle dans le domaine de Laplace pour obtenir une équation beaucoup plus simple à manipuler.
Par exemple, lors de l'étude d'une machine à courant continu :
dans le domaine fréquentiel devient
dans le domaine de Laplace. Ceci n'est valable qu'à conditions initiales nulles : i(0) = 0.
On a utilisé ici des propriétés de la transformation de Laplace, explicitées ci-dessous.
Remarque : la notation « s » (variable de Laplace) est souvent utilisée dans les pays anglo-saxons alors que la notation « p » est utilisée notamment en France et en Allemagne.
On définit aussi, dans les mêmes conditions que ci-dessus, la transformation de Laplace-Carson par[2] :
qui permet d'associer à toute fonction d'une variable une fonction image .
Cette transformation est utilisée par certains ingénieurs car :
- une constante sur [0, +∞[ a pour image la même constante ;
- elle offre dans certains cas une plus grande facilité d'emploi en calcul matriciel et tensoriel.
Inversion
L'inversion de la transformation de Laplace s'effectue par le biais d'une intégrale dans le plan complexe. À l'aide du théorème des résidus, on démontre la formule de Bromwich-Mellin :
où γ est choisi de sorte que :
- l'intégrale soit convergente, ce qui implique que γ soit supérieur à la partie réelle de toute singularité de F(p) ;
- et qu'à l'infini, |F(p)| tende vers 0 au moins aussi rapidement que .
Lorsque cette dernière condition n'est pas satisfaite, la formule ci-dessus est encore utilisable s'il existe un entier n tel que :
- |p–nF(p)| tende vers 0 aussi rapidement que
c'est-à-dire lorsque :
- pour |p| tendant vers l'infini, |F(p)| est majorée par un polynôme en |p|.
En remplaçant F(p) par p–nF(p) dans l'intégrale ci-dessus, on trouve dans le membre de gauche de l'égalité une fonction généralisée à support positif dont la dérivée d'ordre n (au sens des distributions) est la fonction généralisée (elle aussi à support positif) cherchée.
En pratique néanmoins, la formule de Bromwich-Mellin est peu utilisée, et on calcule les inverses des transformées de Laplace à partir des tables de transformées de Laplace.
Propriétés
Linéarité
La transformation de Laplace est linéaire c'est-à-dire que quelles soient les fonctions f, g et deux nombres complexes a et b :
- .
Cette linéarité découle évidemment de celle de l'intégrale.
Continuité
Si est continue et si l'intégrale impropre converge, alors est bien défini pour tout réel et est continue sur . En particulier, .
En effet, la règle d'Abel s'applique ici uniformément par rapport à x[3].
Holomorphie
La transformée de Laplace de est holomorphe et sa dérivée n-ième est (voir infra).
Transformation de Laplace d'une dérivée
Appliquée à la dérivée de f la transformation de Laplace correspond, à une constante additive près, à une multiplication par p de la transformée :
- [4].
De proche en proche ou par récurrence il est possible de montrer pour les dérivations successives[1] :
Cette dernière expression peut s'écrire, avec pour tout entier ,
Notons que, vu la définition donnée plus haut d'une fonction généralisée à support positif (en utilisant la notion de germe), les quantités ne sont pas nulles en général.
Si en revanche f est une fonction usuelle à support positif, 0– est à remplacer partout par 0+.
Plus précisément, écrivons où est l'échelon unité de Heaviside et g est une fonction continûment dérivable (au sens usuel) dans un voisinage de 0. Alors d'après la règle de Leibniz,
- avec
Puisque , , donc .
On a également car .
Maintenant, et . Par définition, car c'est de la transformation monolatérale qu'il s'agit. On obtient donc finalement
En continuant ce raisonnement, on obtient, si g est de classe dans un voisinage de [0, +∞[,
avec pour tout entier .
Exemple
Soit . Alors et . On a et
- . Par conséquent,
Application à la dérivée de la fonction de Heaviside
La fonction de Heaviside vaut 0 pour t < 0, 1 pour t > 0 (sa valeur en 0 n'a aucune importance). Cette fonction étant discontinue, elle n'est pas dérivable au sens habituel. En revanche, sa dérivée au sens des distributions est la « fonction » de Dirac . Il vient
puisque
On notera que si l'on remplaçait, dans la formule de la règle de dérivation, ƒ(0–) par ƒ(0+), on trouverait , ce qui est faux (on va y revenir plus loin). Certaines sources peuvent comporter cette erreur [5].
De même, on voit parfois, la définition suivante de la transformation de Laplace :
avec , voire un manque de précision sur cette limite [6]. Si f est une fonction au sens habituel de ce terme, à support positif, il s'agit d'une intégrale de Lebesgue qui coïncide avec celle correspondant à , puisque est de mesure nulle ; on peut d'ailleurs dans ce cas écrire sans ambiguïté . Il n'en va pas de même si f est une « fonction généralisée », c'est-à-dire une distribution pour Gelfand et Shilov (en), quand celle-ci a une masse non nulle à l'origine. Le prototype est la distribution de Dirac. Au plan algébrique, cette distribution est l'élément neutre dans l'algèbre de convolution des distributions à support positif ; et puisque la transformation de Laplace transforme le produit de convolution en produit ordinaire, il faut donc que ait pour transformée de Laplace . Or, cela ne sera vrai que si . En effet, avec on obtiendrait une transformée de Laplace égale à 0. Cela serait d'autant plus aberrant que la transformation de Laplace ne serait pas injective, puisque .
Multiplication par une puissance de t
La multiplication par dans le domaine temporel correspond, au signe près, à la dérivée n-ième de la transformée :
- .
La formule inverse (pour n = -1) est :
et elle est valide à condition que f soit de la forme où g est une fonction généralisée à support positif. Une manière de démontrer ce résultat est indiquée ci-dessous.
Intégration
La transformation de Laplace d'une intégrale (primitive de f s'annulant en 0) correspond à une multiplication par 1/p :
et si ƒ est une fonction à support positif, continue sur [0, +∞[, on a pour tout a > 0 :
Valeur finale
Supposons f localement intégrable à support positif. Si la limite dans le domaine temporel existe et est finie, alors :
(On notera que c'est la seule propriété où un 0+ apparaît pour la variable .)
Les hypothèses indiquées sont indispensables, comme le montrent les contre-exemples suivants :
- La fonction admet pour limite +∞ lorsque t tend vers +∞. Sa transformée de Laplace est et . Cette dernière limite n'a en réalité aucun sens car l'abscisse de convergence de F est 1, donc 0 n'appartient pas à l'adhérence du domaine de convergence.
- La fonction n'admet pas de limite lorsque t tend vers +∞. Sa transformée de Laplace est , l'abscisse de convergence de F est 0 et (cette dernière limite est pourtant correcte cette fois).
- Si est une fonction rationnelle, existe et est finie si, et seulement si les pôles de appartiennent tous à la réunion du demi-plan gauche ouvert et de l'origine, le pôle en 0, s'il existe, étant simple[8].
Valeur initiale
Si a une abscisse de convergence finie et si la limite dans le domaine temporel existe, alors :
(On notera que c'est la seule propriété où un 0+ apparaît pour la variable [8].)
Transformée de Laplace d'une fonction périodique
Si ƒ est une fonction nulle pour t < 0 et, pour t > 0, périodique de période T, alors pour
Tableau résumé des propriétés de la transformation de Laplace
Domaine temporel | Domaine "p" | Commentaires | |
---|---|---|---|
Linéarité | Résulte des règles de base de l'intégration. | ||
Dérivée de la transformée | est la dérivée première de F. | ||
Dérivées d'ordre n de la transformée | Forme plus générale, dérivée n-ième de F(p). | ||
Dérivée première de la fonction dans le domaine temporel | ƒ est supposée dérivable, et sa dérivée est supposée tendre vers 0 exponentiellement. Peut être obtenue par intégration par parties. | ||
Dérivée seconde | ƒ est supposée deux fois dérivable, avec dérivée seconde convergeant exponentiellement à l'infini. | ||
Dérivée n-ième de ƒ | ƒ est supposée n fois dérivable, avec une dérivée n-ième à convergence exponentielle à l'infini. | ||
Intégration de la transformée de Laplace | |||
Intégration | est la fonction échelon de Heaviside. L'opérateur (u * f)(t) est le produit de convolution de u(t) et ƒ(t). | ||
Dilatation échelle de temps | |||
Décalage sur p | Cette propriété est parfois connue sous le nom de Théorème de l'amortissement (ou Théorème de la modulation) avec . | ||
Décalage domaine temporel | u(t) est la fonction échelon de Heaviside (step function) | ||
Multiplication | L'intégration est effectuée le long de la ligne verticale Re(σ) = c qui est entièrement située dans le rayon de convergence de F[9]. | ||
Produit de convolution | ƒ(t) et g(t) sont étendues sur pour la définition du produit de convolution. | ||
Conjugaison complexe | |||
Fonction de corrélation | |||
Fonction périodique | ƒ(t) est une fonction périodique de période T telle que . Ceci résulte de la propriété de décalage dans le domaine temporel et de la série géométrique. |
Quelques transformées usuelles
La transformée de Laplace monolatérale n'est valide que pour des fonctions (éventuellement généralisées) à support positif. C'est pour cette raison que les fonctions temporelles de cette table sont multiples de (ou composées avec) , fonction échelon unité (Heaviside).
Fonction | Domaine temporel | Transformée de Laplace | Région de convergence | |
---|---|---|---|---|
1 | Distribution de Dirac retardée | |||
1a | Distribution de Dirac | |||
2 | exponentielle-monôme retardée | |||
2a | puissance n-ième | |||
2a.1 | puissance q-ième | |||
2a.2 | échelon unité | |||
2b | échelon retardé | |||
2c | rampe | |||
2d | exponentielle-monôme | |||
2d.1 | exponentielle | |||
3 | approche exponentielle | |||
4 | sinus | |||
5 | cosinus | |||
6 | sinus hyperbolique | |||
7 | cosinus hyperbolique | |||
8 | décroissance exponentielle d'une onde sinusoidale | |||
9 | décroissance exponentielle d'une onde cosinusoidale | |||
10 | n-ième racine | |||
11 | logarithme | |||
12 | fonction de Bessel du premier type, d'ordre n | | ||
13 | fonction de Bessel modifiée du premier type, d'ordre n | | ||
14 | fonction d'erreur | |||
Notes:
|
Exemple d'utilisation de la transformée de Laplace en électricité
On considère un circuit dit « R,C », constitué d'une résistance électrique de valeur R et d'un condensateur de capacité électrique C, placés en série. Dans tous les cas on considère que le circuit n'est placé aux bornes d'un générateur idéal de tension délivrant une tension (en général) variable u(t) qu'à un instant choisi pour origine des dates, et que le condensateur est initialement déchargé.
On a ainsi respectivement pour la charge q(t) du condensateur et l'intensité dans le circuit les conditions initiales suivantes :
Charge d'un condensateur par un échelon de tension
On applique la tension u(t) suivante :
et l'équation différentielle reliant la réponse q(t) à l'entrée u(t) est en appliquant les lois usuelles de l'électricité :
soit encore en posant τ ≡ RC (cette quantité a la dimension d'une durée) et en divisant par R :
On prend la transformée de Laplace membre à membre de cette dernière équation, en notant Q(p) la transformée de q(t), il vient, en prenant en compte le fait que q(0–) = 0 :
ce qui peut aussi s'écrire sous la forme :
- , avec la fonction de transfert du système RC, et la transformée de Laplace de l'entrée[10].
On peut aussitôt inverser cette équation en (on utilise l'entrée numéro 3 de la table ci-dessus avec α = 1/τ) :
L'interprétation physique de cette solution est très simple : il y a superposition d'un régime transitoire
qui décrit la charge progressive du condensateur, la quantité τ ≡ RC donnant l'échelle de temps (c'est un exemple de constante de temps d'un système), à un régime permanent
qui correspond à l'état du condensateur complètement chargé sous la tension continue U0. On montre aisément que le condensateur est à 90 % chargé (q = 0,90 Qm) au bout de la durée T = τ ln(10) ≈ 2,3025 τ.
Le terme (1–e–t/τ) est la fonction de transfert du système dans le domaine temporel.
On voit la facilité d'usage de la transformation de Laplace, qui permet de s'abstraire complètement de la résolution de l'équation différentielle dans l'espace des temps par un passage dans « l'espace p ». Par ailleurs, la prise en compte des conditions initiales est effectuée lors de la transformation.
Notes et références
Notes
- Bourlès 2010 (§12.3.4), Bourlès et Marinescu 2011, § 7.3.4.1.
- Denis-Papin et Kaufmann 1967.
- J.-É. Rombaldi, Exercices et problèmes corrigés pour l'agrégation de mathématiques, De Boeck Supérieur, (lire en ligne), p. 193.
- Bourlès 2010, p. 356.
- (en) Milton Abramowitz et Irene Stegun, Handbook of Mathematical Functions with Formulas, Graphs, and Mathematical Tables [détail de l’édition] (lire en ligne), chap. 29 (« Laplace Transforms »), p. 1020 : 29.2.4. et 29.2.5
- (en) Milton Abramowitz et Irene Stegun, Handbook of Mathematical Functions with Formulas, Graphs, and Mathematical Tables [détail de l’édition] (lire en ligne), chap. 29 (« Laplace Transforms »), p. 1020 : 29.1.1.
- Schwartz 1965, VI,2;2.
- André Desbiens, « Systèmes et commande linéaires GEL-2005. Chapitre 3 : Transformation de Laplace », sur Université Laval, p. 33.
- Bracewell 2000, Table 14.1, p. 385.
- En unité de charge de par la multiplication par C.
Références
- Henri Bourlès, Linear Systems, John Wiley & Sons, , 544 p. (ISBN 978-1-84821-162-9 et 1-84821-162-7, lire en ligne)
- Henri Bourlès et Bogdan Marinescu, Linear Time-Varying Systems: Algebraic-Analytic Approach, Springer, , 638 p. (ISBN 3642197264)
- (en) Ronald N. Bracewell, The Fourier Transform and Its Applications, Boston, McGraw-Hill, , 3e éd. (ISBN 0-07-116043-4).
- M. Denis-Papin et A. Kaufmann, Cours de calcul opérationnel appliqué, Albin Michel, (ASIN B003WR50TY)
- Laurent Schwartz, Méthodes mathématiques pour les sciences physiques, Hermann, (ISBN 2-7056-5213-2)
- (en) D. V. Widder, The Laplace Transform, Dover Publications, , 406 p. (ISBN 978-0-486-47755-8 et 0-486-47755-X)
Voir aussi
Articles connexes
- Transformée bilatérale de Laplace
- Transformée en Z
- Calcul opérationnel
- Transformation de Mellin
- Fonction de transfert
- Théorème de Bernstein sur les fonctions totalement monotones
- Portail de l'analyse