UE MATH 332.

EDP/EDO Suite

La méthode de séparation des variables

1 Introduction
On cherche à déterminer une fonction u = u(x, t) des variables x ∈ R et t ≥ 0 vérifiant le système

 ut = kuxx , 0 < x < l, t > 0
u(x, 0) = f (x), 0 < x < l
u(0, t) = u(l, t) = 0


On rencontre ce système en énergétique:

• u = u(x, t) représente alors la température ( au point d’abscisse x à l’instant t) le long d’une barre homogène de
longueur l.
• L’edp ut = kuxx est appelée équation de la chaleur à une dimension
• La condition initiale u(x, 0) = f (x) représente la température de la barre à l’instant t = 0
• Les relations u(0, t) = u(l, t) = 0 indiquent que la température est nulle aux extrémités x = 0 et x = l de la
barre à tout temps t.
On veut déterminer u sous la forme u(x, t) = X(x)T (t), on sépare donc les variables x et t.

2 Équations vérifiées par X et T

Avec u(x, t) = X(x)T (t) on a les relations:
ut = XT 0 (t) et uxx = X 00 (x)T (t)
où

dT d2 X
T 0 (t) = et X 00 (x) =
dt dx2
La relation ut = kuxx devient alors
T 0 X = kX 00 T
Et en séparant les variables
T0 X 00
=
kT X
A gauche de cette inégalité, on a une fonction de t seulement, et à droite une fonction de x seulement s’il y’a égalité,
nécessairement ces deux quantités sont constantes. On peut donc poser,
T0 X 00
= = −λ, λ∈R
kT X
et en déduire les équations suivantes
T 0 + (λk)T = 0 et X 00 + λX = 0
L’équation T 0 + (λk)T = 0 a pour solution générale
T (t) = Ce−λkt

2.1 X vérifie un problème aux limites

Nous allons voir comment se traduisent les conditions aux limites u(0, t) = u(l, t) = 0.

u(0, t) = 0 ⇔ X(0)T (t) = 0 ⇔ X(0) = 0

De même
u(l, t) = 0 ⇔ X(l)T (t) = 0 ⇔ X(l) = 0
La fonction X vérifie donc le système
X 00 + λX

= 0, 0 < x < l
(P L)
X(0) = X(l) = 0
qui est un problème aux limites.
On va déterminer les valeurs de λ pour lesquelles ce système (problème aux limites) admet des solutions non triviales.
Définition
Les valeurs de λ pour lesquelles le système (P L) possède des solutions non triviales (non identiquement nulles) sont
appelées valeurs propres et les solutions correspondantes des fonctions propres.

2.2 Solutions du problème aux limites

√ √
1. Si λ < 0, on a X(x) = Ae −λx + Be− −λx .
Les relations X(0) = X(l) = 0 conduisent à A = B = 0 et donc X(x) = 0

2. Si λ = 0, on a X(x) = Ax + B
Les relations X(0) = X(l) = 0 conduisent à A = B = 0 et donc X(x) = 0
√ √
3. Si λ > 0, on a X(x) = A cos( λx) + B sin( λx). √
Les relations X(0) = X(l)
√ = 0 conduisent à A = 0 et B sin( λl) = 0. Pour avoir des solutions non triviales; il
est nécessaire que sin( λl) = 0
√
2
4. L’équation sin( λl) = 0 admet une infinité dénombrable de solutions λn = nπ
l , n ∈ N; n ≥ 1
5. Pour chaque valeur de n ≥ 1, on obtient les fonctions
2  nπx 
Tn (t) = e−kλn t = e−k( )
nπ
t
l Xn (x) = sin
l

6. Pour chaque valeur de n ≥ 1, la fonction

un (x, t) = Tn (t)Xn (x)

vérifie le système linéaire 

ut = kuxx , 0 < x < l t > 0
u(0, t) = u(l, t) = 0

7. D’après le principe de superposition toute combinaison linéaire de solutions est encore une solution, par exemple
u1 + u2 + · · · + uN est encore une solution.

