Travaux dirigés de Mathématiques

Ecole Nationale Supérieure de Géologie

1ère année
Didier Maquin
http://w3.cran.univ-lorraine.fr/perso/didier.maquin/fr/
 1ère année E.N.S.G
Travaux Dirigés en Mathématiques
Fonction d’une variable complexe

Exercice 1
Démontrer que :
a) z1 + z2 = z1 + z2
b) |z1 z2 | = |z1 ||z2 |
c) eiθ = ei(θ+2kπ) , k = 0 ± 1, ±2 . . .
d) |z1 + z2 | ≤ |z1 | + |z2 |, analytiquement et graphiquement

Exercice 2
Déterminer et représenter dans le plan complexe les racines cubiques de −1 + i.

Exercice 3
Démontrer l’équivalence des opérateurs suivants :
∂ ∂ ∂
a) = +
∂x ∂z
 ∂ z̄ 
∂ ∂ ∂
b) =i −
∂y ∂z ∂ z̄
si z = x + iy et z̄ = x − iy

Exercice 4
Montrer que :
∂ ∂ ∂
a) ∇ ≡ +i =2
∂x ∂y ∂ z̄
¯ ≡ ∂ ∂ ∂
b) ∇ −i =2
∂x ∂y ∂z

Exercice 5
Montrer que, en coordonnées polaires, les équations de Cauchy-Riemann prennent la forme :

∂u 1 ∂v ∂v 1 ∂u
= et =−
∂r r ∂θ ∂r r ∂θ

Exercice 6
Z
Evaluer z̄dz de z1 = 0 à z2 = 4 + 2i le long de la courbe C si :
C
a) C est définie par z = t2 + it
b) C est formée des segments 0 à 2i et 2i à 4 + 2i

2
 1ère année E.N.S.G
Travaux Dirigés en Mathématiques
Fonction d’une variable complexe

Exercice 1
√ p
Posons z = x + iy avec i = −1, z̄ = x − iy et |z| = x2 + y 2 .
a) Posons : z1 = x1 + iy1 et z2 = x2 + iy2 , on a :
z1 + z2 = x1 + iy1 + x2 + iy2 = x1 + x2 + i(y1 + y2 ) = x1 + x2 − i(y1 + y2 )
z1 + z2 = x1 − iy1 + x2 − iy2 = x1 + iy1 + x2 + iy2 = z1 + z2
b) |z1 z2 | = |(x1 + iy1 )(x2 + iy2 )| = |x1 x2 − y1 y2 + i(x1 y2 + y1 x2 )|
p p
|z1 z2 | = (x1 x2 − y1 y2 )2 + (x1 y2 + y1 x2 )2 = (x1 x2 )2 + (y1 y2 )2 + (x1 y2 )2 + (y1 x2 )2
q
q 2
q 2
x21 + y12 x22 + y22 = x1 + y12 x2 + y22 = |z1 ||z2 |



|z1 z2 | =

c) ei(θ+2kπ) = cos(θ + 2kπ) + i sin(θ + 2kπ) = cos(θ) + i sin(θ) = eiθ

p p p
d) |z1 + z2 | ≤ |z1 | + |z2 | ⇔ (x1 + x2 )2 + (y1 + y2 )2 ≤ x21 + y12 + x22 + y22
En élevant au carré les deux membres de l’égalité :
p
(x1 + x2 )2 + (y1 + y2 )2 ≤ x21 + y12 + 2 (x21 + y12 )(x22 + y22 ) + x22 + y22
p
x1 x2 + y1 y2 ≤ (x21 + y12 )(x22 + y22 )
En élevant de nouveau au carré cette inégalité, on a :
x21 x22 + 2x1 x2 y1 y2 + y12 y22 ≤ x21 x22 + x21 y22 + y12 x22 + y12 y22
2x1 x2 y1 y2 ≤ x21 y22 + y12 x22
Cette dernière inégalité est équivalente à :
(x1 y2 − x2 y2 )2 ≥ 0
qui est toujours vraie. En remontant le raisonnement, on parvient au résultat.

Graphiquement, le résultat provient du fait que |z1 |, |z2 | et |z1 + z2 | représentent les lon-
gueurs des côtés d’un triangle ; l’inégalité triangulaire précise que la somme des longueurs
de deux côtés d’un triangle est supérieure ou égale à la longueur du troisième côté.

3
Exercice 2
p y
Définition : z = x + iy = ρeiθ = ρ(cos(θ) + i sin(θ)) avec ρ = x2 + y 2 et θ = arctan
√ √ x
−1 + i = 2ei arctan(−1) = 2(cos(3π/4 + 2kπ) + i sin(3π/4 + 2kπ))
    
(−1 + i)1/3 = 21/6 cos 3π/4+2kπ
3 + i sin 3π/4+2kπ
3

k = 0 : z1 = 21/6 (cos (π/4) + i sin (π/4))

k = 1 : z2 = 21/6 (cos (11π/12) + i sin (11π/12))
k = 2 : z3 = 21/6 (cos (19π/12) + i sin (19π/12))

Exercice 3
z + z̄ z − z̄ i(z − z̄)
On a x = et y = =− , d’où :
2 2i 2
∂ ∂ ∂z ∂ ∂ z̄ ∂ ∂
= + = +
∂x ∂z ∂x ∂ z̄ ∂x ∂z ∂ z̄
 
∂ ∂ ∂z ∂ ∂ z̄ ∂ ∂ ∂ ∂
= + = i+ (−i) = i −
∂y ∂z ∂y ∂ z̄ ∂y ∂z ∂ z̄ ∂z ∂ z̄

Exercice 4
Par définition, ∇, opérateur différentiel vectoriel est défini par :
∂ ∂
∇≡ +i
∂x ∂y
D’après l’exercice précédent, on a établit que :
 
∂ ∂ ∂ ∂ ∂ ∂
= + et =i −
∂x ∂z ∂ z̄ ∂y ∂z ∂ z̄

d’où :   
∂ ∂ ∂ ∂
∇= + +i i −
∂z ∂ z̄ ∂z ∂ z̄
∂ ∂ ∂ ∂ ∂
∇= + − + =2
∂z ∂ z̄ ∂z ∂ z̄ ∂ z̄

4
 De la même manière :
¯ ≡ ∂ −i ∂
∇
∂x ∂y
  
¯ ∂ ∂ ∂ ∂
∇= + −i i −
∂z ∂ z̄ ∂z ∂ z̄

¯ = ∂ + ∂ + ∂ − ∂ =2 ∂
∇
∂z ∂ z̄ ∂z ∂ z̄ ∂z

Exercice 5
En coordonnées polaires, on a x = r cos(θ) et y = r sin(θ), d’où :

∂u ∂u ∂x ∂u ∂y ∂u ∂u
= + = cos(θ) + sin(θ)
∂r ∂x ∂r ∂y ∂r ∂x ∂y
∂u ∂u ∂x ∂u ∂y ∂u ∂u
= + = −r sin(θ) + r cos(θ)
∂θ ∂x ∂θ ∂y ∂θ ∂x ∂y
∂u ∂u ∂u ∂u
Les dérivées et sont donc solutions d’un système linéaire d’équations en et .
∂r ∂θ ∂x ∂y
On a :
∂u ∂u
∂r sin(θ) ∂r sin(θ)
∂u ∂u
∂u ∂θ r cos(θ) ∂θ r cos(θ) ∂u 1 ∂u
= = = cos(θ) − sin(θ)
∂x cos(θ) sin(θ) r ∂r r ∂θ
−r sin(θ) r cos(θ)

∂u ∂u
cos(θ) ∂r cos(θ) ∂r
∂u ∂u
∂u −r sin(θ) ∂θ −r sin(θ) ∂θ ∂u 1 ∂u
= = = sin(θ) + cos(θ)
∂y cos(θ) sin(θ) r ∂r r ∂θ
−r sin(θ) r cos(θ)
Les conditions de Cauchy-Rieman en coordonnées cartésiennes :
∂u ∂v
=
∂x ∂y
∂u ∂v
=−
∂y ∂x
deviennent en coordonnées polaires :
∂u 1 ∂u ∂v 1 ∂v
cos(θ) − sin(θ) = sin(θ) + cos(θ)
∂r r ∂θ ∂r r ∂θ

∂u 1 ∂u ∂v 1 ∂v
sin(θ) + cos(θ) = − cos(θ) + sin(θ)
∂r r ∂θ ∂r r ∂θ
et ces égalités sont vérifiées si :
∂u 1 ∂v 1 ∂u ∂v
= et =−
∂r r ∂θ r ∂θ ∂r

5
 1ère année E.N.S.G
Travaux Dirigés en Mathématiques
Fonction d’une variable complexe
Intégrale de Cauchy

Exercice 1
d d
Démontrer que a) ez = ez et b) eaz = eaz où a est une constante quelconque. On vérifiera
dz dz
que les fonctions considérées sont bien dérivables à l’aide des équations de Cauchy-Riemann.

Exercice 2
ez
Z
Calculer dz avec z ∈ C si C est un cercle de centre z0 = −1 et de rayon R tel
C z(1 − z)3
que :
a) R < 1 b) 1 ≤ R ≤ 2 c) R ≥ 2

Exercice 3
eiz
Z
Calculer dz si C est le cercle de centre 0 et de rayon R > 1.
C 1 + z2

Exercice 4
eiz
Z
Calculer 2
dz si Γ est le demi-cercle fermé supérieur de rayon R > 1 et de centre 0.
Γ 1+z
Z+∞ +∞
Z
cos(x) sin(x)
En déduire la valeur des intégrales 2
dx et dx.
1+x 1 + x2
−∞ −∞

Exercice 5
sin(πz 2 ) + cos(πz 2 ) e2z
I I
Evaluer dz et dz si C est le cercle |z| = 3.
C (z − 1)(z − 2) C (z − 1)4

Exercice 6
On considère la fonction f (z) = (z 2 + 1)1/2 .
a) Montrer que z = ±i sont des points de branchement de f (z).
b) Montrer qu’un tour complet autour des deux points de branchement à la fois n’entraîne
pas de changement de branche de f (z).
c) Déterminer des coupures rendant uniforme la fonction f (z).

Exercice 7
+∞
xp−1
Z
π
Montrer que dx = pour 0 < p < 1.
1+x sin(pπ)
0

6
 1ère année E.N.S.G
Travaux Dirigés en Mathématiques
Fonction d’une variable complexe
Intégrale de Cauchy

Exercice 1
La fonction de variable complexe f (z) = u(x, y) + iv(x, y) est dérivable si l’on peut écrire
df = ldz avec l = α + iβ ∈ C
I ce qui implique les équations de Cauchy-Riemann :
∂u ∂v
= =α
∂x ∂y
∂v ∂u
=− =β
∂x ∂y
On a ici :
f (z) = ez = ex+iy = ex (cos(y) + i sin(y))
d’où :
u(x, y) = ex cos(y) et v(x, y) = ex sin(y)
On a :
∂u ∂v ∂v ∂u
= ex cos(y) = et = ex sin(y) = −
∂x ∂y ∂x ∂y
La dérivée s’écrit alors :

∂f (z) ∂u ∂v ∂u ∂v
= +i = −i + = ex cos(y) + iex sin(y) = ez
dz ∂x ∂x ∂y ∂y
Le calcul peut être reconduit de la même manière pour f (z) = eaz .

Exercice 2
La fonction possède quatre pôles z1 = 0, pôle simple et z2 = 1, pôle triple.
a) Le cercle C1 de centre z0 = −1 et de rayon R < 1 délimite un domaine D dans lequel la
ez
fonction est holomorphe quelque soit z ∈ D, d’où :
z(1 − z)3
ez
Z
I= 3
dz = 0
C1 z(1 − z)

7
 ez
b) Si le rayon R du cercle C2 est tel que 1 ≤ R < 2, la fonction est holomorphe
(1 − z)3
∀z ∈ D et le point 0 est une singularité.
Z
g(z)
D’après Cauchy, on a : I = dz = 2iπg(0) = 2iπ
C2 z

c) Si le rayon R du cercle C3 est tel que R ≥ 2, on peut écrire :

ez ez ez
Z Z Z
I= 3
dz = 3
dz + 3
dz
C3 z(1 − z) γ0 z(1 − z) γ1 z(1 − z)
| {z }
2iπ

z
Pour l’intégrale le long de γ1 , la fonction g(z) = ez est holomorphe ∀z ∈ D. On a la
relation : Z
(2) 2! g(t)
g (z) = dt
2iπ γ1 (t − z)3
d’où :
et
Z
2iπ (2)
3
dt = g (1)
γ1 t(t − 1) 2
et
Z
On en déduit : 3
dt = −iπg (2) (1)
γ1 t(1 − t)
On a :
ez 0 ez
 z
ez ez 2ez

(2) 0 0 e
g (z) = (g (z)) = − 2 = − 2− 2+ 3
z z z z z z
d’où g (2) (1) = e
et
Z
et 3
dt = −iπe
γ1 t(1 − t)
ez ez ez
Z Z Z
d’où : I = 3
dz = 3
dz + 3
dz = 2iπ − iπe = iπ(2 − e)
C z(1 − z) γ0 z(1 − z) γ1 z(1 − z)

