Transformation de Fourier

Portrait de Joseph Fourier.
Transformée de Fourier
Type
Nom court
(en) FTVoir et modifier les données sur Wikidata
Nommé en référence à
Joseph FourierVoir et modifier les données sur Wikidata
Décrit par
ISO 80000-2:2019 Quantities and units — Part 2: Mathematics (d)Voir et modifier les données sur Wikidata
Aspect de
Analyse de FourierVoir et modifier les données sur Wikidata
Formule
( F f ) ( ω ) = e i ω t f ( t ) d t {\displaystyle \left({\mathcal {F}}f\right)(\omega )=\int \limits _{-\infty }^{\infty }\mathrm {e} ^{-\mathrm {i} \omega t}f(t)\mathrm {d} t} Voir et modifier les données sur Wikidata

modifier - modifier le code - modifier WikidataDocumentation du modèle

En mathématiques, plus précisément en analyse, la transformation de Fourier est une extension, pour les fonctions non périodiques, du développement en série de Fourier des fonctions périodiques. La transformation de Fourier associe à toute fonction intégrable définie sur ℝ et à valeurs réelles ou complexes, une autre fonction sur ℝ appelée transformée de Fourier dont la variable indépendante peut s'interpréter en physique comme la fréquence ou la pulsation.

La transformée de Fourier représente une fonction par la densité spectrale dont elle provient, en tant que moyenne de fonctions trigonométriques de toutes fréquences. La théorie de la mesure ainsi que la théorie des distributions permettent de définir rigoureusement la transformée de Fourier dans toute sa généralité, elle joue un rôle fondamental dans l'analyse harmonique.

Lorsqu'une fonction représente un phénomène physique, comme l'état du champ électromagnétique ou du champ acoustique en un point, on l'appelle signal et sa transformée de Fourier s'appelle son spectre.

Transformation de Fourier pour les fonctions intégrables

Définition

La transformation de Fourier F {\displaystyle {\mathcal {F}}} est une opération qui transforme une fonction intégrable sur ℝ en une autre fonction, décrivant le spectre fréquentiel de cette dernière. Si f est une fonction intégrable sur ℝ, sa transformée de Fourier est la fonction F ( f ) = f ^ {\displaystyle {\mathcal {F}}(f)={\hat {f}}} donnée par la formule[1] :

F ( f ) : ξ f ^ ( ξ ) = + f ( x ) e i ξ x d x {\displaystyle {\mathcal {F}}(f):\xi \mapsto {\hat {f}}(\xi )=\int _{-\infty }^{+\infty }f(x)\,\mathrm {e} ^{-{\rm {i}}\xi x}\,\mathrm {d} x} .

Conventions alternatives

Il est possible de choisir une définition alternative pour la transformation de Fourier. Ce choix est une affaire de convention dont les conséquences ne se manifestent (en général) que par des facteurs multiplicatifs constants. Par exemple, certains scientifiques[Lesquels ?] utilisent ainsi[réf. nécessaire] :

F ( f ) : ν f ^ ( ν ) = + f ( t ) e i 2 π ν t d t {\displaystyle {\mathcal {F}}(f):\nu \mapsto {\hat {f}}(\nu )=\int _{-\infty }^{+\infty }f(t)\,\mathrm {e} ^{-{\rm {i}}2\pi \nu t}\,\mathrm {d} t}

avec t en secondes et ν la fréquence (en hertz).

Certains utilisent (pour des raisons de symétrie avec la transformation de Fourier inverse) la transformation suivante[2] :

F ( f ) : ω f ^ ( ω ) = 1 2 π + f ( t ) e i ω t d t {\displaystyle {\mathcal {F}}(f):\omega \mapsto {\hat {f}}(\omega )={1 \over {\sqrt {2\pi }}}\int _{-\infty }^{+\infty }f(t)\,\mathrm {e} ^{-{\rm {i}}\omega t}\,\mathrm {d} t}

avec t en secondes et ω la pulsation (en radians par seconde).

Cette définition n'est cependant pas adaptée au traitement des produits de convolution : à cause du facteur 1 2 π {\displaystyle {\frac {1}{\sqrt {2\pi }}}} , on a F ( f g ) F ( f ) F ( g ) {\displaystyle {\mathcal {F}}(f*g)\neq {\mathcal {F}}(f)\cdot {\mathcal {F}}(g)} , à moins d'introduire un tel facteur dans la définition du produit de convolution.

L'ensemble de départ est l'ensemble des fonctions intégrables f d'une variable réelle x. L'ensemble d'arrivée est l'ensemble des fonctions d'une variable réelle ξ. Concrètement lorsque cette transformation est utilisée en traitement du signal, on notera volontiers t à la place de x et ω ou ν à la place de ξ qui seront les variables respectives de temps et de pulsation ou de fréquence. On dira alors que f est dans le domaine temporel, et que f ^ {\displaystyle {\hat {f}}} est dans le domaine fréquentiel.

En physique, la transformation de Fourier permet de déterminer le spectre d'un signal. Les phénomènes de diffraction donnent une image de l'espace dual du réseau, ils sont une sorte de « machine à transformation de Fourier » naturelle. Pour ces applications, les physiciens définissent en général la transformation directe avec un facteur 1 2 π {\displaystyle {\frac {1}{\sqrt {2\pi }}}} et la transformation de Fourier inverse avec le même préfacteur.

La notation F ( f ) {\displaystyle {\mathcal {F}}(f)} peut aussi être remplacée par F(ƒ) ou TF(ƒ). Dans cet article, on utilisera exclusivement la première notation.

Il est également d'usage dans certaines communautés scientifiques de noter f(x) pour la fonction de départ et f(p) pour sa transformée, faisant ainsi correspondre à x, y, z les variables duales p, q, r. Cette notation est conforme à l'interprétation physique inspirée par la mécanique quantique : dualité entre position et quantité de mouvement. Cette notation n'est pas retenue ici.

Extension de la transformation de Fourier

Le cadre le plus naturel pour définir les transformations de Fourier est celui des fonctions intégrables. Toutefois, de nombreuses opérations (dérivations, transformation de Fourier inverse) ne peuvent être écrites en toute généralité. On doit à Plancherel[réf. nécessaire] l'introduction de la transformation de Fourier pour les fonctions de carré sommable, pour lesquelles la formule d'inversion est vraie. Puis la théorie des distributions de Schwartz, et plus particulièrement des distributions tempérées permit de trouver un cadre parfaitement adapté.

On peut généraliser la définition de la transformation de Fourier à plusieurs variables, et même sur d'autres groupes que le groupe additif ℝn. Ainsi, on peut la définir sur le groupe additif ℝ/ℤ, c'est-à-dire sur les fonctions de période 1 — on retrouve ainsi les séries de Fourier —, et plus généralement sur des groupes localement compacts, pas nécessairement commutatifs, et en particulier sur des groupes finis. Ces définitions font intervenir les groupes duaux, ainsi que la mesure de Haar.

Propriétés de la transformation de Fourier

