Transformation de Laplace

En mathématiques, la transformation de Laplace est une transformation intégrale qui, à une fonction f — définie sur les réels positifs et à valeurs réelles —, associe une nouvelle fonction F — définie sur les complexes et à valeurs complexes — dite transformée de Laplace de f.

L'intérêt de la transformation de Laplace vient de la conjonction des deux faits suivants :

1. De nombreuses opérations courantes sur la fonction originale f se traduisent par une opération algébrique sur la transformée F. Par exemple, la transformée de la dérivée $f'$ est la fonction $p\mapsto p\,F(p)-f(0).$ De même, la transformée de $t\mapsto f(t-\tau )$ est la fonction $p\mapsto e^{-p\tau }F(p)$ ;
2. Bien que n'étant pas à proprement parler injective, la transformation de Laplace est inversible — au sens qu'il est possible de retrouver (à un ensemble de mesure nulle près) la fonction originale f à partir de sa transformée F.

Ces deux points font de la transformée de Laplace un outil utile à la résolution d'équations différentielles : la méthode consiste à (1) transformer une équation différentielle sur la fonction f en une équation algébrique sur sa transformée F; (2) résoudre cette équation algébrique pour déterminer F; et (3) en déduire f par transformation de Laplace inverse. Dans ce type d'analyse, la transformation de Laplace est souvent interprétée comme un passage du domaine temporel (dans lequel les entrées et sorties sont des fonctions de la variable t représentant le temps) au domaine fréquentiel (la variable $p=Ae^{i\theta }$ pouvant s'interpréter comme une « fréquence complexe » d'amplitude $|p|=A$ et de phase $\mathrm {arg} (p)=\theta$).

La transformation de Laplace généralise la transformation de Fourier, é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.

La transformation de Laplace est également utile en théorie des probabilités, où elle peut être utilisée à la place de la fonction caractéristique pour caractériser la loi d'une variable aléatoire réelle. Cela s'avère particulièrement utile dans l'étude des subordinateurs.

La transformation de Laplace doit son nom à Pierre-Simon de Laplace, qui utilisa une notion proche en 1774 dans le cadre de ses travaux sur la théorie des probabilités.

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 :

$\mathrm {F} (p)={\mathcal {L}}\{f\}(p)=\int _{0^{-}}^{+\infty }{\rm {e}}^{-pt}f(t)\,\mathrm {d} t.$

Plus précisément, cette formule est valide lorsque :

1. Re(p) > α, où α est l'abscisse de convergence (définie plus bas), –∞ ≤ α ≤ +∞ ;
2. 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).

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 β, $f_{\beta }:t\mapsto {\rm {e}}^{-\beta t}f\left(t\right)$. Alors α est la borne inférieure dans ${\overline {\mathbb {R} }}$ 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 :

$e(t)=\mathrm {R} \cdot i(t)+\mathrm {L} {\frac {\mathrm {d} i(t)}{\mathrm {d} t}}$

dans le domaine fréquentiel devient

$\mathrm {E} (p)=\mathrm {R} \cdot \mathrm {I} (p)+p\cdot \mathrm {L} \cdot \mathrm {I} (p)$

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 :

$\phi (p)=p\int _{0^{-}}^{+\infty }{\rm {e}}^{-pt}f(t)\,\mathrm {d} t$

qui permet d'associer à toute fonction d'une variable $t\mapsto f(t)$ une fonction image $p\mapsto \phi (p)$.

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 :

$f(t)={\mathcal {L}}^{-1}\{\mathrm {F} \}(t)={\frac {1}{2\pi {\rm {i}}}}\int _{\gamma -{\rm {i}}\cdot \infty }^{\gamma +{\rm {i}}\cdot \infty }{\rm {e}}^{pt}\mathrm {F} (p)\,{\rm {d}}p,$

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 ${\dfrac {1}{\vert p\vert ^{2}}}$.

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 ${\dfrac {1}{\vert p\vert ^{2}}}$

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 :

${\mathcal {L}}\left\{af+bg\right\}=a\,{\mathcal {L}}\left\{f\right\}+b\,{\mathcal {L}}\left\{g\right\}$.

Cette linéarité découle évidemment de celle de l'intégrale.

Continuité

Si $f:\mathbb {R} _{+}\to \mathbb {C}$ est continue et si l'intégrale impropre $\int _{0}^{\infty }f$ converge, alors ${\mathcal {L}}f(x)$ est bien défini pour tout réel $x\geq 0$ et ${\mathcal {L}}f$ est continue sur $\mathbb {R} _{+}$. En particulier, $\int _{0}^{\infty }f=\lim _{0^{+}}{\mathcal {L}}f$.

En effet, la règle d'Abel s'applique ici uniformément par rapport à x.

Holomorphie

La transformée de Laplace $\mathrm {F} (p)={\mathcal {L}}\{f(t)\}$ de $f$ est holomorphe et sa dérivée n-ième est $\mathrm {F} ^{(n)}(p)=(-1)^{n}{\mathcal {L}}\{t^{n}f(t)\}$ (voir infra).