8
Exercice 3
eiz
La fonction possède deux pôles, z1 = i et z2 = −i. L’intégrale peut être décomposée
1 + z2
en la somme de deux intégrales sur les contours γ1 et γ2 entourant respectivement les deux pôles.

eiz eiz
Z Z
I= 2
dz + 2
dz
γ1 1 + z γ2 1 + z
 iz   iz 
Z e Z e
z+i z−i
I= dz + dz
γ1 z − i γ2 z + i
eiz eiz
I = 2iπ + 2iπ
z=iz+i z−i z=−i
 
1
πe−1 − πe = π −e
e

Exercice 4
eiz
Z
Calculer 2
dz si Γ est le demi-cercle fermé supérieur de rayon R > 1 et de centre 0.
Γ 1+z
Z+∞ +∞
Z
cos(x) sin(x)
En déduire la valeur des intégrales dx et dx.
1 + x2 1 + x2
−∞ −∞
eiz
La fonction dz possède deux pôles z = ±i dont seul z = i se situe à l’intérieur du
1 + z2
domaine d’intégration. D’après Cauchy, on a donc :

eiz eiz
Z
π
I= 2
dz = 2iπ =
Γ 1+z z + i z=i e

On a par ailleurs, avec γ le demi-cercle supérieur de centre 0 et de rayon supérieur à 1 :

ZR
eiz eix eiz
Z Z
I= dz = dz + dz
Γ 1 + z2 1 + x2 γ+ 1 + z2
−R

eiz eiz
Z
dz ≤ sup πR
γ+ 1 + z2 z∈γ + 1 + z
2

En posant z ∈ γ + ⇒ z = Reiθ = R(cos(θ) + i sin(θ)) avec 0 ≤ θ ≤ π, on a :

sup eiR(cos(θ)+i sin(θ)) sup eiR cos(θ) e−R sin(θ))

eiz
Z
0≤θ≤π 0≤θ≤π
dz ≤ πR = πR
γ+ 1 + z2 inf |1 + R2 e2iθ | inf |1 + R2 e2iθ |
0≤θ≤π 0≤θ≤π

sup e−R sin(θ))

eiz
Z
0≤θ≤π
dz ≤ πR
γ+ 1 + z2 inf |1 + R2 e2iθ |
0≤θ≤π

or sup e−R sin(θ)) = 1 (lorsque θ = 0 ou θ = π) et inf |1 + R2 e2iθ | se produit lorsque

0≤θ≤π 0≤θ≤π
π
R e2iθ est
2 le plus négatif possible c’est-à-dire en θ = 2. Dans ce cas R2 e2iθ = R2 eiπ = −R2 et
|1 − R2 | = R2 − 1 car R ≥ 1 et 1 − R2 est négatif.

9
 On obtient donc :
eiz
Z
πR
2
dz ≤ 2
γ+ 1+z R −1
et
eiz
Z
πR
lim dz ∼ lim 2 =0
R→∞ γ+ 1 + z2 R→∞ R − 1

Comme on a :
ZR
eiz eix eiz
Z Z
π
I= dz = dz + dz =
Γ 1 + z2 1 + x2 γ+ 1+z 2 e
−R

on en déduit :
Z∞
eix eiz
Z
π
dx + lim dz =
1 + x2 R→∞ γ + 1 + z 2 e
−∞ | {z }
=0
Z∞ Z∞
cos(x) sin(x) π
dx + i dx =
1 + x2 1 + x2 e
−∞ −∞
En en déduit :
Z∞ Z∞
cos(x) π sin(x)
2
dx = et dx = 0
1+x e 1 + x2
−∞ −∞

Exercice 5
sin(πz 2 ) + cos(πz 2 )
La fonction possède deux pôles simples en z = 1 et z = 2 ; ces deux pôles
(z − 1)(z − 2)
sont à l’intérieur du cercle de rayon 3. On a donc :
sin(πz 2 ) + cos(πz 2 )
I
I= dz
C (z − 1)(z − 2)
sin(πz 2 ) + cos(πz 2 ) sin(πz 2 ) + cos(πz 2 )
       
= 2iπ Rés , 1 + Rés ,2
(z − 2) (z − 1)
sin(πz 2 ) + cos(πz 2 ) sin(πz 2 ) + cos(πz 2 )
    
= 2iπ lim + lim
z→1 (z − 2) z→2 (z − 1)
= 2iπ(1 + 1)
d’où :
sin(πz 2 ) + cos(πz 2 )
I
dz = 4iπ
C (z − 1)(z − 2)
e2z
La fonction possède un pôle z = 1 d’ordre 4. Si C est le cercle |z| = 3 ce pôle est à
(z − 1)4
l’intérieur du contour d’intégration. D’après Cauchy :
Z
(n) n! f (z)
f (a) = dz
2iπ C (z − a)n+1
avec ici f (z) = e2z , a = 1 et n = 3. On a donc :
e2z
Z
2iπ (3)
4
dz = f (1)
C (z − 1) 3!

10
 On a :

f 0 (z) = 2e2z f 00 (z) = 4e2z f (3) (z) = 8e2z d’où f (3) (1) = 8e2

e2z 2iπ 2 8iπe2

Z
dz = 8e =
C (z − 1)4 3! 3

Exercice 7
+∞
xp−1 z p−1
Z Z
Pour calculer dx, on intègre la fonction de variable complexe dz sur le
1+x C 1+z
0
contour suivant :
avec Γ le cercle centré à l’origine, de rayon R et γ le cercle centré à l’origine de rayon ε.
Calculons la valeur de l’intégrale à l’aide de l’intégrale de Cauchy, on a :

z p−1
Z Z
f (z)
dz = dz = 2iπf (−1) = 2iπ(−1)p−1 = 2iπeiπ(p−1)
C 1+z C 1+z

Décomposons le circuit d’intégration en chemins élémentaires ; on a :

z p−1 z p−1 z p−1 z p−1 z p−1

Z Z Z Z Z
dz = dz + dz + dz + dz
C 1+z AB 1 + z Γ+ 1 + z CD 1 + z γ− 1 + z

Evaluons l’intégrale le long de Γ :

z p−1 z p−1
Z
dz ≤ sup 2πR
Γ+ 1 + z z∈Γ 1 + z

En posant z ∈ Γ ⇒ z = Reiθ avec 0 ≤ θ ≤ 2π, on a :

sup Rp−1 eiθ(p−1)

z p−1 Rp−1 eiθ(p−1)
Z
0≤θ≤2π
dz ≤ sup 2πR ≤ 2πR
Γ+ 1+z 0≤θ≤2π 1 + Reiθ inf |1 + Reiθ |
0≤θ≤2π

z p−1 Rp−1
Z
dz ≤ 2πR ∼ Rp−1 quand R tend vers l’infini, avec 0 < p < 1
Γ+ 1+z R−1
z p−1
Z
1
lim dz ∼ lim 1−p = 0
R→∞ Γ+ 1 + z R→∞ R

Evaluons maintenant l’intégrale le long de γ :

z p−1 z p−1
Z
dz ≤ sup 2πε
γ− 1 + z z∈γ 1 + z

En posant z ∈ γ ⇒ z = εeiθ avec 0 ≤ θ ≤ 2π, on a :

sup εp−1 eiθ(p−1)

z p−1 εp−1 eiθ(p−1)
Z
0≤θ≤2π
dz ≤ sup 2πε ≤ 2πε
γ− 1+z 0≤θ≤2π 1 + εeiθ inf |1 + εeiθ |
0≤θ≤2π

z p−1 εp−1
Z
dz ≤ 2πε ∼ εp quand ε tend vers 0, avec 0 < p < 1
γ− 1+z ε−1

11
 z p−1
Z
lim dz ∼ lim εp = 0
ε→0 Γ+ 1+z ε→0

Sur le segment AB, on choisit z = ρeiθ avec θ ≡ 0 et ε ≤ ρ ≤ R ; on a dz = dρ, d’où :

ZR
z p−1 ρp−1
Z
dz = dρ
AB 1+z 1+ρ
ε

Compte tenu du choix précédent, sur le segment CD, on a z = ρe2iπ , avec ε ≤ ρ ≤ R ; on a

dz = dρ d’où :
Zε
p−1 ZR
p−1
z p−1 ρe2iπ ρe2iπ
Z
dz = dρ = − dρ
CD 1 + z 1 + ρe2iπ 1 + ρe2iπ
R ε

En regroupant ces deux intégrales :

ZR ZR
p−1
z p−1 z p−1 ρp−1 ρe2iπ
Z Z
dz + dz = dρ − dρ
AB 1+z CD 1+z 1+ρ 1 + ρe2iπ
ε ε

ZR
z p−1 z p−1 ρp−1
Z Z  
dz + dz = 1 − e2iπ(p−1) dρ = 2iπeiπ(p−1)
AB 1+z CD 1+z 1+ρ
ε

d’où
Z∞
ρp−1 2iπeiπ(p−1) 2iπeiπ(p−1)
dρ = =
1+ρ 1 − e2iπ(p−1) eiπ(p−1) (e−iπ(p−1) − eiπ(p−1) )
0
Z∞
ρp−1 2iπ −π
dρ = −iπ(p−1) iπ(p−1)
=
1+ρ (e −e ) sin((p − 1)π)
0
Z∞
ρp−1 π
dρ =
1+ρ sin(pπ)
0

12
 1ère année E.N.S.G
Travaux Dirigés en Mathématiques
Fonction d’une variable complexe
Séries de Laurent – Théorème des résidus

Exercice 1
Z∞
sin(x) π eiz
Montrer que dx = . La valeur de cette intégrale sera déduite de l’intégrale de
x 2 z
0
sur un contour que l’on précisera.

Exercice 2
√
z ln(z)
On considère la fonction complexe f (z) = . Déterminer la coupure rendant cette
(1 + z)2 Z √
z ln(z)
fonction univalente. Définir un contour C ne franchissant pas cette coupure et calculer 2
dz
C (1 + z)
Z∞ √
x ln(x)
de façon à en déduire la valeur de dx.
(1 + x)2
0

Exercice 3
 
z − ai
Z
dz 1 z  1
Montrer que 2 2
= arctan + c1 = ln + c2
z +a a a 2ai z + ai

Exercice 4
ez
Z
Calculer dz où C est le cercle défini par |z| = 4.
C (z 2 + π 2 )2
a) par l’intégrale de Cauchy b) par le théorème des résidus

Exercice 5
Développer en séries de Laurent les fonctions suivantes :
e2z
 
1
a) en z = 1 b) (z − 3) sin en z = −2
(z − 1)3 z+2
z − sin(z) z
c) en z = 0 d) en z = −2
z3 (z + 1)(z + 2)

Exercice 6
Z ∞
dx
Evaluer .
0 x6 + 1

Exercice 7
Z 2π Z 2π
cos(3θ) π dθ 5π
Montrer que I1 = dθ = et I2 = 2
=
0 5 − 4 cos(θ) 12 0 (5 − 3 sin(θ)) 32

13
 1ère année E.N.S.G
Travaux Dirigés en Mathématiques
Fonction d’une variable complexe
Séries de Laurent – Théorème des résidus

Exercice 1
Z∞
sin(x) eiz
Pour calculer dx, on intègre la fonction de variable complexe sur le contour
x z
0
suivant :

avec Γ le demi-cercle centré à l’origine, de rayon R et γ le demi-cercle centré à l’origine de rayon

ε. En utilisant l’intégrale de Cauchy et en remarquant qu’il n’y a pas de singularités à l’intérieur
du contour d’intégration, on a :

eiz eiz eiz eiz

Z Z Z Z
I= dz + dz + dz + dz = 0
AB z Γ+ z CD z γ− z

Sur le segment AB, on choisit z = ρeiθ avec θ ≡ 0 et ε ≤ ρ ≤ R ; on a dz = dρ, d’où :

ZR
eiz eiρ
Z
dz = dρ
AB z ρ
ε

Compte tenu du choix précédent, sur le segment CD, on a z = ρeiπ = −ρ, avec ε ≤ ρ ≤ R ;
on a dz = −dρ d’où :
Zε Zε ZR
eiz e−iρ e−iρ e−iρ
Z
dz = (−dρ) = dρ = − dρ
CD z −ρ ρ ρ
R R ε

En regroupant ces deux intégrales :

ZR ZR ZR ZR
eiz eiz eiρ e−iρ eiρ − e−iρ
Z Z
sin(ρ)
dz + dz = dρ − dρ = dρ = 2i dρ
AB z CD z ρ ρ ρ ρ
ε ε ε ε

Evaluons maintenant l’intégrale sur Γ (première tentative avec le théorème de l’inégalité

fondamentale) :
eiz eiz
Z
dz ≤ sup πR
Γ+ z z∈Γ z

