title: Calcul opérationnel
text: # Calcul opérationnel
  
  En mathématiques et plus précisément en analyse fonctionnelle, le calcul opérationnel repose essentiellement sur un astucieux changement de variable basé sur la transformée de Laplace permettant l'algébrisation des symboles de dérivation et d'intégration des expressions mathématiques décrivant les phénomènes linéaires. Certains ingénieurs emploient de préférence la transformation de « Laplace-Carson », une constante ayant comme image la même constante.
  
  L'expression : $\phi (p)=p\int _{0}^{+\infty }\mathrm {e} ^{-pt}f(t)\,\mathrm {d} t$ permet d'associer à toute fonction d'une variable $t\mapsto f(t)$ dite « fonction origine » une « fonction image » $p\mapsto \phi (p)$. Ainsi la solution algébrique de l'équation image permet de retrouver, au moyen d'un tableau de correspondance opératoire, la solution de l'équation origine.
  
  La transformation directe est notée :
  
      $\phi (p)$ image de $f(t)$.
  
  La transformation inverse est notée :
  
      $f(t)$ original de $\phi (p)$.
  
  ## Transformations de base
  
  ### Pour une constante « C »
  
  La correspondance entre fonctions originales et fonctions images s'établit comme suit :
  
  $C$
  
  $C$
  
  $Cf(t)$
  
  $C\phi (p)$
  
  $f_{1}(t)+f_{2}(t)$
  
  $\phi _{1}(p)+\phi _{2}(p)$
  
  $C_{1}f_{1}(t)+C_{2}f_{2}(t)$
  
  $C_{1}\phi _{1}(p)+C_{2}\phi _{2}(p)$
  
  ### Image d'une variable « t »
  
      $p\int _{0}^{+\infty }\mathrm {e} ^{-pt}t\,\mathrm {d} t=-\int _{0}^{+\infty }t\,\mathrm {d} (\mathrm {e} ^{-pt})=-\left[t\mathrm {e} ^{-pt}\right]_{0}^{+\infty }+\int _{0}^{+\infty }\mathrm {e} ^{-pt}\,\mathrm {d} t.$
  
  Pour $p>0$, on obtient l'image ${\frac {1}{p}}$.
  
  Ainsi,
  
      $t$ est l'original de ${\frac {1}{p}}$,
      $Ct$ est l'original de ${\frac {C}{p}}$,
      $Ct+C_{1}$ est l'original de ${\frac {C}{p}}+C_{1}$.
  
  D'une manière générale, par récurrence pour tout entier positif n, on obtient : $t^{n}$ original de ${\frac {n!}{p^{n}}}$.
  
  ### Image de l'exponentielle de « at »
  
      $p\int _{0}^{+\infty }\mathrm {e} ^{-pt}\mathrm {e} ^{at}\,\mathrm {d} t={\frac {p}{a-p}}\left[\mathrm {e} ^{(a-p)t}\right]_{0}^{+\infty }.$
  
  Si $a=\alpha +\mathrm {i} \beta$, la parenthèse devient :
  
      $\left[\mathrm {e} ^{(\alpha -p)t}(\cos \beta t+\mathrm {i} \sin \beta t)\right]_{0}^{+\infty }$, expression qui tend vers $-1$ lorsque $p>\alpha$ ; dans ce cas l'image de $\mathrm {e} ^{at}$ est ${\frac {p}{p-a}}$.
  
  Pour $a$ réel, le tableau de correspondance opératoire s'établit comme suit :
  
  | Fonction origine | Fonction image | Fonction origine | Fonction image | Condition(s) |
  | - | - | - | - | - |
  | $\mathrm {e} ^{at}$ | ${\frac {p}{p-a}}$ | $(\mathrm {e} ^{at}-1)$ | ${\frac {a}{p-a}}$ | $p>a,\,a>0$ |
  | $\mathrm {e} ^{-at}$ | ${\frac {p}{p+a}}$ | $-(\mathrm {e} ^{-at}-1)$ | ${\frac {a}{p+a}}$ | $p>-a,\,a<0$ |
  | ${\frac {\mathrm {e} ^{at}-1}{a}}$ | ${\frac {1}{p-a}}$ | $-{\frac {\mathrm {e} ^{-at}-1}{a}}$ | ${\frac {1}{p+a}}$ | - |
  | $\cosh(at)$ | ${\frac {p^{2}}{p^{2}-a^{2}}}$ | $\sinh(at)$ | ${\frac {pa}{p^{2}+a^{2}}}$ | $p>a$ |
  
  Pour $a=\mathrm {i} \omega$
  
  | Fonction origine | Fonction image | Fonction origine | Fonction image | Condition(s) |
  | - | - | - | - | - |
  | $\mathrm {e} ^{\mathrm {i} \omega t}$ | ${\frac {p}{p-\mathrm {i} \omega }}$ | $\mathrm {e} ^{-\mathrm {i} \omega t}$ | ${\frac {p}{p+\mathrm {i} \omega }}$ | - |
  | $\cos(\omega t)$ | ${\frac {p^{2}}{p^{2}+\omega ^{2}}}$ | $\sin(\omega t)$ | ${\frac {p\omega }{p^{2}+\omega ^{2}}}$ | - |
  | ${\frac {1-\cos({\omega t})}{\omega ^{2}}}$ | ${\frac {1}{p^{2}+\omega ^{2}}}$ | - | - | - |
  
  Si $a=\alpha +\mathrm {i} \beta$, l'image de $\mathrm {e} ^{(\alpha +\mathrm {i} \beta )t}$ est : ${\frac {p(p-\alpha +\mathrm {i} \beta )}{(p-\alpha )^{2}+\beta ^{2}}}$
  
  | Fonction origine | Fonction image | Fonction origine | Fonction image | Condition(s) |
  | - | - | - | - | - |
  | $\mathrm {e} ^{\alpha t}\cos \beta t$ | ${\frac {p(p-\alpha )}{(p-\alpha )^{2}+\beta ^{2}}}$ | $\mathrm {e} ^{\alpha t}\sin \beta t$ | ${\frac {p\beta }{(p-\alpha )^{2}+\beta ^{2}}}$ | - |
  
  Si $\alpha <0$, la valeur de $\mathrm {e} ^{(\alpha +\mathrm {i} \beta )t}$ est égale à zéro pour $t={+\infty }$, idem pour la valeur de la fonction image lorsque $p=0$.
  
  ## Hypothèse fondamentale
  
  L'hypothèse fondamentale du calcul opérationnel est que toutes fonctions d'origine f(t) ont une valeur nulle pour toute valeur de t négative. Bien que négligé la plupart du temps dans la pratique, il convient cependant d'écrire les fonctions d'origines comme facteur de la fonction $U(t)$, dite fonction échelon-unité. Exemple : la forme d'origine de ${\frac {p}{p-a}}$ est $U(t)\mathrm {e} ^{at}$.
  
  ### L'échelon unité
  
  La fonction U(t) échelon-unité est nulle pour toute valeur négative de t et égale à 1 pour toute valeur positive de t. Elle est représentée ci-dessous. Son symbole est la lettre grecque Upsilon majuscule et se lit « grand upsilon » de t. Elle se caractérise par son brusque passage de 0 à 1 entre 0- et 0+. Elle admet partout une dérivée nulle sauf en zéro où elle devient infinie.
  
  ### Introduction à la fonction de Dirac (percussion-unité)
  
  Considérons une fonction $h(t)$ telle que représentée ci-dessous. Elle est définie par :
  
  * $h(t)=0$ pour $t<0$ ;
  * $h(t)={\frac {t}{\epsilon }}$ pour $0<t<\epsilon$ ;
  * $h(t)=1$ pour $t>\epsilon$.
  
  La fonction $h(t)$ a pour dérivée $g(t)$, représentée ci-dessous, caractérisée par :
  
  * $g(t)=0$ pour $t<0$ ;
  * $g(t)={\frac {1}{\epsilon }}$ pour $0<t<\epsilon$ ;
  * $g(t)=0$ pour $t>\epsilon$.
  
  Quel que soit $\epsilon$, l'aire du rectangle est égal à l'unité.
  
  **Fonction de Dirac ou percussion-unité**
  
  Si l'on fait tendre $\epsilon$ vers zéro, $h(t)$ tend vers $U(t)$ et $g(t)$ tend vers une fonction notée $U'(t)$ qui est notre fonction de Dirac (ou percussion-unité), caractérisée par deux valeurs :
  
  * $U'(t)=0$ quel que soit t sauf pour $t=0$ où la valeur de $U'(t)$ devient infinie, et
  * $\int _{t_{0}}^{+t}U'(t)\,\mathrm {d} t=1$, quel que soit $t_0$ ≤ 0 et t > 0.
  
  Il vient alors : $\int _{0}^{+t}U'(t)\,\mathrm {d} t=U(t).$
  
  **Image de l'impulsion de Dirac**
  
  L'image de l'impulsion de Dirac est la limite quand $\epsilon$ tend vers zéro de l'expression maintenant bien connue :
  
      $\int _{0}^{+\infty }\mathrm {e} ^{-pt}g(t)\,\mathrm {d} t=\int _{0}^{\epsilon }\mathrm {e} ^{-pt}{\frac {1}{\epsilon }}\,\mathrm {d} t={\frac {1}{\epsilon }}\left[-\mathrm {e} ^{-pt}\right]_{0}^{\epsilon }$
  
  Ce qui est égal à ${\frac {1}{\epsilon }}\left(1-\mathrm {e} ^{-p\epsilon }\right)$ qui tend vers 1 quand $\epsilon$ tend vers zéro (règle de L'Hôpital, par exemple).
  
  L'image de l'impulsion de Dirac d(t) est donc 1.
  
  ## Transformation des dérivées
  
  En dérivant une fonction d'origine : $f'(t)U(t)=[f(t)U(t)]'-f(0)U'(t)$,
  
  on trouve la dérivée de la fonction d'origine : $[f(t)U(t)]'=f'(t)U(t)+f(0)U'(t)$.
  
      $f'(t)U(t)$ est l'original de $p\phi (p)-pf(0)$.
  
  On obtient donc : $[f(t)U(t)]'$ est l'original de $p\phi (p)$.
  
  Dériver une fonction d'origine revient donc à multiplier son image par $p$.
  
  ### Formule de la translation à droite
  
  Soit une fonction $f(t)U(t)$ à laquelle on fait subir une translation de la valeur $T$ à droite et parallèlement à l'axe des $t$ (voir représentations ci-dessus) de telle façon que :
  
  * $f(t)U(t)=0$ pour $t<T$, et
  * $f(t)U(t)=f(t-T)$ pour $t>T$.
  
  La forme d'origine est $F(t-T)U(t-T)$. Son image est :
  
      $p\int _{T}^{+\infty }\mathrm {e} ^{-pt}f(t-T)\,\mathrm {d} t$.
  
  En posant $\tau =t-T$ on obtient :
  
      ${\begin{aligned}p\int _{0}^{+\infty }\mathrm {e} ^{-p(T+\tau )}f(\tau )\,\mathrm {d} \tau &=\mathrm {e} ^{-pT}p\int _{0}^{+\infty }\mathrm {e} ^{-p\tau }f(\tau )\,\mathrm {d} \tau \\&=\mathrm {e} ^{-pT}\phi (p).\end{aligned}}$
  
  Conclusion :
  
  * avec $T>0$ : $f(t-T)$ est la fonction origine de $\mathrm {e} ^{-pT}\phi (p)$, avec $t>T$, et
  * $f(t-T)=0$ pour $t<T$.
  
  ## Application aux équations différentielles linéaires
  
  Ainsi à l'équation différentielle :
  
      $A{\frac {\mathrm {d} ^{3}y}{\mathrm {d} t^{3}}}+B{\frac {\mathrm {d} ^{2}y}{\mathrm {d} t^{2}}}+C{\frac {\mathrm {d} y}{\mathrm {d} t}}+Dy\mapsto f(y,y',y'',y''')$
      avec A, B, C, D étant des constantes et les conditions initiales définies en $y''_{0},y'_{0},y_{0}$,
  
  correspond une équation algébrique image de $f(y,y',y'',y''')$ :
  
      $(Ap^{3}+Bp^{2}+Cp+D)Y=Ap^{3}y_{0}+(Ay'_{0}+By_{0})p^{2}+(Ay''_{0}+By'_{0}+Cy_{0})p+\phi (p).$
  