Fonction Transformée de Fourier
Linéarité a g 1 ( x ) + b g 2 ( x ) {\displaystyle a\cdot g_{1}(x)+b\cdot g_{2}(x)} a g ^ 1 ( ξ ) + b g ^ 2 ( ξ ) {\displaystyle a\cdot {\hat {g}}_{1}(\xi )+b\cdot {\hat {g}}_{2}(\xi )}
Contraction du domaine f ( a x )   {\displaystyle f(a\cdot x)\ } 1 | a | f ^ ( ξ / a ) {\displaystyle {\frac {1}{|a|}}\cdot {\hat {f}}(\xi /a)}
Translation temporelle g ( x + x 0 )   {\displaystyle g(x+x_{0})\ } g ^ ( ξ ) e i ξ x 0 {\displaystyle {\hat {g}}(\xi )\cdot \mathrm {e} ^{\mathrm {i} \xi x_{0}}}
Modulation dans le domaine temporel g ( x ) e i x ξ 0 {\displaystyle g(x)\cdot \mathrm {e} ^{\mathrm {i} x\xi _{0}}} g ^ ( ξ ξ 0 ) {\displaystyle {\hat {g}}(\xi -\xi _{0})}
Produit de convolution ( f g ) ( x ) {\displaystyle (f*g)(x)} f ^ ( ξ ) g ^ ( ξ ) {\displaystyle {\hat {f}}(\xi )\cdot {\hat {g}}(\xi )}
Produit ( f g ) ( x ) {\displaystyle (f\cdot g)(x)} 1 2 π ( f ^ g ^ ) ( ξ ) {\displaystyle {\frac {1}{2\pi }}({\hat {f}}*{\hat {g}})(\xi )}
Dérivation dans le domaine temporel f ( x ) {\displaystyle f'(x)}

(voir conditions ci-dessous)

i ξ f ^ ( ξ ) {\displaystyle {\rm {i}}\xi \cdot {\hat {f}}(\xi )}
Dérivation dans le domaine fréquentiel x f ( x ) {\displaystyle x\cdot f(x)} i f ^ ( ξ ) {\displaystyle \mathrm {i} {\hat {f}}'(\xi )}
Symétrie réelle et paire réelle et paire
réelle paire (à symétrie hermitienne)
réelle et impaire imaginaire pure et impaire
imaginaire pure et paire imaginaire pure et paire
imaginaire pure et impaire réelle et impaire
Forme gaussienne gaussienne
  • La contraction dans un domaine (temporel, spatial ou fréquentiel) implique une dilatation dans l'autre. Un exemple concret de ce phénomène peut être observé par exemple sur un tourne-disque. La lecture d'un 33 tours à 45 tours par minute implique une augmentation de la fréquence du signal audio (a > 1), on contracte le signal audio dans le domaine temporel ce qui le dilate dans le domaine fréquentiel.
  • Si la fonction f est à support borné (c.-à-d. si x 0 R , | x | > x 0 , f ( x ) = 0 {\displaystyle \exists x_{0}\in \mathbb {R} ,\forall |x|>x_{0},f(x)=0} ) alors f ^ {\displaystyle {\hat {f}}} est à support infini. Inversement, si le support spectral de la fonction f ^ {\displaystyle {\hat {f}}} est borné alors f est à support non borné.
  • Si f est une fonction non nulle sur un intervalle borné alors f ^ {\displaystyle {\hat {f}}} est une fonction non nulle sur C {\displaystyle \mathbb {C} } et inversement, si f ^ {\displaystyle {\hat {f}}} est non nulle sur un intervalle borné alors f est une fonction non nulle sur C {\displaystyle \mathbb {C} } .
  • La transformée de Fourier de f est une fonction continue, de limite nulle à l'infini (théorème de Riemann-Lebesgue), notamment bornée par
f ^ f 1 {\displaystyle \|{\hat {f}}\|_{\infty }\leq \|f\|_{1}} .
  • Par changement de variable on trouve des formules intéressantes lorsqu'on effectue une translation, dilatation du graphe de f.
  • Supposons que la fonction g : x i x f ( x ) {\displaystyle g:x\mapsto -\mathrm {i} xf(x)} soit intégrable ; alors on peut dériver la formule de définition sous le signe d'intégration. On constate alors que la dérivée f ^ {\displaystyle {\hat {f}}'} est la transformée de Fourier de g.
  • Si f est localement absolument continue (c.-à-d. dérivable presque partout et égale à « l'intégrale de sa dérivée ») et si f et f ' sont intégrables, alors[3] la transformée de Fourier de la dérivée de f est f ^ ( ξ ) = i ξ f ^ ( ξ ) {\displaystyle {\widehat {f'}}(\xi )=\mathrm {i} \xi {\hat {f}}(\xi )} .

On peut résumer les deux dernières propriétés : notons D l'opération

D f = 1 i f {\displaystyle Df={\frac {1}{\mathrm {i} }}f'}

et M la multiplication par l'argument :

( M f ) ( x ) = x f ( x ) , ( M f ^ ) ( ξ ) = ξ f ^ ( ξ ) {\displaystyle (Mf)(x)=xf(x),\quad (M{\hat {f}})(\xi )=\xi {\hat {f}}(\xi )} .

Alors, si f satisfait des conditions fonctionnelles convenables, D f ^ = + M f ^ {\displaystyle {\widehat {Df}}=+M{\hat {f}}} et M f ^ = D f ^ {\displaystyle {\widehat {Mf}}=-D{\hat {f}}} .

On s'affranchira de ces conditions fonctionnelles en élargissant la classe des objets sur lesquels opère la transformation de Fourier. C'est une des motivations de la définition des distributions.

Transformation de Fourier inverse

Si la transformée de Fourier de f, notée f ^ {\displaystyle {\hat {f}}} , est elle-même une fonction intégrable, la formule dite de transformation de Fourier inverse, opération notée F 1 {\displaystyle {\mathcal {F}}^{-1}} , et appliquée à f ^ {\displaystyle {\hat {f}}} , permet (sous conditions appropriées) de retrouver f à partir des données fréquentielles :

f ( x ) = F 1 ( f ^ ) ( x ) = 1 2 π + f ^ ( ξ ) e + i ξ x d ξ f ^ ( ξ )   = + f ( x ) e i ξ x d x {\displaystyle f(x)={\mathcal {F}}^{-1}({\hat {f}})(x)={1 \over 2\pi }\,\int _{-\infty }^{+\infty }{\hat {f}}(\xi )\,\mathrm {e} ^{+{\rm {i}}\xi x}\,\mathrm {d} \xi \qquad \Leftrightarrow \qquad {\hat {f}}(\xi )\ =\int _{-\infty }^{+\infty }f(x)\,\mathrm {e} ^{-{\rm {i}}\xi x}\,\mathrm {d} x} .

Cette opération de transformation de Fourier inverse a des propriétés analogues à la transformation directe, puisque seuls changent le coefficient multiplicatif et le –i devenu i.

Dans le cas des définitions alternatives, la transformation de Fourier inverse devient :

Définition en fréquence : f ( t ) = + f ^ ( ν ) e + i 2 π ν t d ν f ^ ( ν ) = + f ( t ) e i 2 π ν t d t {\displaystyle f(t)=\int _{-\infty }^{+\infty }{\hat {f}}(\nu )\,\mathrm {e} ^{+{\rm {i}}2\pi \nu t}\,\mathrm {d} \nu \qquad \Leftrightarrow \qquad {\hat {f}}(\nu )=\int _{-\infty }^{+\infty }f(t)\,\mathrm {e} ^{-{\rm {i}}2\pi \nu t}\,\mathrm {d} t} .
Définition en pulsation : f ( t ) = 1 2 π + f ^ ( ω ) e + i ω t d ω f ^ ( ω )   = 1 2 π + f ( t ) e i ω t d t {\displaystyle f(t)={1 \over {\sqrt {2\pi }}}\int _{-\infty }^{+\infty }{\hat {f}}(\omega )\,\mathrm {e} ^{+{\rm {i}}\omega t}\,\mathrm {d} \omega \quad \Leftrightarrow \quad {\hat {f}}(\omega )\ ={1 \over {\sqrt {2\pi }}}\int _{-\infty }^{+\infty }f(t)\,\mathrm {e} ^{-{\rm {i}}\omega t}\,\mathrm {d} t} .

Soit h une fonction complexe définie sur ℝ et deux fois continûment différentiable. On suppose que h vérifie l'estimation

| h ( x ) | C 1 + x 2 {\displaystyle |h(x)|\leq {\frac {C}{1+x^{2}}}}

et que les deux premières dérivées de h sont intégrables sur ℝ. Alors la transformée de Fourier de h vérifie une estimation analogue

| h ^ ( ξ ) | C 1 + ξ 2 {\displaystyle |{\hat {h}}(\xi )|\leq {\frac {C}{1+\xi ^{2}}}} .

Soit y un nombre réel qui, pour le moment, est simplement un paramètre, et notons :

f ( x ) = h ( x ) e i y x {\displaystyle f(x)=h(x)\mathrm {e} ^{-\mathrm {i} yx}} .

On vérifie que f a les mêmes propriétés fonctionnelles que h. Par conséquent, on peut appliquer la formule sommatoire de Poisson à f, avec la période  :

n Z f ( x + 2 π n ) = 1 2 π k Z f ^ ( k ) e i k x {\displaystyle \sum _{n\in \mathbb {Z} }f(x+2\pi n)={\frac {1}{2\pi }}\sum _{k\in \mathbb {Z} }{\hat {f}}(k)\mathrm {e} ^{{\rm {i}}kx}} .

Mais le calcul de f ^ ( k ) {\displaystyle {\hat {f}}(k)} donne :

f ^ ( k ) = R h ( x ) e i ( y + k ) x d x = h ^ ( y + k ) {\displaystyle {\hat {f}}(k)=\int _{\mathbb {R} }h(x)\mathrm {e} ^{-{\rm {i}}(y+k)x}\,\mathrm {d} x={\hat {h}}(y+k)} .

On peut donc réécrire la formule sommatoire de Poisson en termes de h, et il vient :

n Z h ( x + 2 π n ) e i ( x + 2 π n ) y = 1 2 π k Z h ^ ( y + k ) e i k x {\displaystyle \sum _{n\in \mathbb {Z} }h(x+2\pi n)\mathrm {e} ^{-{\rm {i}}(x+2\pi n)y}={\frac {1}{2\pi }}\sum _{k\in \mathbb {Z} }{\hat {h}}(y+k)\mathrm {e} ^{{\rm {i}}kx}} .

On multiplie les deux membres de cette identité par e i x y {\displaystyle \mathrm {e} ^{\mathrm {i} xy}}  :

n Z h ( x + 2 π n ) e 2 i π n y = 1 2 π k Z h ^ ( y + k ) e i ( k + y ) x {\displaystyle \sum _{n\in \mathbb {Z} }h(x+2\pi n)\mathrm {e} ^{-2{\rm {i}}\pi ny}={\frac {1}{2\pi }}\sum _{k\in \mathbb {Z} }{\hat {h}}(y+k)\mathrm {e} ^{{\rm {i}}(k+y)x}} .

On remarque que les séries apparaissant de part et d'autre sont normalement convergentes pour la norme du maximum. On va donc pouvoir échanger la sommation et l'intégration par rapport à y sur l'intervalle [0 ; 1].

À gauche, l'intégration par rapport à y ne laisse subsister qu'un seul terme, celui correspondant à n = 0. À droite, on intègre par rapport à y et l'on effectue dans chaque intégrale le changement de variable y + k = ξ. On obtient ainsi la formule

h ( x ) = 1 2 π R h ^ ( ξ ) e i x ξ d ξ {\displaystyle h(x)={\frac {1}{2\pi }}\int _{\mathbb {R} }{\hat {h}}(\xi )\mathrm {e} ^{{\rm {i}}x\xi }\,\mathrm {d} \xi } .

On passe au cas général de la formule d'inversion de Fourier pour une fonction f intégrable ainsi que sa transformée de Fourier par une méthode de densité. On approche f par une suite de fonctions fp vérifiant les hypothèses fonctionnelles de la présente démonstration. On doit bien sûr supposer que les fp et leurs transformées de Fourier f ^ p {\displaystyle {\hat {f}}_{p}} convergent vers leurs limites respectives f et f ^ {\displaystyle {\hat {f}}} en norme L1(ℝ). On peut construire de telles approximations en tronquant f, c'est-à-dire en le remplaçant par 0 en dehors de l'intervalle [–p, p], et en le régularisant par convolution. Si ϕ est une fonction deux fois continûment différentiable, d'intégrale 1, et à support borné, on pose ϕp(x) = p ϕ(px) et l'on convole la fonction tronquée f|[–p,p] par ϕp. C'est une idée raisonnable d'utiliser ici le même paramètre p.

Soit f une fonction de classe C à support compact. Par le principe de transfert, on peut se contenter d'étudier le cas d'une fonction standard. Dans ce cas, il existe un réel infiniment grand T tel que pour tout réel |x| > T, f(x) = 0. Introduisons une base hilbertienne de L2([–T, T]) donnée par :

e n : x e i n π x / T {\displaystyle e_{n}:x\mapsto {\rm {e}}^{\mathrm {i} n\pi x/T}}

(un calcul immédiat montre qu'elle est bien orthonormée, et le fait qu'elle soit totale se déduit de la densité des fonctions continues et de leur approximation uniforme par des polynômes trigonométriques). Par le lemme de Parseval, on est en mesure d'écrire :

f = n Z c n e n {\displaystyle f=\sum _{n\in \mathbb {Z} }c_{n}e_{n}} c n = 1 2 T T T f ( x ) e i n π x / T d x = 1 2 T f ^ ( n π T ) {\displaystyle c_{n}={\frac {1}{2T}}\int _{-T}^{T}f(x){\rm {e}}^{-\mathrm {i} n\pi x/T}{\rm {d}}x={\frac {1}{2T}}{\widehat {f}}\left({\frac {n\pi }{T}}\right)}

Plus explicitement, pour x standard :

f ( x ) = 1 2 π n Z f ^ ( n π T ) e i n π x / T π T = 1 2 π f ^ ( w ) e i w x d w {\displaystyle f(x)={\frac {1}{2\pi }}\sum _{n\in \mathbb {Z} }{\widehat {f}}\left({\frac {n\pi }{T}}\right){\rm {e}}^{\mathrm {i} n\pi x/T}{\frac {\pi }{T}}={\frac {1}{2\pi }}\int _{-\infty }^{\infty }{\widehat {f}}(w){\rm {e}}^{{\rm {i}}wx}{\rm {d}}w} .

La dernière égalité vient de ce que le membre de gauche est standard, que la somme de Riemann s'effectue sur une partition de longueur infiniment petite (π/T), et donc que le membre de droite est la partie standard du membre intermédiaire. L'égalité recherchée est donc vraie pour toutes les fonctions standard de classe C à support compact et tout x standard. Par le principe de transfert, elle est aussi vérifiée pour toutes les fonctions C à support compact et tout x, puis par densité des fonctions C à support compact dans l'espace des fonctions intégrables, pour toutes les fonctions intégrables dont la transformée est intégrable et pour presque tout x.

Extension à l'espace ℝn

Notons x∙ξ le produit scalaire canonique dans ℝn :

x ξ = j = 1 n x j ξ j {\displaystyle x\cdot \xi =\sum _{j=1}^{n}x_{j}\xi _{j}} .

Si f est une fonction intégrable sur ℝn, sa transformée de Fourier est donnée par la formule :

f ^ ( ξ ) = R n f ( x )   e i x ξ   d x {\displaystyle {\hat {f}}(\xi )=\int _{\mathbb {R} ^{n}}f(x)~{\rm {e}}^{-{\rm {i}}x\cdot \xi }~\mathrm {d} x} .

Si A est une isométrie linéaire directe, f A ^ = f ^ A {\displaystyle {\widehat {f\circ A}}={\hat {f}}\circ A} . Il en résulte que la transformée de Fourier d'une fonction radiale est radiale.

Expression de la transformée de Fourier dans ℝn d'une fonction radiale

Par définition :

f ^ ( t 1 , t 2 , , t n ) = R n f ( x 1 , x 2 , , x n )   e 2 i π ( x 1 t 1 + x 2 t 2 + + x n t n )   d x 1 d x 2 d x n {\displaystyle {\hat {f}}(t_{1},t_{2},\dots ,t_{n})=\int _{\mathbb {R} ^{n}}f(x_{1},x_{2},\dots ,x_{n})~{\rm {e}}^{{\rm {2i\pi }}(x_{1}t_{1}+x_{2}t_{2}+\dots +x_{n}t_{n})}~\mathrm {d} x_{1}\mathrm {d} x_{2}\dots \mathrm {d} x_{n}} .

Si l'on se place dans le cas où f est radiale (ou à symétrie sphérique) alors f ne dépend des variables x1, ..., xn que par l'intermédiaire de la variable ρ = x 1 2 + x 2 2 + + x n 2 {\displaystyle \rho ={\sqrt {x_{1}^{2}+x_{2}^{2}+\dots +x_{n}^{2}}}} . On montre alors que f ^ {\displaystyle {\hat {f}}} ne dépend des variables t1, ..., tn que par l'intermédiaire de la variable τ = t 1 2 + t 2 2 + + t n 2 {\displaystyle \tau ={\sqrt {t_{1}^{2}+t_{2}^{2}+\dots +t_{n}^{2}}}} .

Soit f(x1, ... , xn) = g(ρ).

En notant les vecteurs : ρ = ( x 1 x 2 x n )  et  τ = ( t 1 t 2 t n ) ρ τ = ρ τ cos θ = x 1 t 1 + x 2 t 2 + + x n t n {\displaystyle {\vec {\rho }}=\left({\begin{array}{*{20}{c}}{x_{1}}\\{x_{2}}\\\vdots \\{x_{n}}\end{array}}\right){\text{ et }}{\vec {\tau }}=\left({\begin{array}{*{20}{c}}{t_{1}}\\{t_{2}}\\\vdots \\{t_{n}}\end{array}}\right)\Rightarrow {\vec {\rho }}\cdot {\vec {\tau }}=\rho \tau \cos \theta =x_{1}t_{1}+x_{2}t_{2}+\dots +x_{n}t_{n}}

En passant des coordonnées cartésiennes aux coordonnées polaires dans ℝn :
f ^ ( τ ) = R n g ( ρ )   e 2 i π τ ρ   d n ρ = R n g ( ρ )   e 2 i π τ ρ cos θ   d n ρ {\displaystyle {\hat {f}}({\vec {\tau }})=\int _{\mathbb {R} ^{n}}g(\rho )~{\rm {e}}^{{\rm {2i\pi }}{\vec {\tau }}\cdot {\vec {\rho }}}~\mathrm {d} ^{n}{\vec {\rho }}=\int _{\mathbb {R} ^{n}}g(\rho )~{\rm {e}}^{{\rm {2i\pi }}\tau \rho \cos \theta }~\mathrm {d} ^{n}{\vec {\rho }}}

Considérons la rotation R {\displaystyle {\mathcal {R}}} telle que τ = R ( τ ) f ^ ( τ ) = R n g ( ρ )   e 2 i π τ ρ   d n ρ = R n g ( ρ )   e 2 i π R ( τ ) ρ   d n ρ {\displaystyle {\vec {\tau }}'={\mathcal {R}}({\vec {\tau }})\Rightarrow {\hat {f}}({\vec {\tau }}')=\int _{\mathbb {R} ^{n}}g(\rho )~{\rm {e}}^{{\rm {2i\pi }}{\vec {\tau }}'\cdot {\vec {\rho }}}~\mathrm {d} ^{n}{\vec {\rho }}=\int _{\mathbb {R} ^{n}}g(\rho )~{\rm {e}}^{{\rm {2i\pi }}{\mathcal {R}}({\vec {\tau }}){\vec {\rho }}}~\mathrm {d} ^{n}{\vec {\rho }}}

On ne change pas la valeur de l'intégrale si on remplace ρ {\displaystyle {\vec {\rho }}} par R ( ρ ) {\displaystyle {\mathcal {R}}({\vec {\rho }})} du fait que g est radiale.
f ^ ( τ ) = R n g ( ρ )   e 2 i π R ( τ ) R ( ρ )   d n R ( ρ ) {\displaystyle \Rightarrow {\hat {f}}({\vec {\tau }}')=\int _{\mathbb {R} ^{n}}g(\rho )~{\rm {e}}^{{\rm {2i\pi }}{\mathcal {R}}({\vec {\tau }}){\mathcal {R}}({\vec {\rho }})}~\mathrm {d} ^{n}{\mathcal {R}}({\vec {\rho }})}
Comme R ( τ ) R ( ρ ) = τ ρ = τ ρ cos θ {\displaystyle {\mathcal {R}}({\vec {\tau }})\cdot {\mathcal {R}}({\vec {\rho }})={\vec {\tau }}\cdot {\vec {\rho }}=\tau \rho \cos \theta } et d n R ( ρ ) = d n ρ {\displaystyle \mathrm {d} ^{n}{\mathcal {R}}({\vec {\rho }})=\mathrm {d} ^{n}{\vec {\rho }}}
f ^ ( τ ) = R n g ( ρ )   e 2 i π τ ρ cos θ   d n ρ = f ^ ( τ ) {\displaystyle \Rightarrow {\hat {f}}({\vec {\tau }}')=\int _{\mathbb {R} ^{n}}g(\rho )~{\rm {e}}^{{\rm {2i\pi }}\tau \rho \cos \theta }~\mathrm {d} ^{n}\rho ={\hat {f}}({\vec {\tau }})}

La transformée de Fourier d'une fonction radiale est donc aussi une fonction radiale (qui ne dépend que de τ {\displaystyle \parallel {\vec {\tau }}\parallel } ).

On rappelle la correspondance entre coordonnées sphériques et coordonnées polaires dans ℝn, coordonnées aussi appelées « Coordonnées hypersphériques ».
x i = ρ cos ( φ n i + 1 ) k = 1 n i sin ( φ k ) {\displaystyle x_{i}=\rho \cos(\varphi _{n-i+1})\prod _{k=1}^{n-i}\sin(\varphi _{k})} .
On montre par ailleurs que le jacobien de la transformation des coordonnées cartésiennes en coordonnées hypersphériques est : J = ρ n 1 i = 1 n 2 sin n 1 i ( φ i ) {\displaystyle J=\rho ^{n-1}\prod _{i=1}^{n-2}\sin ^{n-1-i}(\varphi _{i})}
avec φ1 ≤ jn – 2 ∈ [0,π] et φn – 1 ∈ [0,2π].
Il en résulte :
f ^ ( τ ) = R n g ( ρ )   e 2 i π τ ρ cos θ ρ n 1 d ρ d φ n 1 j = 1 n 2 sin n 1 j ( φ j ) d φ j {\displaystyle {\hat {f}}({\vec {\tau }})=\int _{\mathbb {R} ^{n}}g(\rho )~{\rm {e}}^{{\rm {2i\pi }}\tau \rho \cos \theta }\rho ^{n-1}\mathrm {d} \rho \mathrm {d} \varphi _{n-1}\prod _{j=1}^{n-2}\sin ^{n-1-j}(\varphi _{j})\mathrm {d} \varphi _{j}}
Du fait de la symétrie radiale, on ne change rien de l'intégrale si on considère τ {\displaystyle {\vec {\tau }}} parallèle à l'axe x1. Cela revient alors à avoir θ = φ1 (et indépendant des φj ≠ 1)

f ^ ( τ ) = 0 + g ( ρ ) ρ n 1 d ρ ( 0 π sin n 2 ( θ )   e 2 i π τ ρ cos θ d θ ) < 2 > ( j = 2 n 2 φ j = 0 π sin n 1 j ( φ j ) d φ j ) < 1 > ( 0 2 π d φ n 1 ) {\displaystyle \Rightarrow {\hat {f}}({\vec {\tau }})=\int _{0}^{+\infty }g(\rho )\rho ^{n-1}\mathrm {d} \rho \underbrace {\left(\int _{0}^{\pi }\sin ^{n-2}(\theta )~{\rm {e}}^{{\rm {2i\pi }}\tau \rho \cos \theta }\mathrm {d} \theta \right)} _{<2>}\underbrace {\left(\prod _{j=2}^{n-2}\int _{\varphi _{j}=0}^{\pi }\sin ^{n-1-j}(\varphi _{j})\mathrm {d} \varphi _{j}\right)} _{<1>}\left(\int _{0}^{2\pi }\mathrm {d} \varphi _{n-1}\right)}

Calcul de <1>
Posons I j = 0 π sin n 1 j ( φ j ) d φ j {\displaystyle I_{j}=\int _{0}^{\pi }\sin ^{n-1-j}(\varphi _{j})\mathrm {d} \varphi _{j}}
On reconnaît ici la fonction bêta B ( p , q ) = 2 0 π 2 sin 2 p 1 ( α ) cos 2 q 1 ( α ) d α = Γ ( p ) Γ ( q ) Γ ( p + q ) {\displaystyle \mathrm {B} (p,q)=2\int _{0}^{\frac {\pi }{2}}\sin ^{2p-1}(\alpha )\cos ^{2q-1}(\alpha )\mathrm {d} \alpha ={\frac {\Gamma (p)\Gamma (q)}{\Gamma (p+q)}}} avec Γ la fonction gamma et p, q réels positifs.
I j = 0 π sin n 1 j ( φ j ) d φ j = Γ ( n j 2 ) Γ ( 1 2 ) Γ ( n j + 1 2 ) ⇒< 1 >= j = 2 n 2 Γ ( n j 2 ) Γ ( 1 2 ) Γ ( n j + 1 2 ) = π n 3 2 Γ ( n 1 2 ) {\displaystyle \Rightarrow I_{j}=\int _{0}^{\pi }\sin ^{n-1-j}(\varphi _{j})\mathrm {d} \varphi _{j}={\frac {\Gamma ({\frac {n-j}{2}})\Gamma ({\frac {1}{2}})}{\Gamma ({\frac {n-j+1}{2}})}}\Rightarrow <1>=\prod _{j=2}^{n-2}{\frac {\Gamma ({\frac {n-j}{2}})\Gamma ({\frac {1}{2}})}{\Gamma ({\frac {n-j+1}{2}})}}={\frac {\pi ^{\frac {n-3}{2}}}{\Gamma ({\frac {n-1}{2}})}}} avec Γ(1) = 1; Γ(1/2) = π
f ^ ( τ ) = 2 π n 1 2 Γ ( n 1 2 ) 0 + g ( ρ ) ρ n 1 d ρ ( 0 π sin n 2 ( θ )   e 2 i π τ ρ cos θ d θ ) < 2 > {\displaystyle \Rightarrow {\hat {f}}({\vec {\tau }})=2{\frac {\pi ^{\frac {n-1}{2}}}{\Gamma ({\frac {n-1}{2}})}}\int _{0}^{+\infty }g(\rho )\rho ^{n-1}\mathrm {d} \rho \underbrace {\left(\int _{0}^{\pi }\sin ^{n-2}(\theta )~{\rm {e}}^{{\rm {2i\pi }}\tau \rho \cos \theta }\mathrm {d} \theta \right)} _{<2>}}
On notera au passage
2 0 π 2 sin n 2 α d α = Γ ( n 1 2 ) π Γ ( n 2 ) f ^ ( τ ) = 2 π n 2 Γ ( n 2 ) 0 + g ( ρ ) ρ n 1 d ρ ( 0 π sin n 2 ( θ )   e 2 i π τ ρ cos θ d θ ) 0 π sin n 2 α d α {\displaystyle 2\int _{0}^{\frac {\pi }{2}}\sin ^{n-2}\alpha \mathrm {d} \alpha ={\frac {\Gamma ({\frac {n-1}{2}}){\sqrt {\pi }}}{\Gamma ({\frac {n}{2}})}}\Rightarrow {\hat {f}}({\vec {\tau }})=2{\frac {\pi ^{\frac {n}{2}}}{\Gamma ({\frac {n}{2}})}}\int _{0}^{+\infty }g(\rho )\rho ^{n-1}\mathrm {d} \rho {\frac {\left(\int _{0}^{\pi }\sin ^{n-2}(\theta )~{\rm {e}}^{{\rm {2i\pi }}\tau \rho \cos \theta }\mathrm {d} \theta \right)}{\int _{0}^{\pi }\sin ^{n-2}\alpha \mathrm {d} \alpha }}}
Calcul de <2>
Considérons la fonction
L n ( z ) = 0 π e i z cos θ sin n 2 θ d θ 0 π sin n 2 θ d θ {\displaystyle L_{n}(z)={\frac {\int _{0}^{\pi }\mathrm {e} ^{\mathrm {i} z\cos \theta }\sin ^{n-2}\theta \,\mathrm {d} \theta }{\int _{0}^{\pi }\sin ^{n-2}\theta \,\mathrm {d} \theta }}}
On notera Ln(0) = 1
On a alors
d L n ( z ) d z = 0 π i cos θ sin n 2 θ e i z cos θ d θ 0 π sin n 2 θ d θ {\displaystyle {\frac {\mathrm {d} L_{n}(z)}{\mathrm {d} z}}={\frac {\int _{0}^{\pi }i\cos \theta \sin ^{n-2}\theta \mathrm {e} ^{\mathrm {i} z\cos \theta }\mathrm {d} \theta }{\int _{0}^{\pi }\sin ^{n-2}\theta \,\mathrm {d} \theta }}}

En intégrant par parties l'intégrale en numérateur, on établit la relation :

d L n ( z ) d z = z n 1 0 π sin n θ e i z cos θ d θ 0 π sin n 2 θ d θ {\displaystyle -{\frac {\mathrm {d} L_{n}(z)}{\mathrm {d} z}}={\frac {z}{n-1}}{\frac {\int _{0}^{\pi }\sin ^{n}\theta \mathrm {e} ^{\mathrm {i} z\cos \theta }\mathrm {d} \theta }{\int _{0}^{\pi }\sin ^{n-2}\theta \,\mathrm {d} \theta }}}
On notera alors : d L n ( 0 ) d z = 0 {\displaystyle {\frac {\mathrm {d} L_{n}(0)}{\mathrm {d} z}}=0} .
En dérivant une seconde fois :
d 2 L n ( z ) d z 2 = 0 π sin n 2 θ e i z cos θ d θ 0 π sin n 2 θ d θ + 0 π sin n θ e i z cos θ d θ 0 π sin n 2 θ d θ d 2 L n ( z ) d z 2 + n 1 z d L n ( z ) d z + L n ( z ) = 0 {\displaystyle {\frac {d^{2}L_{n}(z)}{dz^{2}}}=-{\frac {\int _{0}^{\pi }\sin ^{n-2}\theta \mathrm {e} ^{\mathrm {i} z\cos \theta }\mathrm {d} \theta }{\int _{0}^{\pi }\sin ^{n-2}\theta \,\mathrm {d} \theta }}+{\frac {\int _{0}^{\pi }\sin ^{n}\theta \mathrm {e} ^{\mathrm {i} z\cos \theta }\mathrm {d} \theta }{\int _{0}^{\pi }\sin ^{n-2}\theta \,\mathrm {d} \theta }}\Rightarrow {\frac {d^{2}L_{n}(z)}{dz^{2}}}+{\frac {n-1}{z}}{\frac {dL_{n}(z)}{dz}}+L_{n}(z)=0} .
On reconnaît ici une équation qui est proche de l'équation différentielle de Bessel. Pour faire disparaître le facteur n – 1 du deuxième terme, posons :
L n ( z ) = a n z m J m ( z ) {\displaystyle L_{n}(z)=a_{n}z^{-m}J_{m}(z)} .
En reportant cette expression dans l'équation différentielle, on arrive à :
d 2 J m ( z ) d z 2 + ( n 1 2 m z ) d J m ( z ) d z + ( m ( m n + 2 ) z 2 + 1 ) J m ( z ) = 0 {\displaystyle {\frac {\mathrm {d} ^{2}J_{m}(z)}{\mathrm {d} z^{2}}}+\left({\frac {n-1-2m}{z}}\right){\frac {\mathrm {d} J_{m}(z)}{\mathrm {d} z}}+\left({\frac {m(m-n+2)}{z^{2}}}+1\right)J_{m}(z)=0} .
Il suffit alors de poser m = n – 2/2 pour arriver à l'équation différentielle de Bessel suivante :
d 2 d z 2 J n 2 2 ( z ) + 1 z d d z J n 2 2 ( z ) + ( 1 ( n 2 2 z ) 2 ) J n 2 2 ( z ) = 0 {\displaystyle {\frac {\mathrm {d} ^{2}}{\mathrm {d} z^{2}}}J_{\frac {n-2}{2}}(z)+{\frac {1}{z}}{\frac {\mathrm {d} }{\mathrm {d} z}}J_{\frac {n-2}{2}}(z)+\left(1-\left({\frac {n-2}{2z}}\right)^{2}\right)J_{\frac {n-2}{2}}(z)=0}
Il s'agit bien d'une équation différentielle de Bessel dont la fonction de Bessel J n 2 2 ( z ) {\displaystyle J_{\frac {n-2}{2}}(z)} est solution. Il en résulte alors la relation suivante, par définition de la fonction de Bessel :
J n 2 2 ( z ) = ( z 2 ) n 2 2 p = 0 + ( 1 ) p p ! Γ ( p + n 2 ) ( z 2 ) 2 p = z n 2 2 a n L n ( z ) L n ( z ) = a n ( 1 2 ) n 2 2 p = 0 + ( 1 ) p p ! Γ ( p + n 2 ) ( z 2 ) 2 p {\displaystyle J_{\frac {n-2}{2}}(z)=\left({\frac {z}{2}}\right)^{\frac {n-2}{2}}\sum _{p=0}^{+\infty }{\frac {(-1)^{p}}{p!\Gamma (p+{\frac {n}{2}})}}\left({\frac {z}{2}}\right)^{2p}={\frac {z^{\frac {n-2}{2}}}{a_{n}}}L_{n}(z)\Rightarrow L_{n}(z)=a_{n}\left({\frac {1}{2}}\right)^{\frac {n-2}{2}}\sum _{p=0}^{+\infty }{\frac {(-1)^{p}}{p!\Gamma (p+{\frac {n}{2}})}}\left({\frac {z}{2}}\right)^{2p}}
Avec
L n ( 0 ) = 1 a n = 2 n 2 2 Γ ( n 2 ) L n ( z ) = Γ ( n 2 ) p = 0 + ( 1 ) p p ! Γ ( p + n 2 ) ( z 2 ) 2 p = 0 π e i z cos θ sin n 2 θ d θ 0 π sin n 2 θ d θ {\displaystyle L_{n}(0)=1\Rightarrow a_{n}=2^{\frac {n-2}{2}}\Gamma \left({\frac {n}{2}}\right)\Rightarrow L_{n}(z)=\Gamma \left({\frac {n}{2}}\right)\sum _{p=0}^{+\infty }{\frac {(-1)^{p}}{p!\Gamma (p+{\frac {n}{2}})}}\left({\frac {z}{2}}\right)^{2p}={\frac {\int _{0}^{\pi }\mathrm {e} ^{\mathrm {i} z\cos \theta }\sin ^{n-2}\theta \,\mathrm {d} \theta }{\int _{0}^{\pi }\sin ^{n-2}\theta \,\mathrm {d} \theta }}}
L n ( 2 π ρ τ ) = Γ ( n 2 ) p = 0 + ( 1 ) p p ! Γ ( p + n 2 ) ( 2 π ρ τ 2 ) 2 p = 0 π e i 2 π ρ τ cos θ sin n 2 θ d θ 0 π sin n 2 θ d θ {\displaystyle \Rightarrow L_{n}(2\pi \rho \tau )=\Gamma \left({\frac {n}{2}}\right)\sum _{p=0}^{+\infty }{\frac {(-1)^{p}}{p!\,\Gamma (p+{\frac {n}{2}})}}\left({\frac {2\pi \rho \tau }{2}}\right)^{2p}={\frac {\int _{0}^{\pi }\mathrm {e} ^{\mathrm {i} 2\pi \rho \tau \cos \theta }\sin ^{n-2}\theta \,\mathrm {d} \theta }{\int _{0}^{\pi }\sin ^{n-2}\theta \,\mathrm {d} \theta }}}
⇒< 2 >= 0 π e 2 i π ρ τ cos θ sin n 2 θ d θ = 2 n 2 2 Γ ( n 2 ) ( 2 π ρ τ ) n 2 2 J n 2 2 ( 2 π ρ τ ) 0 π sin n 2 θ d θ {\displaystyle \Rightarrow <2>=\int _{0}^{\pi }\mathrm {e} ^{2\mathrm {i} \pi \rho \tau \cos \theta }\sin ^{n-2}\theta \,\mathrm {d} \theta =2^{\frac {n-2}{2}}\Gamma \left({\frac {n}{2}}\right)(2\pi \rho \tau )^{-{\frac {n-2}{2}}}J_{\frac {n-2}{2}}(2\pi \rho \tau )\int _{0}^{\pi }\sin ^{n-2}\theta \,\mathrm {d} \theta }
En revenant à l'expression de la transformée de Fourier :
f ^ ( τ ) = 2 π n 2 Γ ( n 2 ) 0 + g ( ρ ) ρ n 1 d ρ ( 0 π sin n 2 ( θ )   e 2 i π τ ρ cos θ d θ ) 0 π sin n 2 α d α = 2 π ( τ ) 2 n 2 0 + g ( ρ ) ( ρ ) n 2 J n 2 2 ( 2 π ρ τ ) d ρ {\displaystyle {\hat {f}}({\vec {\tau }})=2{\frac {\pi ^{\frac {n}{2}}}{\Gamma ({\frac {n}{2}})}}\int _{0}^{+\infty }g(\rho )\rho ^{n-1}\mathrm {d} \rho {\frac {\left(\int _{0}^{\pi }\sin ^{n-2}(\theta )~{\rm {e}}^{{\rm {2i\pi }}\tau \rho \cos \theta }\mathrm {d} \theta \right)}{\int _{0}^{\pi }\sin ^{n-2}\alpha \mathrm {d} \alpha }}=2\pi (\tau )^{\frac {2-n}{2}}\int _{0}^{+\infty }g(\rho )(\rho )^{\frac {n}{2}}J_{\frac {n-2}{2}}(2\pi \rho \tau )\mathrm {d} \rho } .

Si la transformée de Fourier de f est elle-même une fonction intégrable, on a alors la formule d'inversion :

f ( x ) = 1 ( 2 π ) n R n f ^ ( ξ )   e i x ξ   d ξ {\displaystyle f(x)={\frac {1}{(2\pi )^{n}}}\int _{\mathbb {R} ^{n}}{\hat {f}}(\xi )~{\rm {e}}^{{\rm {i}}x\cdot \xi }~\mathrm {d} \xi } .

Par conséquent, la transformation de Fourier de L1 dans C0 est injective (mais pas surjective).

Transformation de Fourier pour les fonctions de carré sommable

Extension de la transformation de L1∩L2 à L2

Le théorème de Plancherel permet de donner un sens à la transformée de Fourier des fonctions de carré sommable sur ℝ.

On commence par un premier résultat préparatoire.

Lemme — Soit h une fonction complexe deux fois continûment dérivable sur ℝ, qui vérifie l'estimation

x R | h ( x ) | C / ( 1 + x 2 ) {\displaystyle \forall x\in \mathbb {R} \quad |h(x)|\leq C/(1+x^{2})} (où C est une constante),

et dont les deux premières dérivées sont intégrables. Ceci implique que la transformée de Fourier h ^ {\displaystyle {\hat {h}}} est bien définie et de carré intégrable. De plus, on a l'identité :

R | h ( x ) | 2 d x = 1 2 π R | h ^ ( ξ ) | 2 d ξ {\displaystyle \int _{\mathbb {R} }|h(x)|^{2}\,{\rm {d}}x={\frac {1}{2\pi }}\int _{\mathbb {R} }|{\hat {h}}(\xi )|^{2}\,{\rm {d}}\xi } .

On reprend la formule établie ci-dessus dans la démonstration de la formule d'inversion de Fourier :

n Z h ( x + 2 π n ) e 2 i π n y = 1 2 π k Z h ^ ( y + k ) e i ( k + y ) x {\displaystyle \sum _{n\in Z}h(x+2\pi n){\rm {e}}^{-2\mathrm {i} \pi ny}={\frac {1}{2\pi }}\sum _{k\in \mathbb {Z} }{\hat {h}}(y+k){\rm {e}}^{\mathrm {i} (k+y)x}} .

On prend le carré du module des deux membres, et l'on intègre sur l'intervalle [0 ; 1] par rapport à y et sur l'intervalle [0 ; 2π] par rapport à x :

0 1 0 2 π m , n Z h ( x + 2 π n ) h ¯ ( x + 2 π m ) e 2 i π ( m n ) y d x d y = 1 4 π 2 0 1 0 2 π j , k Z h ^ ( y + j ) h ^ ¯ ( y + k ) e i ( j k ) x d x d y {\displaystyle \int _{0}^{1}\int _{0}^{2\pi }\sum _{m,n\in \mathbb {Z} }h(x+2\pi n){\bar {h}}(x+2\pi m){\rm {e}}^{2\mathrm {i} \pi (m-n)y}\,{\rm {d}}x\,{\rm {d}}y={\frac {1}{4\pi ^{2}}}\int _{0}^{1}\int _{0}^{2\pi }\sum _{j,k\in \mathbb {Z} }{\hat {h}}(y+j){\bar {\hat {h}}}(y+k){\rm {e}}^{\mathrm {i} (j-k)x}\,{\rm {d}}x\,{\rm {d}}y} .

On peut échanger l'ordre de la sommation et des deux intégrations dans l'expression ci-dessus, parce que les hypothèses faites sur h impliquent que les séries convergent normalement dans l'espace des fonctions continues de x et y, périodiques de période en x et de période 1 en y. L'intégration en y du premier membre ne laisse subsister que les termes pour lesquels m et n sont égaux, et l'intégration en x du deuxième membre ne laisse subsister que les termes pour lesquels j et k sont identiques. Il reste donc :

n Z 0 2 π | h ( x + 2 π n ) | 2 d x = 1 2 π k Z 0 1 | h ^ ( y + k ) | 2 d y {\displaystyle \sum _{n\in \mathbb {Z} }\int _{0}^{2\pi }|h(x+2\pi n)|^{2}\,{\rm {d}}x={\frac {1}{2\pi }}\sum _{k\in \mathbb {Z} }\int _{0}^{1}|{\hat {h}}(y+k)|^{2}\,{\rm {d}}y} .

Il suffit de faire dans le premier membre le changement de variable dans chaque intégrale x + 2π n = x' et dans le second le changement de variable dans chaque intégrale y + k = ξ, et l'on obtient la formule :

R | h ( x ) | 2 d x = 1 2 π R | h ^ ( ξ ) | 2 d ξ {\displaystyle \int _{\mathbb {R} }|h(x')|^{2}\,{\rm {d}}x'={\frac {1}{2\pi }}\int _{\mathbb {R} }|{\hat {h}}(\xi )|^{2}\,{\rm {d}}\xi } .

Après changement de la variable muette x' en x, on obtient la formule annoncée.

Une fois démontrée dans le lemme ci-dessus la formule de Plancherel pour une classe de fonctions suffisamment régulières, on étend par densité la transformation de Fourier à tout L2(ℝ).

Extension de la transformation de Fourier par densité

On adopte encore les mêmes notations que dans la démonstration de la formule d'inversion de Fourier par la formule sommatoire de Poisson, donc ϕ est une fonction deux fois continûment différentiable, à support compact, et d'intégrale 1. On pose ϕp(x) =p ϕ(px).

Soit h une fonction de carré intégrable, et soit p un nombre entier quelconque. On définit

h p = ( h 1 [ p , p ] ) ϕ p {\displaystyle h_{p}=(h1_{[-p,p]})*\phi _{p}}

et l'on peut montrer le résultat suivant :

lim p R | h h p | 2 d x = 0 {\displaystyle \lim _{p\to \infty }\int _{\mathbb {R} }|h-h_{p}|^{2}\,{\rm {d}}x=0} .

La démonstration utilise des techniques classiques d'approximation par régularisation.

D'autre part, les fonctions hp ont les propriétés nécessaires pour appliquer le lemme ci-dessus, et en particulier

R | h p h q | 2 d x = 1 2 π R | h ^ p h ^ q | 2 d ξ {\displaystyle \int _{\mathbb {R} }|h_{p}-h_{q}|^{2}\,{\rm {d}}x={\frac {1}{2\pi }}\int _{\mathbb {R} }|{\hat {h}}_{p}-{\hat {h}}_{q}|^{2}\,{\rm {d}}\xi } .

Comme la suite (hp)p ≥ 1 est de Cauchy dans l'espace L2(ℝ), la suite des transformées de Fourier ( h ^ p ) p 1 {\displaystyle ({\hat {h}}_{p})_{p\geq 1}} est aussi de Cauchy, donc elle converge. Sa limite, que l'on note h ^ {\displaystyle {\hat {h}}} , ne dépend pas du choix de la suite d'approximations. En effet, si gp était une autre suite d'approximations convergeant vers h en moyenne quadratique, et satisfaisant les conditions fonctionnelles sous lesquelles on peut appliquer la formule sommatoire de Poisson, on aurait l'estimation

g p h p 2 g p h 2 + h h p 2 {\displaystyle \|g_{p}-h_{p}\|_{2}\leq \|g_{p}-h\|_{2}+\|h-h_{p}\|_{2}} ,

qui tend vers 0 pour p tendant vers l'infini. Par conséquent, g ^ p h ^ p 2 {\displaystyle \|{\hat {g}}_{p}-{\hat {h}}_{p}\|_{2}} tend aussi vers 0 et l'on conclut que la limite de la suite g ^ p {\displaystyle {\hat {g}}_{p}} est bien h ^ {\displaystyle {\hat {h}}} .

On a ainsi le théorème de Plancherel :

Théorème de Plancherel —  Soit f une fonction complexe sur ℝ et de carré sommable. Alors la transformée de Fourier de f peut être définie comme suit : pour tout p entier, on pose

f p ( x ) = ( f 1 [ p , p ] ) ( x ) = { f ( x ) si  | x | p , 0 sinon. {\displaystyle f_{p}(x)=(f1_{[-p,p]})(x)={\begin{cases}f(x)&{\text{si }}|x|\leq p,\\0&{\text{sinon.}}\end{cases}}}

La suite des transformées de Fourier f ^ p {\displaystyle {\hat {f}}_{p}} converge dans L2(ℝ), et sa limite est la transformée de Fourier f ^ {\displaystyle {\hat {f}}} , c'est-à-dire

lim p R | f ^ ( ξ ) f ^ p ( ξ ) | 2 d ξ = 0 {\displaystyle \lim _{p\to \infty }\int _{\mathbb {R} }|{\hat {f}}(\xi )-{\hat {f}}_{p}(\xi )|^{2}\,{\rm {d}}\xi =0} .

De plus, on a l'identité :

R | f ( x ) | 2 d x = 1 2 π R | f ^ ( ξ ) | 2 d ξ {\displaystyle \int _{\mathbb {R} }|f(x)|^{2}\,{\rm {d}}x={\frac {1}{2\pi }}\int _{\mathbb {R} }|{\hat {f}}(\xi )|^{2}{\rm {d}}\xi } .

De façon similaire, si l'on pose g p ( x ) = p p f ^ ( ξ ) e i x ξ d ξ {\displaystyle g_{p}(x)=\int _{-p}^{p}{\hat {f}}(\xi ){\rm {e}}^{\mathrm {i} x\xi }\,{\rm {d}}\xi } , les gp convergent en moyenne quadratique vers f.

Démonstration du théorème de Plancherel

L'identité suivante résulte du procédé d'extension décrit ci-dessus :

R | h ( x ) | 2 d x = 1 2 π R | h ^ ( ξ ) | 2 d ξ {\displaystyle \int _{\mathbb {R} }|h(x)|^{2}\,{\rm {d}}x={\frac {1}{2\pi }}\int _{\mathbb {R} }|{\hat {h}}(\xi )|^{2}\,{\rm {d}}\xi } .

Considérons alors la suite de fonctions fp = f 1{[–p,p]} . En vertu du théorème de convergence dominée de Lebesgue pour les fonctions de carré sommable, la suite des fp converge en moyenne quadratique vers f et par conséquent, on aura aussi

lim p f ^ f ^ p 2 = 0 {\displaystyle \lim _{p\to \infty }\|{\hat {f}}-{\hat {f}}_{p}\|_{2}=0} .

En d'autres termes, f ^ p {\displaystyle {\hat {f}}_{p}} converge en moyenne quadratique vers f ^ {\displaystyle {\hat {f}}} . La démonstration pour la formule d'inversion est analogue.

Ainsi la transformation de Fourier-Plancherel définit un automorphisme intemporel de l'espace L2, qui est une isométrie, à condition de faire un changement d'échelle si l'on utilise la notation en pulsation

f ^ / 2 π 2 = f 2 π 2 {\displaystyle \|{\hat {f}}/{\sqrt {2\pi }}\|_{2}=\|f{\sqrt {2\pi }}\|_{2}} .

En physique, on interprète le terme | f ^ ( ξ ) / 2 π | 2 {\displaystyle |{\hat {f}}(\xi )/{\sqrt {2\pi }}|^{2}} figurant sous l'intégrale comme une densité spectrale de puissance.

La définition de la transformation de Fourier-Plancherel est compatible avec la définition habituelle de la transformée de Fourier des fonctions intégrables. Sur l'intersection L1(ℝ)} ∩ L2(ℝ) des domaines de définition, on montre à l'aide du théorème de convergence dominée de Lebesgue que les deux définitions coïncident.

La transformation vue comme opérateur de L2(ℝ)

Remarque : ce paragraphe utilise la définition fréquentielle de la transformée de Fourier, pour des raisons d'isométrie.

Nous venons de voir que la transformation de Fourier F {\displaystyle {\mathcal {F}}} induit sur l'espace de Hilbert L2(ℝ) un opérateur linéaire. Nous en récapitulons ici les propriétés :

  • F {\displaystyle {\mathcal {F}}} est un opérateur unitaire de L2. Il s'agit en particulier d'une isométrie. On retrouve le premier fait, connu sous le nom de formule de Parseval, affirmant que pour toutes fonctions f, g ∈ L2(ℝ),

f g = F f F g , i.e.  f ( x ) g ( x ) ¯ d x = f ^ ( ν ) g ^ ( ν ) ¯ d ν {\displaystyle \langle f\mid g\rangle =\langle {\mathcal {F}}f\mid {\mathcal {F}}g\rangle ,\quad {\text{i.e. }}\quad \int _{-\infty }^{\infty }f(x){\overline {g(x)}}\,{\rm {d}}x=\int _{-\infty }^{\infty }{\hat {f}}(\nu ){\overline {{\hat {g}}(\nu )}}\,{\rm {d}}\nu } et en particulier le deuxième fait, connu sous le nom de théorème de Plancherel f 2 = f ^ 2 , i.e.  + | f ( x ) | 2 d x = + | f ^ ( ν ) | 2 d ν {\displaystyle \|f\|_{2}=\|{\hat {f}}\|_{2},\quad {\text{i.e. }}\quad \int _{-\infty }^{+\infty }|f(x)|^{2}\,{\rm {d}}x=\int _{-\infty }^{+\infty }|{\hat {f}}(\nu )|^{2}\,{\rm {d}}\nu }  ;

ϕ n ( x ) = 2 1 / 4 n ! e π x 2 H n ( 2 x π ) {\displaystyle {\phi }_{n}(x)={\frac {2^{1/4}}{\sqrt {n!}}}\,{\rm {e}}^{-\pi x^{2}}H_{n}(2x{\sqrt {\pi }})} ,Hn(x) sont les polynômes d'Hermite « probabilistes », qui s'écrivent H n ( x ) = ( 1 ) n e x 2 2 d n d x n e x 2 2 {\displaystyle H_{n}(x)=(-1)^{n}{\rm {e}}^{\frac {x^{2}}{2}}{\frac {\mathrm {d} ^{n}}{\mathrm {d} x^{n}}}{\rm {e}}^{-{\frac {x^{2}}{2}}}} . Avec ces notations, la formule suivante récapitule la situation ϕ ^ n ( ν ) = ( i ) n ϕ n ( ν ) {\displaystyle {\hat {\phi }}_{n}(\nu )=(-{\rm {i}})^{n}{\phi }_{n}(\nu )} . On retrouve la gaussienne comme première fonction d'Hermite. Ces fonctions appartiennent à la classe de Schwartz S {\displaystyle {\mathcal {S}}} .

Lien avec le produit de convolution

La transformation de Fourier a des propriétés très intéressantes liées au produit de convolution. On rappelle que (d'après l'inégalité de Young pour la convolution) :

  • si f , g L 1 ( R N ) {\displaystyle f,g\in {\rm {L}}^{1}(\mathbb {R} ^{N})} , alors f g L 1 ( R N ) {\displaystyle f*g\in {\rm {L}}^{1}(\mathbb {R} ^{N})} et f g 1 f 1 g 1 {\displaystyle \|f*g\|_{1}\leq \|f\|_{1}\cdot \|g\|_{1}}  ;
  • si f L 1 ( R N ) {\displaystyle f\in {\rm {L}}^{1}(\mathbb {R} ^{N})} et g L 2 ( R N ) {\displaystyle g\in {\rm {L}}^{2}(\mathbb {R} ^{N})} , alors f g L 2 ( R N ) {\displaystyle f*g\in {\rm {L}}^{2}(\mathbb {R} ^{N})} et f g 2 f 1 g 2 {\displaystyle \|f*g\|_{2}\leq \|f\|_{1}\cdot \|g\|_{2}}  ;
  • si f , g L 2 ( R N ) {\displaystyle f,g\in {\rm {L}}^{2}(\mathbb {R} ^{N})} , alors f g L ( R N ) {\displaystyle f*g\in {\rm {L}}^{\infty }(\mathbb {R} ^{N})} et f g f 2 g 2 {\displaystyle \|f*g\|_{\infty }\leq \|f\|_{2}\cdot \|g\|_{2}} .

Ainsi :

  • si f , g L 1 ( R N ) {\displaystyle f,g\in {\rm {L}}^{1}(\mathbb {R} ^{N})} , alors F ( f g ) = F ( f ) F ( g ) {\displaystyle {\mathcal {F}}(f*g)={\mathcal {F}}(f)\,\cdot \,{\mathcal {F}}(g)}  ;
  • par densité, cette égalité tient encore si f L 1 {\displaystyle f\in {\rm {L}}^{1}} et g L 2 {\displaystyle g\in {\rm {L}}^{2}}  ;
  • Si f , g L 2 ( R N ) {\displaystyle f,g\in {\rm {L}}^{2}(\mathbb {R} ^{N})} , alors f g = F 1 [ F ( f ) F ( g ) ] {\displaystyle f\ast g={\mathcal {F}}^{-1}[{\mathcal {F}}(f)\,\cdot \,{\mathcal {F}}(g)]}  ; de plus, l'égalité F ( f g ) = F ( f ) F ( g ) {\displaystyle {\mathcal {F}}(f*g)={\mathcal {F}}(f)\,\cdot \,{\mathcal {F}}(g)} est vraie si f g L 1 {\displaystyle f*g\in {\rm {L}}^{1}} .

Principe d'incertitude

Remarque : ce paragraphe utilise la définition fréquentielle de la transformée de Fourier.

Article détaillé : Principe d'incertitude.

On peut remarquer que les répartitions d'une fonction et de sa transformée de Fourier ont des comportements opposés : plus la masse de f(x) est « concentrée », plus celle de la transformée est étalée, et inversement. Il est en fait impossible de concentrer à la fois la masse d'une fonction et celle de sa transformée.

Ce compromis entre la compaction d'une fonction et celle de sa transformée de Fourier peut se formaliser par un principe d'incertitude en considérant une fonction et sa transformée de Fourier comme des variables conjuguées par la forme symplectique sur le domaine temps-fréquence : par la transformation canonique linéaire, la transformation de Fourier est une rotation de 90° dans le domaine temps–fréquence qui préserve la forme symplectique.

Supposons f intégrable et de carré intégrable. Sans perte de généralité, on supposera f normalisée :

| f ( x ) | 2 d x = 1 {\displaystyle \int _{-\infty }^{\infty }|f(x)|^{2}\,{\rm {d}}x=1} .

Par le théorème de Plancherel, on sait que f ^ ( ν ) {\displaystyle {\hat {f}}(\nu )} est également normalisée.

On peut mesurer la répartition autour d'un point (x = 0 sans perte de généralité) par :

D 0 ( f ) = x 2 | f ( x ) | 2 d x {\displaystyle D_{0}(f)=\int _{-\infty }^{\infty }x^{2}|f(x)|^{2}\,{\rm {d}}x} .

De même pour la fréquence autour du point ν = 0 {\displaystyle \nu =0}  :

D 0 ( f ^ ) = ν 2 | f ^ ( ν ) | 2 d ν {\displaystyle D_{0}({\hat {f}})=\int _{-\infty }^{\infty }\nu ^{2}|{\hat {f}}(\nu )|^{2}\,{\rm {d}}\nu } .

En probabilités, il s'agit des moments d'ordre 2 de |f|2 et de | f ^ | 2 {\displaystyle |{\hat {f}}|^{2}} .

Le principe d'incertitude dit que si f(x) est absolument continue et que les fonctions x·f(x) et f′(x) sont de carrés intégrables, on a alors[4] :

D 0 ( f ) D 0 ( f ^ ) 1 16 π 2 {\displaystyle D_{0}(f)D_{0}({\hat {f}})\geq {\frac {1}{16\pi ^{2}}}} .

Cette inégalité est aussi connue sous le nom d'inégalité de Heisenberg-Gabor ou simplement inégalité de Heisenberg par son utilisation répandue en mécanique quantique.

L'égalité n'est atteinte que pour f ( x ) = C 1 e π x 2 / σ 2 {\displaystyle f(x)=C_{1}\,{\rm {e}}^{{-\pi x^{2}}/{\sigma ^{2}}}} (alors f ^ ( ξ ) = σ C 1 e π σ 2 ξ 2 {\displaystyle {\hat {f}}(\xi )=\sigma C_{1}\,{\rm {e}}^{-\pi \sigma ^{2}\xi ^{2}}} ) pour σ > 0 arbitraire et C1 telle que f est L2–normalisée, soit, si f est une fonction gaussienne (normalisée) centrée en 0 et de variance σ2, et sa transformée de Fourier est une gaussienne de variance σ–2.

Transformation de Fourier sur l'espace de Schwartz

L'espace de Schwartz S ( R n ) {\displaystyle {\mathcal {S}}(\mathbb {R} ^{n})} est l'espace des fonctions f de classe C sur R n {\displaystyle \mathbb {R} ^{n}} , telles que f et toutes ses dérivées soient à décroissance rapide. C'est un sous-espace vectoriel de L1, donc pour lequel la transformée de Fourier est définie. Ces fonctions sont à la fois temporellement et fréquentiellement à décroissance exponentielle. L'intérêt de la classe de Schwartz résulte de la propriété d'échange entre régularité et décroissance à l'infini qu'opère la transformée de Fourier.

  • Toute fonction de Schwartz est de classe C avec des dérivées toutes intégrables. On en déduit que sa transformée de Fourier est à décroissance rapide.
  • Toute fonction de Schwartz est à décroissance rapide. On en déduit que sa transformée de Fourier est de classe C.

Ainsi, on visualise intuitivement pourquoi l'espace de Schwartz est invariant par transformation de Fourier. Cet espace est donc très commode pour l'utilisation de cette dernière. De plus, l'espace de Schwartz est dense dans L1 et dans L2, et pourrait donc servir de base pour la définition de la transformation de Fourier sur ces espaces.

Formule d'inversion de Fourier sur S ( R n ) {\displaystyle {\mathcal {S}}(\mathbb {R} ^{n})}  — 

La transformée de Fourier induit un automorphisme bicontinu de l'espace de Schwartz sur lui-même, dont l'inverse est défini par
( F 1 ϕ ) ( x ) = ( F ϕ ) ( x ) = R n ϕ ( ξ ) e 2 i π x ξ d ξ {\displaystyle ({\mathcal {F}}^{-1}\phi )(x)=({\mathcal {F}}\phi )(-x)=\int _{\mathbb {R} ^{n}}\phi (\xi )\,\mathrm {e} ^{2{\rm {i}}\pi x\cdot \xi }\,\mathrm {d} \xi } .

Remarque : cette formule dépend de la convention choisie pour la transformation de Fourier dans l'espace des fonctions. Elle est valide pour une transformation de Fourier exprimée dans l'espace des fréquences, dont la définition utilise e i 2 π ξ x {\displaystyle {\rm {e}}^{-{\rm {i}}2\pi \xi \cdot x}} .

Démonstration de la formule d'inversion
  • Prouvons d'abord que S {\displaystyle {\mathcal {S}}} est stable par F {\displaystyle {\mathcal {F}}} . Par commodité, nous ne traiterons que le cas n = 1, mais le cas quelconque se traite de manière similaire. Soit donc f S ( R ) {\displaystyle f\in {\mathcal {S}}(\mathbb {R} )} .
  1. D'une part, la décroissance rapide implique que pour tout entier naturel n, x x n f ( x ) {\displaystyle x\mapsto x^{n}f(x)} est intégrable. La fonction F f {\displaystyle {\mathcal {F}}f} est donc définie et C.
  2. D'autre part, pour tout couple d'entiers naturels (n,k), la fonction [ x ( 2 i π x ) n f ( x ) ] ( k ) {\displaystyle [x\mapsto (-2{\rm {i}}\pi x)^{n}f(x)]^{(k)}} est dans S {\displaystyle {\mathcal {S}}} , donc dans L1. Sa transformée de Fourier tend vers 0 à l'infini. Or, en appliquant les propriétés d'échange entre multiplication par un polynôme et dérivation, F ( [ x ( 2 i π x ) n f ( x ) ] ( k ) ) ( ξ ) = ( 2 i π ξ ) k [ F ( x ( 2 i π x ) n f ( x ) ) ] ( ξ ) = ( 2 i π ξ ) k [ F ( f ) ] ( n ) ( ξ ) {\displaystyle {\mathcal {F}}\left([x\mapsto (-2{\rm {i}}\pi x)^{n}f(x)]^{(k)}\right)(\xi )=(2{\rm {i}}\pi \xi )^{k}\left[{\mathcal {F}}\left(x\mapsto (-2{\rm {i}}\pi x)^{n}f(x)\right)\right](\xi )=(2{\rm {i}}\pi \xi )^{k}[{\mathcal {F}}(f)]^{(n)}(\xi )} ,ce qui prouve la décroissance rapide de F f {\displaystyle {\mathcal {F}}f} ainsi que toutes ses dérivées successives. Elle satisfait donc aux conditions d'appartenance à S {\displaystyle {\mathcal {S}}} .
  • Soit f un élément de S {\displaystyle {\mathcal {S}}} donc de L1. D'après le point précédent, F f {\displaystyle {\mathcal {F}}f} appartient aussi à L1. Le théorème d'inversion sur L1 s'applique et donne, en notant ~ {\displaystyle {\tilde {}}} l'opérateur de composition par –Id : f = (   ~ F F ) ( f ) = ( F ~ F ) ( f ) {\displaystyle f=\left(~{\tilde {}}\circ {\mathcal {F}}\circ {\mathcal {F}}\right)(f)=\left({\mathcal {F}}\circ {\tilde {}}\circ {\mathcal {F}}\right)(f)} , ce qui prouve que F : S S {\displaystyle {\mathcal {F}}:{\mathcal {S}}\to {\mathcal {S}}} est bijectif et que son inverse est ~ F {\displaystyle {\tilde {}}\circ {\mathcal {F}}} .

Transformation de Fourier pour les distributions tempérées

On définit la transformée de Fourier d'une distribution tempérée T S ( R n ) {\displaystyle T\in {\mathcal {S}}'(\mathbb {R} ^{n})} comme la distribution définie via son crochet de dualité par

ϕ S ( R n ) F T , ϕ = T , F ϕ {\displaystyle \forall \phi \in {\mathcal {S}}(\mathbb {R} ^{n})\quad \langle {\mathcal {F}}T,\phi \rangle =\langle T,{\mathcal {F}}\phi \rangle } .

De même que sur S {\displaystyle {\mathcal {S}}} , l'opérateur F {\displaystyle {\mathcal {F}}} ainsi défini sur S {\displaystyle {\mathcal {S}}'} est un automorphisme bicontinu.

Les détails et des exemples ne sont pas donnés ici, mais figurent dans l'article relatif aux distributions tempérées.

Remarquons que l'expression de la transformée de Fourier d'une fonction f ressemble au produit scalaire dans L 2 ( C ) , ( f , g ) L 2 := f g ¯ {\displaystyle {\rm {L}}^{2}(\mathbb {C} ),(f,g)_{L^{2}}:=\int f{\bar {g}}} entre f et la conjuguée de e 2 π ξ : x e i 2 π ξ x {\displaystyle e_{2\pi \xi }:x\mapsto {\rm {e}}^{{\rm {i}}2\pi \xi \cdot x}} . Sauf que ( f , e 2 π ξ ) L 2 {\displaystyle (f,e_{2\pi \xi })_{L^{2}}} n'a pas de sens car e2 π ξ n'est pas dans L2. Mais le crochet de dualité des distributions T f , e 2 π ξ {\displaystyle \langle T_{f},e_{2\pi \xi }\rangle } , qui pour les fonctions coïncide avec le produit scalaire de L2, donne sens à cette formulation en tant que produit scalaire.

Cette généralisation va bien plus loin car l'espace des distributions tempérées S ( R n ) {\displaystyle {\mathcal {S}}'(\mathbb {R} ^{n})} englobe les différents objets sur lesquels la transformée de Fourier a été définie : fonctions de R n {\displaystyle \mathbb {R} ^{n}} sommables ou de carré sommable, fonctions de R n {\displaystyle \mathbb {R} ^{n}} périodiques localement sommables ou localement de carré sommable, suites discrètes sommables, suites discrètes périodiques. La transformée de Fourier sur S ( R n ) {\displaystyle {\mathcal {S}}'(\mathbb {R} ^{n})} unifie et généralise les différentes définitions des transformées avec l'unique formalisme des distributions. Nous allons montrer que la transformée de Fourier sur S {\displaystyle {\mathcal {S}}'} généralise les notions d'intégrales de Fourier et de séries de Fourier, en analysant successivement ces espaces.

Compatibilités

Compatibilité avec les espaces de fonctions

Les fonctions intégrables et les fonctions de carré sommable définissent des distributions tempérées. Montrons que les deux notions possibles de transformée de Fourier coïncident dans le cas L1, puis utilisons cette compatibilité pour l'établir dans le cas L2.

Compatibilité avec L1 et L2 — Soit

  • f L 1 ( R n ) {\displaystyle f\in {\rm {L}}^{1}(\mathbb {R} ^{n})} et f ^ C 0 ( R n ) {\displaystyle {\hat {f}}\in C_{0}(\mathbb {R} ^{n})} sa transformée de Fourier,

ou bien

  • f L 2 ( R n ) {\displaystyle f\in {\rm {L}}^{2}(\mathbb {R} ^{n})} et f ^ L 2 ( R n ) {\displaystyle {\hat {f}}\in {\rm {L}}^{2}(\mathbb {R} ^{n})} sa transformée de Fourier.
Dans ces deux cas, f ^ {\displaystyle {\hat {f}}} définit une distribution tempérée égale à la transformée de Fourier de Tf, c'est-à-dire
F T f = T f ^ {\displaystyle {\mathcal {F}}T_{f}=T_{\hat {f}}} .
Démonstration
  • Soit f L 1 {\displaystyle f\in {\rm {L}}^{1}} . Il s'agit de vérifier que pour tout ϕ S {\displaystyle \phi \in {\mathcal {S}}} ,
R n f ( x ) ϕ ^ ( x ) d x = R n f ^ ( y ) ϕ ( y ) d y {\displaystyle \int _{\mathbb {R} ^{n}}f(x){\hat {\phi }}(x){\rm {d}}x=\int _{\mathbb {R} ^{n}}{\hat {f}}(y)\phi (y){\rm {d}}y} ,

c'est-à-dire

R n f ( x ) ( R n ϕ ( y ) e i x y d y ) d x = R n ( R n f ( x ) e i x y d x ) ϕ ( y ) d y {\displaystyle \int _{\mathbb {R} ^{n}}f(x)\left(\int _{\mathbb {R} ^{n}}\phi (y){\rm {e}}^{-{\rm {i}}x\cdot y}{\rm {d}}y\right){\rm {d}}x=\int _{\mathbb {R} ^{n}}\left(\int _{\mathbb {R} ^{n}}f(x){\rm {e}}^{-{\rm {i}}x\cdot y}{\rm {d}}x\right)\phi (y){\rm {d}}y} .

Cela résulte simplement du théorème de Fubini, appliqué à la fonction intégrable ( x , y ) f ( x ) ϕ ( y ) e i x y {\displaystyle (x,y)\mapsto f(x)\phi (y){\rm {e}}^{-{\rm {i}}x\cdot y}} .

  • Les deux applications continues f F T f {\displaystyle f\mapsto {\mathcal {F}}T_{f}} et f T f ^ {\displaystyle f\mapsto T_{\hat {f}}} , de L2 dans l'espace séparé S {\displaystyle {\mathcal {S}}'} , sont égales car elles coïncident, d'après le point précédent, sur le sous-espace dense L1 ∩ L2.

Enfin, les fonctions périodiques intégrables sur une période sont exactement les fonctions à la fois périodiques et localement intégrables, et donc définissent des distributions régulières.

Compatibilité avec L1per — La transformée de Fourier d'une distribution régulière Tf définie par une fonction T-périodique f L 1 ( [ 0 , T [ ) {\displaystyle f\in {\rm {L}}^{1}([0,T[)} , est la distribution à support discret correspondant à la suite de ses coefficients de Fourier :

F T f = n Z c n ( f ) δ n a v e c c n ( f ) = 0 T f ( x ) e i 2 π T n x d x {\displaystyle {\mathcal {F}}T_{f}=\sum _{n\in \mathbb {Z} }c_{n}(f)\delta _{n}\quad {\rm {avec}}\quad c_{n}(f)=\int _{0}^{T}f(x){\rm {e}}^{-{\rm {i}}{\frac {2\pi }{T}}nx}\,\mathrm {d} x} .

Le résultat énoncé ne concerne que les fonctions périodiques de la variable réelle mais s'étend facilement aux fonctions périodiques sur un réseau de ℝN. Comme la transformation de Fourier[Laquelle ?] possède une réciproque définie sur le même domaine, elle est de ce fait bijective, alors la démonstration de ce résultat sera une conséquence du théorème sur les distributions périodiques[Lequel ?].

Compatibilité avec les espaces de suites

Les suites, c'est-à-dire les signaux discrets, peuvent parfois s'exprimer comme des distributions sur ℝ à support dans ℤ. À une suite donnée a := ( a n ) n Z {\displaystyle a:=(a_{n})_{n\in \mathbb {Z} }} correspond en effet de manière unique une série de masses de Dirac T a := k Z a k δ k {\displaystyle T_{a}:=\sum _{k\in \mathbb {Z} }a_{k}\delta _{k}} . Lorsque cette suite est sommable, cette série de masses de Dirac a un sens en tant que distribution tempérée d'ordre 0.

Compatibilité de F {\displaystyle {\mathcal {F}}} avec l1 — Soit une suite sommable à valeurs complexes notée a := ( a n ) n Z {\displaystyle a:=(a_{n})_{n\in \mathbb {Z} }} . Sa transformée de Fourier à temps discret est une fonction 1-périodique qui coïncide avec la transformée de Fourier de la série de masses de Dirac associée à a.

T F T D [ ( a n ) n Z ] = F ( k Z a k δ k ) {\displaystyle \mathbf {TFTD} [(a_{n})_{n\in \mathbb {Z} }]={\mathcal {F}}\left(\sum _{k\in \mathbb {Z} }a_{k}\delta _{k}\right)} .
Démonstration de la compatibilité de F {\displaystyle {\mathcal {F}}} avec l1

Lorsque a est sommable, la somme T a := k Z a k δ k {\displaystyle T_{a}:=\sum _{k\in \mathbb {Z} }a_{k}\delta _{k}} définit bien une distribution d'ordre 0. En effet, pour une fonction test ϕ S ( R ) {\displaystyle \phi \in {\mathcal {S}}(\mathbb {R} )} , k Z a k δ k , ϕ = k Z ( a k ϕ ( k ) ) | k Z a k δ k , ϕ | ( k Z | a k | ) ϕ . {\displaystyle {\begin{array}{lcl}\langle \sum _{k\in \mathbb {Z} }a_{k}\delta _{k},\phi \rangle &=&\sum _{k\in \mathbb {Z} }(a_{k}\phi (k))\\\left|\langle \sum _{k\in \mathbb {Z} }a_{k}\delta _{k},\phi \rangle \right|&\leq &\left(\sum _{k\in \mathbb {Z} }|a_{k}|\right)\|\phi \|_{\infty }.\end{array}}}

Par continuité de la transformation de Fourier et formule de la transformée du Dirac F ( δ k ) = e 2 π k {\displaystyle {\mathcal {F}}(\delta _{k})=e_{-2\pi k}} ,

F T a = ξ k Z a k e i 2 π k ξ = T F T D [ a ] {\displaystyle {\mathcal {F}}T_{a}=\xi \mapsto \sum _{k\in \mathbb {Z} }a_{k}\mathrm {e} ^{-\mathrm {i} 2\pi k\xi }=\mathbf {TFTD} [a]} .

On retrouve bien la transformée de Fourier en temps discret.

Par densité, la démonstration s'étend aux séries de carré sommable. Notons en outre que la transformation de Fourier des distributions périodiques donne une définition de la transformée de Fourier discrète de suites non nécessairement sommables : les suites à croissance polynomiale.

En particulier, la transformée de Fourier discrète (TFD) s'interprète également comme la transformée d'une distribution tempérée. En effet, une suite finie de N points { x k } k = 0 N 1 {\displaystyle \lbrace x_{k}\rbrace _{k=0}^{N-1}} s'identifie de manière unique avec une suite N-périodique obtenue par périodisation, c'est-à-dire convolution avec un peigne de Dirac.

Compatibilité de F {\displaystyle {\mathcal {F}}} avec la TFD — La TFD d'une suite x(•) à l'ordre N est la transformée de Fourier de la distribution à support dans ℤ obtenue par périodisation de x(•) à la période N, c'est-à-dire convolution par un peigne de Dirac WN :

T F D N [ x ( . ) ] ( k ) = F [ x ~ ( . ) ] ( k ) {\displaystyle \mathbf {TFD} _{N}[x(.)](k)={\mathcal {F}}[{\tilde {x}}(.)](k)} avec x ~ = W N ( n = 0 N 1 x ( n ) δ n ) {\displaystyle {\tilde {x}}=W_{N}\ast \left(\sum _{n=0}^{N-1}x(n)\delta _{n}\right)} .

Signaux discrets et signaux périodiques

Nous pouvons retenir que formellement, la transformée de Fourier échange discrétisation et périodisation.

  • Le spectre d'un signal discret x[•] obtenu par échantillonnage à la période T présente un spectre périodique, résultant de la périodisation du spectre du signal continu :
T F T D ( x [ . ] ) = F ( x ( . ) ) W 2 π T {\displaystyle \mathbf {TFTD} (x[.])={\mathcal {F}}(x(.))\ast W_{\frac {2\pi }{T}}} .

Si la multiplication n'est pas définie entre distribution, on donne dans le cas du peigne un sens à x [ . ] = x ( . ) W T {\displaystyle x[.]=x(.)\cdot W_{T}} , et la formulation de convolution est encore vérifiée : F ( x ( . ) W T ) = F ( x ( . ) ) F ( W T ) {\displaystyle {\mathcal {F}}(x(.)\cdot W_{T})={\mathcal {F}}(x(.))\ast {\mathcal {F}}(W_{T})} .

  • Le spectre d'un signal T-périodique xT(•), c'est-à-dire la somme de sa série de Fourier, est celui obtenu par discrétisation du spectre du signal tronqué sur une seule période.
F ( x T ( . ) ) = F ( x ( . ) ) W T {\displaystyle {\mathcal {F}}(x_{T}(.))={\mathcal {F}}(x(.))\cdot W_{T}} avec x = x T .1 [ 0 , T ] {\displaystyle x=x_{T}.1_{[0,T]}} .

Liens avec d'autres transformations

Lien avec les transformations de Laplace

La transformée de Fourier d'une fonction f est un cas particulier de la transformée bilatérale de Laplace de cette même fonction définie par : L b i l { f } ( p ) = + f ( t ) e p t d t {\displaystyle {\mathcal {L}}_{bil}\{f\}(p)=\int _{-\infty }^{+\infty }f(t)\,{\rm {e}}^{-pt}{\rm {d}}t} avec p C {\displaystyle p\in \mathbb {C} } .

On constate alors que F { f } ( ξ ) = L b i l { f } ( i ξ ) {\displaystyle {\mathcal {F}}\{f\}(\xi )={\mathcal {L}}_{bil}\{f\}({\rm {i}}\xi )} .

On peut également écrire ce lien en utilisant la transformée de Laplace « usuelle » par :

F { f } ( ξ ) = L { f + } ( + i ξ ) + L { f } ( i ξ ) {\displaystyle {\mathcal {F}}\{f\}(\xi )={\mathcal {L}}\{f^{+}\}(+{\rm {i}}\xi )+{\mathcal {L}}\{f^{-}\}(-{\rm {i}}\xi )} [réf. nécessaire]

où les fonctions f + et f sont définies par :

f + ( t ) = f ( + t ) ,    si  t 0 , 0   sinon . {\displaystyle f^{+}(t)=f(+t),\ {\text{ si }}t\geq 0,0\ {\text{sinon}}.}
f ( t ) = f ( t ) ,    si  t 0 , 0   sinon . {\displaystyle f^{-}(t)=f(-t),\ {\text{ si }}t\leq 0,0\ {\text{sinon}}.}

Lien avec les séries de Fourier

Parallèle formel

La transformée de Fourier est définie de façon semblable : la variable d'intégration x est remplacée par nΔx, n étant l'indice de sommation, et l'intégrale par la somme. On a alors

f ^ ( k ) = Δ t n = f ( n ) e i 2 π k n Δ t {\displaystyle {\hat {f}}(k)=\Delta t\sum _{n=-\infty }^{\infty }f(n){\rm {e}}^{-{\rm {i}}2\pi kn\Delta t}} .

On trouvera quelques remarques à ce sujet dans Analyse spectrale.

Lien direct

Cependant, comme indiqué par l'étude théorique dans la section précédente, un lien direct entre séries et transformées de Fourier est possible par la théorie des distributions. En reprenant de façon plus pratique l'exposé précédent, la transformée de Fourier (définition fréquentielle) d'une fonction périodique f de période T est un peigne de Dirac de période fréquentielle ν T = 1 / T {\displaystyle \nu _{T}=1/T} , modulé par des coefficients complexes cn :

f ^ ( ν ) = n = + c n δ ( ν n ν T ) , {\displaystyle {\hat {f}}(\nu )=\sum _{n=-\infty }^{+\infty }c_{n}\delta \left(\nu -n\nu _{T}\right),}

où les cn sont précisément les coefficients de la série de Fourier (complexe) de f. Pour le voir, il suffit de vérifier que la formule de transformation inverse de f ^ ( ν ) {\displaystyle {\hat {f}}(\nu )} (définition en fréquences) donne précisément la série de Fourier de f, et donc qu'elle est égale à f presque partout (en supposant que la série de Fourier de f converge).

Cela permet d'unifier le formalisme des séries de Fourier avec celui de la transformation de Fourier.

Avec la définition standard de la transformée de Fourier, il faut remplacer la formule précédente par :

f ^ ( ξ ) = 2 π n = + c n δ ( ξ 2 π n T ) . {\displaystyle {\hat {f}}(\xi )=2\pi \sum _{n=-\infty }^{+\infty }c_{n}\delta \left({\xi \over 2\pi }-{\frac {n}{T}}\right).}

Avec la définition pulsatoire, et en notant la pulsation de f {\displaystyle f} par ω T = 2 π / T {\displaystyle \omega _{T}=2\pi /T} , elle devient

f ^ ( ω ) = 2 π n = + c n δ ( ω n ω T 2 π ) . {\displaystyle {\hat {f}}(\omega )={\sqrt {2\pi }}\sum _{n=-\infty }^{+\infty }c_{n}\delta \left({\omega -n\omega _{T} \over 2\pi }\right).}

Par exemple, après quelques manipulations, on a les transformées de Fourier fréquentielles suivantes :

  • F e 2 i π ν T t = δ ( ν ν T ) {\displaystyle {\cal {F}}e^{2i\pi \nu _{T}t}=\delta (\nu -\nu _{T})} (Dirac décalé) ;
  • F cos ( 2 π ν T t ) = 1 2 [ δ ( ν ν T ) + δ ( ν + ν T ) ] {\displaystyle {\cal {F}}\cos(2\pi \nu _{T}t)={1 \over 2}\left[\delta (\nu -\nu _{T})+\delta (\nu +\nu _{T})\right]}
  • F sin ( 2 π ν T t ) = 1 2 i [ δ ( ν ν T ) δ ( ν + ν T ) ] {\displaystyle {\cal {F}}\sin(2\pi \nu _{T}t)={1 \over 2i}[\delta (\nu -\nu _{T})-\delta (\nu +\nu _{T})]}

Il y a encore une formule utile qui donne les coefficients c n {\displaystyle c_{n}} de la série de Fourier d'une fonction périodique f dès que l'on connait la transformation de Fourier de sa « restriction » g = 1 τ f {\displaystyle g=1_{\tau }f} à la période centrale τ = [ T / 2 ,   T / 2 ] {\displaystyle \tau =[-T/2,\ T/2]} ( g ^ {\displaystyle {\hat {g}}} existe nécessairement si f est localement intégrable puisque τ {\displaystyle \tau } est compacte). En effet, par comparaison de la formule des coefficients de la série de Fourier de f {\displaystyle f} avec celle donnant la transformée de Fourier inverse de g {\displaystyle g} , on obtient facilement, pour la définition fréquentielle, que

c n = 1 T g ^ ( n ν T ) = ν T g ^ ( n ν T ) . {\displaystyle c_{n}={\frac {1}{T}}{\hat {g}}(n\nu _{T})=\nu _{T}{\hat {g}}(n\nu _{T}).}

Pour la définition standard de la transformée de Fourier, cette formule devient c n = 1 T g ^ ( n ω T ) , {\displaystyle c_{n}={1 \over T}{\hat {g}}(n\omega _{T}),} avec ω T = 2 π / T {\displaystyle \omega _{T}=2\pi /T} , et pour la définition pulsatoire, elle devient c n = 1 2 π ω T g ^ ( n ω T ) . {\displaystyle c_{n}={1 \over {\sqrt {2\pi }}}\omega _{T}{\hat {g}}(n\omega _{T}).}

Cette formule permet l'utilisation de l'imposante machinerie disponible pour la transformation de Fourier (convolution, décalage, produit, distributions, tables, etc.) pour le calcul des coefficients de Fourier d'une fonction périodique. On peut ainsi facilement obtenir la série de Fourier de trains d'ondes pulsées de forme carrée, triangulaire, demi-sinusoïdale, etc.

Par exemple, quelle est la série de Fourier correspondant à un train de pulses étroits, de masse 1 et de période T grande relativement à la durée des pulses ? Approximons chaque pulse par un Dirac δ {\displaystyle \delta } . La transformée de Fourier fréquentielle de δ {\displaystyle \delta } est la fonction identiquement égale à 1 (voir table ci-dessous). Donc la formule précédente donne c n = 1 T . {\displaystyle c_{n}={1 \over T}.} Ainsi la série de Fourier du train de pulses est

f ( x ) = n = + 1 T e i 2 π n T x = 1 T + n = 1 2 T cos ( 2 π n T x ) {\displaystyle f(x)=\sum _{n=-\infty }^{+\infty }{1 \over T}{\rm {e}}^{{\rm {i}}2\pi {\tfrac {n}{T}}x}={1 \over T}+\sum _{n=1}^{\infty }{2 \over T}\cos(2\pi {\tfrac {n}{T}}x)}

(au sens des distributions).

Autre interprétation

Comme on l'a vu plus haut, il est d'autre part possible d'interpréter l'intégrale de la transformée de Fourier comme une somme finie de n oscillateurs harmoniques, où n est un entier non standard[5] ; cela revient à identifier (en un sens différent) la transformation de Fourier aux coefficients d'une série de Fourier.

Transformée

On utilise les variables normalisées suivantes : F = f f e = f Δ t = f | Δ t = 1 {\displaystyle F={f \over f_{e}}=f\Delta t=f|_{\Delta t=1}} , Ω = 2 π F = 2 π f Δ t = ω Δ t = ω | Δ t = 1 {\displaystyle \Omega =2\pi F=2\pi f\Delta t=\omega \Delta t=\omega |_{\Delta t=1}} .

Transformation de Fourier (analyse) Transformation inverse (synthèse)
X ( f ) = Δ t n = x ( n ) e i 2 π f n Δ t {\displaystyle X(f)=\Delta t\sum _{n=-\infty }^{\infty }x(n){\rm {e}}^{-{\rm {i}}2\pi fn\Delta t}} x ( n ) = f e X ( f ) e i 2 π f n Δ t d f {\displaystyle x(n)=\int _{f_{e}}X(f){\rm {e}}^{{\rm {i}}2\pi fn\Delta t}{\rm {d}}f}
X ( w ) = Δ t n = x ( n ) e i ω n Δ t {\displaystyle X(w)=\Delta t\sum _{n=-\infty }^{\infty }x(n){\rm {e}}^{-{\rm {i}}\omega n\Delta t}} x ( n ) = 1 2 π ω 2 = 2 π f e X ( w ) e i w n Δ t d w {\displaystyle x(n)={1 \over 2\pi }\int _{\omega _{2}=2\pi f_{e}}X(w){\rm {e}}^{{\rm {i}}wn\Delta t}{\rm {d}}w}
X ( F ) = n = x ( n ) e i 2 π n F {\displaystyle X(F)=\sum _{n=-\infty }^{\infty }x(n){\rm {e}}^{-{\rm {i}}2\pi nF}} x ( n ) = 1 X ( f ) e i 2 π n F d F {\displaystyle x(n)=\int _{1}X(f){\rm {e}}^{{\rm {i}}2\pi nF}{\rm {d}}F\,\!}
X ( Ω ) = n = x ( n ) e i n Ω {\displaystyle X(\Omega )=\sum _{n=-\infty }^{\infty }x(n){\rm {e}}^{-{\rm {i}}n\Omega }} x ( n ) = 1 2 π 2 π X ( Ω ) e i n Ω d Ω {\displaystyle x(n)={1 \over 2\pi }\int _{2\pi }X(\Omega ){\rm {e}}^{{\rm {i}}n\Omega }{\rm {d}}\Omega }

Généralisation

La transformée de Fourier se généralise pratiquement telle quelle aux groupes abéliens localement compacts, grâce à la dualité de Pontryagin.

En traitement d'images, on effectue des transformations de Fourier à deux dimensions : si f est une fonction de ℝ2 dans ℝ, sa transformée de Fourier est définie par :

f ^ ( u , v ) = f ( x , y ) e i ( u x + v y ) d x d y {\displaystyle {\hat {f}}(u,v)=\int _{-\infty }^{\infty }\int _{-\infty }^{\infty }f(x,y)\mathrm {e} ^{-{\rm {i}}(ux+vy)}\,\mathrm {d} x\mathrm {d} y} .

Tables des principales transformées de Fourier

Les tableaux suivants présentent les transformations de Fourier de certaines fonctions. Les transformées de Fourier de f (x), g(x) et h(x) sont notées respectivement , ĝ et ĥ. N'apparaissent que les trois conventions les plus courantes. Il peut être utile de noter que l'entrée sur la dualité indique une relation entre la transformée de Fourier d'une fonction et la fonction d'origine, ce qui peut être considéré comme une relation entre la transformation de Fourier et son inverse.

Relations entre fonctions à une variable

Les transformées de Fourier de ce tableau sont traitées dans (en) Arthur Erdélyi, Tables of Integral Transforms, Vol. 1, McGraw-Hill, ou (en) David Kammler, A First Course in Fourier Analysis, USA, Prentice Hall, .

Fonction Transformée de Fourier

ξ est la fréquence

Transformée de Fourier
ω = 2πξ est la pulsation ou fréquence angulaire
Transformée de Fourier

définition alternative

Remarques
f ( x ) {\displaystyle f(x)\,} f ^ ( ξ ) = f ( x ) e 2 π i x ξ d x {\displaystyle {\begin{aligned}&{\hat {f}}(\xi )\\&=\int _{-\infty }^{\infty }f(x){\rm {e}}^{-2\pi {\rm {i}}x\xi }\,\mathrm {d} x\end{aligned}}} f ^ ( ω ) = 1 2 π f ( x ) e i ω x d x {\displaystyle {\begin{aligned}&{\hat {f}}(\omega )\\&={\frac {1}{\sqrt {2\pi }}}\int _{-\infty }^{\infty }f(x){\rm {e}}^{-{\rm {i}}\omega x}\,\mathrm {d} x\end{aligned}}} f ^ ( ν ) = f ( x ) e i ν x d x {\displaystyle {\begin{aligned}&{\hat {f}}(\nu )\\&=\int _{-\infty }^{\infty }f(x){\rm {e}}^{-{\rm {i}}\nu x}\,\mathrm {d} x\end{aligned}}} Définition
a f ( x ) + b g ( x ) {\displaystyle a\cdot f(x)+b\cdot g(x)\,} a f ^ ( ξ ) + b g ^ ( ξ ) {\displaystyle a\cdot {\hat {f}}(\xi )+b\cdot {\hat {g}}(\xi )\,} a f ^ ( ω ) + b g ^ ( ω ) {\displaystyle a\cdot {\hat {f}}(\omega )+b\cdot {\hat {g}}(\omega )\,} a f ^ ( ν ) + b g ^ ( ν ) {\displaystyle a\cdot {\hat {f}}(\nu )+b\cdot {\hat {g}}(\nu )\,} Linéarité
f ( x a ) {\displaystyle f(x-a)\,} e 2 π i a ξ f ^ ( ξ ) {\displaystyle {\rm {e}}^{-2\pi {\rm {i}}a\xi }{\hat {f}}(\xi )\,} e i a ω f ^ ( ω ) {\displaystyle {\rm {e}}^{-{\rm {i}}a\omega }{\hat {f}}(\omega )\,} e i a ν f ^ ( ν ) {\displaystyle {\rm {e}}^{-{\rm {i}}a\nu }{\hat {f}}(\nu )\,} Décalage dans le temps
f ( x ) e i a x {\displaystyle f(x){\rm {e}}^{{\rm {i}}ax}\,} f ^ ( ξ a 2 π ) {\displaystyle {\hat {f}}\left(\xi -{\frac {a}{2\pi }}\right)\,} f ^ ( ω a ) {\displaystyle {\hat {f}}(\omega -a)\,} f ^ ( ν a ) {\displaystyle {\hat {f}}(\nu -a)\,} Décalage dans le domaine des fréquences, relation duale de la formule du décalage dans le temps
f ( a x ) {\displaystyle f(ax)\,} 1 | a | f ^ ( ξ a ) {\displaystyle {\frac {1}{|a|}}{\hat {f}}\left({\frac {\xi }{a}}\right)\,} 1 | a | f ^ ( ω a ) {\displaystyle {\frac {1}{|a|}}{\hat {f}}\left({\frac {\omega }{a}}\right)\,} 1 | a | f ^ ( ν a ) {\displaystyle {\frac {1}{|a|}}{\hat {f}}\left({\frac {\nu }{a}}\right)\,} Changement d'échelle des temps.
Si | a | {\displaystyle |a|} est grand, alors f (ax) est resserré autour de 0 et 1 | a | f ^ ( ω a ) {\displaystyle {\frac {1}{|a|}}{\hat {f}}\left({\frac {\omega }{a}}\right)\,} s'étale et s’aplatit.
f ^ ( x ) {\displaystyle {\hat {f}}(x)\,} f ( ξ ) {\displaystyle f(-\xi )\,} f ( ω ) {\displaystyle f(-\omega )\,} 2 π f ( ν ) {\displaystyle 2\pi f(-\nu )\,} Dualité.
Ici doit être calculée en utilisant la même formule que dans la colonne transformation de Fourier. Cela résulte d'un changement de la variable "muette", de x à ξ ou ω ou ν.
d n f ( x ) d x n {\displaystyle {\frac {\mathrm {d} ^{n}f(x)}{\mathrm {d} x^{n}}}\,} ( 2 π i ξ ) n f ^ ( ξ ) {\displaystyle (2\pi {\rm {i}}\xi )^{n}{\hat {f}}(\xi )\,} ( i ω ) n f ^ ( ω ) {\displaystyle ({\rm {i}}\omega )^{n}{\hat {f}}(\omega )\,} ( i ν ) n f ^ ( ν ) {\displaystyle ({\rm {i}}\nu )^{n}{\hat {f}}(\nu )\,}
x n f ( x ) {\displaystyle x^{n}f(x)\,} ( i 2 π ) n d n f ^ ( ξ ) d ξ n {\displaystyle \left({\frac {\rm {i}}{2\pi }}\right)^{n}{\frac {\mathrm {d} ^{n}{\hat {f}}(\xi )}{\mathrm {d} \xi ^{n}}}\,} i n d n f ^ ( ω ) d ω n {\displaystyle {\rm {i}}^{n}{\frac {\mathrm {d} ^{n}{\hat {f}}(\omega )}{\mathrm {d} \omega ^{n}}}} i n d n f ^ ( ν ) d ν n {\displaystyle {\rm {i}}^{n}{\frac {\mathrm {d} ^{n}{\hat {f}}(\nu )}{\mathrm {d} \nu ^{n}}}} Relation duale de la précédente
( f g ) ( x ) {\displaystyle (f*g)(x)\,} f ^ ( ξ ) g ^ ( ξ ) {\displaystyle {\hat {f}}(\xi ){\hat {g}}(\xi )\,} 2 π f ^ ( ω ) g ^ ( ω ) {\displaystyle {\sqrt {2\pi }}{\hat {f}}(\omega ){\hat {g}}(\omega )\,} f ^ ( ν ) g ^ ( ν ) {\displaystyle {\hat {f}}(\nu ){\hat {g}}(\nu )\,} La notation fg signifie le produit de convolution de f et g — cette règle est le théorème de Fubini
f ( x ) g ( x ) {\displaystyle f(x)g(x)\,} ( f ^ g ^ ) ( ξ ) {\displaystyle \left({\hat {f}}*{\hat {g}}\right)(\xi )\,} 1 2 π ( f ^ g ^ ) ( ω ) {\displaystyle {\frac {1}{\sqrt {2\pi }}}\left({\hat {f}}*{\hat {g}}\right)(\omega )\,} 1 2 π ( f ^ g ^ ) ( ν ) {\displaystyle {\frac {1}{2\pi }}\left({\hat {f}}*{\hat {g}}\right)(\nu )\,} Relation duale du théorème de Fubini
Si f (x) est purement réelle f ^ ( ξ ) = f ^ ( ξ ) ¯ {\displaystyle {\hat {f}}(-\xi )={\overline {{\hat {f}}(\xi )}}\,} f ^ ( ω ) = f ^ ( ω ) ¯ {\displaystyle {\hat {f}}(-\omega )={\overline {{\hat {f}}(\omega )}}\,} f ^ ( ν ) = f ^ ( ν ) ¯ {\displaystyle {\hat {f}}(-\nu )={\overline {{\hat {f}}(\nu )}}\,} Symétrie hermitienne.
z est la notation du complexe conjugué de z.
Si f (x) est purement réelle et paire  (ξ),  (ω) et  (ν) sont purement réelles et paires.
Si f (x) est purement réelle et impaire  (ξ),  (ω) et  (ν) sont purement imaginaires et impaires.
Si f (x) est imaginaire pur f ^ ( ξ ) = f ^ ( ξ ) ¯ {\displaystyle {\hat {f}}(-\xi )=-{\overline {{\hat {f}}(\xi )}}\,} f ^ ( ω ) = f ^ ( ω ) ¯ {\displaystyle {\hat {f}}(-\omega )=-{\overline {{\hat {f}}(\omega )}}\,} f ^ ( ν ) = f ^ ( ν ) ¯ {\displaystyle {\hat {f}}(-\nu )=-{\overline {{\hat {f}}(\nu )}}\,} z est le complexe conjugué de z.
f ( x ) ¯ {\displaystyle {\overline {f(x)}}} f ^ ( ξ ) ¯ {\displaystyle {\overline {{\hat {f}}(-\xi )}}} f ^ ( ω ) ¯ {\displaystyle {\overline {{\hat {f}}(-\omega )}}} f ^ ( ν ) ¯ {\displaystyle {\overline {{\hat {f}}(-\nu )}}} Conjugaison complexe
f ( x ) cos ( a x ) {\displaystyle f(x)\cos(ax)} f ^ ( ξ a 2 π ) + f ^ ( ξ + a 2 π ) 2 {\displaystyle {\frac {{\hat {f}}\left(\xi -{\frac {a}{2\pi }}\right)+{\hat {f}}\left(\xi +{\frac {a}{2\pi }}\right)}{2}}} f ^ ( ω a ) + f ^ ( ω + a ) 2 {\displaystyle {\frac {{\hat {f}}(\omega -a)+{\hat {f}}(\omega +a)}{2}}\,} f ^ ( ν a ) + f ^ ( ν + a ) 2 {\displaystyle {\frac {{\hat {f}}(\nu -a)+{\hat {f}}(\nu +a)}{2}}} Peut se déduire de la formule d'Euler : cos ( a x ) = e i a x + e i a x 2 {\displaystyle \cos(ax)={\frac {{\rm {e}}^{{\rm {i}}ax}+{\rm {e}}^{-{\rm {i}}ax}}{2}}}
f ( x ) sin ( a x ) {\displaystyle f(x)\sin(ax)} f ^ ( ξ a 2 π ) f ^ ( ξ + a 2 π ) 2 i {\displaystyle {\frac {{\hat {f}}\left(\xi -{\frac {a}{2\pi }}\right)-{\hat {f}}\left(\xi +{\frac {a}{2\pi }}\right)}{2{\rm {i}}}}} f ^ ( ω a ) f ^ ( ω + a ) 2 i {\displaystyle {\frac {{\hat {f}}(\omega -a)-{\hat {f}}(\omega +a)}{2{\rm {i}}}}} f ^ ( ν a ) f ^ ( ν + a ) 2 i {\displaystyle {\frac {{\hat {f}}(\nu -a)-{\hat {f}}(\nu +a)}{2{\rm {i}}}}} Peut se déduire de la formule d'Euler : sin ( a x ) = e i a x e i a x 2 i {\displaystyle \sin(ax)={\frac {{\rm {e}}^{{\rm {i}}ax}-{\rm {e}}^{-{\rm {i}}ax}}{2{\rm {i}}}}}

Fonctions de carré intégrable à une variable

Les transformées de Fourier de ce tableau peuvent être trouvées dans les deux références précédentes ou dans (en) George Campbell ; Ronald Foster, Fourier Integrals for Practical Applications, New York, USA, D. Van Nostrand Company, Inc, .

Fonction Transformée de Fourier

ξ est la fréquence

Transformée de Fourier
ω = 2πξ est la pulsation ou fréquence angulaire
Transformée de Fourier

définition alternative

Remarques
f ( x ) {\displaystyle f(x)\,} f ^ ( ξ ) = f ( x ) e 2 π i x ξ d x {\displaystyle {\begin{aligned}&{\hat {f}}(\xi )\\&=\int _{-\infty }^{\infty }f(x){\rm {e}}^{-2\pi ix\xi }\,\mathrm {d} x\end{aligned}}} f ^ ( ω ) = 1 2 π f ( x ) e i ω x d x {\displaystyle {\begin{aligned}&{\hat {f}}(\omega )\\&={\frac {1}{\sqrt {2\pi }}}\int _{-\infty }^{\infty }f(x){\rm {e}}^{-{\rm {i}}\omega x}\,\mathrm {d} x\end{aligned}}} f ^ ( ν ) = f ( x ) e i ν x d x {\displaystyle {\begin{aligned}&{\hat {f}}(\nu )\\&=\int _{-\infty }^{\infty }f(x){\rm {e}}^{-{\rm {i}}\nu x}\,\mathrm {d} x\end{aligned}}}
rect ( a x ) {\displaystyle \operatorname {rect} (ax)\,} 1 | a | sinc ( ξ a ) {\displaystyle {\frac {1}{|a|}}\cdot \operatorname {sinc} \left({\frac {\xi }{a}}\right)} 1 2 π a 2 sinc ( ω 2 π a ) {\displaystyle {\frac {1}{\sqrt {2\pi a^{2}}}}\cdot \operatorname {sinc} \left({\frac {\omega }{2\pi a}}\right)} 1 | a | sinc ( ν 2 π a ) {\displaystyle {\frac {1}{|a|}}\cdot \operatorname {sinc} \left({\frac {\nu }{2\pi a}}\right)} Pour la fonction rectangulaire voir fonction porte ; la fonction sinus cardinal normalisé est définie ici par sinc(x) = sin(πx)/πx
sinc ( a x ) {\displaystyle \operatorname {sinc} (ax)\,} 1 | a | rect ( ξ a ) {\displaystyle {\frac {1}{|a|}}\cdot \operatorname {rect} \left({\frac {\xi }{a}}\right)\,} 1 2 π a 2 rect ( ω 2 π a ) {\displaystyle {\frac {1}{\sqrt {2\pi a^{2}}}}\cdot \operatorname {rect} \left({\frac {\omega }{2\pi a}}\right)} 1 | a | rect ( ν 2 π a ) {\displaystyle {\frac {1}{|a|}}\cdot \operatorname {rect} \left({\frac {\nu }{2\pi a}}\right)} Relation duale de la précédente.
La fonction porte est un filtre passe-bas idéal, et la fonction sinus cardinal est la réponse impulsionnelle non causale d'un tel filtre. La fonction sinc est la fonction sinus cardinal normalisée : sinc(x) = sin(πx)/πx
sinc 2 ( a x ) {\displaystyle \operatorname {sinc} ^{2}(ax)} 1 | a | tri ( ξ a ) {\displaystyle {\frac {1}{|a|}}\cdot \operatorname {tri} \left({\frac {\xi }{a}}\right)} 1 2 π a 2 tri ( ω 2 π a ) {\displaystyle {\frac {1}{\sqrt {2\pi a^{2}}}}\cdot \operatorname {tri} \left({\frac {\omega }{2\pi a}}\right)} 1 | a | tri ( ν 2 π a ) {\displaystyle {\frac {1}{|a|}}\cdot \operatorname {tri} \left({\frac {\nu }{2\pi a}}\right)} La fonction tri(x) est la fonction triangulaire.
tri ( a x ) {\displaystyle \operatorname {tri} (ax)} 1 | a | sinc 2 ( ξ a ) {\displaystyle {\frac {1}{|a|}}\cdot \operatorname {sinc} ^{2}\left({\frac {\xi }{a}}\right)\,} 1 2 π a 2 sinc 2 ( ω 2 π a ) {\displaystyle {\frac {1}{\sqrt {2\pi a^{2}}}}\cdot \operatorname {sinc} ^{2}\left({\frac {\omega }{2\pi a}}\right)} 1 | a | sinc 2 ( ν 2 π a ) {\displaystyle {\frac {1}{|a|}}\cdot \operatorname {sinc} ^{2}\left({\frac {\nu }{2\pi a}}\right)} Relation duale de la précédente
e a x u ( x ) {\displaystyle {\rm {e}}^{-ax}u(x)\,} 1 a + 2 π i ξ {\displaystyle {\frac {1}{a+2\pi {\rm {i}}\xi }}} 1 2 π ( a + i ω ) {\displaystyle {\frac {1}{{\sqrt {2\pi }}(a+{\rm {i}}\omega )}}} 1 a + i ν {\displaystyle {\frac {1}{a+{\rm {i}}\nu }}} La fonction u(x) est la fonction marche de Heaviside et a > 0.
e α x 2 {\displaystyle {\rm {e}}^{-\alpha x^{2}}\,} π α e ( π ξ ) 2 α {\displaystyle {\sqrt {\frac {\pi }{\alpha }}}\cdot {\rm {e}}^{-{\frac {(\pi \xi )^{2}}{\alpha }}}} 1 2 α e ω 2 4 α {\displaystyle {\frac {1}{\sqrt {2\alpha }}}\cdot {\rm {e}}^{-{\frac {\omega ^{2}}{4\alpha }}}} π α e ν 2 4 α {\displaystyle {\sqrt {\frac {\pi }{\alpha }}}\cdot {\rm {e}}^{-{\frac {\nu ^{2}}{4\alpha }}}} Nota : pour les deux premières transformations de Fourier, la fonction gaussienne eαx2 est, pour un choix judicieux de α, sa propre transformée de Fourier. Pour qu'elle soit intégrable on doit avoir Re(α) > 0.
e a | x | {\displaystyle \operatorname {e} ^{-a|x|}\,} 2 a a 2 + 4 π 2 ξ 2 {\displaystyle {\frac {2a}{a^{2}+4\pi ^{2}\xi ^{2}}}} 2 π a a 2 + ω 2 {\displaystyle {\sqrt {\frac {2}{\pi }}}\cdot {\frac {a}{a^{2}+\omega ^{2}}}} 2 a a 2 + ν 2 {\displaystyle {\frac {2a}{a^{2}+\nu ^{2}}}} Pour Re(a) > 0. Ceci signifie que la transformée de Fourier d'une densité de probabilité d'une distribution de Laplace est une densité de probabilité d'une loi de Cauchy.
sech ( a x ) {\displaystyle \operatorname {sech} (ax)\,} π a sech ( π 2 a ξ ) {\displaystyle {\frac {\pi }{a}}\operatorname {sech} \left({\frac {\pi ^{2}}{a}}\xi \right)} 1 a π 2 sech ( π 2 a ω ) {\displaystyle {\frac {1}{a}}{\sqrt {\frac {\pi }{2}}}\operatorname {sech} \left({\frac {\pi }{2a}}\omega \right)} π a sech ( π 2 a ν ) {\displaystyle {\frac {\pi }{a}}\operatorname {sech} \left({\frac {\pi }{2a}}\nu \right)} La sécante hyperbolique est sa propre transformée de Fourier.
e a 2 x 2 2 H n ( a x ) {\displaystyle {\rm {e}}^{-{\frac {a^{2}x^{2}}{2}}}H_{n}(ax)\,} 2 π ( i ) n a e 2 π 2 ξ 2 a 2 H n ( 2 π ξ a ) {\displaystyle {\frac {{\sqrt {2\pi }}(-{\rm {i}})^{n}}{a}}{\rm {e}}^{-{\frac {2\pi ^{2}\xi ^{2}}{a^{2}}}}H_{n}\left({\frac {2\pi \xi }{a}}\right)} ( i ) n a e ω 2 2 a 2 H n ( ω a ) {\displaystyle {\frac {(-{\rm {i}})^{n}}{a}}{\rm {e}}^{-{\frac {\omega ^{2}}{2a^{2}}}}H_{n}\left({\frac {\omega }{a}}\right)} ( i ) n 2 π a e ν 2 2 a 2 H n ( ν a ) {\displaystyle {\frac {(-{\rm {i}})^{n}{\sqrt {2\pi }}}{a}}{\rm {e}}^{-{\frac {\nu ^{2}}{2a^{2}}}}H_{n}\left({\frac {\nu }{a}}\right)} Hn est le ne polynôme d'Hermite. Si a = 1 alors les fonctions d'Hermite-Gauss sont des fonctions propres de l'opérateur transformée de Fourier.

Distributions à une variable

Les transformées de Fourier de ce tableau sont traitées dans (en) Arthur Erdélyi, Tables of Integral Transforms, Vol. 1, McGraw-Hill, ou dans (en) David Kammler, A First Course in Fourier Analysis, USA, Prentice Hall, .

Fonction Transformée de Fourier

ξ est la fréquence

Transformée de Fourier
ω = 2πξ est la pulsation ou fréquence angulaire
Transformée de Fourier

définition alternative

Remarques
f ( x ) {\displaystyle f(x)\,} f ^ ( ξ ) = f ( x ) e 2 π i x ξ d x {\displaystyle {\begin{aligned}&{\hat {f}}(\xi )=\int _{-\infty }^{\infty }f(x){\rm {e}}^{-2\pi {\rm {i}}x\xi }\,{\rm {d}}x\end{aligned}}} f ^ ( ω ) = 1 2 π f ( x ) e i ω x d x {\displaystyle {\begin{aligned}&{\hat {f}}(\omega )={\frac {1}{\sqrt {2\pi }}}\int _{-\infty }^{\infty }f(x){\rm {e}}^{-{\rm {i}}\omega x}\,{\rm {d}}x\end{aligned}}} f ^ ( ν ) = f ( x ) e i ν x d x {\displaystyle {\begin{aligned}&{\hat {f}}(\nu )=\int _{-\infty }^{\infty }f(x){\rm {e}}^{-{\rm {i}}\nu x}\,{\rm {d}}x\end{aligned}}}
1 {\displaystyle 1} δ ( ξ ) {\displaystyle \delta (\xi )} 2 π δ ( ω ) {\displaystyle {\sqrt {2\pi }}\cdot \delta (\omega )} 2 π δ ( ν ) {\displaystyle 2\pi \delta (\nu )} δ(ξ) est la distribution de Dirac.
δ ( x ) {\displaystyle \delta (x)\,} 1 {\displaystyle 1} 1 2 π {\displaystyle {\frac {1}{\sqrt {2\pi }}}\,} 1 {\displaystyle 1} Relation duale de la précédente.
e i a x {\displaystyle {\rm {e}}^{{\rm {i}}ax}} δ ( ξ a 2 π ) {\displaystyle \delta \left(\xi -{\frac {a}{2\pi }}\right)} 2 π δ ( ω a ) {\displaystyle {\sqrt {2\pi }}\cdot \delta (\omega -a)} 2 π δ ( ν a ) {\displaystyle 2\pi \delta (\nu -a)}
cos ( a x ) {\displaystyle \cos(ax)} δ ( ξ a 2 π ) + δ ( ξ + a 2 π ) 2 {\displaystyle {\frac {\delta \left(\xi -{\frac {a}{2\pi }}\right)+\delta \left(\xi +{\frac {a}{2\pi }}\right)}{2}}} 2 π δ ( ω a ) + δ ( ω + a ) 2 {\displaystyle {\sqrt {2\pi }}\cdot {\frac {\delta (\omega -a)+\delta (\omega +a)}{2}}\,} π ( δ ( ν a ) + δ ( ν + a ) ) {\displaystyle \pi \left(\delta (\nu -a)+\delta (\nu +a)\right)} Résulte de la formule d'Euler.
sin ( a x ) {\displaystyle \sin(ax)} δ ( ξ a 2 π ) δ ( ξ + a 2 π ) 2 i {\displaystyle {\frac {\delta \left(\xi -{\frac {a}{2\pi }}\right)-\delta \left(\xi +{\frac {a}{2\pi }}\right)}{2i}}} 2 π δ ( ω a ) δ ( ω + a ) 2 i {\displaystyle {\sqrt {2\pi }}\cdot {\frac {\delta (\omega -a)-\delta (\omega +a)}{2i}}} i π ( δ ( ν a ) δ ( ν + a ) ) {\displaystyle -{\rm {i}}\pi \left(\delta (\nu -a)-\delta (\nu +a)\right)} Résulte de la formule d'Euler.
cos ( a x 2 ) {\displaystyle \cos \left(ax^{2}\right)} π a cos ( π 2 ξ 2 a π 4 ) {\displaystyle {\sqrt {\frac {\pi }{a}}}\cos \left({\frac {\pi ^{2}\xi ^{2}}{a}}-{\frac {\pi }{4}}\right)} 1 2 a cos ( ω 2 4 a π 4 ) {\displaystyle {\frac {1}{\sqrt {2a}}}\cos \left({\frac {\omega ^{2}}{4a}}-{\frac {\pi }{4}}\right)} π a cos ( ν 2 4 a π 4 ) {\displaystyle {\sqrt {\frac {\pi }{a}}}\cos \left({\frac {\nu ^{2}}{4a}}-{\frac {\pi }{4}}\right)}
sin ( a x 2 ) {\displaystyle \sin \left(ax^{2}\right)\,} π a sin ( π 2 ξ 2 a π 4 ) {\displaystyle -{\sqrt {\frac {\pi }{a}}}\sin \left({\frac {\pi ^{2}\xi ^{2}}{a}}-{\frac {\pi }{4}}\right)} 1 2 a sin ( ω 2 4 a π 4 ) {\displaystyle {\frac {-1}{\sqrt {2a}}}\sin \left({\frac {\omega ^{2}}{4a}}-{\frac {\pi }{4}}\right)} π a sin ( ν 2 4 a π 4 ) {\displaystyle -{\sqrt {\frac {\pi }{a}}}\sin \left({\frac {\nu ^{2}}{4a}}-{\frac {\pi }{4}}\right)}
x n {\displaystyle x^{n}\,} ( i 2 π ) n δ ( n ) ( ξ ) {\displaystyle \left({\frac {\rm {i}}{2\pi }}\right)^{n}\delta ^{(n)}(\xi )\,} i n 2 π δ ( n ) ( ω ) {\displaystyle {\rm {i}}^{n}{\sqrt {2\pi }}\delta ^{(n)}(\omega )\,} 2 π i n δ ( n ) ( ν ) {\displaystyle 2\pi {\rm {i}}^{n}\delta ^{(n)}(\nu )\,} Ici n est un entier naturel et δ ( n ) ( ξ ) {\displaystyle \delta ^{(n)}(\xi )\,}
est la nième dérivée (au sens des distributions) de la distribution de Dirac.
On peut en déduire les transformées de tous les polynômes.
δ ( n ) ( x ) {\displaystyle \delta ^{(n)}(x)\,} ( 2 π i ξ ) n {\displaystyle (2\pi {\rm {i}}\xi )^{n}\,} ( i ω ) n 2 π {\displaystyle {\frac {({\rm {i}}\omega )^{n}}{\sqrt {2\pi }}}\,} ( i ν ) n {\displaystyle ({\rm {i}}\nu )^{n}\,} δ ( n ) ( ξ ) {\displaystyle \delta ^{(n)}(\xi )\,} est la nième dérivée (au sens des distributions) de la distribution de Dirac.
1 x {\displaystyle {\frac {1}{x}}} i π sgn ( ξ ) {\displaystyle -{\rm {i}}\pi \operatorname {sgn}(\xi )} i π 2 sgn ( ω ) {\displaystyle -{\rm {i}}{\sqrt {\frac {\pi }{2}}}\operatorname {sgn}(\omega )} i π sgn ( ν ) {\displaystyle -{\rm {i}}\pi \operatorname {sgn}(\nu )} Ici sgn(ξ) est la fonction signe.
On notera que 1/x n'est pas une distribution.
On doit utiliser la valeur principale de Cauchy pour étudier les fonctions de Schwartz.
Cette règle est utile quand on étudie la transformation de Hilbert.
1 x n := ( 1 ) n 1 ( n 1 ) ! d n d x n log | x | {\displaystyle {\begin{aligned}&{\frac {1}{x^{n}}}\\&:={\frac {(-1)^{n-1}}{(n-1)!}}{\frac {{\rm {d}}^{n}}{{\rm {d}}x^{n}}}\log |x|\end{aligned}}} i π ( 2 π i ξ ) n 1 ( n 1 ) ! sgn ( ξ ) {\displaystyle -{\rm {i}}\pi {\frac {(-2\pi {\rm {i}}\xi )^{n-1}}{(n-1)!}}\operatorname {sgn}(\xi )} i π 2 ( i ω ) n 1 ( n 1 ) ! sgn ( ω ) {\displaystyle -{\rm {i}}{\sqrt {\frac {\pi }{2}}}\cdot {\frac {(-{\rm {i}}\omega )^{n-1}}{(n-1)!}}\operatorname {sgn}(\omega )} i π ( i ν ) n 1 ( n 1 ) ! sgn ( ν ) {\displaystyle -{\rm {i}}\pi {\frac {(-{\rm {i}}\nu )^{n-1}}{(n-1)!}}\operatorname {sgn}(\nu )} 1/xn est la distribution homogène (en)
définie par la dérivée de ( 1 ) n 1 ( n 1 ) ! d n d x n log | x | {\displaystyle {\frac {(-1)^{n-1}}{(n-1)!}}{\frac {{\rm {d}}^{n}}{{\rm {d}}x^{n}}}\log |x|}
| x | α {\displaystyle |x|^{\alpha }\,} 2 sin ( π α 2 ) Γ ( α + 1 ) | 2 π ξ | α + 1 {\displaystyle -{\frac {2\sin \left({\frac {\pi \alpha }{2}}\right)\Gamma (\alpha +1)}{|2\pi \xi |^{\alpha +1}}}} 2 2 π sin ( π α 2 ) Γ ( α + 1 ) | ω | α + 1 {\displaystyle {\frac {-2}{\sqrt {2\pi }}}\cdot {\frac {\sin \left({\frac {\pi \alpha }{2}}\right)\Gamma (\alpha +1)}{|\omega |^{\alpha +1}}}} 2 sin ( π α 2 ) Γ ( α + 1 ) | ν | α + 1 {\displaystyle -{\frac {2\sin \left({\frac {\pi \alpha }{2}}\right)\Gamma (\alpha +1)}{|\nu |^{\alpha +1}}}} Formule valide pour α réel avec -1 < α < 0.
Si α complexe avec Re (α) > −1, alors | x | α {\displaystyle |x|^{\alpha }\,} est une fonction localement intégrable et est donc une distribution tempérée.
La fonction α ↦ |x|α est une fonction holomorphe du demi-plan complexe réel dans l'espace des distributions tempérées. Elle admet une unique extension méromorphe qui est une distribution tempérée également notée | x | α {\displaystyle |x|^{\alpha }\,} uniquement pour α ≠ −2, −4,....
1 | x | {\displaystyle {\frac {1}{\sqrt {|x|}}}\,} 1 | ξ | {\displaystyle {\frac {1}{\sqrt {|\xi |}}}} 1 | ω | {\displaystyle {\frac {1}{\sqrt {|\omega |}}}} 2 π | ν | {\displaystyle {\frac {\sqrt {2\pi }}{\sqrt {|\nu |}}}} Cas particulier de la précédente, pour α = –1/2.
sgn ( x ) {\displaystyle \operatorname {sgn}(x)} 1 i π ξ {\displaystyle {\frac {1}{{\rm {i}}\pi \xi }}} 2 π 1 i ω {\displaystyle {\sqrt {\frac {2}{\pi }}}{\frac {1}{{\rm {i}}\omega }}} 2 i ν {\displaystyle {\frac {2}{{\rm {i}}\nu }}} La transformation de Fourier doit ici être prise comme la valeur principale de Cauchy.
u ( x ) {\displaystyle u(x)} 1 2 ( 1 i π ξ + δ ( ξ ) ) {\displaystyle {\frac {1}{2}}\left({\frac {1}{{\rm {i}}\pi \xi }}+\delta (\xi )\right)} π 2 ( 1 i π ω + δ ( ω ) ) {\displaystyle {\sqrt {\frac {\pi }{2}}}\left({\frac {1}{{\rm {i}}\pi \omega }}+\delta (\omega )\right)} π ( 1 i π ν + δ ( ν ) ) {\displaystyle \pi \left({\frac {1}{{\rm {i}}\pi \nu }}+\delta (\nu )\right)} La fonction u(x) est la fonction de Heaviside.
n = δ ( x n T ) {\displaystyle \sum _{n=-\infty }^{\infty }\delta (x-nT)} 1 T k = δ ( ξ k T ) {\displaystyle {\frac {1}{T}}\sum _{k=-\infty }^{\infty }\delta \left(\xi -{\frac {k}{T}}\right)} 2 π T k = δ ( ω 2 π k T ) {\displaystyle {\frac {\sqrt {2\pi }}{T}}\sum _{k=-\infty }^{\infty }\delta \left(\omega -{\frac {2\pi k}{T}}\right)} 2 π T k = δ ( ν 2 π k T ) {\displaystyle {\frac {2\pi }{T}}\sum _{k=-\infty }^{\infty }\delta \left(\nu -{\frac {2\pi k}{T}}\right)} Transformée de Fourier du peigne de Dirac. On utilise aussi le fait que
n = e i n x = 2 π k = δ ( x + 2 π k ) {\displaystyle \sum _{n=-\infty }^{\infty }{\rm {e}}^{{\rm {i}}nx}=2\pi \sum _{k=-\infty }^{\infty }\delta (x+2\pi k)} sont considérées comme des distributions.
J 0 ( x ) {\displaystyle J_{0}(x)} 2 rect ( π ξ ) 1 4 π 2 ξ 2 {\displaystyle {\frac {2\,\operatorname {rect} (\pi \xi )}{\sqrt {1-4\pi ^{2}\xi ^{2}}}}} 2 π rect ( ω 2 ) 1 ω 2 {\displaystyle {\sqrt {\frac {2}{\pi }}}\cdot {\frac {\operatorname {rect} \left({\frac {\omega }{2}}\right)}{\sqrt {1-\omega ^{2}}}}} 2 rect ( ν 2 ) 1 ν 2 {\displaystyle {\frac {2\,\operatorname {rect} \left({\frac {\nu }{2}}\right)}{\sqrt {1-\nu ^{2}}}}} La fonction J0(x) est la fonction de Bessel d'ordre zéro de la 1re espèce.
J n ( x ) {\displaystyle J_{n}(x)} 2 ( i ) n T n ( 2 π ξ ) rect ( π ξ ) 1 4 π 2 ξ 2 {\displaystyle {\frac {2(-{\rm {i}})^{n}T_{n}(2\pi \xi )\operatorname {rect} (\pi \xi )}{\sqrt {1-4\pi ^{2}\xi ^{2}}}}} 2 π ( i ) n T n ( ω ) rect ( ω 2 ) 1 ω 2 {\displaystyle {\sqrt {\frac {2}{\pi }}}{\frac {(-{\rm {i}})^{n}T_{n}(\omega )\operatorname {rect} \left({\frac {\omega }{2}}\right)}{\sqrt {1-\omega ^{2}}}}} 2 ( i ) n T n ( ν ) rect ( ν 2 ) 1 ν 2 {\displaystyle {\frac {2(-{\rm {i}})^{n}T_{n}(\nu )\operatorname {rect} \left({\frac {\nu }{2}}\right)}{\sqrt {1-\nu ^{2}}}}} Généralisation de 315. La fonction Jn(x) est la fonction de Bessel d'ordre n de la 1re espèce. La fonction Tn(x) est le polynôme de Tchebychev de 1re espèce.
log | x | {\displaystyle \log \left|x\right|} 1 2 1 | ξ | γ δ ( ξ ) {\displaystyle -{\frac {1}{2}}{\frac {1}{\left|\xi \right|}}-\gamma \delta (\xi )} π 2 | ω | 2 π γ δ ( ω ) {\displaystyle -{\frac {\sqrt {\frac {\pi }{2}}}{\left|\omega \right|}}-{\sqrt {2\pi }}\gamma \delta (\omega )} π | ν | 2 π γ δ ( ν ) {\displaystyle -{\frac {\pi }{\left|\nu \right|}}-2\pi \gamma \delta (\nu )} γ est la constante d'Euler-Mascheroni.
( i x ) α {\displaystyle \left(\mp {\rm {i}}x\right)^{-\alpha }} ( 2 π ) α Γ ( α ) u ( ± ξ ) ( ± ξ ) α 1 {\displaystyle {\frac {\left(2\pi \right)^{\alpha }}{\Gamma (\alpha )}}u\left(\pm \xi \right)\left(\pm \xi \right)^{\alpha -1}} 2 π Γ ( α ) u ( ± ω ) ( ± ω ) α 1 {\displaystyle {\frac {\sqrt {2\pi }}{\Gamma (\alpha )}}u(\pm \omega )(\pm \omega )^{\alpha -1}} 2 π Γ ( α ) u ( ± ν ) ( ± ν ) α 1 {\displaystyle {\frac {2\pi }{\Gamma \left(\alpha \right)}}u(\pm \nu )(\pm \nu )^{\alpha -1}} Formule correcte pour 0 < α < 1. La formule de dérivation permet de déduire la formule pour des exposants plus élevés. u est la fonction de Heaviside.

Fonctions de deux variables

Fonction Transformée de Fourier

ξ est la fréquence

Transformée de Fourier
ω = 2πξ est la pulsation ou fréquence angulaire
Transformée de Fourier

définition alternative

Remarques
f ( x , y ) {\displaystyle f(x,y)} f ^ ( ξ x , ξ y ) = f ( x , y ) e 2 π i ( ξ x x + ξ y y ) d x d y {\displaystyle {\begin{aligned}&{\hat {f}}(\xi _{x},\xi _{y})\\&=\iint f(x,y){\rm {e}}^{-2\pi {\rm {i}}(\xi _{x}x+\xi _{y}y)}\,{\rm {d}}x\,{\rm {d}}y\end{aligned}}} f ^ ( ω x , ω y ) = 1 2 π f ( x , y ) e i ( ω x x + ω y y ) d x d y {\displaystyle {\begin{aligned}&{\hat {f}}(\omega _{x},\omega _{y})\\&={\frac {1}{2\pi }}\iint f(x,y){\rm {e}}^{-{\rm {i}}(\omega _{x}x+\omega _{y}y)}\,{\rm {d}}x\,{\rm {d}}y\end{aligned}}} f ^ ( ν x , ν y ) = f ( x , y ) e i ( ν x x + ν y y ) d x d y {\displaystyle {\begin{aligned}&{\hat {f}}(\nu _{x},\nu _{y})\\&=\iint f(x,y){\rm {e}}^{-{\rm {i}}(\nu _{x}x+\nu _{y}y)}\,{\rm {d}}x\,{\rm {d}}y\end{aligned}}} Les variables ξx, ξy, ωx, ωy, νx, νy sont réelles. Les intégrales couvrent tout le plan complexe.
e π ( a 2 x 2 + b 2 y 2 ) {\displaystyle {\rm {e}}^{-\pi \left(a^{2}x^{2}+b^{2}y^{2}\right)}} 1 | a b | e π ( ξ x 2 a 2 + ξ y 2 b 2 ) {\displaystyle {\frac {1}{|ab|}}{\rm {e}}^{-\pi \left({\frac {\xi _{x}^{2}}{a^{2}}}+{\frac {\xi _{y}^{2}}{b^{2}}}\right)}} 1 2 π | a b | e 1 4 π ( ω x 2 a 2 + ω y 2 b 2 ) {\displaystyle {\frac {1}{2\pi \cdot |ab|}}{\rm {e}}^{-{\frac {1}{4\pi }}\left({\frac {\omega _{x}^{2}}{a^{2}}}+{\frac {\omega _{y}^{2}}{b^{2}}}\right)}} 1 | a b | e 1 4 π ( ν x 2 a 2 + ν y 2 b 2 ) {\displaystyle {\frac {1}{|ab|}}{\rm {e}}^{-{\frac {1}{4\pi }}\left({\frac {\nu _{x}^{2}}{a^{2}}}+{\frac {\nu _{y}^{2}}{b^{2}}}\right)}} La fonction et ses transformées sont toutes des gaussiennes.
circ ( x 2 + y 2 ) {\displaystyle \operatorname {circ} \left({\sqrt {x^{2}+y^{2}}}\right)} J 1 ( 2 π ξ x 2 + ξ y 2 ) ξ x 2 + ξ y 2 {\displaystyle {\frac {J_{1}\left(2\pi {\sqrt {\xi _{x}^{2}+\xi _{y}^{2}}}\right)}{\sqrt {\xi _{x}^{2}+\xi _{y}^{2}}}}} J 1 ( ω x 2 + ω y 2 ) ω x 2 + ω y 2 {\displaystyle {\frac {J_{1}\left({\sqrt {\omega _{x}^{2}+\omega _{y}^{2}}}\right)}{\sqrt {\omega _{x}^{2}+\omega _{y}^{2}}}}} 2 π J 1 ( ν x 2 + ν y 2 ) ν x 2 + ν y 2 {\displaystyle {\frac {2\pi J_{1}\left({\sqrt {\nu _{x}^{2}+\nu _{y}^{2}}}\right)}{\sqrt {\nu _{x}^{2}+\nu _{y}^{2}}}}} La fonction est définie par circ(r) = 1 sur 0 ≤ r ≤ 1, et est nulle partout ailleurs. Le résultat est la distribution de l'amplitude de la tache d'Airy. J1 est la fonction de Bessel de première espèce d'ordre 1[6].
i x + i y {\displaystyle {\frac {\rm {i}}{x+{\rm {i}}y}}} 1 ξ x + i ξ y {\displaystyle {\frac {1}{\xi _{x}+{\rm {i}}\xi _{y}}}} 1 ω x + i ω y {\displaystyle {\frac {1}{\omega _{x}+{\rm {i}}\omega _{y}}}} 2 π ν x + i ν y {\displaystyle {\frac {2\pi }{\nu _{x}+{\rm {i}}\nu _{y}}}}

Notes et références

(en) Cet article est partiellement ou en totalité issu de l’article de Wikipédia en anglais intitulé « Fourier transform » (voir la liste des auteurs).
  1. Commission électrotechnique internationale, « Mathématiques - Fonctions : Transformations integrales », dans IEC 60050 Vocabulaire électrotechnique international, 1987/1994 (lire en ligne), p. 103-04-01
  2. Alain Bouvier, Michel George et François Le Lionnais, Dictionnaire des mathématiques, Presses universitaires de France, (1re éd. 1979), p. 361.
  3. Walter Rudin, Analyse réelle et complexe [détail des éditions], p. 174 de l'édition de 1975-77.
  4. (en) Mark Pinsky (en), Introduction to Fourier Analysis and Wavelets, Brooks/Cole, (ISBN 978-0-8218-7198-0, lire en ligne), p. 131.
  5. Ou plus précisément à l’ombre de cette somme[réf. nécessaire].
  6. Stein et Weiss 1971, Thm. IV.3.3.

Voir aussi

Sur les autres projets Wikimedia :

  • Série et transformée de Fourier en physique, sur Wikiversity

Articles connexes

Bibliographie

Liens externes

  • Alain Yger, Espaces de Hilbert et analyse de Fourier (2008), cours de 3e année de licence, université Bordeaux I
  • (en) « FTL-SE », programme éducatif sur les transformées de Fourier d'images
  • J. Fourier, Théorie analytique de la chaleur (1822), chap. III (fondements de la transformée de Fourier), en ligne et commenté sur le site BibNum
  • Œil et physiologie de la vision : « les signaux électrophysiologiques ».
  • icône décorative Portail de l'analyse