14
 En posant z ∈ Γ ⇒ z = Reiθ = R(cos(θ) + i sin(θ)) avec 0 ≤ θ ≤ π, on a :

sup eiR(cos(θ)+i sin(θ))

eiz
Z
0≤θ≤π
dz ≤ πR = sup eiR cos(θ) e−R sin(θ)) π
Γ+ z R 0≤θ≤π

eiz
Z
dz ≤ sup e−R sin(θ)) π
Γ+ z 0≤θ≤π
R 
eiz
or sup e−R sin(θ)) = 1 (lorsque θ = 0 ou θ = π) et l’on a lim Γ+ z dz ≤ π.
0≤θ≤π R→∞

Ce résultat est trop pauvre pour conclure, la majoration est trop importante. Calculons
explicitement l’intégrale (seconde tentative) ; on a dz = iReiθ dθ, d’où :
π π
eiz eiR(cos(θ)+i sin(θ))
Z Z Z
dz = iθ
iReiθ dθ = i eiR(cos(θ)+i sin(θ)) dθ = I1
Γ+ z 0 Re 0
Z π Z π
|I1 | ≤ eiR(cos(θ) e−R sin(θ) dθ = e−R sin(θ) dθ = I2
0 0
Or, pour 0 ≤ θ ≤ π/2, on a sin(θ) = sin(π − θ). On peut donc écrire :
Z π/2 Z π
−R sin(θ)
I2 = e dθ + e−R sin(θ) dθ
0 π/2

En posant, dans la seconde intégrale, u = π − θ, du = −dθ et avec u = π/2 quand θ = π/2

et u = 0 quand θ = π, on obtient :
Z π/2 Z 0 Z π/2
I2 = e−R sin(θ) dθ + e−R sin(u) (−du) = 2 e−R sin(θ) dθ
0 π/2 0

Sur l’intervalle [0, π/2], la fonction que l’on intègre est monotone. L’intégrale I2 est un ma-
jorant de |I1 | ; cherchons un majorant de I2 ; on a :

2θ
sin(θ) ≥ pour 0 ≤ θ ≤ π/2
π

Zπ/2  π h iπ/2 π
I2 ≤ 2 e−2Rθ/π dθ = 2 − e−2Rθ/π = (1 − e−R )
2R 0 R
0

15
 On a donc : π 
lim I1 ≤ lim I2 ≤ lim (1 − e−R ) = 0
R→∞ R→∞ R→∞ R
d’où :
eiz
Z 
lim dz =0
R→∞ Γ+ z
Evaluons maintenant l’intégrale sur γ (première tentative : théorème de l’inégalité fondamen-
tale) ; compte tenu du calcul précédent, on s’aperçoit que l’on peut directement écrire :

eiz
Z 
lim dz ≤ π
ε→∞ γ− z

Là encore, ce résultat est trop pauvre pour pouvoir conclure. Explicitons le calcul :

Z0 Z0
eiz eiε(cos(θ)+i sin(θ)) iθ
Z
dz = iεe dθ = i eiε cos(θ) e−ε sin(θ) dθ
γ− z εeiθ
π π

On cherche à évaluer :
 0 
eiz
Z  Z
lim dz = i lim  eiε cos(θ) e−ε sin(θ) dθ
ε→0 γ− z ε→0
π

Comme la fonction que l’on intègre est bornée, l’intégrale est convergente et l’on peut per-
muter limite et intégrale :
 
0 Z0
eiz
Z  Z   
iε cos(θ) −ε sin(θ)

lim dz = i  lim e e  = i dθ = −iπ
dθ
ε→0 γ− z
 ε→0
π | {z } π
=1

En rassemblant les différents résultats, on obtient :

ZR
sin(ρ)
lim I = 0 = 2i lim dρ − iπ
R→∞ R→∞ ρ
ε→0 ε→0 ε

d’où :
Z∞
sin(ρ) π
dρ =
ρ 2
0

Exercice 2
√
z ln(z) p
La fonction de variable complexe f (z) = 2
est multivalente parce que (z) et ln(z)
(1 + z)
le sont. Il faut donc définir une coupure du plan complexe permettant de rendre univalentes ces
fonctions tout en permettant l’intégration sur le demi-axe réel positif (intégrale de 0 à ∞). On
considère donc le plan complexe coupé par ce demi-axe, origine comprise, et on définit le contour
C d’intégration dans le plan complexe suivant :
Examinons alors la valeur des√ fonctions sur les différentes parties constituant le contour d’in-
z ln(z)
tégration. La fonction f (z) = possède un point singulier en z = 1. Pour z ∈ AB, on
(1 + z)2

16
choisit de poser z = xeiθ avec θ = 0, soit z = x, dans ce cas, on a, pour z ∈ DC, z = xe2iπ .
√ √ √
En√ce qui concerne la fonction z, on a, pour z ∈ AB, z = x et pour z ∈ DC,
√ √ √
z = xe2iπ = xeiπ = − x.
√ √ √
Pour le point particulier z = −1, on a z= −1 = eiπ = eiπ/2 = i.

En ce qui concerne la fonction ln(z), on a, pour z ∈ AB, ln(z) = ln(|z|) + i arg z = ln(x) et
pour z ∈ DC, ln(z) = ln(|z|) + i arg z = ln(x) + 2iπ.

Pour le point particulier z = −1, on a ln(−1) = ln(eiπ ) = ln(1) + iπ = iπ.

z ∈ AB z ∈ DC z = −1
√ √ √
z x − x i
ln(z) ln(x) ln(x) + 2iπ iπ

On peut donc écrire :

Z √ ZR √ Z √ Zε √ Z √
z ln(z) x ln(x) z ln(z) x(ln(x) + 2iπ) z ln(z)
dz = dx + dz − dx + dz
C (1 + z)2 (1 + x)2 Γ+ (1 + z)2 (1 + x)2 γ− (1 + z)2
ε R

Le théorème des résidus permet de calculer cette intégrale. La seule singularité à l’intérieur
du contour défini par C est le point z = −1 ; on a donc :
Z √
z ln(z) √
2
dz = 2iπRés( z ln(z), −1)
C (1 + z)

Le point z = −1 est un pôle double, on a donc :

√ 
√ √


0 ln(z) z
Rés( z ln(z), −1) = lim z ln(z) = lim √ +
z→−1 z→−1 2 z z
√ ln(−1) √
 
π
Rés( z ln(z), −1) = √ − −1 = − i
2 −1 2
On obtient donc : Z √
z ln(z) π 
dz = 2iπ − i = iπ 2 + 2π
C (1 + z)2 2
De plus, on a, d’après le théorème de l’inégalité fondamentale :
√ √ q 2
Z √ sup| z||ln(z)| ε ln (ε) + 4π 2
z ln(z) z∈Γ
2
dz ≤ 2πR = 2πε
Γ+ (1 + z) inf |(1 + z)2 | (R − 1)2
z∈Γ

√ √ q 2
Z
z ln(z) R ln (R)
ln(R)
2
dz ∼ R 2
∼ √
Γ+ (1 + z) R R
Z √
z ln(z) ln(R)
lim 2
dz = lim √ =0
R→∞ Γ+ (1 + z) R→∞ R
de plus : √
√ sup| z||ln(z)| √ q 2
Z
z ln(z) z∈γ ε ln (ε) + 4π 2
dz ≤ 2πε = 2πε
γ− (1 + z)2 inf |(1 + z)2 | (ε − 1)2
z∈γ

17
 √
√
Z
z ln(z)
2
dz ∼ −ε ε ln(ε)
γ− (1 + z)
q q
car ln2 (ε) + 4π 2 ∼ ln2 (ε) = − ln(ε) et (ε − 1)2 ∼ 1.
√
√
Z
z ln(z)
lim 2
dz = lim −ε ε ln(ε) = 0
ε→0 γ− (1 + z) ε→0

On obtient donc :
Z √ Zε √
 R 
x ln(x) x(ln(x) + 2iπ) 
lim  dx − dx = iπ 2 + 2π
R→∞ (1 + x)2 (1 + x)2
ε→0 ε R

Z √ ZR √
 R 
x ln(x) x(ln(x) + 2iπ) 
lim  2
dx + dx = iπ 2 + 2π
R→∞ (1 + x) (1 + x)2
ε→0 ε ε

ZR √ ZR √
 
x ln(x) x
lim 2 dx + 2iπ dx = 2π + iπ 2
R→∞ (1 + x)2 (1 + x)2
ε→0 ε ε

d’où :
Z∞ √ Z∞ √
x ln(x) x π
dx = π on a aussi dx =
(1 + x)2 (1 + x)2 2
0 0

On aurait pu effectuer un autre choix :

z = xe2iπ ∈ AB z = xe4iπ ∈ DC z = e3iπ = −1

√ √ √
z − x x −i
ln(z) ln(x) + 2iπ ln(x) + 4iπ 3iπ

On a toujours :
√ 
√ √


0 ln(z) z
Rés( z ln(z), −1) = lim z ln(z) = lim √ + avec z = e3iπ = −1
z→−1 z→−1 2 z z
 √ 
ln(z) z 3iπ −i 3π
√ + = + =− +i
2 z z z=e3iπ 2(−i) −1 2
d’où : Z √  
z ln(z) 3π
2
dz = 2iπ − + i = −3iπ 2 − 2π
C (1 + z) 2

Z √ Zε √
 R 
x ln(x) + 2iπ x(ln(x) + 4iπ) 
lim − dx + dx = −3iπ 2 − 2π
R→∞ (1 + x)2 (1 + x)2
ε→0 ε R

ZR √ ZR √
 
x ln(x) x
lim −2 2
dx − 6iπ dx = −3iπ 2 − 2π
R→∞ (1 + x) (1 + x)2
ε→0 ε ε

d’où les résultats obtenus précédemment.

18
Exercice 3
Comme il s’agit ici d’une intégrale indéfinie, il faut trouver une primitive.

Première méthode (en se ramenant à une forme connue) :

Z Z Z
dz dz 1 1/a
2 2
= z 2 = 2 dz
z +a a2 a2 + 1 a 1 + az 2
du
La fonction à intégrer est du type 1+u2
dont la primitive est arctan(u). On obtient donc :
Z
dz 1 z 
= arctan + c1
z 2 + a2 a a

Seconde méthode (par décomposition en éléments simples) :

 
1 1 1 1 1
f (z) = 2 = = −
z + a2 (z − ia)(z + ia) 2ai z − ai z + ai
Z Z Z 
1 1 1
f (z) = −
2ai z − ai z + ai
Z
1
f (z) = (ln(z − ai) − ln(z + ai)) + c2
2ai
 
z − ai
Z
dz 1
= ln + c2
z 2 + a2 2ai z + ai

Exercice 6
1
La fonction f (z) = est paire. On souhaite évaluer l’intégrale dur l’axe réel positif. On
z6 + 1
choisit le contour d’intégration ci-dessous :

π kπ
Cherchons les pôles de la fonction. z 6 + 1 = 0 ⇔ z 6 = −1 d’où z = ei( 6 + 3
)
, k = [0, 5]. Les
pôles à l’intérieur du domaine d’intégration sont ceux pour k = 1, 2, 3.

On a donc : I
dz π π 5π
I= = 2iπ(Rés(ei 6 ) + Rés(ei 2 ) + Rés(ei 6 ))
C z6+1
1 g(z) g(a)
Chaque pôle est simple. Avec f (z) = = , on a Rés(a) = 0 , h0 (z) = 6z 5 .
z6 +1 h(z) h (a)

π 1 1 5π π 1 1 π 5π 1 1 π
Rés(ei 6 ) = = e−i 6 Rés(ei 2 ) = = e−i 2 Rés(ei 6 )= = e−i 6
i 5π 6 i 5π 6 i 25π 6
6e 6 6e 2 6e 6

19
 √ √
iπ 3 i 3 i 2π
I = (− − |{z}
−i + − )=
3 | 2{z 2} π | 2 {z 2} 3
−i
5π e 2 π
ei 6 e−i 6

ZR Z
dx dz
I1 = 6
+
z +1 γ z6 +1
−R

En posant z ∈ γ ⇒ z = Reiθ = R(cos(θ) + i sin(θ)) avec 0 ≤ θ ≤ π, on a :

Z
dz πR πR
6
≤ 6
≤ 6
γ z +1 inf |z + 1| R −1
z∈γ
Z
dz πR
lim = lim 6 =0
R→∞ γ z6 +1 R→∞ R − 1

Or la fonction f (x) est paire donc on a :

ZR ZR
dx dx 2π
lim 6
= lim 2 =
R→∞ x + 1 R→∞ x6 +1 3
−R 0

d’où :
Z∞
dx π
=
x6 + 1 3
0

Exercice 7
Z 2π
cos(3θ)
I1 = dθ
0 5 − 4 cos(θ)
z + z̄ z 3 + z̄ 3 1
On pose z = eiθ , c’est-à-dire cos(θ) = et cos(3θ) = avec z̄ = et dz = ieiθ dθ =
2 2 z
izdθ.

1
z 3 + z̄ 3 z3 + z6 + 1
Z Z Z
dz dz
z3 dz
I1 = = 1 = 2 2
2(5 − 2(z + z̄)) iz 2(5 − 2(z + z )) iz 3
2z (5 − z (z + 1)) iz
C(0,1) C(0,1) C(0,1)

z6 + 1 z6 + 1
Z Z
1 1
I1 = 3 2
dz = dz
2i z (5z − 2z − 2) 2i z 3 (z − 2)(−2z + 1)
C(0,1) C(0,1)

z6 + 1
Z
1
I1 = − dz
4i z 3 (z − 2)(z − 21 )
C(0,1)

On a donc z = 0, pôle triple, z = 2 et z = 1/2, pôles simples. Seuls les pôles z = 0 et z = 1/2
sont à l’intérieur du cercle unité. On a donc :
 
1
I1 = 2iπ − (Rés(0) + Rés(1/2))
4i

20
 1 6 1 65


2 +1 64 +1 65
Rés(1/2) = = = − 64
3 =−
1 3 1 − 18 23 12



2 2 −2 16
!00
1 z6 + 1 1
Rés(0) = lim z 3 3 = lim (f (z))00
z (z − 2) z − 12


2! z→0 2! z→0

6z 5 z6 + 1 z6 + 1
(f (z))0 = − −
(z − 2) z − 12 (z − 2)2 z − 21



2
(z − 2) z − 12
30z 4 6z 5 6z 5
(f (z))00 = − −
(z − 2) z − 12 (z − 2)2 z − 12



2
(z − 2) z − 12
6z 5 z6 + 1 z6 + 1
− − 2 +
(z − 2)2 z − 21 (z − 2)3 z − 12



2
(z − 2)2 z − 21
6z 5 z6 + 1 z6 + 1
−
+
+2
1 2 1 2 1 3


(z − 2) z − 2 (z − 2)2 z − 2 (z − 2) z − 2
2 1 1 2
f 00 (0) =
+ + +
− 12

2
2
3
(−2)3 (−2)2 − 12 (−2)2 − 21 (−2) − 12
1 21
f 00 (0) =
+1+1+8=
2 2
    
1 65 1 21 π −65 + 63 π
I1 = 2iπ − − + =− =
4i 12 2 2 2 12 12
Z 2π
dθ
I2 =
0 (5 − 3 sin(θ))2
z − z̄ 1
On pose z = eiθ , c’est-à-dire sin(θ) = avec z̄ = et dz = ieiθ dθ = izdθ.
2i z

(2iz)2 dz
Z Z Z
dz dz
I2 = 3 = 3 =
iz(5 − 2i (z − z̄))2 iz(5 − 2iz (z 2 − 1))2 iz(10iz − 3(z 2 − 1))2
C(0,1) C(0,1) C(0,1)
Z Z Z
4 zdz 4 zdz 4 zdz
I2 = − 2 2
=− 2 2
=−
i (10iz − 3(z − 1)) i (3z − 10iz − 3) i (z − 3i)2 (3z − i)2
C(0,1) C(0,1) C(0,1)
Z
4 zdz
I2 = −
9i (z − 3i)2 (z − 3i )2
C(0,1)
i
Seul le pôle z = 3 est à l’intérieur du cercle unité.

   2 !0  0
i i z z
Rés = lim z− = lim = lim (f (z))0
3 z→ 3i 3 (z − 3i)2 (z − 3i )2 z→ 3i (z − 3i)2 z→ 3i

1 2z
(f (z))0 = −
(z − 3i)2 (z − 3i)3

21
 2 3i
 
i 1 9 (2i)27
f0 = i − = −
3 ( 3 − 3i)2 ( 3i − 3i)3 (−8i)2 3(−8i)3
 
0 i 9 (2i)27 9 9 45
f =− − =− − =−
3 64 3.512i 64 256 256
   
4 i 8π 45 5π
I2 = − 2iπRés =− − =
9i 3 9 256 32

22
 1ère année E.N.S.G
Travaux Dirigés en Mathématiques
Equations aux dérivées partielles

Exercice 1
Résoudre le problème suivant à l’aide de la méthode de séparation des variables.

∂u(x, y) ∂u(x, y)
=4 avec u(x = 0, y) = 8e−3y
∂x ∂y

Exercice 2
Résoudre le problème suivant à l’aide de la méthode de séparation des variables et par appli-
cation du principe de superposition des solutions.

∂u(x, y) ∂u(x, y)
=4 avec u(x = 0, y) = 8e−3y + 4e−5y
∂x ∂y

Exercice 3
Résoudre le problème suivant (équation hyperbolique avec conditions initiales et aux limites)
à l’aide de la méthode de séparation des variables et de la méthode de superposition des solutions.

∂ 2 y(x, t) ∂ 2 y(x, t)
= 16 0 < x2, t > 0
∂t2 ∂x2
Conditions aux limites : y(x = 0, t > 0) = 0
y(x = 2, t > 0) = 0

Conditions initiales : y(x, t = 0+ ) = 6 sin(πx) − 3 sin(4πx)

∂y
(x, t = 0+ ) = 0
∂t

Exercice 4
Résoudre le problème suivant (équation parabolique avec conditions initiales et aux limites) à
l’aide de la méthode de séparation des variables et de la méthode de superposition des solutions.

∂p(x, t) ∂ 2 p(x, t)
=
∂t ∂x2
∂
Conditions aux limites : (p(x = 0, t > 0) = 0
∂x
p(x = 1, t > 0) = 0

Condition initiale : p(x, t = 0+ ) = 1 (0 ≤ x ≤ 1)

23
Exercice 5
Résoudre le problème suivant (équation parabolique avec conditions initiales et aux limites) à
l’aide de la méthode de séparation des variables et de la méthode de superposition des solutions.
∂p(x, t) ∂ 2 p(x, t)
=
∂t ∂x2
Conditions aux limites : p(x = 2, t > 0) = 0

p(x = 0, t > 0) = 0

Condition initiale : p(x, t = 0+ ) = 1 (0 ≤ x ≤ 2)

Exercice 6
Résoudre le problème suivant (équation parabolique avec conditions initiales et aux limites) à
l’aide de la méthode de séparation des variables et de la méthode de superposition des solutions.
∂p(x, t) ∂ 2 p(x, t)
=
∂t ∂x2
Conditions aux limites : p(x = 2, t > 0) = 0

p(x = 0, t > 0) = 0

Condition initiale : p(x, t = 0+ ) = kx (0 ≤ x ≤ 2)

Exercice 7
Résoudre le problème suivant (équation hyperbolique avec conditions initiales et aux limites)
à l’aide de la méthode de séparation des variables, de la méthode de superposition des solutions
et de la méthode de décomposition en séries de Fourier.

Une corde de longueur L est tendue entre les points (0, 0) et (L, 0) sur l’axe des x. Au temps
t = 0+ , sa forme est donnée par la fonction f (x), 0 < x < L et elle est libérée. Trouver le
déplacement de la corde à tout instant ultérieur.

L’équation de la corde vibrante est :

∂2y 2
2∂ y
= a 0 < x < L, t > 0
∂t2 ∂x2
où y(x, t) désigne l’écartement par rapport à l’axe des x au temps t.

Conditions aux limites (extrémités de la corde fixées en x = 0 et x = L) :

y(x = 0, t > 0) = 0
y(x = L, t > 0) = 0

Conditions initiales :

Forme initiale : y(x, t = 0+ ) = f (x)

∂y
Vitesse initiale : yt (x, t = 0+ ) = (x, t = 0+ ) = 0
∂t

24
Exercice 8
Résoudre le problème suivant (équation elliptique avec conditions initiales et aux limites) à
l’aide de la méthode de séparation des variables, de la méthode de superposition des solutions et
de la méthode de décomposition en séries de Fourier.

Soit une plaque carrée dont les côtés ont la longueur unité. Trois faces de cette plaque sont
isolées et ses côtés maintenus à 0◦ C tandis que le quatrième côté est maintenu à la température
θ1 . Déterminer la température d’état stationnaire(ou de régime permanent) en tout point de la
plaque.

Les conditions aux limites du problème sont les suivantes :

θ(x = 0, 0 ≤ y ≤ 1) = θ(x = 1, 0 ≤ y ≤ 1) = 0
θ(0 ≤ x ≤ 1, y = 0) = 0
θ(0 ≤ x ≤ 1, y = 1) = θ1

Le problème thermique transitoire (ou instationnaire) est régi par une équation aux dérivées
partielles de type parabolique :

∂ ∂2 ∂2
θ(x, y, t) = θ(x, y, t) + θ(x, y, t)
∂t ∂x2 ∂y 2
L’état stationnaire ou le régime permanent, est caractérisé par l’indépendance de la tempé-
rature par rapport au temps : ∂θ∂t = 0.

Le problème thermique à résoudre se ramène alors à l’équation de Laplace en deux dimensions

(problème elliptique) :
∂2 ∂2
θ(x, y) + θ(x, y) = 0
∂x2 ∂y 2

25
 1ère année E.N.S.G
Travaux Dirigés en Mathématiques
Equations aux dérivées partielles

Application de la méthode de Fourier

— méthode de la séparation des variables, pour obtenir des solutions à variables séparées ;
— application du principe de superposition pour construire la solution qui vérifiera les condi-
tions imposées (souvent la condition initiale).

Exercice 1
Méthode de séparation des variables ; on pose u(x, y) = f (x)g(y) d’où :

∂f (x) ∂g(y)
f 0 (x)g(y) = 4f (x)g 0 (y) avec f 0 (x) = et g 0 (y) =
dx dy
soit :
f 0 x)g(y) 4f (x)g 0 (y)
=
4f (x)g(y) 4f (x)g(y)
f 0 (x) g 0 (y)
=
4f (x) g(y)
Puisque f ne dépend que de x, g ne dépend que de y et que x et y sont des variables
indépendantes, les deux rapports précédents doivent être égaux à une constante notée λ :

f 0 (x) g 0 (y)
=λ =λ
4f (x) g(y)

dont les solutions sont de la forme :

f (x) = Ae4λx g(y) = Beλy

La solution de l’équation est donc de la forme :

u(x, y) = Ce4λx eλy

Il reste à vérifier la condition initiale :

u(x = 0, y) = 8e−3y

soit :
u(x = 0, y) = Ceλy
d’où :
C=8 et λ = −3
La solution du problème complet est donc :

u(x, y) = 8e−12x−3y

26
Exercice 2
Comme pour l’exercice 1,

u(x, y) = Ce4λx eλy = Ceλ(4x+y)

est solution de l’équation. Deα(4x+y) et Eeβ(4x+y) sont donc également des solutions. Par appli-
cation du principe de superposition, leur somme est également une solution :

u(x, y) = Deα(4x+y) + Eeβ(4x+y)

D’après la condition aux limites :

u(x = 0, y) = 8e−3y + 4e−5y = Deαy + Eeβy

on identifie les constantes D, α, E, β :

D = 8; α = −3; E = 4; β = −5

La solution cherchée est donc :

u(x, y) = 8e−12x−3y − 4e−20x−5y

Exercice 3
Avec y(x, t) = f (x)g(t), l’équation :

∂ 2 y(x, t) ∂ 2 y(x, t)
= 16
∂t2 ∂x2
s’écrit :
∂ 2 g(t) ∂ 2 f (t)
f (x)g 00 (t) = 16f 00 (x)g(t) avec g 00 (t) = et f 00 (t) =
dt2 dx2
soit :
f (x)g 00 (t) 16f 00 (x)g(t)
=
16f (x)g(t) 16f (x)g(t)
ou :
g 00 (t) f 00 (x)
=
16g(t) f (x)
d’où :
g 00 (t) f 00 (x)
= = −λ
16g(t) f (x)
On pose, par commodité d’écriture λ = ω 2 soit :

g 00 (t)
= −ω 2
16g(t)

f 00 (x)
= −ω 2
f (x)
La solution des deux équations précédentes est de la forme :

f (x) = A cos(ωx) + B sin(ωx)

27
 g(t) = C cos(4ωt) + D sin(4ωt)
La solution de l’équation est donc de la forme :

y(x, t) = (A cos(ωx) + B sin(ωx))(C cos(4ωt) + D sin(4ωt))

d’où :
∂y(x, t)
= 4ω(A cos(ωx) + B sin(ωx))(−C sin(4ωt) + D cos(4ωt))
dt
La prise en compte des conditions aux limites :

y(x = 0, t > 0) = 0

y(x = 2, t > 0) = 0
conduit à :
y(x = 0, t > 0) = A(C cos(4ωt) + D sin(4ωt)) = 0
d’où :
A=0
y(x = 2, t > 0) = B sin(2ω)(C cos(4ωt) + D sin(4ωt)) = 0
qui implique :
nπ
sin(2ω) = 0 d’où 2ω = nπ soit ω = , n = 0, 1, 2, . . .
2
Il reste à vérifier les conditions initiales :

y(x, t = 0+ ) = 6 sin(πx) − 3 sin(4πx)

∂y
yt (x, t = 0+ ) =
(x, t = 0+ ) = 0
∂t
On exploite tout d’abord la seconde équation :
∂y
= 4ωB sin(ωx)(−C sin(4ωt) + D cos(4ωt))
∂t
∂y
(x, t = 0+ ) = 4ωB sin(ωx)D cos(4ωt)) = 0
∂t
d’où :
D=0
Il reste donc :
y(x, t) = B sin(ωx)C cos(4ωt)
soit :
y(x, t) = E sin(ωx) cos(4ωt)
ou :  nπx 
y(x, t) = En sin cos(2nπt)
2
On applique le principe de superposition pour vérifier la condition initiale supplémentaire :

