Fonction de Heaviside

C'est une primitive de la distribution de Dirac en théorie des distributions. La valeur de H(0) a très peu d'importance, puisque la fonction est le plus souvent utilisée dans une intégrale. Certains auteurs donnent H(0) = 0, d'autres H(0) = 1. La valeur H(0) = 0,5 est souvent utilisée, parce que la fonction obtenue est ainsi symétrique. La définition est alors :

${\displaystyle \forall x\in \mathbb {R} ,\ H(x)=\left\{{\begin{matrix}0&\mathrm {si} &x<0\\{\frac {1}{2}}&\mathrm {si} &x=0\\1&\mathrm {si} &x>0.\end{matrix}}\right.}$

La valeur de la fonction en 0 est parfois notée avec un indice : la fonction H_a satisfait l'égalité H_a(0) = a pour a un réel quelconque.

La fonction est utilisée dans les mathématiques du traitement du signal pour représenter un signal obtenu en fermant un interrupteur à un instant donné et en le maintenant fermé indéfiniment.

La dérivée au sens des distributions de la fonction de Heaviside est la distribution de Dirac : ${\displaystyle \ H'=\delta }$ .

En effet, en partant tout d'abord de l'expression de la dérivation au sens des distributions :

${\displaystyle \langle H',\phi \rangle =-\langle H,\phi '\rangle }$

En appliquant ceci à l'échelon de Heaviside, nous obtenons :

${\displaystyle \langle H',\phi \rangle =-\int _{-\infty }^{+\infty }H(x)\phi '(x)\mathrm {d} x}$

Une primitive de ${\displaystyle \phi '(x)}$ est ${\displaystyle \phi (x)}$.

Nous avons alors :

${\displaystyle \langle H',\phi \rangle =-\int _{0}^{+\infty }\phi '(x)\mathrm {d} x=-\lim _{x\to \infty }\phi (x)+\phi (0)}$

Or, ${\displaystyle \phi \in {\mathcal {D}}}$, l'espace des fonctions test sur ${\displaystyle \mathbb {R} }$, donc ${\displaystyle \lim _{x\to \infty }\phi (x)=0}$.

D'où on déduit l'expression de la dérivée de l'échelon de Heaviside (au sens des distributions) :

${\displaystyle \langle H',\phi \rangle =\phi (0)=\langle \delta ,\phi \rangle ,}$

par définition de l'impulsion de Dirac, ${\displaystyle \delta }$.

Une primitive (au sens des distributions) de la fonction de Heaviside est la fonction rampe ${\displaystyle R:=XH:x\mapsto xH(x)}$. En effet,

${\displaystyle R'=(XH)'=X'H+XH'=H+X\delta =H.}$

La fonction de Heaviside est parfois utilisée pour modéliser des phénomènes variant rapidement. Toutefois, un phénomène est rarement discontinu et l'introduction d'une fonction de Heaviside dans les équations de comportement donne parfois des résultats aberrants. Par exemple, si l'on considère le démarrage d'une machine ou d'un véhicule, on considère souvent que l'accélération est nulle avant le démarrage et a une valeur fixe en phase de démarrage, la « fonction accélération » a(t) est donc modélisée par une fonction marche. Cependant, l'accélération est causée par une action mécanique associée à une déformation de la matière ; la matière ne peut pas passer d'un état « repos » à un état « déformé » instantanément, donc dans la réalité, la transition est plus « douce ». Par ailleurs, une variation instantanée de l'accélération correspondrait à un à-coup infini (fonction de Dirac) vis-à-vis des équations de mouvement, ce qui n'est pas possible.

La première solution consiste à remplacer la fonction de Heaviside par une fonction rampe c'est-à-dire une fonction linéaire passant de y = 0 à y = 1 lorsque x passe de 0 à une valeur définie δx :

${\displaystyle \forall x\in \mathbb {R} ,\ \mathrm {rampe} (x)={\begin{cases}0\qquad {\text{si }}x<0\\{\frac {x}{\delta x}}\qquad {\text{si }}0\leqslant x<\delta x\\1\qquad {\text{si }}x\geqslant \delta x\end{cases}}}$

Cette fonction est continue mais n'est pas dérivable en 0 et en δx.

Pour avoir une fonction dérivable, on utilise fréquemment une fonction polynomiale de degré 3 ; elle est dérivable une fois et sa dérivée est discontinue en début et fin de transition :

si la transition se fait sur un intervalle δx (constante réelle), alors on définit la fonction de progression

${\displaystyle \Delta (x)={\frac {x}{\delta x}}}$

et

${\displaystyle \forall x\in \mathbb {R} ,\ \mathrm {marche} (x)={\begin{cases}0\qquad {\text{si }}x<0\\\Delta ^{2}(x)\cdot (3-2\Delta (x))\qquad {\text{si }}0\leqslant x<\delta x\\1\qquad {\text{si }}x\geqslant \delta x\end{cases}}}$

De manière générale, si la fonction passe de y = h₀ à y = h₁ lorsque x passe de x₀ à x₁, on a :

${\displaystyle \Delta (x)={\frac {x-x_{0}}{x_{1}-x_{0}}}}$

${\displaystyle a=h_{1}-h_{0}}$

${\displaystyle \forall x\in \mathbb {R} ,\ \mathrm {marche} (x)={\begin{cases}h_{0}\qquad {\text{si }}x<x_{0}\\h_{0}+a\cdot \Delta ^{2}(x)\cdot (3-2\Delta (x))\qquad {\text{si }}x_{0}\leqslant x<x_{1}\\h_{1}\qquad {\text{si }}x\geqslant x_{1}\end{cases}}}$

On peut par exemple utiliser un polynôme de degré 5 sur une très courte durée (fonction souvent appelée step5, littéralement « marche5 ») ; la transition est continue, dérivable deux fois mais la dérivée troisième est discontinue en début et en fin de transition : avec les mêmes notations, pour une transition entre 0 et 1 pour x allant de 0 à δx

${\displaystyle \forall x\in \mathbb {R} ,\ \mathrm {marche5} (x)={\begin{cases}0\qquad {\text{si }}x<0\\\Delta ^{3}(x)\cdot (10-15\Delta (x)+6\Delta ^{2}(x))\qquad {\text{si }}0\leqslant x<\delta x\\1\qquad {\text{si }}x\geqslant \delta x\end{cases}}}$

et de manière plus générale :

${\displaystyle \forall x\in \mathbb {R} ,\ \mathrm {marche5} (x)={\begin{cases}h_{0}\qquad {\text{si }}x<x_{0}\\h_{0}+a\cdot \Delta ^{3}(x)\cdot (10-15\Delta (x)+6\Delta ^{2}(x))\qquad {\text{si }}x_{0}\leqslant x<x_{1}\\h_{1}\qquad {\text{si }}x\geqslant x_{1}\end{cases}}}$

• Fonction porte