Transformée de Fourier

Joseph Fourier Carl Frederich Gauss Leonhard Euler

1768-1830 1777-1855 1707-1783

November 29, 2021 1 / 21

Plan

1 Transformée de Fourier continue

2 Transformée de Fourier discrète
3 Transformée de Fourier IRL

November 29, 2021 2 / 21

Transformée de Fourier continue

November 29, 2021 3 / 21

Transformée de Fourier continue

Série de Fourier
∞
X
f ∈ L2 ([0, 1]) ⇐⇒ f (x) = fˆn e 2iπnx
n=−∞

where Z 1
fˆn = f (x)e −2iπnx dx
0

fˆn sont les coefficients de Fourier (espace spectral)

November 29, 2021 4 / 21

Propriétés

f réelle ⇐⇒ ∀k, fˆ−k = fˆk ∗

November 29, 2021 5 / 21

Propriétés

f réelle ⇐⇒ ∀k, fˆ−k = fˆk ∗

(f\+ g )n = fˆn + ĝn

November 29, 2021 5 / 21

Propriétés

f réelle ⇐⇒ ∀k, fˆ−k = fˆk ∗

(f\+ g )n = fˆn + ĝn
c n = λfˆn
λf

November 29, 2021 5 / 21

Propriétés

f réelle ⇐⇒ ∀k, fˆ−k = fˆk ∗

(f\+ g )n = fˆn + ĝn
c n = λfˆn
λf
g (x) = f (x + δx) ⇐⇒ ĝn = fˆn e 2inπδx → translation de l’espace
physique ⇐⇒ rotation dans l’espace spectral
R1
(g ∗ f )(x) = 0 g (x − t)f (t)dt

November 29, 2021 5 / 21

Propriétés

f réelle ⇐⇒ ∀k, fˆ−k = fˆk ∗

(f\+ g )n = fˆn + ĝn
c n = λfˆn
λf
g (x) = f (x + δx) ⇐⇒ ĝn = fˆn e 2inπδx → translation de l’espace
physique ⇐⇒ rotation dans l’espace spectral
R1
(g ∗ f )(x) = 0 g (x − t)f (t)dt
h = g ∗ f ⇐⇒ ĥn = fˆn × ĝn
convolution physique ⇐⇒ produit espace spectral

November 29, 2021 5 / 21

Propriétés

Orthogonalité Z 1
e 2iπkx e −2iπpx dx = δkp
0

November 29, 2021 6 / 21

Propriétés

Orthogonalité Z 1
e 2iπkx e −2iπpx dx = δkp
0

Identité de Parseval
Z 1 ∞
X
|f (x)|2 dx = |fˆn |2
0 n=−∞

Projection de f sur la base des polynômes trigonométriques {e 2inπx }

November 29, 2021 6 / 21

Fonctions de base

1 2inπx 1 −2inπx
f (x) = cos(2πnx) = e + e
2 2
i 2inπx i
g (x) = sin(2πnx) = − e + e −2inπx
2 2
fˆ1 = fˆ−1 = 1/2, fˆk = 0 pour kkk =
6 1
ĝ1 = −ĝ−1 = −i/2, ĝk = 0 pour kkk =
6 1

November 29, 2021 7 / 21

Coefficients de Fourier d’une fonction réelle

∞
X
f (x) = fˆn e 2iπnx
n=−∞

Euler e 2iπnx
= cos(2πnx) + isin(2πnx)
X∞
=⇒ f (x) = an cos(2πnx) + bn sin(2πnx)
n=0
where an = 2Re[fˆn ], bn = −2Im[fˆn ]

November 29, 2021 8 / 21

Exemple

1 1
5
0.8 0.8

4 0.6 0.6

0.4 0.4
3
0.2 0.2
signal

1
= 1 0

-0.2
-0.2 0

-0.2

-0.4 -0.4
0
-0.6 -0.6

-1 -0.8 -0.8

-1 -1
-2 0 0.2 0.4 0.6 0.8 1 0 0.2 0.4 0.6 0.8 1
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1
time

0.8

0.6

0.4

0.2

+1 0

-0.2

-0.4

-0.6

-0.8

-1
0 0.2 0.4 0.6 0.8 1

1 1

0.8 0.8

0.6 0.6

0.4 0.4

0.2 0.2

-2 0

-0.2
+1 0

-0.2

-0.4 -0.4

-0.6 -0.6