y(x, t = 0+ ) = 6 sin(πx) − 3 sin(4πx)

soit :  n πx   n πx 
1 2
y(x, t) = En1 sin cos(2n1 πt) + En2 sin cos(2n2 πt)
2 2

28
  n πx   n πx 
1 2
y(x, t = 0+ ) = En1 sin + En2 sin
2 2
y(x, t = 0+ ) = 6 sin(πx) − 3 sin(4πx)
Par identification, on obtient alors :

En1 = 6; En2 = −3; n1 = 2; n2 = 8

soit la solution recherchée :

y(x, t) = 6 sin(πx) cos(4πt) − 3 sin(4πx) cos(16πt)

Interprétation physique : ce problème peut s’interpréter comme un problème de corde vi-

brante. Les extrémités de la corde sont fixées en x = 0 et x = 2 et la forme initiale donnée de la
corde est :
f (x) = 6 sin(πx) − 3 sin(4πx)
La corde est alors lâchée avec une vitesse initiale nulle. L’équation :

y(x, t) = 6 sin(πx) cos(4πt) − 3 sin(4πx) cos(16πt)

donne alors, dans ces conditions, le déplacement de tout point x à un instant ultérieur.

Programme Matlab

x = 0:0.01:2;
y = 6*sin(pi*x)-3*sin(4*pi*x);
whitebg
plot(x,y)
axis([0 2 -10 10]);
axis(axis)
M = moviein(100);
vidObj = VideoWriter(’corde.mp4’,’MPEG-4’);
open(vidObj);
set(gca,’nextplot’,’replacechildren’);
for j = 0:99
t = j*0.01;
y = 6*sin(pi*x)*cos(4*pi*t)-3*sin(4*pi*x)*cos(16*pi*t);
plot(x,y,’r’,’Linewidth’,2)
M(:,j+1) = getframe;
writeVideo(vidObj,getframe);
end
close(vidObj)
movie(M)
map=colormap;

