TD I : Manipulation des transformées de Fourier

et Laplace en physique

Promotion 140
Mardi 5 Octobre 2021

L'objectif de ce TD est de vous montrer que les transformations de Fourier et de Laplace sont des outils très
utiles en physique. On se propose de reprendre des exemples d'équations dans diérents domaines an de montrer
l'universalité de leur usage.

Mécanique du point

Dans cet exercice, nous allons étudier le système masse-ressort suivant :

On s'intéresse à la trajectoire r(t) d'une particule chargée (de charge q et de masse m) dans un champ électro-
magnétique statique (E = E0 z et B = B0 y) avec comme condition initiale r(0) = 0 et ṙ(0) = 0. On doit donc
résoudre l'équation :

mr̈ = q(E + ṙ × B) (1)

1. Établir les 3 équations du mouvement vériées par x(t), y(t) et z(t).

Solution :


 mẍ + qB0 ż = 0
mÿ = 0
mz̈ − qB0 ẋ = qE0


2. Montrer que la particule n'aura pas de mouvement suivant y.

Solution :
La deuxième équation précédente portant sur y(t) avec les conditions initiales y(0) = 0 et ẏ(0) = 0 nous
donne bien y(t) = 0 ∀t.

3. Prendre la transformée de Laplace du système restant. On posera Ω = qB0

m et α = − qE
m pour alléger les
0

notations.
Solution :

mp2 x̃ + qB0 pz̃ = 0


mp2 z̃ − qB0 px̃ = qEp
0

En introduisant les notations de l'énoncé ce système se simplie en :


px̃ + Ωpz̃ = 0
α
Ωpz̃ − p2 z̃ = p

4. Résoudre le système précédent.

Solution :
Après quelques étapes de calcul, on arrive facilement à :
(
ωα
x̃ = p2 (p2 +Ω2 )
α
z̃ = p(p2 +Ω2 )

1
5. Décomposer ces transformées de Laplace en éléments simples pour nalement trouver la solution dans l'espace
temporel.
Solution :
La technique pour décomposer une fraction rationnelle en éléments simples est tout à fait classique. Elle
consiste dans un premier temps à écrire la décomposition avec des coecients :
ωα A B
= 2+ 2
p2 (p2 2
+Ω ) p p + Ω2

Puis, pour trouver A on multiple les 2 côtés de l'équation par le dénominateur de A et on prend la limite
de l'équation quand p est égale à la racine du dénominateur de A. Soit :
ωα A Bp2
lim = lim
p→0 p2 + Ω2 p→0 + p2 + Ω2

D'où on tire A = Ω2 .
ωα
On réitère cette technique pour chacun des pôles de la fraction rationnelle. Soit ici :
ωα A(p2 + Ω2 )
lim 2
= 2 lim 2 +B
p →−Ω2
2 p p →−Ω p2
Et cette fois, on obtient B = Ω2 .
− ωα Au nal, on a donc :
  
α 1 1
 x̃ =
Ω2 2 − p2 +Ω2
p
α 1 p
 z̃ = Ω2 p − p2 +Ω2

D'où on tire facilement : α

x(t) = (Ωt − sin(Ωt))
Ω2
et : α
z(t) = (1 − cos(Ωt))
Ω2

Équation d'onde avec conditions initiales

Dans cet exercice, nous allons utiliser conjointement la transformée de Fourier et la transformée de Laplace. Nous
nous intéressons à la résolution de l'équation d'onde à une dimension (par exemple, la corde vibrante) avec des
conditions initiales. Il s'agit donc de résoudre l'équation :
∂ 2 u(x, t) ∂ 2 u(x, t)
2
− c2 =0 (2)
∂t ∂x2
avec les conditions initiales : u(x, 0) = f (x) et ∂t u(x, 0) = g(x).
6. Établir la transformée de Laplace par rapport au temps de l'équation (2).
Solution :
Comme on cherche une solution temporelle suite à des conditions initiales c'est la transformée de Laplace
qui est l'outil le plus approprié. Le calcul (directe application du cours) donne simplement :
∂ 2 û(x, p)
p2 û(x, p) − pf (x) − g(x) − c2 =0
∂x2
7. Prendre la transformée de Fourier par rapport à la variable x de l'équation précédente.
Solution :
Cette fois-ci nous n'imposons pas de conditions spéciques en x = 0, l'outil le plus adéquat pour la variable
x reste la transformée de Fourier. Le calcul (directe application du cours) donne simplement :

˜ p) = p 1
û(k, f˜(k) + 2 2 g̃(k)
c2 k 2 + p2 c k + p2

2
8. Prendre la transformée de Laplace inverse de l'équation précédente.

Solution :
Soit on connaît  par c÷ur la TL des cosinus et sinus, soit on décompose en fraction rationnelle élémentaire
pour retrouver des exponentielles. Le calcul (directe application du cours encore) donne simplement :
sin(ckt)
ũ(k, t) = f˜(k) cos(ckt) + g̃(k)
ck

9. Obtenir nalement la solution de l'équation d'onde avec conditions initiales.

Solution :
Calculons les transformées de Fourier inverses de chaque terme séparément :
h i h i
F −1 f˜(k) cos(ckt) = F −1 f˜ ∗ F −1 [cos(ckt)]
 
1 1
= f (x) ∗ δ(x − ct) + δ(x + ct)
2 2
= f (x − ct) + f (x + ct)

   
sin(ckt) 1 sin(ctk)
F −1 g̃(k) = F −1 2ct ∗ F −1 [g̃(k)]
ck 2c ctk
1
= Π2ct ∗ g(x)
2c
Z x+ct
1
= g(x0 )dx0
2c x−ct

Et on retrouve le résultat bien connu :

Z x+ct
1 1 1
u(x, t) = f (x − ct) + f (x + ct) + g(x0 )dx0
2 2 2c x−ct

Équation de la chaleur

On cherche à déterminer la solution u(x, tde l'équation de la chaleur :

∂T ∂2T
=a 2 (3)
∂t ∂x
avec la condition initiale : T (x, 0) = T0 (x) où T0 (x) est la température à l'instant t = 0.
10. Prendre la transformée de Fourier par rapport à la variable x de l'équation puis résoudre cette équation
diérentielle.
Solution :

∂ T̃
= −ak 2 T̃
∂t
D'où on tire : 2 2
T̃ (k, t) = C e−ak t

Avec la condition initiale sur T (x, t) on a nalement T̃ (k, t) = T˜0 (k)e−ak

2 2
t

11. Revenir dans l'espace réel pour retrouver la solution spatio-temporelle.

3
Solution :
Il faut se rappeler la propriété (qui vient de la convention de la transformation de Fourier choisie) :
F [f ∗ g] = F [f ] F [g]

Avec dans notre cas f = T( x) et g telle que F [g] = e−ak t . Par transformation de Fourier inverse on
2 2

montre :
1 x2
g(x, t) = √ e− 4at
4πat
Finalement, on obtient la solution :
Z ∞
1 (x−x0 )2
T (x, t) = √ T0 (x0 )e− 4at ddx0
4πat −∞