-0.8 -0.8

-1 -1
0 0.2 0.4 0.6 0.8 1 0 0.2 0.4 0.6 0.8 1

November 29, 2021 9 / 21

Transformée de Fourier d’un créneau

Z 0.75
1
Ĥk = e −2iπkx dx = [e −2iπkx ]0.25
0.75
0.25 −2iπk
1 1 3 iπk 1 sin(k π2 ) 1 π
= − [e 2 − e 2 iπk ] = = sinc(k )
πk 2i kπ 2 2
(−1)k+1
Ĥk = si k impair
πk
Ĥk = 0 si k pair

November 29, 2021 10 / 21

Transformée de Fourier en deux dimensions
XX
f (x, y ) = fˆnm e 2iπ(nx+my )
n m
avec Z 1Z 1
fˆnm = f (x, y )e −2iπ(nx+my ) dxdy
0 0

November 29, 2021 11 / 21

Coefficients de Fourier

fˆ0 : valeur moyenne (constante)

November 29, 2021 12 / 21

Coefficients de Fourier

fˆ0 : valeur moyenne (constante)

Définition: Fonction à bande limitée
∃N > 0, fˆp = 0 pour p > N

November 29, 2021 12 / 21

Coefficients de Fourier

fˆ0 : valeur moyenne (constante)

Définition: Fonction à bande limitée
∃N > 0, fˆp = 0 pour p > N
Comment représenter fˆn ?
Définition: Densité d’énergie spectrale

Sn = |fˆn |2

November 29, 2021 12 / 21

Coefficients de Fourier

fˆ0 : valeur moyenne (constante)

Définition: Fonction à bande limitée
∃N > 0, fˆp = 0 pour p > N
Comment représenter fˆn ?
Définition: Densité d’énergie spectrale

Sn = |fˆn |2

Comportement asymptotique de fˆn quand |n| → ∞?

November 29, 2021 12 / 21

Décroissance des coefficients

Z 1
fˆn = f (x)e −2iπnx dx
0

November 29, 2021 13 / 21

Décroissance des coefficients

Z 1
fˆn = f (x)e −2iπnx dx
0

Si f est dérivable, alors en intégrant par parties u 0 = e −2iπnx , v = f (x)

Z Z
u 0 v = [uv ] − uv 0

November 29, 2021 13 / 21

Décroissance des coefficients

Z 1
fˆn = f (x)e −2iπnx dx
0

Si f est dérivable, alors en intégrant par parties u 0 = e −2iπnx , v = f (x)

Z Z
u 0 v = [uv ] − uv 0

Z 1
1 1
fˆn = [ e −2iπnx f (x)]10 − e −2iπnx f 0 (x)dx
−2iπn 0 −2iπn
Le terme de bord s’annule

November 29, 2021 13 / 21

Décroissance des coefficients
Z 1
1
fˆn = e −2iπnx f 0 (x)dx
2iπn 0

November 29, 2021 14 / 21

Décroissance des coefficients
Z 1
1
fˆn = e −2iπnx f 0 (x)dx
2iπn 0
Si f est deux fois dérivable et si f0 est périodique, alors en intégrant à
nouveau par parties,
Z 1
1 1
fˆn = [ e −2iπnx f 0 (x)]10 + e −2iπnx f 00 (x)dx
(−2iπn)2 0 (−2iπn)2

November 29, 2021 14 / 21

Décroissance des coefficients
Z 1
1
fˆn = e −2iπnx f 0 (x)dx
2iπn 0
Si f est deux fois dérivable et si f0 est périodique, alors en intégrant à
nouveau par parties,
Z 1
1 1
fˆn = [ e −2iπnx f 0 (x)]10 + e −2iπnx f 00 (x)dx
(−2iπn)2 0 (−2iπn)2

Si f est p fois dérivable et si f (k) est périodique pour k ≤ p alors

Z 1
1
ˆ
fn = p
e −2iπnx f (p) (x)dx
0 (−2iπn)

November 29, 2021 14 / 21

Décroissance des coefficients
Z 1
1
fˆn = e −2iπnx f 0 (x)dx
2iπn 0
Si f est deux fois dérivable et si f0 est périodique, alors en intégrant à
nouveau par parties,
Z 1
1 1
fˆn = [ e −2iπnx f 0 (x)]10 + e −2iπnx f 00 (x)dx
(−2iπn)2 0 (−2iπn)2