29
 Figure 1 – Animation de l’évolution de la corde (nécessite Adobe Reader)

Exercice 4
On résout ce problème en deux étapes :
a) On recherche dans un premier temps la forme générale de la solution du problème de
diffusion suivant (équation d’évolution et conditions aux limites) :

∂p(x, t) ∂ 2 p(x, t)
=
∂t ∂x2
∂
p(x = 0, t > 0) = 0
∂x
p(x = 1, t > 0) = 0
en appliquant la méthode de séparation des variables :

p(x, t) = f (x)g(t)

b) A partir de la solution précédente, on construit la solution du problème complet (équation

d’évolution conditions aux limites et condition initiale) sous la forme d’une série des solu-
tions de a) (application du principe de superposition des solutions).

Problème intermédiaire a) :

En utilisant :
∂p(x, t) dg(t) ∂ 2 p(x, t) d2 f (x)
= f (x) et = g(t)
∂t dt ∂x2 dx2

30
l’équation de diffusion prend la forme :

dg(t) d2 f (x)
f (x) = g(t)
dt dx2
soit :
1 d2 f (x) 1 dg(t)
2
= = constante
f (x) dx g(t) dt
On cherche la solution de l’équation précédente sous la forme :

1 dg(t)
= −η 2
g(t) dt
soit :
g(t) = A exp(−η 2 t)
Le choix d’une constante réelle négative −η 2 est guidé par la nécessité de construire une
solution définie quelque soit l’instant t positif (en particulier pour t → ∞).

La solution de l’équation :
1 d2 f (x)
= −η 2
f (x) dx2
est de la forme :
f (x) = B1 cos(ηx) + B2 sin(ηx)
On a donc un solution générale de la forme :

p(x, t) = (B1 cos(ηx) + B2 sin(ηx))A exp(−η 2 t)

donc :
∂p(x, t)
= η(−B1 sin(ηx) + B2 cos(ηx))A exp(−η 2 t)
∂x
En prenant en compte la condition initiale :
∂
(p(x = 0, t > 0) = ηB2 A exp(−η 2 t) = 0
∂x
donc : B2 = 0 soit :
p(x, t) = B1 cos(ηx)A exp(−η 2 t)
La condition imposée en x = 1 s’écrit :

p(x = 1, t > 0) = B1 cos(η)A exp(−η 2 t) = 0

elle implique :
cos(η) = 0
d’où
π
η = (1 + 2n)
2
π2
η 2 = (1 + 2n)2
4
On en déduit la forme de la solution du problème a) :
2
 
 π  2π
p(x, t) = An cos (1 + 2n) x exp −(1 + 2n) t
2 4

31
Construction de la solution du problème complet sous forme d’une série

La famille de solution :
π2
 
 π 
p(x, t) = An cos (1 + 2n) x exp −(1 + 2n)2 t
2 4
ne satisfait pas la condition initiale

p(x, t = 0+ ) = 1 (0 ≤ x ≤ 1)

On va rechercher la solution du problème complet sous forme d’une série :

∞  2

X  π  2π
p(x, t) = An cos (1 + 2n) x exp −(1 + 2n) t
2 4
n=0

on en déduit :
∞
+
X  π 
p(x, t = 0 ) = An cos (1 + 2n) x
2
n=0
La condition initiale impose :
∞
+
X  π 
p(x, t = 0 ) = An cos (1 + 2n) x (0 ≤ x ≤ 1)
2
n=0

Le problème posé est donc de trouver la forme des coefficients An tels que :
∞
X  π 
An cos (1 + 2n) x = 1
2
n=0

Le principe de la méthode suivante est “d’extraire” le terme An de la sommation infinie. On

va utiliser pour cela une propriété dite d’orthogonalité de la base des fonctions cos (1 + 2n) π2 x
vis-à-vis du produit scalaire sur 0, 1.

On multiplie tout d’abord les deux membres de l’égalité :

∞
X  π 
An cos (1 + 2n) x = 1
2
n=0

par le facteur cos 1 + 2n0 ) π2 x , soit :

∞
 π X
0
 π   π 
cos 1 + 2n ) x An cos (1 + 2n) x = cos (1 + 2n0 ) x
2 2 2
n=0

ou encore :
∞
X  π   π   π 
An cos (1 + 2n0 ) x cos (1 + 2n) x = cos (1 + 2n0 ) x
2 2 2
n=0

On intègre ensuite les deux membres de l’égalité précédente entre x = 0 et x = 1 :

Z1 X
∞ Z1
 π   π   π 
An cos (1 + 2n ) x cos (1 + 2n) x dx = cos (1 + 2n0 ) x dx
0
2 2 2
0 n=0 0

32
 On réécrit le premier membre en inversant les opérations d’intégration et de sommation et
en extrayant de l’intégrale le terme constant An :

∞ Z1 Z1
X  π   π   π 
An 0
cos (1 + 2n ) x cos (1 + 2n) x dx = cos (1 + 2n0 ) x dx
2 2 2
n=0 0 0

La propriété d’orthogonalité de la base des fonctions cos 1 + 2n) π2 x va permettre de réduire

la sommation infinie à un seul terme. En effet, en utilisant les résultats suivants :

1
cos(a) cos(b) = (cos(a + b) + cos(a − b))
2
 π   π  1
cos 1 + 2n0 ) x cos 1 + 2n) x = (cos((1 + n + n0 )πx) + cos((n − n0 )πx)))
2 2 2
puis :

Z1 Z1
 π  
0 π  1
I= cos 1 + 2n ) x cos 1 + 2n) x dx = (cos((1 + n + n0 )πx) + cos((n − n0 )πx)))
2 2 2
0 0

Si n 6= n0 :
 
1 1 1 1 1
sin((1 + n + n0 )πx) 0 + sin((n − n0 )πx) 0
 
I= 0 0
2 (1 + n + n )π (n − n )π

I=0
Si n = n0 :
Z1  
1 1 1
I= (cos((1 + 2n)πx + 1)dx = [sin((1 + 2n)πx)]10 + [x]10
2 2 (1 + 2n)π
0
 
1 1 1
I= sin((1 + 2n)π) + 1 =
2 (1 + 2n)π 2
On obtient donc :
∞ Z1
X  π   π 
I1 = An cos (1 + 2n0 ) x cos (1 + 2n) x dx
2 2
n=0 0

Z1  π  0

0 π
 1
I1 = A n0 cos 1 + 2n ) x cos 1 + 2n ) x dx = An0
2 2 2
0

D’autre part :

Z1  π  2 h 
0 π
i1 2 
0 π

cos (1 + 2n0 ) x dx = sin (1 + 2n ) x = sin (1 + 2n )
2 (1 + 2n0 )π 2 0 (1 + 2n0 )π 2
0

Z1 0
 π  2(−1)n
cos (1 + 2n0 ) x dx =
2 (1 + 2n0 )π
0

33
 L’égalité :

∞ Z1 Z1
X  π   π   π 
An 0
cos (1 + 2n ) x cos (1 + 2n) x dx = cos (1 + 2n0 ) x dx
2 2 2
n=0 0 0

se réduit donc à : 0 0
1 2(−1)n 4(−1)n
An0 = d’où An0 =
2 (1 + 2n0 )π (1 + 2n0 )π
On obtient donc la valeur recherchée du coefficient An :
4(−1)n
An =
(1 + 2n)π

On a bien finalement :
∞
X 4(−1)n  π 
p(x, t = 0+ ) = cos (1 + 2n) x = 1 (0 ≤ x ≤ 1)
(1 + 2n)π 2
n=0

La solution complète s’écrit donc :

∞
4(−1)n 2
 
X  π  2π
p(x, t) = cos (1 + 2n) x exp −(1 + 2n) t
(1 + 2n)π 2 4
n=0

Interprétation physique

Le problème étudié représente le problème physique dit de “la consolidation d’une couche de
sol sous surcharge”. On s’intéresse à l’évolution de la pression de l’eau (pression interstitielle)
contenue dans cette couche de sol. La base de la couche de sol (x = 0) repose sur un socle rigide
et imperméable, cette dernière condition se traduit par une condition de débit nul :
∂
(x = 0, t > 0) = 0
∂x
La surface (x = 1) de la couche de sol est drainée, ce qui s’écrit mathématiquement comme
une condition de pression constante (nulle) :

p(x = 1, t > 0) = 0

La condition initiale du calcul :

p(x, t = 0+) = 1 (0 ≤ x ≤ 1)

correspond à la perturbation provoquée par l’application instantanée de la surcharge mécanique

sur la surface.

L’étude de la solution de ce problème transitoire montre que le champ de pression interstitielle

dans la couche de sol tend à retrouver son état d’équilibre (pression nulle dans la couche) et que
la “perturbation”, c’est-à-dire l’écart par rapport à l’équilibre, tend à s’annuler.

34
Exercice 5
On résout ce problème en deux étapes :
a) On recherche dans un premier temps la forme générale de la solution du problème de
diffusion suivant (équation d’évolution et conditions aux limites) :

∂p(x, t) ∂ 2 p(x, t)
=
∂t ∂x2
p(x = 2, t > 0) = 0
p(x = 0, t > 0) = 0
en appliquant la méthode de séparation des variables :

p(x, t) = f (x)g(t)

b) A partir de la solution précédente, on construit la solution du problème complet (équation

d’évolution conditions aux limites et condition initiale) sous la forme d’une série des solu-
tions de a) (application du principe de superposition des solutions).

Problème intermédiaire a) :

En utilisant :
∂p(x, t) dg(t) ∂ 2 p(x, t) d2 f (x)
= f (x) et = g(t)
∂t dt ∂x2 dx2
l’équation de diffusion prend la forme :

dg(t) d2 f (x)
f (x) = g(t)
dt dx2
soit :
1 d2 f (x) 1 dg(t)
2
= = constante
f (x) dx g(t) dt
On cherche la solution de l’équation précédente sous la forme :

1 dg(t)
= −η 2
g(t) dt

soit :
g(t) = A exp(−η 2 t)
Le choix d’une constante réelle négative −η 2 est guidé par la nécessité de construire une
solution définie quelque soit l’instant t positif (en particulier pour t → ∞).

La solution de l’équation :
1 d2 f (x)
= −η 2
f (x) dx2
est de la forme :
f (x) = B1 cos(ηx) + B2 sin(ηx)
On a donc un solution générale de la forme :

p(x, t) = (D cos(ηx) + E sin(ηx)) exp(−η 2 t)

35
 Les conditions aux limites :
