F

MENUCours Outils et Méthodes pour la Physique

Jusqu'à maintenant nous n'avons envisagé que des signaux périodiques dont nous avons mis en évidence la décomposition en série de Fourier, et où le poids de chaque harmonique se distribue selon un histogramme que l'on appelle spectre du signal.
Hélas, certains signaux ne sont pas périodiques. On peut penser aux impulsions lumineuses, au potentiel d'action d'une fibre nerveuse, au signal acoustique du langage etc. Les notions d'harmoniques et de spectre ont-ils encore un sens ? Ce chapitre répond à cette question et montre le rôle unificateur que joue la transformée de Fourier dans différents domaines de la physique.

Transformation de Fourier

Du discret au continu

Nous avons établi dans le chapitre précédent que tout signal périodique physique peut s'écrire en notation complexe \[ f_T(t)=\sum_{n=-\infty}^\infty \underline{c_n}\,\mathrm{e}^{i n\,2\pi\nu\,t} \quad\text{avec}\quad \nu=\frac{1}{T} \] où \({\underline{c_n}}\) est un ensemble discret de coefficients complexes donnés par \[ \underline{c_n}=\frac{1}{T}\int_{-T/2}^{T/2}f(t) \, \mathrm{e}^{-i n\,2\pi\nu\, t}\,\mathrm{d}t \quad\text{avec}\quad n\in \mathbb{Z} \]