8. On pose alors:
) t sin nπx
2  
An e−k(
X X nπ
u(x, t) = un (x, t) = l
l
9. La fonction u ainsi définie vérifie le système

ut = kuxx , 0 < x < l t > 0
u(0, t) = u(l, t) = 0

10. Pour que cette fonction vérifie la condition initiale u(x, 0) = f (x), il est nécessaire que:
X  nπx 
An sin = f (x)
l

11. Cette dernière relation est en fait le développement en série de sinus de la fonction f et par conséquent; on a
Z l
2  nπx 
An = f (x) sin dx
l 0 l

12. On peut donc conclure que la fonction u définie par la relation

Z l
) t sin nπx 2
2    nπx 
An e−k(
X nπ
u(x, t) = l avec An = f (x) sin dx
l l 0 l

vérifie le système 
 ut = kuxx , 0 < x < l t > 0
u(x, 0) = f (x), 0 < x < l
u(0, t) = u(l, t) = 0

2.3 Exemples
Exemple 1

 ut = kuxx , 0 < x < 1 t > 0
u(x, 0) = sin(πx), 0 < x < 1
u(0, t) = u(1, t) = 0


Ici, l = 1 et f (x) = sin(πx)

On sait que X 2
u(x, t) = An e−k(nπ) t sin (nπx)
On détermine An par la condition initiale
X
u(x, 0) = f (x) ⇔ An sin(nπx) = sin(πx)
n≥1

On en déduit par identification que An = 1 et An = 0 si n ≥ 2.

Par conséquent
2
u(x, t) = e−k(π) t sin (πx)
On peut observer que
u(x, t) → 0 quand t → +∞
On remarque dans cet exemple que l’on n’a pas besoin d’utiliser la formule définissant An comme coefficient du
développement en série de sinus de f (x); car f (x) = sin(πx) est déjà un développement en série de sinus.

Exemple 2

 ut = kuxx , 0 < x < π t > 0
u(x, 0) = 5 sin(3x) + 7 sin(10x), 0 < x < π
u(0, t) = u(π, t) = 0


Ici, l = π et f (x) = 5 sin(3x) + 7 sin(10x)

On sait que X 2
u(x, t) = An e−kn t sin (nx)
On détermine An par la condition initiale
X
u(x, 0) = f (x) ⇔ An sin(nx) = 5 sin(3x) + 7 sin(10x)
n≥1

On en déduit par identification que A3 = 5; A10 = 7 et An = 0 si n 6= 3 et n 6= 10.

En conclusion
2 2
u(x, t) = A3 e−3 kt
sin(3x) + A10 e−k10 t sin(10x) = 5e−9kt sin(3x) + 7e−100kt sin(10x)

Exemple 3

ut = 3uxx
 0 < x < π, t > 0
u(0, t) = 0 et u(l, t) = 0

u(x, 0) = f (x)


Ici, k = 3, l = π, f (x) = 3 sin 2x − 6 sin 5x La solution est alors donnée par

∞
X 2
u(x, t) = An e−3n t sin nx
0

Les coefficients An sont déterminés par la condition initiale u(x, 0) = f (x) c’est-à-dire
∞
X
u(x, t) = An sin nx = 3 sin 2x − 6 sin 5x
0

On obtient par identification A2 = 3, A5 = −6 et An = 0 lorsque n 6= 2, 5. La solution est alors donnée par

2 2
u(x, t) = A2 e−3.2 2t
sin(2x) + A5 e−3.5 t sin(5x) = 3e−12t sin(2x) − 6e−75t sin(5x).
Exemple 4

 ut = uxx , 0 < x < 1, t > 0
u(x, 0) = x, 0<x<1
u(0, t) = u(1, t) = 0


Ici l = 1, k = 1 et f (x) = x
On sait que X 2
u(x, t) = An e−(nπ) t sin (nπx)
On détermine An par la condition initiale
X
u(x, 0) = f (x) ⇔ An sin(nπx) = x
n≥1

On déduit
1
2(−1)n+1
Z
An = 2 x sin (nπx) dx =
0 nπ
Par conséquent
2 X (−1)n+1 −n2 π2 t
u(x, t) = e sin (nπx)
π n

3 Autres problèmes
1. Il existe d’autres types de conditions aux limites.
• ux (0, t) = 0 et ux (l, t) = 0
• u(0, t) = 0 et ux (l, t) = 0
• u(0, t) = u(l, t) et ux (0, t) = ux (l, t)
2. La méthode ci-dessus est inopérante pour résoudre un problème avec conditions aux limites non homogènes
c’est-à-dire par exemple 

 ut = kuxx , 0 < x < l, t > 0
u(x, 0) = f (x), 0 < x < l


 u(0, t) = g(t)
u(l, t) = h(t)


lorsque g et h sont non nulles. On doit donc commencer par opérer un changement de fonction inconnue pour
se ramener à une EDP avec conditions aux limites homogènes. Ce changement de variables consiste à poser
u = v + w où v est la nouvelle fonction inconnue et w une fonction à déterminer. Lorsque les fonctions g et h
sont constantes, w est la solution dite à l’état stable c’est-à-dire indépendante du temps. Par exemple, pour le
système 

 ut = kuxx , 0 < x < l, t > 0
u(x, 0) = f (x), 0 < x < l

(S1 )

 u(0, t) = A
u(l, t) = B


w vérifie le système
 wxx = w00

= 0, 0 < x < l, t > 0
w(0) = A
u(l) = B


4 Travaux dirigés
4.1 Solution d’équilibre
Déterminer la solution à l’état stable (équilibre) de l’équation de la chaleur
ut = kuxx vérifiant les conditions aux limites données.
1. u(0, t) = 200, u(10, t) = 100
2. u(0, t) = 100, ux (10, t) = 50
3. ux (0, t) = 8, u(10, t) = 100
4. ux (0, t) = 30, ux (10, t) = 10
 5. u(0, t) + ux (0, t) = 10, u(10, t) = 100
6. u(0, t) + ux (0, t) = 0, u(10, t) − ux (10, t) = 200
7. u(0, t) − 10ux (0, t) = 30, u(10, t) − 5ux (10, t) = 0

solutions
Il s’agit de déterminer une solution w de l’équation de la chaleur indépendante du temps (wt = 0)et vérifiant les
conditions aux limites.
1. wxx = 0, w(0) = 200, w(10) = 100, on trouve : w(x) = −10x + 200
2. wxx = 0, w(0) = 100, w0 (10) = 50, on trouve : w(x) = 50x + 100
3. wxx = 0, w0 (0) = 8, w(10) = 100, on trouve : w(x) = 8x + 20
4. w n’existe pas
5. wxx = 0, w(0) + w0 (0) = 10, w(10) = 100, on trouve : w(x) = 10x
6. wxx = 0, w(0) + w0 (0) = 0, w(10) − w0 (10) = 200, on trouve : w(x) = 25x − 25
7. wxx = 0, w(0) − 10w0 (0) = 30, w(10) − 5w0 (10) = 200, on trouve : w(x) = −2x + 10

4.2 Résolution explicite

1. 2.1 Résoudre l’équation de la chaleur ci-dessous pour la condition initiale donnée
(
ut = 2uxx 0 < x < 6, t > 0
ux (6, t) = 0 et ux (6, t) = 0

(a) u(x, 0) = π + 3 cos(πx) − 4 cos( 3πx

2 ) − cos(3πx)
(b) u(x, 0) = 4
(c) u(x, 0) = x2
(d) u(x, 0) = 0
2.2 Pour chacune des solutions déterminée en 1.1, calculer

lim u(x, t)
t→+∞

solutions
(a) La solution est alors donnée par
∞
X 2
u(x, t) = An e−3n t sin nx
0

Les coefficients An sont déterminés par la condition initiale u(x, 0) = f (x) c’est-à-dire
∞
X
u(x, t) = An sin nx = 3 sin 2x − 6 sin 5x
0

On obtient par identification A2 = 3, A5 = −6 et An = 0 lorsque n 6= 2, 5. La solution est alors donnée

par
2 2
u(x, t) = A2 e−3.2 2t sin(2x) + A5 e−3.5 t sin(5x) = 3e−12t sin(2x) − 6e−75t sin(5x).
(b) La méthode de séparation de variables: u(x, t) = X(x)T (t) conduit au problème aux limites

X 00 + λX = 0, X 0 (0) = X 0 (6) = 0

et à l’edo T 0 + 2λT = 0.
nπ 2
n ∈ N , les fonctions propres Xn (x) = cos nπx


On déduit les valeurs propres λn = 6 6 ainsi que
2 2
− 2n36π t
Tn (t) = e
La solution est alors donnée par
∞
X 2n2 π 2 t nπx
u(x, t) = An e − 36 cos
0
6
 Les coefficients An sont déterminés par la condition initiale u(x, 0) = f (x) c’est-à-dire
∞
X nπx
An cos = f (x) = dévellopement en série de cosinus de f(x)
0
6

donc: Z 6 Z 6
1 2 nπx
A0 = f (x)dx An = f (x) cos dx
6 0 6 0 6
i.
2 9π 2 t 3πx 2
u(x, t) = π + 3e−2π t cos(πx) − 4e− 2 cos( ) − e−18π t cos(3πx)
2
ii. u(x, t) = 4
iii.
iv. u(x, t) = 0
(c) 2.1 Résoudre l’équation de la chaleur ci-dessous pour la condition initiale donnée après avoir transformé
les conditions aux limites en données nulles.
(
ut = 2uxx 0 < x < 6, t > 0
u(0, t) = 40 et u(6, t) = 10

(a) u(x, 0) = −5x + 40 + 5 sin(2πx) − 2 sin( 5πx

2 )
(b) u(x, 0) = 0
2.2 Pour chacune de solutions trouvée en 2.1, déterminer

lim u(4, t)
t→+∞

Solutions
On détermine d’abord la solution d’équilibre w(x). Ici, w(x) = −5x + 40.
On pose ensuite u(x, t) = v(x, t) + w(x). v vérifie l’équation de la chaleur avec conditions aux limites limites et
condition initiale v(x, 0) = u(x, 0) − w(x)
(a)
2 25π 2 t 5πx
u(x, t) = −5x + 40 + 5e−8π t sin(2πx) − 2e− 2 sin( )
2
(b)
2. On considère le problème de la conduction de la chaleur.

ut = 9uxx
 0 < x < 4, t > 0
u(0, t) = 32 et u(4, t) = 32

u(x, 0) = f (x)


a. Transformer les conditions aux limites en données nulles.

b. Déterminer l’expression de la solution générale
c. Déterminer la limite de u(x, t) lorsque t → +∞

solutions
a. On détermine w vérifiant w00 = 0, w(0) = 32, w(4) = 32 et on pose u = v + w
b. 2
 nπx 
π 2 t/16
X
u(x, t) = 32 + An e−n sin
4
1
R4 nπx


avec An = 2 0
(f (x) − 32) sin 4 dx
c. lim u(x, t) = 32

3. On considère le problème de la conduction de la chaleur, avec les conditions aux limites ci-dessous.

ut = 9uxx , 0 < x < 10, t > 0

 a) u(0, t) = 0, et u(10, t) = 0
b) ux (0, t) = 0, et ux (10, t) = 0
c) u(0, t) = 0, et u(10, t) = 100
d) u(0, t) = 100, et u(10, t) = 50
Pour une condition initiale commune u(x, 0) = 300, déterminer la température au milieu de la barre après un
temps très long. Quel ensemble de conditions initiales permet d’obtenir la température la plus élevée en ce point.
4. On veut résoudre le système suivant:

ut = uxx − cu 0 < x < l, t > 0

ux (0, t) = 0 et u(l, t) = 0 t > 0

u(x, 0) = u0 (x) 0 ≤ x ≤ l


où c ∈ R est une constante réelle et u0 une fonction de classe C 1

a. Déterminer µ ∈ R de sorte que, par le changement de fonction inconnue v(x, t) = eµt u(x, t), on se ramène à
l’équation vt = vxx .
b. Écrire le nouveau système pour v.
c. Trouver la solution du nouveau système par la méthode de séparation des variables. En déduire la solution
u(x, t) du système initial.

5 Méthode de séparation des variables pour l’équation des ondes

5.1 Position du problème
On considère le problème suivant qui peut modéliser les vibrations transversales d’une corde de longueur l fixée à ses
extrémités. 

 ut = c2 uxx ; 0 < x < l, t > 0
u(x, 0) = f (x); 0 ≤ x ≤ l


 u t (x, 0) = g(x); 0 ≤ x ≤ l
u(0, t) = u(l, t) = 0 t > 0


5.2 Problème aux limites associé

1. On pose u(x, t) = X(x)T (t), les relations suivantes sont vérifiées:
T 00 X 00
= = −λ λ constante
c2 T X
2. La fonction T vérifie
T 00 + λc2 T = 0
et la fonction X vérifie le problème aux limites
X 00 + λX = 0