p(x = 2, t > 0) = 0
p(x = 0, t > 0) = 0
imposent :
p(x = 0, t > 0) = D exp(−η 2 t) = 0
d’où D = 0 et :
p(x, t) = E sin(ηx)) exp(−η 2 t)
p(x = 2, t > 0) = E sin(2η) exp(−η 2 t)
d’où :
sin(2η) = 0 ⇒ 2η = (1 + n)π
soit :
π
η = (1 + n) n = 0, 1, 2, . . .
2
π2
η 2 = (1 + n)2
4
La solution du problème intermédiaire est donc de la forme :
2
 
 π  2π
p(x, t) = En sin (1 + n) x exp −(1 + n) t
2 4

Construction de la solution du problème complet sous forme d’une série

La famille de solution :
2
 
 π  2π
p(x, t) = En sin (1 + n) x exp −(1 + n) t
2 4

ne satisfait pas la condition initiale

p(x, t = 0+ ) = 1 0 ≤ x ≤ 2

On va rechercher la solution du problème complet sous forme d’une série :

∞  2

X  π  2π
p(x, t) = En sin (1 + n) x exp −(1 + n) t
2 4
n=0

on en déduit :
∞
X  π 
p(x, t = 0+ ) = En sin (1 + n) x
2
n=0

La condition initiale impose :

∞
X  π 
p(x, t = 0+ ) = En sin (1 + n) x = 1 0 ≤ x ≤ 2
2
n=0

Le problème posé est donc de trouver la forme des coefficients En telle que :
∞
X  π 
En sin (1 + n) x = 1
2
n=0

36
 Le principe de la méthode suivante est “d’extraire” le terme En de la sommation infinie. On

va utiliser pour cela une propriété dite d’orthogonalité de la base des fonctions sin (1 + n) π2 x
vis-à-vis du produit scalaire sur 0, 2.

On multiplie tout d’abord les deux membres de l’égalité :

∞
X  π 
En sin (1 + n) x = 1
2
n=0

par le facteur :  π 
sin (1 + n0 ) x
2
soit :
∞
 π X0
 π   π 
sin (1 + n ) x En sin (1 + n) x = sin (1 + n0 ) x
2 2 2
n=0
ou encore :
∞
X  π   π   π 
En sin (1 + n0 ) x sin (1 + n) x = sin (1 + n0 ) x
2 2 2
n=0
On intègre ensuite les deux membres de l’égalité précédente entre x = 0 et x = 2 :
Z2 X
∞ Z2
 π   π   π 
En sin (1 + n0 ) x sin (1 + n) x dx = sin (1 + n0 ) x dx
2 2 2
0 n=0 0

On réécrit le premier membre en inversant les opérations d’intégration et de sommation et

en extrayant de l’intégrale le terme constant En :
∞ Z2 Z2
X  π 
0
 π   π 
En sin (1 + n ) x sin (1 + n) x dx = sin (1 + n0 ) x dx
2 2 2
n=0 0 0

En utilisant le résultat suivant :

1
sin(a) sin(b) = (cos(a − b) − cos(a + b))
2
d’où :
 π   π  1  π   π 
sin (1 + n0 ) x sin (1 + n) x = cos (n − n0 ) x − cos (2 + n + n0 ) x
2 2 2 2 2
puis :
Z2 Z2 
 π   π  1  π   π 
I= sin (1 + n0 ) x sin (1 + n) x dx = cos (n − n0 ) x − cos (2 + n + n0 ) x dx
2 2 2 2 2
0 0

Si n 6= n0 :
 i2 
1 2 h 
0 π
i2 2 h 
0 π
I= sin (n − n ) x − sin (2 + n + n ) x =0
2 (n − n0 )π 2 0 (2 + n + n0 )π 2 0

Si n = n0 :
Z2  
1 1 1
I= (1 − cos ((2 + 2n)πx)) dx = [x]20 + [sin((2 + 2n)πx]20 =1
2 2 (2 + 2n)π
0

37
 Le membre gauche de l’égalité initiale s’écrit :

Z2 X
∞  π   π 
I1 = En sin (1 + n0 ) x sin (1 + n) x dx
2 2
0 n=0

Z2  π   π 
I1 = En0 sin (1 + n0 ) x sin (1 + n0 ) x dx = En0
2 2
0

D’autre part :

Z2  π  2 h 
0 π
i2
sin (1 + n0 ) x = − cos (1 + n ) x
2 (1 + n0 )π 2 0
0

Z2  π  2
sin (1 + n0 ) x = − cos((1 + n0 )π) − 1


2 0
(1 + n )π
0

Z2  π  2
sin (1 + n0 ) x = 1 − cos((1 + n0 )π)


2 0
(1 + n )π
0

L’égalité :

Z2 X
∞ Z2
 π  
0 π   π 
En sin (1 + n ) x sin (1 + n) x = sin (1 + n0 ) x
2 2 2
0 n=0 0

se réduit donc à :
2
1 − cos((1 + n0 )π)


En0 = 0
(1 + n )π
On obtient donc la valeur recherchée du coefficient En
2 (1 − cos((1 + n)π))
En =
(1 + n)π

On a finalement :
∞
X 2 (1 − cos((1 + n)π))  π 
p(x, t = 0+ ) = sin (1 + n) x = 1
(1 + n)π 2
n=0

En observant que :

(1 − cos(1 + n)π)
=1 si n = 0, 2, 4 . . . n pair
2
(1 − cos(1 + n)π)
=0 si n = 1, 3, 5 . . . n impair
2
La solution complète s’écrit :
∞  2

4X 1  π  2π
p(x, t) = sin (1 + 2n) x exp −(1 + 2n) t
π (1 + 2n)π 2 4
n=0

38
Exercice 6
On applique la même méthode que pour l’exercice précédent. On aboutit à :
2
 
 π  2π
p(x, t) = En sin (1 + n) x exp −(1 + n) t
2 4
Construction de la solution du problème complet sous forme d’une série

La famille de solution :
2
 
 π  2π
p(x, t) = En sin (1 + n) x exp −(1 + n) t
2 4
ne satisfait pas la condition initiale

p(x, t = 0+ ) = kx 0 ≤ x ≤ 2

On va rechercher la solution du problème complet sous forme d’une série :

∞  2

X  π  2π
p(x, t) = En sin (1 + n) x exp −(1 + n) t
2 4
n=0

on en déduit :
∞
X  π 
p(x, t = 0+ ) = En sin (1 + n) x
2
n=0
La condition initiale impose :
∞
X  π 
p(x, t = 0+ ) = En sin (1 + n) x = kx 0 ≤ x ≤ 2
2
n=0

Le problème posé est donc de trouver la forme des coefficients En telle que :
∞
X  π 
En sin (1 + n) x = kx
2
n=0

Le principe de la méthode suivante est “d’extraire” le terme En de la sommation infinie. On

va utiliser pour cela une propriété dite d’orthogonalité de la base des fonctions sin (1 + n) π2 x
vis-à-vis du produit scalaire sur 0, 2.

On multiplie tout d’abord les deux membres de l’égalité :

∞
X  π 
En sin (1 + n) x = kx
2
n=0

par le facteur :  π 
sin (1 + n0 ) x
2
soit :
∞
 π X
0
 π   π 
sin (1 + n ) x En sin (1 + n) x = kx sin (1 + n0 ) x
2 2 2
n=0
ou encore :
∞
X  π   π   π 
En sin (1 + n0 ) x sin (1 + n) x = kx sin (1 + n0 ) x
2 2 2
n=0

39
 On intègre ensuite les deux membres de l’égalité précédente entre x = 0 et x = 2 :

Z2 X
∞ Z2
 π  
0 π   π 
En sin (1 + n ) x sin (1 + n) x dx = k x sin (1 + n0 ) x dx
2 2 2
0 n=0 0

On réécrit le premier membre en inversant les opérations d’intégration et de sommation et

en extrayant de l’intégrale le terme constant En :

∞ Z2 Z2
X  π  0
 π   π 
En sin (1 + n ) x sin (1 + n) x dx = k x sin (1 + n0 ) x dx
2 2 2
n=0 0 0

Le calcul du membre de gauche de l’égalité a été effectué lors de l’exercice 5 :

∞ Z2
X  π   π 
En sin (1 + n0 ) x sin (1 + n) x dx = En0
2 2
n=0 0

Il faut maintenant calculer le membre de droite :

Z2  π 
I=k x sin (1 + n0 ) x dx
2
0

par partie avec u = x, du = dx, dv = sin (1 + n0 ) π2 x et

En procédant par intégration
2 0 π


v = − (1+n 0 )π cos (1 + n ) 2 x :

Z2   2 Z2
 π  0 2 
0 π 2 
0 π

x sin (1 + n ) x dx = − cos (1 + n ) x x + cos (1 + n ) x dx
2 (1 + n0 )π 2 0
0 (1 + n )π 2
0 0

Z2  π  4 4 h 
0 π
i2
x sin (1 + n0 ) x dx = − 0


cos (1 + n )π + sin (1 + n ) x dx
2 (1 + n0 )π (1 + n0 )2 π 2 2 0
0

Z2  π  4 4
x sin (1 + n0 ) x dx = − cos (1 + n0 )π + sin (1 + n0 )π



2 0
(1 + n )π | {z 0 2 2
} (1 + n ) π | {z }
0 =(−1)n0 +1 =0

d’où : !
0 0
4(−1)n +1 4k (−1)n
I=k − =
(1 + n0 )π π 1 + n0
d’où la valeur de En :
4k (−1)n
En =
π 1+n
et l’expression de la solution complète :
∞
4k X (−1)n 2
 
 π  2π
p(x, t) = sin (1 + n) x exp −(1 + n) t
π 1+n 2 4
n=0

40
Exercice 7
Application de la méthode de séparation des variables :

y(x, t) = f (x)g(t)

L’équation :
∂2y 2
2∂ y
= a
∂t2 ∂x2
s’écrit :
∂ 2 g(t) ∂ 2 f (t)
f (x)g 00 (t) = a2 f 00 (x)g(t) avec g 00 (t) = et f 00 (t) =
dt2 dx2
soit :
f (x)g 00 (t) a2 f 00 (x)g(t)
=
a2 f (x)g(t) a2 f (x)g(t)
ou :
g 00 (t) f 00 (x)
=
a2 g(t) f (x)
d’où :
g 00 (t) f 00 (x)
= = −λ
a2 g(t) f (x)
On pose, par commodité d’écriture λ = ω 2 soit :
g 00 (t)
= −ω 2
a2 g(t)

f 00 (x)
= −ω 2
f (x)
La solution des deux équations précédentes est de la forme :

f (x) = A cos(ωx) + B sin(ωx)

g(t) = C cos(aωt) + D sin(aωt)

La solution de l’équation est donc de la forme :

y(x, t) = (A cos(ωx) + B sin(ωx))(C cos(aωt) + D sin(aωt))

d’où :
∂y(x, t)
= aω(A cos(ωx) + B sin(ωx))(−C sin(aωt) + D cos(aωt))
dt
La prise en compte des conditions aux limites :

y(x = 0, t > 0) = 0

y(x = L, t > 0) = 0
conduit à :
y(x = 0, t > 0) = A(C cos(aωt) + D sin(aωt)) = 0
d’où :
A=0
y(x = L, t > 0) = B sin(Lω)(C cos(aωt) + D sin(aωt)) = 0

41
qui implique :
nπ
sin(Lω) = 0 d’où Lω = nπ soit ω = , n = 0, 1, 2, . . .
L
Il reste à vérifier les conditions initiales :

y(x, t = 0+ ) = f (x)

∂y
(x, t = 0+ ) = 0
yt (x, t = 0+ ) =
∂t
On exploite tout d’abord la seconde équation :
∂y
= aωB sin(ωx)(−C sin(aωt) + D cos(aωt))
∂t
∂y
(x, t = 0+ ) = 4ωB sin(ωx)D cos(aωt)) = 0
∂t
d’où :
D=0
Il reste donc :
y(x, t) = B sin(ωx)C cos(aωt)
soit :
y(x, t) = E sin(ωx) cos(aωt)
ou :  nπx   
nπat
y(x, t) = En sin cos
L L
On applique le principe de superposition pour vérifier la condition initiale supplémentaire :

y(x, t = 0+ ) = f (x)

et construire la solution à l’aide d’une série infinie :

∞  nπx   
X nπat
y(x, t) = En sin cos
L L
n=1

en vérifiant la condition initiale :

∞
X  nπx 
y(x, t = 0+ ) = En sin = f (x)
L
n=1

Pour résoudre ce problème, on peut utiliser la technique du prolongement en série de Fourier

de la fonction f (x) définie sur l’intervalle [0, L]. On choisit évidemment ici, pour le prolongement,
une fonction g(x) impaire, de période 2L, égale à f pour 0 < x < L, de manière à obtenir une
série de sinus.

Dans le cas général, soit une fonction g(x) de période 2L, définie sur [−L, L], le développement
en série de Fourier de la fonction g(x) s’écrit :
∞ 
X  nπx   nπx 
g(x) = a0 + an cos + bn sin
L L
n=1