Si f est p fois dérivable et si f (k) est périodique pour k ≤ p alors

Z 1
1
ˆ
fn = p
e −2iπnx f (p) (x)dx
0 (−2iπn)

fn décroit comme n−p lorsque n → ∞

November 29, 2021 14 / 21

Décroissance des coefficients
Z 1
1
fˆn = e −2iπnx f 0 (x)dx
2iπn 0
Si f est deux fois dérivable et si f0 est périodique, alors en intégrant à
nouveau par parties,
Z 1
1 1
fˆn = [ e −2iπnx f 0 (x)]10 + e −2iπnx f 00 (x)dx
(−2iπn)2 0 (−2iπn)2

Si f est p fois dérivable et si f (k) est périodique pour k ≤ p alors

Z 1
1
ˆ
fn = p
e −2iπnx f (p) (x)dx
0 (−2iπn)

fn décroit comme n−p lorsque n → ∞

f est C ∞ ⇐⇒ fˆn décroit plus vite que toute puissance algébrique de n
=⇒ convergence spectrale.
November 29, 2021 14 / 21
Décroissance des coefficients

Régularité ⇐⇒ Décroissance à l’infini

f =1 ⇐⇒ fˆn = δn

⇐⇒

November 29, 2021 15 / 21

Décroissance des coefficients

Décroissance à l’infini ⇐⇒ Régularité

δ∗f =f ⇐⇒ δ̂n fˆn = fˆn

⇐⇒

November 29, 2021 16 / 21

La transformée de Fourier d’une Gaussienne est une
Gaussienne
f ∈ L2 (R) Z ∞
f (t) = fˆk e 2iπkx dk
−∞
Z ∞
fˆk = f (x)e −2iπkx dx
−∞

2 π 1/2 −π 2 k 2 /a
f (x) = e −ax fˆk = a e
November 29, 2021 17 / 21
Dérivation

∞
X
f (x) = fˆn e 2iπnx
n=−∞
∞
X
=⇒ f 0 (x) = (2inπ)fˆn e 2iπnx
n=−∞

Dérivation physique ⇐⇒ multiplication spectrale

=⇒ résolution spectrale d’EDP devient problème algébrique

November 29, 2021 18 / 21

Résolution spectrale d’équations - I

d 2f
= g sur ]0, 1[(E )
dx 2
avec g(x) connue, de carré intégrable
et des conditions aux limites périodiques pour f
On cherche f sous la forme d’une série de Fourier:
∞
X
f (x) = fˆn e 2iπnx , fˆn inconnus
n=−∞

∞
X
g (x) = ĝn e 2iπnx , ĝn connus
n=−∞
Z 1
ĝn = g (x)e −2iπnx dx
0

November 29, 2021 19 / 21

Résolution spectrale d’équations - II

∞
df X
= (2iπn)2 fˆn e 2iπnx
dx 2 n=−∞

November 29, 2021 20 / 21

Résolution spectrale d’équations - II

∞
df X
= (2iπn)2 fˆn e 2iπnx
dx 2 n=−∞
R1
Pour p entier, on calcule 0 (E )e −2iπpx dx

November 29, 2021 20 / 21

Résolution spectrale d’équations - II

∞
df X
= (2iπn)2 fˆn e 2iπnx
dx 2 n=−∞
R1
Pour p entier, on calcule 0 (E )e −2iπpx dx
Z 1
g (x)e −2iπpx dx = ĝp
0
1
d 2 f −2iπpx
Z
e dx = (2iπp)2 fˆp
0 dx 2
ĝp
=⇒ fˆp = − 2 2
4π p

November 29, 2021 20 / 21

Résolution spectrale d’équations - II

∞
df X
= (2iπn)2 fˆn e 2iπnx
dx 2 n=−∞
R1
Pour p entier, on calcule 0 (E )e −2iπpx dx
Z 1
g (x)e −2iπpx dx = ĝp
0
1
d 2 f −2iπpx
Z
e dx = (2iπp)2 fˆp
0 dx 2
ĝp
=⇒ fˆp = − 2 2
4π p
∞
X ĝn 2iπnx
=⇒ f (x) = − 2 n2
e
n=−∞
4π

November 29, 2021 20 / 21

Conclusion

La Transformée de Fourier est un outil pour

l’analyse du signal
la résolution d’équations
Comment la mettre en oeuvre en pratique ?

November 29, 2021 21 / 21