X(0) = X(l) = 0

3. Les valeurs propres et les fonctions propres du problème aux limites ci-dessus sont données par
 nπ 2  nπx 
λn = et Xn (x) = sin n≥1
l l
4. Pour n ≥ 1 fixé, on a comme solution générale
 nπx     
nπct  nπx  nπct
un (x, t) = An sin cos + Bn sin sin
l l l l

5. En utilisant le principe de superposition à l’extrême

X X  nπx  
nπct
  nπx  
nπct

u(x, t) = un (x, t) = An sin cos + Bn sin sin
l l l l
n≥1

6. En utilisant les conditions initiales u(x, 0) = f (x) et ut (x, 0) = g(x), on obtient

2 l
Z  nπx  Z l
2  nπx 
An = f (x) sin dx et Bn = g(x) sin dx
l 0 l nπc 0 l
6 Travaux dirigés
6.1
Résoudre l’équation des cordes vibrantes pour les données initiales indiquées
(
utt = 4uxx , 0 < x < π, t > 0
(1)
u(0, t) = 0, u(π, t) = 0

1. u(x, 0) = 12 sin(2x) − 16 sin(5x) + 24 sin(6x)

ut (x, 0) = 0
2. u(x, 0) = 0
ut (x, 0) = 6

3. u(x, 0) = 0
ut (x, 0) = 12 sin(2x) − 16 sin(5x) + 24 sin(6x)
4. u(x, 0) =
ut (x, 0) = 12 sin(2x) − 16 sin(5x) + 24 sin(6x)

6.2
Résoudre l’équation des cordes vibrantes


 utt = 100uxx , 0 < x < 2, t > 0

 u(0, t) = 0, u(2, t) = 0


u(x, 0) = 32 sin(πx) + e2 sin(3πx) + 25 sin(6πx) (2)


 ut (x, 0) = 6 sin(2πx) − 16 sin( 5πx )



2

6.3
Résoudre l’équation des cordes vibrantes

utt = 25uxx , 0 < x < 2, t > 0




 u(0, t) = 0, u(2, t) = 0

(3)


 u(x, 0) = x − x2

ut (x, 0) = π

6.4
Utiliser la méthode de séparation des variables pour résoudre le problème des cordes pour une corde non fixée à ses
extrémités.
utt = α2 uxx , 0 < x < l, t > 0




 u (0, t) = 0, u (l, t) = 0
x x
(4)


 u(x, 0) = f (x)

ut (x, 0) = g(x)
Si vous avez des questions, contacter moi via l’adresse mail suivante: [email protected]