42
 ZL
1
a0 = g(x)dx
2L
−L

ZL ZL
2  nπx  1  nπx 
an = g(x) cos dx = g(x) cos dx
2L L L L
−L −L

ZL ZL
2  nπx  1  nπx 
bn = g(x) sin dx = g(x) sin dx
2L L L L
−L −L

Si g(x) est une fonction paire, on a g(−x) = g(x) alors :

 nπx   nπx 
g(x) cos est paire et g(x) sin est impaire
L L
d’où :
ZL ZL
1 2  nπx 
bn = 0 a0 = g(x)dx an = g(x) cos dx
2L L L
−L 0

Si g(x) est une fonction impaire, on a g(−x) = −g(x) alors :

 nπx   nπx 
g(x) cos est impaire et g(x) sin est paire
L L
d’où :
ZL
2  nπx 
a0 = an = 0 bn = g(x) sin dx
L L
0

Pour l’exercice, c’est ce second cas qui nous intéresse. On considère une fonction g(x) impaire
telle que :

g(x) = f (x) pour 0 < x < L

g(−x) = −g(x) = −f (x) pour −L < x < 0

On a :
∞
X  nπx 
g(x) = bn sin
L
n=1
avec :
ZL ZL
2  nπx  2  nπx 
bn = g(x) sin dx = f (x) sin dx
L L L L
0 0

On identifie donc le coefficient En :

∞
X  nπx 
+
y(x, t = 0 ) = En sin = f (x)
L
n=1

ZL
2  nπx 
En = f (x) sin dx
L L
0

43
et la solution du problème s’écrit :

ZL
 
∞  
X 2  nπu   nπx  nπat
y(x, t) =  f (u) sin du sin
 cos
L L L L
n=1 0

Les termes de cette série représentent les modes naturels ou normaux
de vibration. La fré-
quence fn du nième mode normal s’obtient à partir du terme en cos nπat
L et a pour expression :
nπa na
2πfn = ou fn =
L 2L

Exercice 8
Application de la méthode de séparation des variables :

θ(x, y) = f (x)g(y)

L’équation :
∂2 ∂2
θ(x, y) + θ(x, y) = 0
∂x2 ∂x2
s’écrit :
∂ 2 g(y) ∂ 2 f (x)
f (x)g 00 (y) + f 00 (x)g(y) = 0 avec g 00 (y) = et f 00 (t) =
dy 2 dx2
soit :
f 00 (x)g(y) f (x)g 00 (y)
=−
f (x)g(y) f (x)g(y)
ou :
f 00 (x) g 00 (y)
=− =
f (x) g(y)
d’où :
g 00 (t) f 00 (x)
2
= = −λ2
a g(t) f (x)
soit :
f 00 (x)
= −λ2
f (x)
g 00 (y)
= λ2
g(y)
La solution des deux équations précédentes est de la forme :

f (x) = A cos(λx) + B sin(λx)

g(y) = C cosh(λt) + D sinh(λy)

La solution de l’équation est donc de la forme :

θ(x, y) = (A cos(λx) + B sin(λx))(C cosh(λy) + D sinh(λy))

La prise en compte des conditions aux limites :

θ(x = 0, 0 ≤ y ≤ 1) = θ(x = 1, 0 ≤ y ≤ 1) = 0
θ(0 ≤ x ≤ 1, y = 0) = 0

44
 θ(0 ≤ x ≤ 1, y = 1) = θ1

conduit à
θ(x = 0, 0 ≤ y ≤ 1) = A(C cosh(λy) + D sinh(λy)) = 0
d’où :
A=0
θ(0 ≤ x ≤ 1, y = 0) = 0 = C(A cos(λx) + B sin(λx)) = 0
d’où :
C=0
soit :
θ(x, y) = E sin(λx)) sinh(λy))
θ(x = 1, 0 ≤ y ≤ 1) = E sin(λx)) sinh(λy)) = 0
soit :
λ = (1 + n)π n = 0, 1, 2, . . .
d’où :
θ(x, y) = E sin((1 + n)πx)) sinh((1 + n)πy))
Pour satisfaire la dernière conditions aux limites θ(0 ≤ x ≤ 1, y = 1) = θ1 , on applique le
théorème de superposition et l’on construit une solution à l’aide d’une série infinie :
∞
X
θ(x, y) = En sin((1 + n)πx)) sinh((1 + n)πy))
n=0

∞
X
θ(0 ≤ x ≤ 1, y = 1) = θ1 = En sin((1 + n)πx)) sinh((1 + n)π))
n=0
soit :
∞ 
X 
En sinh((1 + n)π)) sin((1 + n)πx)) = θ1
n=0

On retrouve ci-dessus l’expression d’une série de Fourier en sinus. On peut donc considérer,
pout l’obtenir, un prolongement de la fonction sur l’intervalle [−1, 1], fonction de période 2, de
la fonction f (x) = θ1 :

g(x) = f (x) = θ1 pour 0 < x < 1

g(−x) = −f (x) = −θ1 pour −1 < x < 0

Au point de discontinuité, on remplace f (x) par 21 (f (0 + ε) + f (0 − ε)) = 0.

On a :
∞
X  nπx 
g(x) = bn sin
L
n=1
avec :
ZL ZL
2  nπx  2  nπx 
bn = g(x) sin dx = f (x) sin dx
L L L L
0 0
avec :
L=1 f (x) = θ1 0<x<1

45
donc :
Z1 Z1
bn = 2 f (x) sin(nπx)dx = 2θ1 sin(nπx)dx
0 0
 
1 2θ1
bn = 2θ1 − [cos(nπx)]10 = −
(cos(nπ) − 1)
nπ nπ
 
4θ1 1 − cos(nπ)
bn =
nπ 2
1 − cos(nπ)
=1 pour n = 1, 3, 5, . . .
2
avec : 1 − cos(nπ)
=0 pour n = 2, 4, 6, . . .
2

46
 1ère année E.N.S.G
Travaux Dirigés en Mathématiques
Séries de Fourier

Exercice 1
Chercher les coefficients de Fourier correspondant à la fonction :

0 −5 < x < 0
f (x) =
3 0<x<5

Ecrire la série de Fourier correspondante. Comment doit-on définir f (x) aux points −5, 0 et
5 pour que la série de Fourier converge vers f (x) pour −5 < x < 5.

Exercice 2
Développer en série de Fourier la fonction périodique, de période 2π, f (x) = x2 , 0 < x < 2π.
Démontrer alors l’égalité :
1 1 1 π2
1+ + + + ... =
4 9 16 6

Exercice 3
Montrer qu’une fonction paire ne peut avoir aucun terme en sinus dans son développement
en série de Fourier cosinus.

Exercice 4
Développer f (x) = x, 0 < x < 2 en série de Fourier sinus puis en série de Fourier cosinus.

Exercice 5
Ecrire l’identité de Parseval correspondant à la série de Fourier cosinus de l’exercice précédent.
En déduire la somme S de la série :
1 1 1
1+ 4
+ 4 + 4 + ...
2 3 4

Exercice 6
On considère un signal périodique f (t), de période T et sa décomposition en série de Fourier :
+∞ Z
X 2π 1
cn einωt avec ω = et cn = T /2f (t)e−inωt
n=−∞
T T
−T /2

Montrer que si f (t) est paire, on a :

Z T /2
2
cn = f (t) cos(nωt)dt
T 0

47
 T T
Montrer que, si de surcroît, f (t) possède un centre de symétrie pour t = + k , on a :
4 2
cn = 0 pour n pair
ZT /4
4
cn = f (t) cos(nωt)dt pour n impair
T
0

On considère le signal f(t) de période T représenté ci-dessous :

Ce signal est pair et possède un centre de symétrie en T /4. Calculer la valeur de t1 pour
obtenir c3 = 0. Calculer les harmoniques 1, 5 et 7. Quels sont les autres harmoniques nuls ?

48
 1ère année E.N.S.G
Travaux Dirigés en Mathématiques
Séries de Fourier

Soit f (x) une fonction définie sur (−L, L) (période 2L). Le développement en série de Fourier
de f (x) est donné par :
∞
a0 X   nπx   nπx 
f (x) = + an cos + bn sin
2 L L
n=1

ZL
1  nπx 




 an = f (x) cos dx


 L L
−L


n = 0, 1, 2, . . .
ZL


1  nπx 




 bn = f (x) sin dx


 L L
−L
En notation complexe :
+∞ ZL
X inπx 1 inπx
f (x) = cn e L avec cn = f (x)e− L dx
n=−∞
2L
−L

Exercice 1

La période de cette fonction vaut 2L = 10 d’où L = 5 :

Z5
1  nπx 
an = f (x) cos dx
5 5
−5

Z0 Z5
1  nπx  1  nπx 
an = 0 cos dx + 3 cos dx
5 5 5 5
−5 0

Z5
1  nπx 
an = 3 cos dx
5 5
0
3 5 h  nπx i5
an = sin = 0 si n 6= 0
5 nπ 5 0
Z5   Z5
3 0πx 3
Si n = 0, a0 = cos dx = dx = 3
5 5 5
0 0

49
 Z5
1  nπx 
bn = f (x) sin dx
5 5
−5

Z0 Z5
1  nπx  1  nπx 
bn = 0 sin dx + 3 sin dx
5 5 5 5
−5 0

Z5
3  nπx 
bn = sin dx
5 5
0
  nπx i5 
3 5 h 3
bn = − cos = (1 − cos(nπ))
5 nπ 5 0 nπ
La série de Fourier correspondante s’écrit :
∞  
3 X 3  nπx 
f (x) = + (1 − cos(nπ)) sin
2 nπ 5
n=1

soit :        
3 6 πx 1 3πx 1 5πx
f (x) = + sin + sin + sin + ...
2 π 5 3 5 5 5
Puisque f (x) satisfait les conditions de Dirichlet (f (x) est définie et univoque sauf peut être
en un nombre fini de points de l’intervalle (−L, L), f (x) est périodique en dehors de l’intervalle
(−L, L) avec une période 2L, f (x) et f 0 (x) sont continues par morceaux dans (−L, L)), on peut
+ (x− )
dire que la série converge vers f (x) en tout point de continuité et vers f (x )+f
2 aux points de
discontinuité. Aux points x = −5, 0 et 5 qui sont des points de discontinuité, la série converge
vers 3+0 3
2 = 2 . Si l’on redéfinit f (x) comme suit :


 3/2 x = −5
 0 −5 < x < 0


f (x) = 3/2 x = 0
3 0<x<5




3/2 x = 5


alors la série converge vers f (x) pour −5 < x < 5.

Exercice 2
La période de la fonction vaut 2L = 2π d’où L = π.

Z2π Z2π
1  nπx  1  nπx 
an = f (x) cos dx = x2 cos dx
π π π π
0 0

uv 0 = [uv] − u0 v
R R
En intégrant par partie :

Z2π 2π Z2π

 

1  nπx  1 sin(nx) sin(nx) 
an = x2 cos dx =  x2 − 2x dx
π π π n 0 n
0 0

50
  2π !
1 sin(nx) cos(nx) sin(nx) 4
an = x2 + 2x −2 = pour n 6= 0
π n n2 n3 0 n2

Z2π
1 8π 2
Si n = 0, a0 = x2 dx =
π 3
0

Z2π  2π !
1  nπx  1 cos(nx) sin(nx) cos(nx)
bn = x2 sin dx = −x2 + 2x +2
π π π n n2 n3 0
0
4π
bn = −
n
On a donc :
∞
4π 2 X
 
4 4π
f (x) = x2 = + cos(nx) − sin(nx)
3 n2 n
n=1

Code Matlab

%--------Laboratoire d’Automatique et de Recherche Appliquée--------%

% %
% Développement en série de Fourier de la fonction f(x) = x*x %
% avec 0 < x < 2*pi %
% %
%----------------------------DM/06/10/92----------------------------%
%
n = input(’Nombre de termes de la série : ’);
pas = 2*pi/100;
k=0;
s = zeros(200,1);
for x=-2*pi:pas:2*pi
k = k + 1;
for n=1:n
s(k) = s(k) + (4/(n*n))*cos(n*x) - (4*pi/n)*sin(n*x);
end
s(k) = s(k) + (4*pi*pi)/3;
end
x = -2*pi:pas:2*pi;
x1 = 0:pas:2*pi;
x2 = ( x1.*x1 x1.*x1);
plot(x,s,x,x2)
title(( ’Développement pour n = ’,num2str(n)))

51
 ∞
4π 2 X 4
Au point x = 0, la série de Fourier se réduit à + . D’après les conditions de
3 n2
n=1

52
Dirichlet, la série converge, au point x = 0, vers 12 (0 + 4π 2 ) = 2π 2 . On a donc l’égalité :
∞ ∞
4π 2 X 4 2
X 1 π2
+ = 2π d’où =
3 n2 n2 6
n=1 n=1

Exercice 3
Si l’on pose :
∞
a0 X   nπx   nπx 
f (x) = + an cos + bn sin
2 L L
n=1