Transformation de Laplace d'une dérivée

Appliquée à la dérivée $f'$ de f, la transformation de Laplace correspond, à une constante additive près, à une multiplication par p de la transformée :

${\mathcal {L}}\{f'\}=p{\mathcal {L}}\{f\}-f(0^{-})$.

Démonstration

Soit à calculer :

${\mathcal {L}}\{f'\}=\int _{0^{-}}^{\infty }{\rm {e}}^{-pt}f'(t)\,{\rm {d}}t.$

En intégrant par parties, on obtient :

${\mathcal {L}}\{f'\}=\left[{\rm {e}}^{-pt}f(t)\right]_{0^{-}}^{\infty }+p\int _{0^{-}}^{\infty }{\rm {e}}^{-pt}f(t)\,{\rm {d}}t,$

soit finalement : ${\mathcal {L}}\{f'\}=p{\mathcal {L}}\{f\}-f(0^{-}).$

De proche en proche ou par récurrence, il est possible de montrer pour les dérivations successives :

${\mathcal {L}}\{f'\}=p{\mathcal {L}}\{f\}-f(0^{-})$ ;

${\mathcal {L}}\{f''\}=p^{2}{\mathcal {L}}\{f\}-pf(0^{-})-f'(0^{-})$ ;

${\mathcal {L}}\left\{f^{(n)}\right\}=p^{n}{\mathcal {L}}\{f\}-p^{n-1}f(0^{-})-\cdots -f^{(n-1)}(0^{-})$.

Cette dernière expression peut s'écrire, avec $\partial _{0}^{i}f\left(0^{-}\right):=f^{\left(i\right)}\left(0^{-}\right)$ pour tout entier $i\geq 0$,

${\mathcal {L}}\left(f^{\left(n\right)}\right)=p^{n}{\mathcal {L}}\left(f\right)-{\frac {p^{n}-\partial _{0}^{n}}{p-\partial _{0}}}f\left(0^{-}\right).$

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 $f(0^{-}),...,f^{(n-1)}(0^{-})$ 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 $f=g\Upsilon$ où $\Upsilon$ 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,

$f'=g'\Upsilon +g\Upsilon '$ avec $\Upsilon '=\delta .$

Puisque $g\delta =g(0)\delta$, $f'=g'+g(0)\delta$, donc ${\mathcal {L}}\{f'\}={\mathcal {L}}\{g'\}+g(0)$.

On a également ${\mathcal {L}}\{f'\}=p{\mathcal {L}}\{f\}-f(0^{-})=p{\mathcal {L}}\{f\}$ car $f(0^{-})=0$.

Maintenant, ${\mathcal {L}}\{f\}={\mathcal {L}}\{g\Upsilon \}$ et $g(0)=(g\Upsilon )(0^{+})$. Par définition, ${\mathcal {L}}\{g'\}={\mathcal {L}}\{g'\Upsilon \}$ car c'est de la transformation monolatérale qu'il s'agit. On obtient donc finalement

${\mathcal {L}}\{g'\Upsilon \}=p{\mathcal {L}}\{g\Upsilon \}-(g\Upsilon )(0^{+}).$

En continuant ce raisonnement, on obtient, si g est de classe $C^{n}$ dans un voisinage de [0, +∞[,

${\mathcal {L}}\left(g^{\left(n\right)}\Upsilon \right)=p^{n}{\mathcal {L}}\left(g\Upsilon \right)-{\frac {p^{n}-\partial _{0}^{n}}{p-\partial _{0}}}(g\Upsilon )\left(0^{+}\right)$

avec $\partial _{0}^{i}(g\Upsilon )\left(0^{+}\right):=(g^{\left(i\right)}\Upsilon )\left(0^{+}\right)$ pour tout entier $i\geq 0$.

Exemple

Soit $g(t)=\cos(\omega t)$. Alors ${\mathcal {L}}\{g\Upsilon \}(p)={\frac {p}{p^{2}+\omega ^{2}}}$ et $g\Upsilon (0^{+})=1$. On a $g'(t)=-\omega \sin(\omega t)$ et

$p{\mathcal {L}}\{g\Upsilon \}(p)-g\Upsilon (0^{+})={\frac {p^{2}}{p^{2}+\omega ^{2}}}-1=-{\frac {\omega ^{2}}{p^{2}+\omega ^{2}}}$. Par conséquent,

${\mathcal {L}}\{t\mapsto \sin(\omega t)\Upsilon (t)\}(p)={\frac {\omega }{p^{2}+\omega ^{2}}}.$

Application à la dérivée de la fonction de Heaviside

La fonction de Heaviside $\Upsilon$ 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 $\delta$. Il vient

${\mathcal {L}}\left(\delta \right)=p{\mathcal {L}}\left(\Upsilon \right)-\Upsilon \left(0^{-}\right)=1-0=1,$

puisque

${\mathcal {L}}\left(\Upsilon \right)={\dfrac {1}{p}},\Re (p)>0.$

Si l'on remplaçait, dans la formule de la règle de dérivation, ƒ(0^–) par ƒ(0⁺), on trouverait ${\mathcal {L}}\left(\delta \right)=0$, ce qui est faux (voir plus loin). Certaines sources peuvent comporter cette erreur.

De même, on voit parfois la définition suivante de la transformation de Laplace :

$F\left(p\right)=\int _{\alpha }^{+\infty }{\rm {e}}^{-pt}f\left(t\right)~{\rm {d}}t$

avec $\alpha =0^{+}$, voire un manque de précision sur cette limite. 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 à $\alpha =0^{-}$, puisque $\left\{0\right\}$ est de mesure nulle ; on peut d'ailleurs dans ce cas écrire sans ambiguïté $\alpha =0$. 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, quand celle-ci a une masse non nulle à l'origine. Le prototype est la distribution de Dirac. Au plan algébrique, cette distribution $\delta$ est l'élément neutre dans l'algèbre de convolution ${\mathcal {D}}_{+}^{\prime }$ des distributions à support positif ; et puisque la transformation de Laplace transforme le produit de convolution en produit ordinaire, il faut donc que $\delta$ ait pour transformée de Laplace ${\mathcal {L}}\left(\delta \right)=1$. Or, cela ne sera vrai que si $\alpha =0^{-}$. En effet, avec $\alpha =0^{+}$ 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 $\delta \neq 0$.

Multiplication par une puissance de t

La multiplication par $t^{n},n\in \mathbb {N}$ dans le domaine temporel correspond, au signe près, à la dérivée n-ième de la transformée :

${\mathcal {L}}\left\{t^{n}f\left(t\right)\right\}=\left(-1\right)^{n}{\dfrac {\mathrm {d} ^{n}{\mathcal {L}}\left\{f\right\}}{\mathrm {d} p^{n}}}$.

Démonstration

(1) Supposons f localement intégrable à support positif. La transformée de Laplace de f est donc définie pour $\Re (p)>\alpha$, où $\alpha$ est l'abscisse de convergence, par

${\mathcal {L}}\left\{f\right\}(p)=\int _{0}^{+\infty }f(t){\rm {e}}^{-pt}~{\rm {d}}t$.

La fonction $p\mapsto f(t){\rm {e}}^{-pt}$ est holomorphe. Soit $\beta >\alpha$ et $\Re (p)>\beta$. Alors $\left\vert tf\left(t\right){\rm {e}}^{-pt}\right\vert \leq \left\vert tf\left(t\right){\rm {e}}^{-\beta t}\right\vert$ et par croissances comparées, la fonction $t\mapsto \left\vert tf\left(t\right){\rm {e}}^{-\beta t}\right\vert$ est intégrable sur [0, +∞[. La fonction $p\mapsto \int _{0}^{+\infty }f(t){\rm {e}}^{-pt}~{\rm {d}}t$ est donc holomorphe, et sa dérivée s'obtient en dérivant sous le signe somme :

${\frac {{\rm {d}}{\mathcal {L}}\left\{f\right\}}{{\rm {d}}p}}\left(p\right)=\int _{0}^{+\infty }f\left(t\right)\left(-t{\rm {e}}^{-pt}\right)~{\rm {d}}t=-{\mathcal {L}}\left\{tf\left(t\right)\right\}\left(p\right)$.

Ceci prouve le résultat dans le cas n = 1. Le cas général s'ensuit, par récurrence.

(2) Ce résultat est encore valide lorsque f est une distribution à support positif.

La formule inverse (pour n = –1) est :

${\mathcal {L}}\left\{{\frac {f(t)}{t}}\right\}\left(p\right)=\int _{p}^{\infty }\mathrm {F} (\sigma )\,\mathrm {d} \sigma$

et elle est valide à condition que f soit de la forme $t\mapsto tg(t)$ 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.

Démonstration

$\int _{p}^{\infty }\mathrm {F} (\sigma )\,\mathrm {d} \sigma =\int _{p}^{\infty }\int _{0}^{\infty }\mathrm {e} ^{-\sigma t}f(t)\,\mathrm {d} t\,\mathrm {d} \sigma =\int _{0}^{\infty }f(t)\int _{p}^{\infty }\mathrm {e} ^{-\sigma t}\,\mathrm {d} \sigma \,\mathrm {d} t=\int _{0}^{\infty }f(t){\frac {1}{t}}\mathrm {e} ^{-pt}\,\mathrm {d} t={\mathcal {L}}\left\{{\frac {f(t)}{t}}\right\}\left(p\right)$.

Intégration

La transformation de Laplace d'une intégrale (primitive de f s'annulant en 0) correspond à une multiplication par 1/p :

${\mathcal {L}}\left\{\int _{0^{-}}^{t}f(\tau )\mathrm {d} \tau \right\}={\frac {1}{p}}{\mathcal {L}}\{f\}$

et si ƒ est une fonction à support positif, continue sur [0, +∞[, on a pour tout a > 0 :

${\mathcal {L}}\left\{\int _{a}^{t}f(\tau )~{\rm {d}}\tau \right\}={\frac {1}{p}}{\mathcal {L}}\{f\}+{1 \over p}\int _{a}^{0}f(\tau )\mathrm {d} \tau .$

Valeur finale

Supposons f localement intégrable à support positif. Si la limite dans le domaine temporel existe et est finie, alors :

$\lim _{t\to +\infty }f(t)=\lim _{p\in \mathbb {R} ,p\to 0^{+}}p\mathrm {F} (p).$

(On notera que c'est la seule propriété où un 0⁺ apparaît pour la variable $p$.)

Démonstration

Soit $l=\lim \limits _{t\rightarrow +\infty }f\left(t\right)$. L'existence de cette limite finie implique que l'abscisse de convergence de la transformée de Laplace $F(p)$ est $\leq 0$.

On a $\lim _{t\to +\infty }\Upsilon (t)=1$ ; la transformée de Laplace de $\Upsilon$ est ${\frac {1}{p}}$, et $\lim _{p\in \mathbb {R} ,p\to 0^{+}}p{\frac {1}{p}}=1$. En soustrayant $l\Upsilon (t)$ de $f(t)$, on est donc ramené au cas d'une fonction, de nouveau notée f, telle que $\lim \limits _{t\rightarrow +\infty }f\left(t\right)=0$.

Alors, pour tout $\varepsilon >0$, il existe $A>0$ tel que pour tout $t>A$, $\vert f(t)\vert \leq \varepsilon$. On a

$pF(p)=p\int _{0}^{A}f(t){\rm {e}}^{-pt}{\rm {d}}t+p\int _{A}^{+\infty }f(t){\rm {e}}^{-pt}~{\rm {d}}t.$

Prenons $p\in \mathbb {R} ,p>0$. On a

$\left\vert \int _{0}^{A}f(t){\rm {e}}^{-pt}~{\rm {d}}t\right\vert \leq \int _{0}^{A}\vert f(t)\vert ~{\rm {d}}t<+\infty$

et par suite

$\lim \limits _{p\in \mathbb {R} ,p\rightarrow 0^{+}}p\int _{0}^{A}f(t){\rm {e}}^{-pt}~{\rm {d}}t=0.$

Par conséquent, il existe un réel $\alpha >0$ tel que pour $p\in \mathbb {R}$ et $0<p<\alpha$

$\left\vert p\int _{0}^{A}f(t){\rm {e}}^{-pt}~{\rm {d}}t\right\vert <\varepsilon .$

D'autre part,

$\left\vert \int _{A}^{+\infty }f(t){\rm {e}}^{-pt}~{\rm {d}}t\right\vert \leq \int _{0}^{+\infty }\varepsilon {\rm {e}}^{-pt}~{\rm {d}}t={\frac {\varepsilon }{p}}$

donc il existe $\beta >0$ tel que pour $p\in \mathbb {R}$ et $0<p<\beta$

$\left\vert p\int _{A}^{+\infty }f(t){\rm {e}}^{-pt}~{\rm {d}}t\right\vert \leq \varepsilon .$

Par conséquent, si $p\in \mathbb {R}$ et $0<p<\min(\alpha ,\beta )$

$\left\vert p\int _{0}^{+\infty }f(t){\rm {e}}^{-pt}~{\rm {d}}t\right\vert \leq 2\varepsilon$

ce qui entraîne que $p\int _{0}^{+\infty }f(t){\rm {e}}^{-pt}~{\rm {d}}t\rightarrow 0$ lorsque $p\in \mathbb {R}$ tend vers 0⁺.

Les hypothèses indiquées sont indispensables, comme le montrent les contre-exemples suivants :

• la fonction $f:t\mapsto {\rm {e}}^{t}\Upsilon (t)$ admet pour limite +∞ lorsque t tend vers +∞. Sa transformée de Laplace est $F\left(p\right)={\frac {1}{p-1}}$ et $\lim \limits _{p\rightarrow 0}pF\left(p\right)=0$. 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 $f:t\mapsto \sin t\Upsilon (t)$ n'admet pas de limite lorsque t tend vers +∞. Sa transformée de Laplace est $F\left(p\right)={\frac {1}{p^{2}+1}}$, l'abscisse de convergence de F est 0 et $\lim \limits _{p\in \mathbb {R} ,p\rightarrow 0^{+}}pF\left(p\right)=0$ (cette dernière limite est pourtant correcte cette fois) ;
• si $F(p)$ est une fonction rationnelle, $\lim _{t\to +\infty }f(t)$ existe et est finie si et seulement si les pôles de $F(p)$ appartiennent tous à la réunion du demi-plan gauche ouvert et de l'origine, le pôle en 0, s'il existe, étant simple.

Valeur initiale

Si $F(p)$ a une abscisse de convergence finie et si la limite dans le domaine temporel existe, alors :

$\lim _{t\to 0^{+}}f(t)=\lim _{p\in \mathbb {R} ,p\to +\infty }p\mathrm {F} (p)$.

(On notera que c'est la seule propriété où un 0⁺ apparaît pour la variable $t$.)

Démonstration

Soit $l=\lim \limits _{t\rightarrow 0^{+}}f\left(t\right)$. On a $\lim _{t\to 0^{+}}\Upsilon (t)=1$ ; la transformée de Laplace de $\Upsilon$ est ${\frac {1}{p}}$, et $\lim _{p\in \mathbb {R} ,p\to +\infty }p{\frac {1}{p}}=1$. En soustrayant $l\Upsilon (t)$ de $f(t)$, on est donc ramené au cas d'une fonction, de nouveau notée f, telle que $\lim \limits _{t\rightarrow 0^{+}}f\left(t\right)=0$.

Soit $\varepsilon >0$. Il existe par hypothèse $\eta >0$ tel que pour tout t tel que $0<t<\eta$, on a $\vert f(t)\vert \leq \varepsilon$. D'autre part,

$p\int _{0}^{+\infty }f(t){\rm {e}}^{-pt}~{\rm {d}}t=I_{1}+I_{2}$

avec

$I_{1}=p\int _{0}^{\eta }f(t){\rm {e}}^{-pt}~{\rm {d}}t,~I_{2}=p\int _{\eta }^{+\infty }f(t){\rm {e}}^{-pt}~{\rm {d}}t.$

Soit $\alpha$ un réel strictement supérieur à l'abscisse de convergence de $F(p)$ et $p\in \mathbb {R} ,~p>\alpha$. On a

$\left\vert I_{2}\right\vert =\left\vert p\int _{\eta }^{+\infty }f(t){\rm {e}}^{-\left(p-\alpha \right)t}{\rm {e}}^{-\alpha t}~{\rm {d}}t\right\vert \leq p{\rm {e}}^{-\left(p-\alpha \right)\eta }\int _{0}^{+\infty }\left\vert f\left(t\right)\right\vert {\rm {e}}^{-\alpha t}~{\rm {d}}t$

où l'intégrale de droite est convergente, donc $I_{2}\rightarrow 0$ lorsque $p\rightarrow +\infty$. Par conséquent, il existe un réel $A>0$ tel que $\left\vert I_{2}\right\vert \leq \varepsilon$ dès que $p\in \mathbb {R}$ et $p>A$.

D'autre part,

$\left\vert I_{1}\right\vert \leq p\varepsilon \int _{0}^{\eta }{\rm {e}}^{-pt}~{\rm {d}}t=\varepsilon \left(1-{\rm {e}}^{-p\eta }\right)$

et ce terme tend vers $\varepsilon$ lorsque $p\rightarrow +\infty$, donc il existe un réel $B>0$ tel que $\left\vert I_{1}\right\vert \leq 2\varepsilon$ dès que $p\in \mathbb {R}$ et $p>B$. Finalement, pour $p\in \mathbb {R}$ et $p>\max(A,B)$ on a

$\left\vert pF\left(p\right)\right\vert \leq 3\varepsilon .$

Or, $\varepsilon >0$ est arbitrairement petit, donc ce terme tend vers 0 lorsque $p\in \mathbb {R}$ et $p\rightarrow +\infty$.

Convolution

La transformation de Laplace change le produit de convolution en produit :

${\mathcal {L}}\{f*g\}={\mathcal {L}}\{f\}{\mathcal {L}}\{g\}$.

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 $Re\left(p\right)>0$

${\mathcal {L}}\{f\}(p)={\frac {1}{1-{\rm {e}}^{-Tp}}}\int _{0}^{T}{\rm {e}}^{-pt}f(t)\,\mathrm {d} t.$

Démonstration

On utilise la relation de Chasles pour décomposer l'intégrale sur chaque période :

$\int _{0}^{\infty }{\rm {e}}^{-pt}f(t)\,\mathrm {d} t=\int _{0}^{T}{\rm {e}}^{-pt}f(t)\,\mathrm {d} t+\int _{T}^{2T}{\rm {e}}^{-pt}f(t)\,\mathrm {d} t+\int _{2T}^{3T}{\rm {e}}^{-pt}f(t)\,\mathrm {d} t+\ldots$

On fait un changement de variables pour ramener les intégrales sur [0,T]

$\int _{0}^{\infty }{\rm {e}}^{-pt}f(t)\,{\rm {d}}t=\int _{0}^{T}{\rm {e}}^{-pu}f(u)\,\mathrm {d} u+\int _{0}^{T}{\rm {e}}^{-p(u+T)}f(u+T)\,\mathrm {d} u+\int _{0}^{T}{\rm {e}}^{-p(u+2T)}f(u+2T)\,\mathrm {d} u+\ldots$

Comme ƒ est périodique, on peut simplifier les intégrales en

$\int _{0}^{\infty }{\rm {e}}^{-pt}f(t)\,\mathrm {d} t=\int _{0}^{T}{\rm {e}}^{-pu}f(u)\,\mathrm {d} u+{\rm {e}}^{-pT}\int _{0}^{T}{\rm {e}}^{-pu}f(u)\,\mathrm {d} u+{\rm {e}}^{-2pT}\int _{0}^{T}{\rm {e}}^{-pu}f(u)\,\mathrm {d} u+\ldots$

On regroupe les termes :

$\int _{0}^{\infty }{\rm {e}}^{-pt}f(t)\,\mathrm {d} t=\left(1+{\rm {e}}^{-pT}+{\rm {e}}^{-2pT}+\ldots \right)\int _{0}^{T}{\rm {e}}^{-pu}f(u)\,\mathrm {d} u.$

Cette série géométrique converge (car e^–pT < 1). Il vient alors

${\mathcal {L}}\{f\}(p)={\frac {1}{1-{\rm {e}}^{-Tp}}}\int _{0}^{T}{\rm {e}}^{-pt}f(t)\,\mathrm {d} t.$

Développement de Heaviside

Si la valeur finale est une fraction rationnelle F(p)/G(p) telle que le degré de F est inférieur au degré de G, alors pour retrouver la valeur initiale, on décompose en éléments simples et on applique la linéarité.

Supposons que

$G(p)=(p-a_{1})^{k_{1}}\cdots (p-a_{n})^{k_{n}}$

Alors on peut écrire

${\frac {F(p)}{G(p)}}={\frac {c_{1,1}}{(p-a_{1})}}+\cdots +{\frac {c_{1,k_{1}}}{(p-a_{1})^{k_{1}}}}+\cdots +{\frac {c_{n,k_{n}}}{(p-a_{n})^{k_{n}}}}$

et alors

${\begin{aligned}{\mathcal {L}}^{-1}\left({\frac {F(p)}{G(p)}}\right)(t)&={\mathcal {L}}^{-1}\left({\frac {c_{1,1}}{(p-a_{1})}}\right)(t)+\cdots +{\mathcal {L}}^{-1}\left({\frac {c_{1,k_{1}}}{(p-a_{1})^{k_{1}}}}\right)(t)+\cdots +{\mathcal {L}}^{-1}\left({\frac {c_{n,k_{n}}}{(p-a_{n})^{k_{n}}}}\right)(t)\\&=c_{1,1}\mathrm {e} ^{a_{1}t}+\cdots +c_{1,k_{1}}{\frac {t^{k_{1}-1}}{(k_{1}-1)!}}\mathrm {e} ^{a_{1}t}+\cdots +c_{n,k_{n}}{\frac {t^{k_{n}-1}}{(k_{n}-1)!}}\mathrm {e} ^{a_{n}t}\end{aligned}}$

Tableau résumé des propriétés de la transformation de Laplace

Propriétés de la transformée de Laplace unilatérale

	Domaine temporel	Domaine "p"	Commentaires
Linéarité	$af(t)+bg(t)$	$a\mathrm {F} (p)+b\mathrm {G} (p)$	Résulte des règles de base de l'intégration.
Dérivée de la transformée	$tf(t)$	$-\mathrm {F} '(p)$	$\mathrm {F} '$ est la dérivée première de F.
Dérivées d'ordre n de la transformée	$t^{n}f(t)$	$(-1)^{n}\mathrm {F} ^{(n)}(p)$	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	$f'(t)$	$p\mathrm {F} (p)-f\left(0^{-}\right)$	ƒ 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	$f''(t)$	$p^{2}\mathrm {F} (p)-pf\left(0^{-}\right)-f'\left(0^{-}\right)$	ƒ est supposée deux fois dérivable, avec dérivée seconde convergeant exponentiellement à l'infini.
Dérivée n-ième de ƒ	$f^{(n)}(t)$	$p^{n}\mathrm {F} (p)-p^{n-1}f\left(0^{-}\right)-\cdots -f^{(n-1)}\left(0^{-}\right)$	ƒ 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	${\frac {f(t)}{t}}$	$\int _{p}^{\infty }\mathrm {F} (\sigma )\,\mathrm {d} \sigma$
Intégration	$\int _{0}^{t}f(\tau )\,\mathrm {d} \tau =(u*f)(t)$	${1 \over p}\mathrm {F} (p)$	$u(t)$ 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	$f(at)\$	${\frac {1}{|a|}}\mathrm {F} \left({p \over a}\right)$
Décalage sur p	${\rm {e}}^{at}f(t)$	$\mathrm {F} (p-a)$	Cette propriété est parfois connue sous le nom de Théorème de l'amortissement (ou Théorème de la modulation) avec $a<0$.
Décalage domaine temporel	$f(t-a)u(t-a)$	${\rm {e}}^{-ap}\mathrm {F} (p)$	u(t) est la fonction échelon de Heaviside (step function).
Multiplication	$f(t)g(t)$	${\frac {1}{2\pi {\rm {i}}}}\lim _{\mathrm {T} \to \infty }\int _{c-{\rm {i}}\mathrm {T} }^{c+{\rm {i}}\mathrm {T} }\mathrm {F} (\sigma )G(p-\sigma )\,\mathrm {d} \sigma$	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.
Produit de convolution	$(f*g)(t)=\int _{0}^{t}f(\tau )g(t-\tau )\,{\rm {d}}\tau$	$\mathrm {F} (p)\cdot \mathrm {G} (p)$	ƒ(t) et g(t) sont étendues sur $\mathbb {R}$ pour la définition du produit de convolution.
Conjugaison complexe	$f^{*}(t)$	$\mathrm {F} ^{*}(p^{*})$
Fonction de corrélation	$f(t)\star g(t)$	$\mathrm {F} ^{*}(-p^{*})\cdot \mathrm {G} (p)$
Fonction périodique	$f(t)$	${\frac {1}{1-{\rm {e}}^{-\mathrm {T} p}}}\int _{0}^{\mathrm {T} }{\rm {e}}^{-pt}f(t)\,\mathrm {d} t$	ƒ(t) est une fonction périodique de période T telle que $f(t)=f(t+\mathrm {T} ),\;\forall t\geq 0$. 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) $\Upsilon$, fonction échelon unité (Heaviside).

Table des transformées de Laplace usuelles

	Fonction	Domaine temporel $x(t)={\mathcal {L}}^{-1}\left\{\mathrm {X} (p)\right\}$	Transformée de Laplace $\mathrm {X} (p)={\mathcal {L}}\left\{x(t)\right\}$	Région de convergence
1	Distribution de Dirac retardée	$\delta (t-\tau )\$	${\rm {e}}^{-\tau p}\$	$\forall \ p\,$
1a	Distribution de Dirac	$\delta (t)\$	$1\$	$\forall \ p\,$
2	exponentielle-monôme retardée	${\frac {(t-\tau )^{n}}{n!}}{\rm {e}}^{-\alpha (t-\tau )}\cdot \Upsilon (t-\tau )$	${\frac {{\rm {e}}^{-\tau p}}{(p+\alpha )^{n+1}}}$	$\operatorname {Re} (p)>0\,$
2a	puissance n-ième	${t^{n} \over n!}\cdot \Upsilon (t)$	${1 \over p^{n+1}}$	$\operatorname {Re} (p)>0\,$
2a.1	puissance q-ième	${t^{q} \over \Gamma (q+1)}\cdot \Upsilon (t)$	${1 \over p^{q+1}}$	$\operatorname {Re} (p)>0\,$
2a.2	échelon unité	$\Upsilon (t)\$	${1 \over p}$	$\operatorname {Re} (p)>0\,$
2b	échelon retardé	$\Upsilon (t-\tau )\$	${{\rm {e}}^{-\tau p} \over p}$	$\operatorname {Re} (p)>0\,$
2c	rampe	$t\cdot \Upsilon (t)\$	${\frac {1}{p^{2}}}$	$\operatorname {Re} (p)>0\,$
2d	exponentielle-monôme	${\frac {t^{n}}{n!}}{\rm {e}}^{-\alpha t}\cdot \Upsilon (t)$	${\frac {1}{(p+\alpha )^{n+1}}}$	$\operatorname {Re} (p)>-\alpha \,$
2d.1	exponentielle	${\rm {e}}^{-\alpha t}\cdot \Upsilon (t)\$	${1 \over p+\alpha }$	$\operatorname {Re} (p)>-\alpha \$
3	approche exponentielle	$(1-{\rm {e}}^{-\alpha t})\cdot \Upsilon (t)\$	${\frac {\alpha }{p(p+\alpha )}}$	$\operatorname {Re} (p)>0\$
4	sinus	$\sin(\omega t)\cdot \Upsilon (t)\$	${\omega \over p^{2}+\omega ^{2}}$	$\operatorname {Re} (p)>0\$
5	cosinus	$\cos(\omega t)\cdot \Upsilon (t)\$	${p \over p^{2}+\omega ^{2}}$	$\operatorname {Re} (p)>0\$
6	sinus hyperbolique	$\sinh(\alpha t)\cdot \Upsilon (t)\$	${\alpha \over p^{2}-\alpha ^{2}}$	$\operatorname {Re} (p)>|\alpha |\$
7	cosinus hyperbolique	$\cosh(\alpha t)\cdot \Upsilon (t)\$	${p \over p^{2}-\alpha ^{2}}$	$\operatorname {Re} (p)>|\alpha |\$
8	décroissance exponentielle d'une onde sinusoidale	${\rm {e}}^{-\alpha t}\sin(\omega t)\cdot \Upsilon (t)\$	${\omega \over (p+\alpha )^{2}+\omega ^{2}}$	$\operatorname {Re} (p)>-\alpha \$
9	décroissance exponentielle d'une onde cosinusoidale	${\rm {e}}^{-\alpha t}\cos(\omega t)\cdot \Upsilon (t)\$	${p+\alpha \over (p+\alpha )^{2}+\omega ^{2}}$	$\operatorname {Re} (p)>-\alpha \$
10	n-ième racine	${\sqrt[{n}]{t}}\cdot \Upsilon (t)$	$p^{-(n+1)/n}\cdot \Gamma \left(1+{\frac {1}{n}}\right)$	$\operatorname {Re} (p)>0\,$
11	logarithme	$\ln \left({t \over t_{0}}\right)\cdot \Upsilon (t)$	$-{t_{0} \over p}\ [\ \ln(t_{0}p)+\gamma \ ]$	$\operatorname {Re} (p)>0\,$
12	fonction de Bessel du premier type, d'ordre n	$J_{n}(\omega t)\cdot \Upsilon (t)$	${\frac {\omega ^{n}\left(p+{\sqrt {p^{2}+\omega ^{2}}}\right)^{-n}}{\sqrt {p^{2}+\omega ^{2}}}}$	$\operatorname {Re} (p)>0\,$ $(n>-1)\,$
13	fonction de Bessel modifiée du premier type, d'ordre n	$\mathrm {I} _{n}(\omega t)\cdot \Upsilon (t)$	${\frac {\omega ^{n}\left(p+{\sqrt {p^{2}-\omega ^{2}}}\right)^{-n}}{\sqrt {p^{2}-\omega ^{2}}}}$	$\operatorname {Re} (p)>|\omega |\,$ $(n>-1)\,$
14	fonction d'erreur	$\mathrm {erf} (t)\cdot \Upsilon (t)$	${{\rm {e}}^{p^{2}/4}\operatorname {erfc} \left(p/2\right) \over p}$	$\operatorname {Re} (p)>0\,$
Notes: $\Upsilon (t)\,$ représente la fonction de Heaviside. $\delta (t)\,$ représente la fonction de Dirac. $\Gamma (z)\,$ est la fonction Gamma. $\gamma \,$ est la constante d'Euler-Mascheroni. $t\,$, est un nombre réel, il représente typiquement le temps, mais peut désigner n'importe quelle autre quantité. $p\,$ est un nombre complexe. $q\,$ est un nombre réel ($q+1>0$). $\alpha \,$, $\beta \,$, $\tau \,$, et $\omega \,$ sont des nombres réels. $n\,$ est un entier.

 

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 $i\left(t\right)\equiv {\frac {{\rm {d}}q}{{\rm {d}}t}}$ les conditions initiales suivantes :

$q\left(0^{-}\right)=0,i\left(0^{-}\right)=0.$

Charge d'un condensateur par un échelon de tension

On applique la tension u(t) suivante :

$u(t)={\begin{cases}0,{\text{ si }}t<0\\\mathrm {U} _{0}=cte,{\text{ si }}t\geq 0\end{cases}},$

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é :

$\mathrm {U} _{0}\Upsilon (t)=\mathrm {R} {\frac {\mathrm {d} q}{\mathrm {d} t}}+{\frac {q(t)}{\mathrm {C} }},$

soit encore en posant τ ≡ RC (cette quantité a la dimension d'une durée) et en divisant par R :

${\frac {\mathrm {C} \mathrm {U} _{0}}{\tau }}\Upsilon (t)={\frac {q(t)}{\tau }}+{\frac {\mathrm {d} q}{\mathrm {d} t}}.$

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 :

$\mathrm {Q} (p)=\mathrm {C} \mathrm {U} _{0}{\frac {\frac {1}{\tau }}{p\left({\frac {1}{\tau }}+p\right)}},$

ce qui peut aussi s'écrire sous la forme :

$\mathrm {Q} (p)=\mathrm {H} (p)\mathrm {U} (p)$, avec $\mathrm {H} (p)\equiv {\frac {\tau ^{-1}}{\tau ^{-1}+p}},$ la fonction de transfert du système RC, et $\mathrm {U} (p)={\frac {\mathrm {U} _{0}\mathrm {C} }{p}}$ la transformée de Laplace de l'entrée.

On peut aussitôt inverser cette équation en (on utilise l'entrée numéro 3 de la table ci-dessus avec α = 1/τ) :

$q(t)=\mathrm {U} _{0}\mathrm {C} \left(1-{\rm {e}}^{-t/\tau }\right)\Upsilon (t).$

L'interprétation physique de cette solution est très simple : il y a superposition d'un régime transitoire

$q_{\mathrm {trans} }\left(t\right)=-\mathrm {U} _{0}\mathrm {C} \,{\rm {e}}^{-t/\tau },$

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

$\mathrm {Q_{perm}} =\mathrm {U} _{0}\mathrm {C} \equiv \mathrm {Q_{m}}$

qui correspond à l'état du condensateur complètement chargé sous la tension continue U₀. On montre aisément que le condensateur est à 90 % chargé (q = 0,90 Q_m) 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

Références

• Henri Bourlès, Linear Systems, John Wiley & Sons, 2010, 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, 2011, 638 p. (ISBN 3642197264)
• (en) Ronald N. Bracewell, The Fourier Transform and Its Applications, Boston, McGraw-Hill, 2000, 3^e éd. (ISBN 0-07-116043-4).
• M. Denis-Papin et A. Kaufmann, Cours de calcul opérationnel appliqué, Albin Michel, 1967 (ASIN B003WR50TY)
• Laurent Schwartz, Méthodes mathématiques pour les sciences physiques, Hermann, 1965 (ISBN 2-7056-5213-2)
• (en) D. V. Widder, The Laplace Transform, Dover Publications, 2011, 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