Phénomène de Gibbs

Un article de Wikipédia, l'encyclopédie libre.
Aller à : navigation, rechercher

Lors de l'étude des séries de Fourier et des transformées de Fourier, il apparaît parfois une déformation du signal, connue sous le nom de phénomène de Gibbs. Ce phénomène est un effet de bord qui se produit à proximité d'une discontinuité, lors de l'analyse d'une fonction dérivable par morceaux.

Histoire[modifier | modifier le code]

Le phénomène fut mis pour la première fois en évidence en 1848 par Henry Wilbraham (en)[1], mais cette découverte ne connut guère d'écho.

En 1898, Albert Michelson développa un système mécanique capable de calculer et sommer la série de Fourier d'un signal donné en entrée. Il observa alors un effet d'amplification des discontinuités, qui persistait malgré l'augmentation du nombre de coefficients calculés.

Alors que Michelson soupçonnait un défaut dans la fabrication de son engin, Josiah Willard Gibbs montra que le phénomène était d'origine mathématique et se produisait dans des conditions très générales[2]. En 1906, Maxime Bôcher donna la première interprétation satisfaisante du phénomène[3] auquel il donna le nom de phénomène de Gibbs.

Description du phénomène[modifier | modifier le code]

Phénomène de Gibbs
Approximation du créneau à l'ordre 10
Approximation du créneau à l'ordre 50
Approximation du créneau à l'ordre 250

Le phénomène de Gibbs est, en quelque sorte, un « défaut d'approximation » pour une fonction continue de classe C1 par morceaux. Pour une telle fonction f, le théorème de Dirichlet assure que la série de Fourier de f converge simplement vers la fonction f sur l'intervalle où f est C1 par morceaux. En tout point x de continuité la somme de la série de Fourier est f(x).

Le polynôme trigonométrique S_N(f), N-ième somme partielle de la série de Fourier, est une fonction continue ; il est donc normal qu'il ne puisse approcher uniformément la fonction au voisinage des points de discontinuité. Inversement, sur un segment sur lequel f est dérivable, on observe une convergence uniforme, conformément au théorème de Weierstrass trigonométrique (c'est le cas des zones de « plateau » dans l'exemple de la fonction créneau).

Au point de discontinuité x, S_N(f) subit une forte oscillation, une sorte de « ressaut » qui se mesure en comparant les valeurs en x-\frac{\pi}N et x+\frac{\pi}N. En effet, toujours d'après le théorème de Dirichlet, la série de Fourier de f converge aussi simplement aux points de discontinuités mais vers la régularisée de Dirichlet, i.e. la demi-somme des valeurs de f de part et d'autre du point de discontinuité. Lorsque N devient grand, l'amplitude de ces oscillations tend vers une limite strictement plus grande que l'amplitude de la discontinuité, alors que la largeur de la zone d'oscillation tend vers 0.

Il est remarquable que le phénomène s'exprime quantitativement de façon indépendante de la fonction considérée. Si la fonction a une discontinuité d'amplitude \Delta y, alors S_N(f), tout en restant continue, connaîtra un « ressaut » en ordonnée valant de l'ordre de 9 % de plus au voisinage de la discontinuité.

Synthesis square.gif

Explication pour le signal carré[modifier | modifier le code]

Animation des harmoniques d'un signal carré. Plus le nombre d’harmoniques est grand, plus le phénomène de Gibbs sera visible.

Considérons un signal carré impair de période L=2\pi. On se place sur la discontinuité en 0, et le saut est fixé à \pi/2. On se place dans le cas d'une somme S_N(f) pour N pair, sans perte de généralité. On a alors

S_N f(x) = \sin(x) + \frac{1}{3} \sin(3x) + \cdots + \frac{1}{N-1} \sin((N-1)x).

En particulier, le théorème de Dirichlet nous donne en x=0

S_N f(0)=0=\frac{-\frac{\pi}4+\frac{\pi}4}2=\frac{f(0^-)+f(0^+)}2.

De plus,

S_N f\left(\frac{2\pi}{2N}\right)=\sin\left(\frac{\pi}N\right)+\frac13\sin\left(\frac{3\pi}N\right)+\cdots+\frac1{N-1}\sin\left(\frac{(N-1)\pi}N\right)

soit, en utilisant la fonction sinus cardinal normalisée \operatorname{sinc}(x)=\tfrac{\sin(\pi x)}{\pi x} :

S_Nf\left(\frac{2\pi}{2N}\right)=\frac{\pi}2\left[\frac2N\operatorname{sinc}\left(\frac1N\right)+\frac2N\operatorname{sinc}\left(\frac3N\right)+\cdots+\frac2N\operatorname{sinc}\left(\frac{(N-1)}N\right)\right].

Or le terme de droite n'est autre qu'une approximation de l'intégrale \frac{\pi}2\int_0^1 \operatorname{sinc}(x)\, \mathrm{d}x par une méthode du point médian de pas 2/N. Comme le sinus cardinal est continu, la somme converge vers la valeur exacte de l'intégrale pour N \to \infty :


\lim_{N\to+\infty}S_Nf\left(\frac{2\pi}{2N}\right)=\frac{\pi}2\int_0^1\operatorname{sinc}(x)\, \mathrm{d}x=\frac12\int_0^\pi\frac{\sin(t)}{t}\, \mathrm{d}t=\frac{\pi}4+\frac{\pi}2\times0{,}089490\dots

De même, on a :

\lim_{N\to+\infty} S_N f\left(-\frac{2\pi}{2N}\right)=-\frac{\pi}2\int_0^1\operatorname{sinc}(x)\,\mathrm{d}x=-\frac{\pi}4-\frac{\pi}2\times0{,}089490\dots

La quantité

\int_0^\pi\frac{\sin t}t\,\mathrm{d}t=\frac{\pi}2+\pi\times0{,}089392\dots=1{,}851937052\dots

est appelée constante de Gibbs-Wilbraham. C'est l'apparition de cette constante qui explique la surélévation du signal près de la discontinuité.

Liens avec la causalité[modifier | modifier le code]

Les appareils de mesure se comportent comme des filtres passe-bas : ils sont incapables de réagir à des signaux de trop haute fréquence. Supposons que l'on mesure un échelon de tension : il sera déformé et fera apparaître des oscillations avant et après la discontinuité du signal réel.

Cependant, cela signifierait que les oscillations se sont produites avant que l'échelon ne soit mesuré, ce qui briserait la causalité. En réalité, cela peut avoir deux interprétations : soit les appareils ne sont pas des filtres parfaits (ce qui est compréhensible), soit l'échelon de tension ne peut être parfaitement discontinu.

Si on considère l'une ou l'autre de ces corrections, la transformée de Fourier est simplement décalée : les oscillations débutent au moment où l'échelon est envoyé, l'échelon n'étant observé qu'après un léger retard.

Notes et références[modifier | modifier le code]

(en) Cet article est partiellement ou en totalité issu de l’article de Wikipédia en anglais intitulé « Gibbs phenomenon » (voir la liste des auteurs).

  1. (en) H. Wilbraham, « On a certain periodic function », Cambridge Dublin Math. J., vol. 3,‎ , p. 198-201.
  2. (en) J. W. Gibbs, « Fourier Series », Nature, vol. 59,‎ , p. 200 et 1899, p. 606.
  3. (en) M. Bôcher, « Introduction to the theory of Fourier's series », Ann. Math., vol. 7,‎ , p. 81-152.

Voir aussi[modifier | modifier le code]

Articles connexes[modifier | modifier le code]

Bibliographie[modifier | modifier le code]

(en) Antoni Zygmund, Trigonometrical series, Dover Publications, 1955.