alors
∞
a0 X   nπx   nπx 
f (−x) = + an cos − bn sin
2 L L
n=1

Si f (x) est paire alors on a f (x) = f (−x), d’où :

∞ ∞
a0 X   nπx   nπx  a
0
X  nπx   nπx 
+ an cos + bn sin = + an cos − bn sin
2 L L 2 L L
n=1 n=1

c’est-à-dire :
∞
X  nπx 
2 bn sin =0 ⇒ aucun terme en sinus
L
n=1

Exercice 4
Pour être développée en série de Fourier sinus, f (x) doit être une fonction impaire, d’où la
figure suivante :

On a alors une période 2L = 4, d’où L = 2. On a donc an = 0 et :

ZL ZL
1  nπx  2  nπx 
bn = f (x) sin dx = f (x) sin dx
L L L L
−L 0

ZL    nπx    nπx 2

2  nπx  −2 −4
bn = x sin dx = (x) cos − (1) sin
L L nπ 2 n2 π 2 2 0
0
4
bn = − cos(nπ)
nπ

53
 On obtient donc :
∞   nπx  4   πx  1     
X 4 2πx 1 3πx
f (x) = − cos(nπ) sin = sin − sin + sin − ...
nπ 2 π 2 2 2 3 2
n=1

Pour être développée en série de Fourier cosinus, f (x) doit être une fonction paire, d’où la
figure suivante :

La demi-période L est égale à 2. On a donc bn = 0 et :

ZL ZL
1  nπx  2  nπx 
an = f (x) cos dx = f (x) cos dx
L L L L
−L 0

ZL    nπx    nπx 2

2  nπx  2 −4
an = x cos dx = (x) sin − (1) cos
L L nπ 2 n2 π 2 2 0
0
4
bn = (cos(nπ) − 1) si n 6= 0
n2 π 2
Z2
Si n = 0, a0 = xdx = 2
0
On obtient alors :
∞
X4  nπx 
f (x) = 1 + (cos(nπ) − 1) cos
n2 π 2 2
n=1
     
8  πx  1 3πx 1 5πx
1 + 2 cos + 2 cos + 2 cos + ...
π 2 3 2 5 2

Exercice 5
En résumé, on a :
4
L = 2, a0 = 2, an = (cos(nπ) − 1) pour n 6= 0, bn = 0
n2 π 2
L’identité de Parseval s’écrit :
ZL ∞
1 2 a2 X 2
(f (x)) dx = 0 + (an + b2n )
L 2
−L n=1

Z2 ∞
1 2
X 16
x dx = 2 + (cos(nπ) − 1)2
2 n4 π 4
−2 n=1

54
  
8 64 1 1 1
=2+ 4 4
+ 4 + 4 + ...
3 π 1 3 5
On obtient donc :
1 1 1 π4
+ + + . . . =
14 34 54 96
On a, d’autre part :
1 1 1
S =1+ 4
+ 4 + 4 + ...
2 3 4
   
1 1 1 1 1 1
S= + + + ... + + + + ...
14 34 54 24 44 64
   
1 1 1 1 1 1 1
S= + + + ... + 4 + + + ...
14 34 54 2 14 24 34
π4 1
S= + 4S
96 2
d’où :
15 π4 16 4 π 4
S= ⇒ π =
16 96 15.96 90

55
 1ère année E.N.S.G
Travaux Dirigés en Mathématiques
Transformation de Laplace

Exercice 1
En utilisant la définition de la transformée de Laplace, démontrer les propriétés suivantes :
1 s
a) L(f (at)) = F b) L(e−at f (t)) = F (s + a)
a a  
df (t)
c) L(f (t − τ )) = e−sτ F (s) d) L = sF (s) − f (0+ )
∞  dt
Z
e) L  f1 (t − τ )f2 (τ )dτ  = F1 (s)F2 (s)
0

Exercice 2
Calculer les transformées de Laplace des fonctions suivantes :

a) tn (avec n entier) b) eat t2 c) e−at sin(ωt) d) e−at cos(ωt)

Exercice 3
Résoudre l’équation différentielle suivante à l’aide de la transformée de Laplace :

d2 y(t) dy(t)
2
+ y(t) = t avec les conditions initiales y(0) = 1 et = −2
dt dt t=0

Exercice 4
Calculer la transformée de Laplace de la fonction périodique, de période 2T , suivante :

Exercice 5
Résoudre l’équation intégrale (équation de Volterra) :

Zt
y(u)f (t − u)du = g(t) avec f (t) = e−t et g(t) = sin(t)
0

Exercice 6
Résoudre le système d’équations différentielles suivant :

dy1 (t) dy2 (t)
 3 dt + dt + 2y1 (t) = 1



avec y1 (0) = 0 et y2 (0) = 0

 dy1 (t) dy2 (t)
+4 + 3y2 (t) = 0


dt dt

56
Exercice 7
Calculer les transformées inverses de Laplace suivantes :
   
−1 1 −1 1 − 2s
L , L
1 + s + s2 s(1 + s)2

Exercice 8
Calculer à l’aide de l’intégrale de convolution :
 
−1 s
L
(s + 1)(s2 + 1)

57
 1ère année E.N.S.G
Travaux Dirigés en Mathématiques
Transformation de Laplace

Exercice 1
Z∞ Z∞
−st
a) L(f (at)) = f (at)e dt = f (at)e−sat/a dt
0 0
En posant u = at, on a du = adt d’où dt = du/a, donc :
Z∞ Z∞
du 1 s 1 s
f (u)e−su/a = e− a u du = F
a a a a
0 0

Z∞ Z∞
−at −at −st
b) L(e f (t)) = e f (t)e dt = e−(s+a)t f (t)dt = F (s + a)
0 0
Z∞ Z∞
c) L(f (t − τ )) = f (t − τ )e−st dt = f (t − τ )e−s(t−τ +τ ) dt
0 0
En posant u = t − τ , on a du = dt, d’où :
Z∞ Z∞ Z∞
−su −sτ −su −sτ −sτ
L(f (t − τ )) = f (u)e e du = f (u)e e du = e f (u)e−su du
−τ 0 0

L(f (t − τ )) = e−sτ F (s)

  Z∞ Z∞
df (t) df (t) −st ∞
e dt = f (t)e−st 0 − −sf (t)e−st dt = sF (s) − f (0)

d) L =
dt dt
0 0
Z∞
e) Posons f (t) = f1 (t − τ )f2 (τ )dτ
 0
Z∞ Z∞ Z∞
 ∞ 
Z
L(f (t)) =  f1 (t − τ )f2 (τ )dτ  e−st dt = f2 (τ )  f1 (t − τ )e−st dt dτ
0 0 0 0
En posant u = t − τ , du = dt dans l’intégrale en t (τ est constant dans cette intégrale) et en
utilisant le théorème de Fubini (fonctions intégrables au sens de Lebesgue), on obtient :

Z∞
∞
Z∞
 ∞ 
Z Z
L(f (t)) = f2 (τ )  f1 (u)e−s(u+τ ) du dτ = f2 (τ )e−sτ  f1 (u)e−su du dτ
0 0 0 0
Z∞
  ∞
Z
L(f (t)) =  f2 (τ )e−sτ dτ   f1 (u)e−su du = L(f1 (t))L(f2 (t)) = F1 (s)F2 (s)
0 0

58
Exercice 2
Z∞  n −st ∞ Z∞  −st  Z∞
n −st t e e n
a) L(t ) =n
t e dt = − − nt n−1
− dt = tn−1 e−st dt
s 0 s s
0 0 0
n
L(t ) = L(tn−1 )
n
s
Sans faire une démonstration par récurrence rigoureuse, on peut partir des transformées
1 1
connues : L(t0 ) = L(f (t) ≡ 1) = et L(t) = 2 dont on peut vérifier qu’elle sont bien liées par
s s
la récurrence que l’on a établie.

Calculons L(t2 ) et L(t3 ) avec la formule de récurrence :

2 2 3 3×2
L(t2 ) = L(t) = 3 L(t3 ) = L(t2 ) =
s s s s4
On perçoit la construction du résultat :

n!
L(tn ) =
sn+1
2
b) On a L(t2 ) = 3 et la multiplication par eat correspond à une translation dans le domaine
s
fréquentiel, d’où :
2
L(eat t2 ) =
(s − a)3
ω at correspond à une translation dans le
c) On a L(sin(ωt) = s2 +ω 2 et la multiplication par e
domaine fréquentiel, d’où :
ω
L(e−at sin(ωt)) =
(s + a)2 + ω 2
s
d) On a L(cos(ωt) = s2 +ω 2
d’où :
s+a
L(e−at cos(ωt)) =
(s + a)2 + ω 2

Exercice 3
En transformée de Laplace :
1
s2 Y (s) − sy(0+ ) − y 0 (0+ ) + Y (s) =
s2
1
(s2 + 1)Y (s) − s + 2 =
s2
1 s−2 1 1 s 2
Y (s) = + = 2− 2 + 2 − 2
s2 (s2 + 1) 2
s +1 s s +1 s +1 s +1
1 s 3
Y (s) = 2
+ 2 − 2
s s +1 s +1
d’où :
y(t) = t + cos(t) − 3 sin(t)

59
Exercice 4
On a :
Zt
y(u)et−u du = sin(t)
0
soit :
y(t) ∗ e−t = sin(t)
En prenant la transformée de Laplace de l’égalité, on obtient :
1 1
Y (s) = 2
s+1 s +1

s 1
Y (s) = +
s2 + 1 s2 + 1
et en revenant aux fonctions temporelles :

y(t) = cos(t) + sin(t)

Exercice 5
En transformées de Laplace :

 3(sY1 (s) − y1 (0+ )) + sY2 (s) − y2 (0+ ) + 2Y1 (s) = 1


s
+ +
sY1 (s) − y1 (0 ) + 4(sY2 (s) − y2 (0 )) + 3Y2 (s) = 0


 (3(s)3s + 2)Y1 (s) + sY2 (s) = 1


a)
s
sY1 (s) + (4s + 3)Y2 (s) = 0 b)


−sY1 (s)
b)→ Y2 (s) =
4s + 3
s2 Y1 (s) 1
a)→ (3s + 2)Y1 (s) − =
4s + 3 s
(3s + 2)(4s + 3) − s 2 1
Y1 (s) =
4s + 3 s
d’où :
4s + 3
Y1 (s) =
s(11s2 + 17s + 6
et :
−sY1 (s) −1
Y2 (s) = = 2
4s + 3 11s + 17s + 6
On décompose ensuite en éléments simples afin d’établir les originaux temporels.

4s + 3 4s + 3
Y1 (s) = 2
= 6
s(11s + 17s + 6 11s(s + 11 )(s + 1)
A B C 1 3 1
Y1 (s) = + 6 + = − 6 −
s s + 11 s+1 2s 10(s + 11 )
5(s + 1)

60
 1 3 1
y1 (t) = − e−6t/11 − e−t
2 10 5
−1 −1 A B
Y2 (s) = = 6 = 6 + s+1
11s2 + 17s + 6 11(s + 11 )(s + 1) s + 11
−1 1
Y2 (s) = 6 + 5(s + 1)
5(s + 11 )
1
y2 (t) = (e−t − e6t/11 )
5

Exercice 6
On a :
 
1 1 −1 1 1 −at
=  √ 2 or on a L = e sin(ωt)
1 + s + s2 1 2

3 (s + a)2 + ω 2 ω
s+ 2 + 2

d’où : √ !
1 2 3
2
= √ e−t/2 sin t
1+s+s 3 2
Puis :
1 − 2s A B C 1 1 3
2
= + + 2
= − −
s(1 + s) s 1 + s (1 + s) s 1 + s (1 + s)2
 
−1 1 − 2s
L = 1 − e−t − 3te−t
s(1 + s)2

Exercice 7
    
−1 s −1 1 s
L =L
(s + 1)(s2 + 1) s+1 2
s +1
or, on a :    
−1 1 −t −1 s
L =e et L 2
= cos(t)
(s + 1 s +1
On peut donc écrire :
  Zt Zt
−1 s −(t−τ ) −t
L = e cos(τ )dτ = e eτ cos(τ )dτ
(s + 1)(s2 + 1)
0 0

Zt
Posons I = eτ cos(τ )dτ , en intégrant par parties, on obtient :
0

Zt
I = [e τ
sin(τ )]t0 − eτ sin(τ )dτ
0

61
 Zt
 

I = [eτ sin(τ )]t0 − [−eτ cos(τ )]t0 + eτ cos(τ )dτ 

0

I = [e τ
sin(τ )]t0 − [−e cos(τ )]t0 − I
τ

2I = et sin(t) + et cos(t) − 1
On obtient donc :
  Zt
−1 s −t 1
L =e eτ cos(τ )dτ = e−t (et sin(t) + et cos(t) − 1)
(s + 1)(s2 + 1) 2
0
 
−1 s 1
L = (sin(t) + cos(t) − e−t )
(s + 1)(s2 + 1) 2

62