Considérons dorénavant un signal \(f(t)\) non périodique. On peut toujours considérer que le signal représente un unique motif d'un signal périodique dont la période tendrait vers l'infini. En d'autres termes, on peut considérer que le signal est de fréquence \(\nu\to 0\) et que les harmoniques balayent tout l'espace des réels ; le spectre discret devenant alors continu. Voyons comment cela se traduit mathématiquement. D'après les relations précédentes, un signal périodique peut s'écrire \[ f_T(t)=\sum_{n=-\infty}^\infty \left[ \nu\int_{-T/2}^{T/2}f(t')\,\mathrm{e}^{-i n\,2\pi\nu\,t'}\,\mathrm{d}t'\right]\,\mathrm{e}^{i n\,2\pi\nu\,t} \] Faisons tendre \(T\to \infty\) c'est-à-dire \(\nu\to 0\) : \[ f(t)=\lim_{\nu\to 0} \sum_{n=-\infty}^\infty \left[ \mathrm{e}^{i n\,2\pi\nu\,t} \int_{-\infty}^{\infty}f(t')\,\mathrm{e}^{-i n\,2\pi\nu\, t'}\,\mathrm{d}t'\right]\, \nu \]

On reconnaît ici une intégrale au sens de Riemann. Ainsi, un signal non périodique peut formellement s'écrire \[ f(t)=\int_{-\infty}^\infty \mathrm{e}^{i2\pi\nu\, t} \left[\int_{-\infty}^{\infty}f(t')\,\mathrm{e}^{-i2\pi\nu\,t'}\,\mathrm{d}t'\right]\, \mathrm{d}\nu \] Le terme entre crochets est appelé transformée de Fourier.

Intégrale de Fourier

Tout signal physique peut se décomposer en une intégrale de Fourier, de la forme

\begin{equation} f(t)=\int_{-\infty}^\infty \widehat{f}(\nu)\,\mathrm{e}^{i2\pi\,\nu\, t}\, \mathrm{d}\nu \tag{1.a} \end{equation}

où \(\widehat{f}(\nu)\) désigne la transformée de Fourier du signal. \(\widehat{f}(\nu)\) est une fonction continue à valeurs complexes, définie par

\begin{equation} \widehat{f}(\nu)=\int_{-\infty}^{\infty}f(t)\,\mathrm{e}^{-i2\pi\,\nu\,t}\,\mathrm{d}t \tag{1.b} \end{equation}

Ce théorème concerne les signaux de carré sommable pour lesquelles \(\int_\mathbb{R}|f(t)|^2\, \mathrm{d}t\) est finie, c'est-à-dire des signaux qui transportent une énergie finie comme c'est le cas en physique.

La transformée de Fourier \(\widehat{f}(\nu)\) est une fonction continue de la fréquence et joue le même rôle que les coefficients de Fourier du signal périodique. La relation (1.a) peut s'interpréter comme une somme continue d'harmoniques complexes, la fonction \(\widehat{f}(\nu)\) traduisant le poids relatif des diverses fréquences.

Exercice

Fonction porte

Calculer la transformée de Fourier d'une fonction porte de largeur temporelle \(\tau\) définie par : \[ \Pi_\tau(t)=\begin{cases} 1& \text{ si }|t|<\frac12 \tau\\ 0&\text{sinon} \end{cases} \]

Rép. — \(\widehat{\Pi}_\tau(\nu)=\tau\, \text{sinc}(\pi\nu\tau)\)

Spectre en puissance d'une fonction porte
Spectre en puissance d'une fonction porte de largeur temporelle \(\tau\).

Le graphe de \(|\widehat{f}|\) en fonction de la fréquence est appelé spectre d'amplitude du signal, alors que le graphe de \(|\widehat{f}|^2\) est son spectre en énergie. Par exemple, la densité spectrale de puissance d'une fonction porte à la forme d'un sinus cardinal au carré (Fig. 2). On voit notamment qu'un tel signal contient des harmoniques dont les valeurs s'étendent jusqu'à l'infini. Toutefois, l'essentiel des harmoniques se situe dans un bande de fréquence définie par \[ |\nu|<\frac{1}{\tau} \]

En d'autres termes, plus brève est l'impulsion, plus vaste est le domaine des harmoniques qui constitue le signal. De façon équivalente, plus courte est la résolution temporelle avec laquelle on veut analyser un signal plus large sera la bande passante de l'analyseur employée pour traiter le signal.

La relation de Parseval \(\overline{f^2}=\displaystyle\sum_{n=-\infty}^\infty |\underline{c_n}|^2\) rencontrée dans le chapitre précédent prend également une forme intégrale lors du passage au continu. On obtient

Relation de Parseval-Plancherel

\begin{equation} \int_\mathbb{R}|f(t)|^2\, \mathrm{d}t=\int_\mathbb{R}|\widehat{f}(\nu)|^2\, \mathrm{d}\nu \tag{2} \end{equation}

Autrement dit, l'énergie totale d'un signal ne dépend pas de la représentation choisie : fréquentielle ou temporelle.

Exercice

Vérifier la relation (2) sur l'exemple de la fonction porte sachant que \[ \int_\mathbb{R} \left(\frac{\sin x}{x}\right)^2\, \mathrm{d}x=\pi \]

Remarque

Certains préfèrent une décomposition spectrale en pulsation (\(\omega\)) plutôt qu'en fréquence (\(\nu\)). Dans ce cas les transformations de Fourier s'écrivent \[ f(t)=\frac{1}{\sqrt{2\pi}}\int_\mathbb{R} \widehat{f}(\omega)\,\mathrm{e}^{i\,\omega t}\, \mathrm{d}\omega \quad\text{avec}\quad \widehat{f}(\omega)=\frac{1}{\sqrt{2\pi}}\int_\mathbb{R}f(t)\,\mathrm{e}^{-i\,\omega t}\,\mathrm{d}t \]

Quelques propriétés de la Transformée de Fourier

On dit que \(f(t)\) et \(\widehat{f}(\nu)\) forment une paire de transformées de Fourier. On passe de l'une à l'autre par une transformation de Fourier (TF) ou transformation de Fourier inverse (TF\(^{-1}\)) : \[ f(t) \quad \rightleftharpoons\quad \widehat{f}(\nu) \]

quelques propriétés de la transformée de Fourier des signaux réels :

Théorème de modulation

Pour capturer ou transmettre un signal utile \(f(t)\), on procède parfois à une modulation d'amplitude. Cela consiste à utiliser le signal \(f(t)\) comme facteur modulant l'amplitude d'un signal sinusoïdal de haute fréquence \(\nu_0\). Le signal ainsi constitué s'écrit \[ s(t)=f(t)\cdot \cos(2\pi \nu_0 t) \] Le signal \(s(t)\) est appelé signal modulé et l'harmonique est la porteuse.

Exemple

Les premières émissions radiophoniques furent transmises par modulation d'amplitude d'une onde électromagnétique (la porteuse) de fréquence située dans la bande [150 kHz—281 kHz].

Effet d'une modulation dans l'espace de Fourier
Effet d'une modulation dans l'espace de Fourier.

Cherchons quelle est l'allure du spectre d'un tel signal modulé. Appelons \(\widehat{f}(\nu)\) la transformée de Fourier du signal utile puis calculons celle de \(s(t)\) en utilisant l'identité \(\cos(x)=\frac12(\mathrm{e}^{ix}+\mathrm{e}^{-ix})\) :

\[ \begin{array}{rcl} \widehat{s}(\nu) &= &\displaystyle \int_\mathbb{R}\frac12 f(t)\left(\mathrm{e}^{i2\pi \nu_0 t}+\mathrm{e}^{-i2\pi\nu_0 t}\right)\mathrm{e}^{-i2\pi \nu t}\, \mathrm{d}t\\[4mm] &= &\displaystyle\frac12\int_\mathbb{R}f(t)\,\mathrm{e}^{-i2\pi(\nu-\nu_0)t}\, \mathrm{d}t+\frac12\int_\mathbb{R}f(t)\,\mathrm{e}^{-i2\pi(\nu+\nu_0)t}\, \mathrm{d}t\\[4mm] \widehat{s}(\nu) &= &\dfrac12 \left[\widehat{f}(\nu-\nu_0)+\widehat{f}(\nu+\nu_0)\right] \end{array} \]

Théorème de modulation

Après modulation d'une porteuse de fréquence \(\nu_0\), le spectre est simplement translaté de \(\pm \nu_0\) (au facteur 1/2 près).

Cette translation dans l'espace des fréquences peut présenter plusieurs intérêts :

Relation temps-fréquence

Comme on l'a vu, toute dilatation de l'échelle des temps conduit à une contraction inverse de l'échelle des fréquences et réciproquement. Mathématiquement, on a

\[ f(t/a) \quad\rightleftharpoons\quad a\widehat{f}(a\,\nu) \]

Remarque

Le facteur multiplicatif \(a\) figurant devant \(\widehat{f}\) est sans importance puisqu'il ne change pas le poids relatif des diverses composantes spectrales. Il faut le voir comme une constante de normalisation permettant de respecter la relation de Parseval.

Fonction porte et son spectre en puissance
Fonction porte et son spectre en puissance.

Illustrons cette propriété sur l'exemple de la fonction porte. Il est naturel de définir la largeur temporelle d'une fonction porte par \(\Delta t=\tau\). De même, on peut introduire une largeur spectrale du spectre via \(\Delta \nu=2/\tau\). Ainsi les largeurs spectrale et temporelle sont liées via \[ \Delta \nu \times \Delta t=2 \] Autrement dit, augmenter la durée du signal d'un facteur \(a\) entraîne une contraction du même facteur de la largeur spectrale. Cette propriété est généralisable. Le produit de l'extension temporelle \(\Delta t\) d'un signal par l'extension en fréquence \(\Delta \nu\) est de l'ordre de l'unité

Dualité temps-fréquence

\begin{equation} \Delta \nu \times \Delta t \sim 1 \tag{3} \end{equation} où le symbole \(\sim 1 \) signifie que le produit donne une constante de l'ordre de l'unité qui dépend de la façon dont sont définies les extensions temporelle et spectrale.

Donnons un exemple concret : la transmission analogique d'une émission de télévision dans le standard PAL nécessite d'émettre 25 fois par seconde 625\(\times\)720 (lignes\(\times\)pixels) signaux individuels. Ainsi, la durée de chaque signal qui règle la teinte d'un point de l'écran à un instant donné ne peut dépasser \[ \Delta t=\frac{1}{25\times 625\times 720}=\frac{1}{11{,}25\cdot 10^6}\,\mathrm{s} \] Dans ces conditions, (3) montre que le signal occupe dans l'espace de Fourier une bande de largeur \(\Delta\nu\sim 10\,\mathrm{MHz}\). Les signaux de télévision étant véhiculés par des ondes hertziennes (autour de 200 MHz), on comprend pourquoi on ne peut recevoir qu'un nombre relativement limité de canaux.

Impulsion de Dirac

Définition

L'impulsion de Dirac a été inventée pour étendre le concept de transformée de Fourier aux fonctions périodiques. Mathématiquement, l'impulsion de Dirac relève de la théorie des distributions qui sort largement du cadre de ce cours. On se contentera donc d'une approche heuristique.

Représentation d'une impulsion de Dirac
Représentation d'une impulsion de Dirac.

On appelle impulsion de Dirac, un signal, noté \(\delta(t)\), qui est nul partout sauf en \(t=0\) avec la condition \[ \int_\mathbb{R}\delta(t)\, \mathrm{d}t=1 \] On peut interpréter l'impulsion de Dirac comme une fonction porte dont ont fait tendre la largeur vers 0 tout en maintenant une aire sous la courbe égale à 1 : \[ \delta(t)=\lim_{\tau\to 0}\frac{1}{\tau}\Pi_\tau(t) \] On la représente par une flèche située en \(t=0\).

Nous avons déjà calculé la transformée de Fourier d'une porte. En faisant tendre \(\tau\to 0\) on obtient la transformée de Fourier d'une impulsion de Dirac : \[ \text{TF}[\delta(t)]=\frac{1}{\tau}\lim_{\tau\to 0}\tau\, \text{sinc}(\pi\nu\tau)=1 \] En vertu de la propriété de dualité, on en déduit donc \[ \delta(t)\quad\rightleftharpoons\quad 1 \qquad\text{et}\qquad 1\quad\rightleftharpoons\quad \delta(\nu) \] L'impulsion est donc constituée d'une infinité d'harmoniques de même poids. Inversement, la fonction constante a pour transformée de Fourier un Dirac centré en \(\nu=0\).

La fonction de Dirac joue un rôle intéressant lorsqu'elle intervient dans une intégration. En effet, puisque \(\delta(t)\) est nulle partout sauf en \(t=0\), on peut écrire \(f(t)\times \delta(t)=f(0)\times \delta(t)\). Par conséquent,

Définition intégrale de l'impulsion de Dirac

\begin{equation} \int_\mathbb{R} f(t)\,\delta(t)\, \mathrm{d}t=f(0) \tag{4} \end{equation}

Exercice

Calculer la transformée de Fourier d'une impulsion de Dirac à l'aide de la définition (4).

Spectre d'un sinus

Une fois définie l'impulsion de Dirac, on obtient aisément la transformée de Fourier des fonctions circulaires. Pour cela il suffit d'utiliser le théorème de modulation : \[ \begin{array}{lcr} f(t)\cos(2\pi\nu_0\, t) &\rightleftharpoons& \dfrac12 \left[\widehat{f}(\nu-\nu_0)+\widehat{f}(\nu+\nu_0)\right]\\[3mm] f(t)\sin(2\pi\nu_0\, t) &\rightleftharpoons& \dfrac{1}{2i}\left[\widehat{f}(\nu-\nu_0)-\widehat{f}(\nu+\nu_0)\right] \end{array} \] En choisissant \(f(t)=1\) on en déduit

TF des fonctions circulaires

\[ \begin{array}{lcr} \cos(2\pi\nu_0\, t) &\rightleftharpoons & \dfrac12 \left[\delta(\nu-\nu_0)+\delta(\nu+\nu_0)\right]\\[3mm] \sin(2\pi\nu_0\, t) &\rightleftharpoons & \dfrac{1}{2i}\left[\delta(\nu-\nu_0)-\delta(\nu+\nu_0)\right] \end{array} \]

Transformation de FOurier des fonctions trigonométriques
Représentation de \(f(t)\) et \(\widehat{f}(\nu)\) pour les fonctions circulaires.

Relation entre série de Fourier et transformée de Fourier

Considérons tout d'abord un signal périodique \(f(t)\) de période \(T_0\) qui peut se décomposer en série de Fourier. On a donc \[ f(t)=a_{0}+\sum_{n=1}^{\infty}a_{n}\cos(n\, 2\pi \nu_0\, t)+b_{n}\sin(n\, 2\pi \nu_0 \, t) \] Ce signal n'est pas de et ne présente pas de transformée de Fourier au sens classique du terme. Cependant, on peut définir une transformée de Fourier d'un tel signal au sens des distributions. En effet, nous venons de rencontrer les paires de transformées suivantes : \[ \begin{array}{lcr} \cos(2\pi\nu_0\, t) &\rightleftharpoons& \dfrac12 \left[\delta(\nu-\nu_0)+\delta(\nu+\nu_0)\right]\\[3mm] \sin(2\pi\nu_0\, t) &\rightleftharpoons& \dfrac{1}{2i}\left[\delta(\nu-\nu_0)-\delta(\nu+\nu_0)\right] \end{array} \] On peut donc écrire par linéarité \[ \begin{array}{rcl} \widehat{f}(\nu)&=&\displaystyle a_{0}\,\delta(\nu)+\sum_{n=1}^{\infty}\dfrac{a_{n}}{2}\left[\delta(\nu-n\nu_0)+\delta(\nu+n\nu_0)\right]+ \dfrac{b_{n}}{2i}\left[\delta(\nu-n\nu_0)-\delta(\nu+n\nu_0)\right]\\[3mm] &=&\displaystyle a_{0}\,\delta(\nu)+\sum_{n=1}^{\infty}\dfrac{a_{n}-ib_n}{2} \delta(\nu-n\nu_0)+\dfrac{a_n+ib_n}{2}\delta(\nu+n\nu_0)\\[3mm] \widehat{f}(\nu)&=&\displaystyle\sum_{n=-\infty}^{\infty}\underline{c_n}\, \delta(\nu-n\nu_0) \end{array} \] On obtient ainsi un ensemble d'impulsions de Dirac situées tous les multiples de \(\nu_0\) et dont le poids est le coefficient de Fourier associé à la fréquence \(n \nu_0\).

Transformée de Fourier d'un signal périodique

Au sens des distributions, la transformée de Fourier d'un signal périodique de fréquence \(\nu_0\) est un peigne de Dirac de pas \(\nu_0\), modulé par les coefficients de Fourier.

Relation entre série de Fourier et transformée de Fourier
Relation entre série de Fourier et transformée de Fourier.

Autrement dit, la périodicité d'un signal se traduit dans l'espace de Fourier par l'existence d'un peigne de Dirac dont la modulation est intimement liée à la forme du motif de base ; précisément à sa transformée de Fourier comme nous allons le voir. En effet, si l'on note \(m(t)\) le motif de base défini par \[ m(t)=\begin{cases} f(t)&\text{ si }t\in [0,T_0[ \\ 0 &\text{ sinon} \end{cases} \] on peut écrire le signal \(f(t)\) sous la forme \[ f(t)=\sum_{n=-\infty}^{+\infty}m(t-nT_0) \] Il s'agit bien d'un signal périodique de période \(T_0\). En vertu du théorème de Fourier, on peut décomposer \(f(t)\) en série de Fourier : \[ f(t)=\sum_{n=-\infty}^{+\infty}\underline{c_n}\, \mathrm{e}^{i\, n\, 2\pi\nu_0 t} \quad\text{avec}\quad \underline{c_n}=\frac{1}{T_0}\int_{0}^{T_0}f(t) \, \mathrm{e}^{-i\, n\, 2\pi\nu_0 t}\,\mathrm{d}t \] Or, dans l'intervalle [0,\(T_0\)], on peut remplacer \(f(t)\) par \(m(t)\). Et comme \(m(t)\) est nulle en dehors cet intervalle, on a \[ \underline{c_n}=\frac{1}{T_0}\int_\mathbb{R}m(t) \, \mathrm{e}^{-i\, 2\pi(n\nu_0) t}\,\mathrm{d}t=\frac{1}{T_0}\widehat{m}(n\nu_0) \]

Finalement, la transformée de Fourier d'un signal périodique de période \(T_0\) est un peigne de Dirac de pas \(1/T_0\), modulé par la transformée de Fourier du motif de base.

La transformée de Fourier d'un signal périodique est un peigne de Dirac modulé par la transformée de Fourier du motif de base
La transformée de Fourier d'un signal périodique de période \(T_0\) est un peigne de Dirac de pas \(\nu_0\), modulé par la transformée de Fourier du motif de base, ici un sinus cardinal.

Exercice

Peigne de Dirac temporel

Montrer qu'un peigne de Dirac de pas \(T_0\) s'écrit \[ f(t)=\frac{1}{T_0}\sum_{n=-\infty}^{+\infty}\mathrm{e}^{i\frac{2\pi n t}{T_0}} \] En déduire la TF d'un peigne de Dirac de pas \(T_0\).

La transformée de Fourier en sciences-physiques

L'analyse spectrale est un outil dont le champ d'application est très large. Voici quelques exemples tirés des sciences-physiques.

Réponse d'un filtre

Revenons aux filtres analogiques étudiés dans le chapitre précédent. La relation entre la sortie et l'entrée est modélisée par une équation différentielle linéaire à coefficients constants du type \[ \alpha_0\,e(t)+\alpha_1 \frac{\mathrm{d}e(t)}{\mathrm{d}t}+\ldots+\alpha_n \frac{\mathrm{d}^ne(t)}{\mathrm{d}t^n}= \beta_0\,s(t)+\beta_1 \frac{\mathrm{d}s(t)}{\mathrm{d}t}+\ldots+\beta_m \frac{\mathrm{d}^m s(t)}{\mathrm{d}t^m} \] Prenons maintenant la TF de cette équation en utilisant la propriété \[ \frac{\mathrm{d}f(t)}{\mathrm{d}t} \quad\rightleftharpoons\quad i2\pi\nu\widehat{f}(\nu)=i\omega \widehat{f}(\nu) \] On obtient \[ \left[\alpha_0+i\omega\alpha_1+\ldots+(i\omega)^n\alpha_n \right]\widehat{e}(\nu)= \left[\beta_0+i\omega\beta_1+\ldots+(i\omega)^m\beta_m \right]\widehat{s}(\nu) \] Si l'on pose la fonction de transfert comme le rapport \[ \underline{H}(\nu)= \frac{\alpha_0+\mathrm{j}\omega\, \alpha_1+\ldots+(\mathrm{j}\omega)^n\, \alpha_n}{\beta_0+\mathrm{j}\omega\, \beta_1+\ldots+(\mathrm{j}\omega)^m\, \beta_m} \quad\text{avec}\quad \omega=2\pi\nu \] on trouve que dans l'espace de Fourier, l'action d'un filtre se résume à une simple multiplication : \[ \widehat{s}(\nu)=\underline{H}(\nu)\times \widehat{e}(\nu) \] Comme on l'a déjà vu pour les signaux périodiques, chaque composante spectrale du signal d'entrée est multipliée par la fonction de transfert. On retrouve alors le signal de sortie temporel par une simple transformée de Fourier inverse : \[ s(t)=\int_\mathbb{R} \left[\underline{H}(\nu)\widehat{e}(\nu)\right]\, \mathrm{e}^{-i2\pi\nu\,t}\, \mathrm{d}\nu \]

Imaginons maintenant que l'on excite un filtre à l'aide d'une impulsion de très courte durée assimilable à une impulsion de Dirac : \(e(t)=K\delta(t)\). Le signal de sortie est alors appelé réponse impulsionnelle. Puisque \(\widehat{e}(\nu)=K\), la réponse impulsionnelle s'écrit \[ s(t)=K\int_\mathbb{R} \underline{H}(\nu)\, \mathrm{e}^{-i2\pi\nu\,t}\, \mathrm{d}\nu \] Autrement dit, \(\underline{H}\) est la transformée de Fourier de la réponse impulsionnelle. On peut donc caractériser complètement un filtre en analysant sa réponse impulsionnelle. En pratique, on échantillonne la réponse impulsionnelle puis on utilise un de calcul rapide de transformée de Fourier.

Profil spectral d'une raie

Dans une lampe à décharge ou un tube fluorescent, la lumière est produite par une décharge électrique dans une ampoule contenant un gaz. Un système dispersif, tel un prisme ou un réseau de diffraction, permet d'effectuer une analyse spectrale de la lumière émise par une telle source. Le spectre ainsi obtenu est constitué de raies caractéristiques du gaz et représente son spectre d'émission.

Champ électrique rayonné lors d'une émission atomique
Champ électrique rayonné lors d'une émission atomique.

À l'échelle microscopique, la décharge électrique excite les électrons périphériques des atomes constitutifs du gaz en les portant à des niveaux d'énergie instables. C'est lors de la désexcitation qu'une émission d'onde électromagnétique se produit. Dans une approche classique de l'interaction entre matière et rayonnement on peut montrer que le champ électrique rayonné prend la forme d'un oscillateur amorti : \[ E(t)=E_0 \cos(2\pi\nu_0\, t)\cdot g(t) \quad\text{avec}\quad g(t)=\begin{cases} \mathrm{e}^{-\lambda t}&\text{si }t\geq 0\\ 0&\text{sinon} \end{cases} \] où \(\lambda\) est un coefficient d'amortissement lié à la durée de vie de l'état excité, et \(\nu_0\) la fréquence de la transition observée.

L'intensité du rayonnement est proportionnelle à \(E^2(t)\), et d'après le théorème de Parseval on a \[ \int_\mathbb{R}E^2(t)\, \mathrm{d}t=\int_\mathbb{R}|\widehat{E}(\nu)|^2\, \mathrm{d}\nu \] où \(|\widehat{E}(\nu)|^2\) représente la densité spectrale d'énergie du signal. C'est précisément ce que l'on mesure lors d'une analyse spectroscopique, à une multiplicative près. Cherchons donc quelle est l'allure de la raie spectrale associée à une telle émission en déterminant la transformé de Fourier de \(E(t)\). Pour cela, commençons d'abord par déterminer la TF de \(g(t)\) : \[ \begin{array}{rcl} \widehat{g}(\nu)&=&\displaystyle\int_\mathbb{R}g(t)\, \mathrm{e}^{-i2\pi\nu\, t}\, \mathrm{d}t\\[3mm] &=&\displaystyle\int_0^\infty \mathrm{e}^{-\lambda t}\mathrm{e}^{-i2\pi\nu\, t}\, \mathrm{d}t\\[3mm] &=&\left[-\dfrac{\mathrm{e}^{-i2\pi\nu \,t-\lambda t}}{\lambda+i2\pi\nu}\right]_0^\infty\\[3mm] \widehat{g}(\nu)&=&\dfrac{1}{\lambda+i2\pi\nu} \end{array} \]

Courbe de Lorentz
Courbe de Lorentz.

Le spectre en énergie de \(g(t)\) s'écrit \(|\widehat{g}(\nu)|^2=\dfrac{1}{\lambda^2+4\pi^2\nu^2}\). Son graphe a la forme d'une courbe en cloche appelée courbe de Lorentz ou lorentzienne.

Profil d'une raie atomique
Profil d'une raie atomique.

Utilisons maintenant le théorème de modulation : \[ E(t)=E_0 \cos(2\pi\nu_0\, t)\times g(t) \quad\rightleftharpoons\quad \frac12E_0[\widehat{g}(\nu-\nu_0)+\widehat{g}(\nu+\nu_0)] \] Si l'on se restreint aux fréquences positives, le profil spectral est une lorentzienne centrée sur la fréquence \(\nu_0\). \[ |\widehat{E}(\nu)|^2=\dfrac{{E_0}^2/4}{\lambda^2+\left[2\pi(\nu-\nu_0)\right]^2} \] Sa largeur spectrale est liée au temps de vie de l'état excité \(\tau=1/\lambda\) via la relation \[ \Delta\nu_{1/2}\times \tau=\frac{1}{\pi} \] On retrouve la dualité temps-fréquence déjà discutée. D'ailleurs en multipliant la relation précédente par la constante de Planck \(h\), on obtient la relation d'indétermination d'Heisenberg \[ \Delta E\times \tau\sim \hbar \]

Forme d'un train d'ondes quasi-harmoniques
Forme d'un train d'ondes quasi-harmoniques. Les flèches indiquent les instants de désexcitations aléatoires.

En réalité, une source est constituée d'un grand nombre d'atomes se désexcitant de façon imprévisible. Par ailleurs, les collisions inter-atomiques ont pour effet de provoquer des désexcitations supplémentaires et de réduire la durée de vie du niveau excité. Le rayonnement produit a l'allure d'un train d'ondes dont le spectre reste lorentzien, mais s'élargit du fait des .

RMN

La résonance magnétique nucléaire (RMN) est une technique d'analyse moléculaire qui consiste à plonger un échantillon dans un champ magnétique \(\overrightarrow{B_0}\) intense. Dans ce cas, certains noyaux atomiques se comportent comme des oscillateurs dont la fréquence dépend du champ \(B_0\) et de l'environnement électronique du noyau.

En envoyant un pulse électromagnétique sur l'échantillon, on excite tous ces . En calculant la TF de la réponse de l'échantillon on accède au spectre RMN qui contient autant de pics qu'il y a d'oscillateurs différents.

RMN
Schéma d'une chaine de mesure RMN.

Spectroscopie par transformée de Fourier

Principe de l'interféromètre de Michelson
Principe de l'interféromètre de Michelson.

La spectroscopie par transformée de Fourier est une importante technique qui permet d'accéder au spectre d'émission d'une source ou au spectre d'absorption d'un échantillon. Le principe repose sur l'utilisation d'un interféromètre de Michelson réglé en lame d'air. Il s'agit de mesurer le signal d'interférence en fonction du décalage optique introduit par le déplacement du miroir mobile de l'interféromètre. On accède alors à un signal — dit interférogramme — \(I(\tau)\) où \(\tau\) est un temps qui correspond au retard optique \(\tau=2x/c\) introduit par le déplacement du miroir. Ce signal permet de remonter au spectre de la source par le calcul d'une transformée de Fourier.

Interféromgramme pour une source monochromatique Interféromgramme pour une raie gaussienne
Relation de dualité entre spectre et interférogramme.
FTIR
Schéma d'un spectrophotomètre

Cette technique est implémentée dans de nombreux spectrophotomètres notamment , car elle offre une meilleure résolution que les spectromètres à réseau.

Techniquement, on montre que l'interférogramme s'écrit \[ I(\tau)=2I_0\,[1+\underbrace{\int_\mathbb{R}\widehat{s}(\nu)\cos(2\pi\nu\, \tau)\, \mathrm{d}\nu}_{\text{terme d'interférence }f(\tau)}] \quad\text{avec}\quad \int_\mathbb{R}\widehat{s}(\nu)\, \mathrm{d}\nu=1 \] où \(\widehat{s}(\nu)\) est le spectre normalisé de la source. Comme \(\widehat{s}(\nu)\) est une quantité réelle, le terme d'interférence est la transformée de Fourier inverse du spectre ; réciproquement le spectre est la transformée de Fourier du terme d'interférence : \[ f(\tau)=\int_\mathbb{R}\widehat{s}(\nu)\,\mathrm{e}^{i2\pi \nu \, \tau}\, \mathrm{d}\nu \quad\rightleftharpoons\quad \widehat{s}(\nu)=\int_\mathbb{R}f(\tau)\,\mathrm{e}^{-i2\pi \nu \, \tau}\, \mathrm{d}\tau \]

La diffraction, une transformée de Fourier naturelle

Éclairons une pupille diffractante par une onde plane lumineuse en incidence normale, puis observons la lumière diffractée sur un capteur placé à . On montre que l'onde diffractée dans la direction donnée par le vecteur d'onde \(\overrightarrow{k}\) s'écrit [7] \[ \underline{\psi}(\overrightarrow{k})\propto \iint_{(S)}\underline{t}(x,y)\,\mathrm{e}^{i\, \overrightarrow{k}\cdot \overrightarrow{r}}\, \mathrm{d}x \mathrm{d}y \quad\text{avec}\quad k=\frac{2\pi}{\lambda} \] où \(\overrightarrow{r}=(x,y)\) repère un point P de la pupille et \(\underline{t}(x,y)\) désigne la transparence de la pupille diffractante au point P.

Paramétrage du problème de diffraction en champ lointain
Paramétrage du problème de diffraction en champ lointain.

On peut exprimer \(\overrightarrow{k}\) en fonction des angles de diffraction \(\theta_x\) et \(\theta_y\) : \[ \overrightarrow{k}=\frac{2\pi}{\lambda}\sin\theta_x\,\overrightarrow{u_x}+ \frac{2\pi}{\lambda}\sin\theta_y\,\overrightarrow{u_y}+ k_z\,\overrightarrow{u_z} \] de sorte que le terme de phase s'écrit \[ \overrightarrow{k}\cdot \overrightarrow{r}=\frac{2\pi}{\lambda}(x\sin\theta_x+y\sin\theta_y)=2\pi(x\nu_x+y\nu_y) \] où l'on fait intervenir les variables duales \((\nu_x,\nu_y)\) à \((x,y)\). \(\nu_x\) et \(\nu_y\) sont appelées fréquences spatiales. Finalement, l'onde diffractée dans la direction \(\overrightarrow{k}\) s'écrit \[ \underline{\psi}(\overrightarrow{k})\propto \iint_{(S)}\underline{t}(x,y)\,\mathrm{e}^{i\, 2\pi(x\nu_x+y\nu_y)}\, \mathrm{d}x \mathrm{d}y=\mathrm{TF_{2d}}[t(x,y)] \] où \(\mathrm{TF_{2d}}\) est la généralisation de la transformée de Fourier à deux dimensions.

Par exemple, prenons comme pupille diffractante un écran percé par une ouverture rectangulaire de longueur \(2a\) et de largeur \(2b\). Sa fonction de transparence est nulle en dehors du trou et vaut 1 pour tout point situé dans le trou ; en d'autres termes elle s'écrit comme un produit de deux fonctions portes : \[ t(x,y)=\Pi_{2a}(x)\times \Pi_{2b}(y) \] Le calcul de la \(\mathrm{TF_{2d}}\) est grandement simplifié car le domaine d'intégration est un rectangle : \[ \underline{\psi}(\overrightarrow{k})\propto \mathrm{TF}[\Pi_{2a}(x)]\times \mathrm{TF}[\Pi_{2b}(y)]= (4ab)\mathrm{sinc}(2\pi a\nu_x)\times \mathrm{sinc}(2\pi b\nu_y) \] L'intensité lumineuse détectée étant proportionnelle au carré de l'amplitude diffractée, on trouve finalement \[ I(\overrightarrow{k})=I_\text{max}\left[\mathrm{sinc}(2\pi a\nu_x)\times \mathrm{sinc}(2\pi b\nu_y)\right]^2 \] La répartition de l'intensité est représenté Fig. 17. On retrouve la propriété de dualité selon laquelle une petite dimension spatiale entraîne un étalement dans l'espace de Fourier, c'est-à-dire un élargissement des taches de diffraction. Notez que la symétrie de l'obstacle se retrouve dans la tache de diffraction.

Pupille diffractante (à gauche) et tache de diffraction (à droite)
Pupille diffractante (à gauche) et tache de diffraction (à droite). Simulation réalisée pour \(\lambda=700\,\mathrm{nm}\), \(a=0{,}56~\mathrm{mm}\), \(b=3a\) et \(z=7~\mathrm{m}\).

De la même manière, on peut obtenir des informations sur la structure des cristaux en observant la diffraction des par ces cristaux. On montre que les rayons X sont essentiellement diffusés par les électrons de sorte que la fonction de transparence d'un échantillon cristallin est proportionnelle à la densité électronique \(\rho\). Le cliché obtenu après diffraction des rayons X par un cristal est donnée par \(|\mathrm{TF_{3d}}[\rho(x,y,z)]|^2\).

Pour en savoir plus...

  1. R. Renaux-petel L'analyse de Fourier en physique2015. Disponible sur fermedesetoiles.fr
  2. B. Houchmandzadeh Mathématiques pour la Physique2010. Disponible sur hal.archives-ouvertes.fr
  3. J. Peatross et al. Physics of Light and Optics2015. Disponible sur optics.byu.edu
  4. H. Gie & J-P. Sarmant Électromagnétisme, Volume 2Paris : Technique et documentation, Lavoisier, 1985.
  5. J-P. Lecardonnel A propos de sons périodiques sans fondamental BUP№767, 1994.
  6. R. Bracewell L'analyse de Fourier Pour la science №142, p.74-80, 1989.
  7. J. Roussel Diffraction de Fraunhofer[en ligne], 2018. Disponible sur femto-physique.fr