22EdO2 20-21
22EdO2 20-21
22EdO2 20-21
Tiaret
Département de Mathématiques
Le problème de Sturm-Liouville
Niveau M1
Ziane Mohamed
1 Introduction 3
2 Généralités 7
2.1 Espaces fonctionnels . . . . . . . . . . . . . . . . . . . . . . . 7
2.2 Le problème de Cauchy . . . . . . . . . . . . . . . . . . . . . . 8
3 Le problème de Sturm-Liouville 9
3.1 Propriétés . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 12
3.2 Alternative de Fredholm . . . . . . . . . . . . . . . . . . . . . 16
4 La fonction de Green 17
4.1 Fonction de Green . . . . . . . . . . . . . . . . . . . . . . . . 18
5 Travaux Dirigés 23
Chapitre 1
Introduction
Ly = f, (1.1)
1. Pour des raisons de simplicité, nous allons restreindre l'étude à l'intervalle fermé
borné [a, b] (a < b) bien que la théorie ne se limite nullement dans ce cadre.
4
où S est le sous espace vectoriel de fonctions C 2 ([a, b]) qui vérient les données
aux bords homogènes, λ est un paramètre, l'équation (1.3) est ainsi écrite
sous la forme (L − λI)y = 0 (I : désigne l'application identité).
et
S = {y ∈ C([a, b], R)/y(0) = y(π) = 0} .
5
et
Z π Z π
hy1 , Ly2 i = y1 (t)Ly2 (t) dt = − y1 (t)y200 (t) dt
0 0
π Z π
= −y1 (t)y20 (t) + y10 (t)y20 (t) dt
0 0
Z π
0 0
= −(y1 (π)y2 (π) − y1 (0)y2 (0)) + y10 (t)y20 (t) dt
0
Z π
= y10 (t)y20 (t) dt.
0
où λ est une constante arbitraire, a2 (t) > 0 et µ(t) > 0 pour tut t ∈ [a, b].
Posons :
Z t
a1 (t) a0 (t)p(t) µ(t)p(t)
p(t) = exp − dt , q(t) = , ω(t) =
a1 (t) a2 (t) a2 (t)
p(t)g(t)
et f (t) = .
a2 (t)
1
Multiplions l'équation (1.5) par p(t), on obtient une équation de la
a2 (t)
forme
d dy
− p(t) + (q(t) − λω(t))y(t) = f (t) (1.6)
dt dt
qui sera appelée équation de Sturm-Liouville.
3
On ajoute des conditions aux limites (séparés) du type
Généralités
p
on a donc kf kω = hf, f iω .
Théorème 2.1 Supposons que les fonctions a1 (t), a2 (t) et a3 (t) sont conti-
nues sur [a, b]. Alors le problème de Cauchy (2.1)-(2.2) admet une solution
unique.
0
y 0 1 y 0
= + .
v −a2 (t) −a1 (t) v a3 (t)
0 1 0
où A(t) = et b(t) = .
−a2 (t) −a1 (t) a3 (t)
y
Posons x = , il vient que
v
0 dx
x (t) = = A(t)x(t) + g(t), (2.4)
dt
y(t 0 )
x(t0 ) =
v(t0 )
Le problème de Sturm-Liouville
[a, b], où p ∈ C 1 ([a, b], R), p(t) > 0 dans [a, b], ω ∈ C([a, b], R),
sur l'intervalle
ω(t) > 0 pour tout t ∈ [a, b] et q ∈ C([a, b], R), avec les conditions aux limites
(
(Sλ ) Lλ [y](t) = f (t), pour t ∈ [a, b],
(P )
(F ) F1 [y] = 0, F2 [y] = 0.
Remarque 3.1 Les exemples les plus courants de conditions aux limites
sont :
y(a) = y(b) = 0 et y 0 (a) = y 0 (b) = 0.
1. de sorte que F1,2 [y] ne se réduit pas à la condition 0 = 0!
10
y 00 (t) + λy(t) = 0,
y 0 (0) = y 0 (π) = 0.
pour λ > 0 le problème admet seulement la solution trivial (y(t) ≡ 0).
Pour λ < 0 la solution générale de l'équation y00 (t) + λy(t) = 0 est de la
forme √ √
y(t) = C1 sin( λt) + C2 cos( λt).
Tenu compte les conditions aux limites y0 (0) = 0 et y0 (π) = 0, on obtient
√
C1 λ = 0,
√ √ √ √
C1 λ cos( λπ) − C2 λ sin( λπ) = 0.
√
Par conséquant C√1 = 0 et C2 sin( λπ) √ = 0. Si C2 = 0 alors y(t) = 0. Pour2
C2 6= 0 donc sin( λπ) = 0. Il vient λπ = kπ, k = 1, 2, · · · D'où, λk = k
sont les valeurs propres et yk (t) = cos(kt) sont les fonctions
Z propres. π
On munit l'espace C([0, π]) par le produit scalaire hf, gi = f (t)g(t) dt, les
0
fonctions propres forment un sous espace orthogonal
Z π
π
hyk , yk i = cos2 (kt) dt = ,
0 2
et pour k 6= m, on a :
Zπ
hyk , ym i = cos(kt) cos(mt) dt
0
Zπ
1
= [cos(k + m)t + cos(k − m)t] dt
2
"0 π #
1 sin(k + m)t sin(k − m)t
= + = 0.
2 k+m k−m
0
3.1 Propriétés 12
r
y (t) 2
Les fonctions ϕk (t) = k = cos(kt) (k = 1, 2, · · · ) forment un es-
kyk (t)k π
pace orthonormal dans C([0, π]).
Exemple 3.2 Soit le problème de Sturm-Liouville suivant :
y 00 (t) + λy(t) = 0, pour t ∈ [0, π],
y(0) = 0,
y(π) + y 0 (π) = 0.
Exemple 3.3 On considère l'équation d'Euler t2 y00 (t) + ty0 (t) + λ2 y(t) = 0
sur l'intervalle t ∈ [0, e], et les conditions y(1) = 0 et y(e) = 0.
d 0 λ2
ty (t) + y(t) = 0.
dt t
3.1 Propriétés
Z b
(vL0 [u] − uL0 [v]) dt
aZ
b
= − (v(pu0 )0 − u(pv 0 )0 ) dt
a b Z
b
0 0
= −(v(pu ) − u(pv )) + [v 0 (pu0 ) − u0 (pv 0 )] dt
a
a
= −p(b)[u0 (b)v(b) − u(b)v 0 (b)] + p(a)[u0 (a)v(a) − u(a)v 0 (a)]
= M [u, v](b) − M [u, v](a).
Théorème 3.2 Soient y1 , y2 ∈ C 2 ([a, b]) satisfont les conditions aux limites
(F ). Alors
Démonstration:
1. Comme y1 et y2 vérient la condition [F1 ], alors le système
(
α0 y1 (a) + α1 y10 (a) = 0
α0 y2 (a) + α1 y20 (a) = 0
possède une solution non triviale (α0 , α1 ) 6= (0, 0). Par conséquent
Zb Z b
hL0 [y1 ], y2 i−hy1 , L0 [y2 ]i = y2 (t)L0 [y1 ](t) dt− y1 (t)L0 [y2 ](t) dt = 0,
a
a
il resulte que hL0 [y1 ], y2 i = hy1 , L0 [y2 ]i. L'opérateur L0 est donc sy-
métrique sur le sous-espace S de fonctions dans C 2 ([a, b]) qui satisfait
aux conditions aux limites (F ).
1 i
y(t) = (y(t) + y(t)) + (iy(t) − iy(t)) = y1 (t) + iy2 (t),
2 2
1
L0 [y1 ] = L0 [(y + y)/2] = (λωy + λωy)
2
1
= λω(y + y)
2
= λωy1 ,
3.2 Alternative de Fredholm 16
de même
(
(N H) y 00 (t) + p(t)y 0 (t) + q(t)y(t) = f (t), pour t ∈ [a, b],
(P )nh
(F )nh F1 [y] = γ, F2 [y] = δ.
La fonction de Green
Lu = f.
2
Nous allons maintenant examiner un exemple simple dont nous expli-
querons l'existence et la détermination de la fonction de Green de façon
intuitive.
Exemple 4.1 Pour f ∈ C([a, b], R). Soit le problème à deux points suivants
y 00 (t) = f (t), pour t ∈ [0, 1],
y(0) = y(1) = 0.
1. George Green (1793-1841)
où
t(x − 1), si 0 ≤ t < x,
G(x, t) =
x(t − 1), si x ≤ t ≤ 1.
yp0 (t) = v10 (t)y1 (t) + v1 (t)y10 (t) + v20 (t) + y2 (t) + v2 (t) + y20 (t).
alors,
f (t)
v10 (t)y10 (t) + v20 (t)y20 (t) = .
p(t)
le système 0
v1 (t)y1 (t) + v20 (t)y2 (t) = 0,
f (t)
v1 (t)y10 (t) + v20 (t)y20 (t) = ,
p(t)
On a
f (t)y2 (t)
v10 (t) = − , (4.6)
p(t)W [y1 , y2 ](t)
f (t)y1 (t)
v20 (t) = . (4.7)
p(t)W [y1 , y2 ](t)
Supposons que,
a1 y1 (t0 ) + a2 y2 (t0 ) = 0
a1 y10 (t0 ) + a2 y20 (t0 ) = 0
4.1 Fonction de Green 21
Z x
v1 (x) = − Cy2 (t)f (t) dt
Z xb
v2 (x) = Cy1 (t)f (t) dt,
a
On peut écrire
Z x Z b Z b
yp (x) = G1 (x, t)f (t) dt + G2 (x, t)f (t) dt = G(x, t)f (t) dt.
a x a
4.1 Fonction de Green 22
∂G2 ∂G1 1
lim+ (t − , t) − lim+ (t + , t) = .
→0 ∂x →0 ∂x p(t)
Chapitre 5
Travaux Dirigés
Exercice 5.1 Résoudre les problèmes aux limites linéaires homogènes sui-
vants 00
y (t) + y(t) = 0, t ∈ [0, 1]
(5.1)
y(0) = y(1) = 0,
y 00 (t) + π 2 y(t) = 0, t ∈ [0, 1]
(5.2)
y(0) = y(1) = 0,
00
y (t) + π 2 y(t) = 0, t ∈ [0, 1]
y(0) + y(1) = 0, (5.3)
0
y (0) + y 0 (1) = 0,
00
y (t) + λy(t) = 0, t ∈ [a, b]
(5.4)
y(a) = y(b) = 0,
Solution 5.1 :
( √ √
k1 cos(a λ) + k2 sin(a λ) = 0,
√ √
k1 cos(b λ) + k2 sin(b λ) = 0.
√
Si d 6= 0, ce qui est possible si ∀k ∈ Z : (b − a) λ 6= kπ , alors
y ≡ 0. √
S'il existe k ∈ Z : (b − a) λ = kπ alors, d=0 et λ = λn =
n2 π 2
(b−a)2
, n ∈ Z.
Dans ce cas, le problème (5.4) admet un innité de solutions
et l'on a les fonctions propres.
√ √
Si sin(a λ) = 0 alors, k1 = 0 et sin(b λ) = 0 ; donc y(t) =
nπ
k sin(t (b−a) ).
√ √
Si cos(a λ) = 0 alors, k2 = 0 et cos(b λ) = 0 ; donc y(t) =
nπ
k cos(t (b−a) ).
√ √ √
Si sin(a λ) 6= 0 et cos(a λ) 6= 0 alors y(t) = k[sin(t λ) −
√ √
tan(a λ) cos(t λ)], k ∈ R∗ et donc, à chaque fois, le problème
(5.4) admet une innité (et non une double innité comme pour
le problème (5.3)) de solutions.
et
S = {y ∈ C([a, b], R)/y(0) = y(π) = 0} .
Si y1 , y 2 ∈ S , on a : y1 (0) = y1 (π) = y2 (0) = y2 (π) = 0, nous obtenons en
intégrant par parties :
Z π Z π
hLy1 , y2 i = Ly1 (t)y2 (t) dt = − y100 (t)y2 dt
0 0
π Z π
= −y10 (t)y2 (t) + y10 (t)y20 (t) dt
0 0
Z π
0 0
= −(y1 (π)y2 (π) − y1 (0)y2 (0)) + y10 (t)y20 (t) dt
0
Z π
= y10 (t)y20 (t) dt,
0
26
et
Z π Z π
hy1 , Ly2 i = y1 (t)Ly2 (t) dt = − y1 (t)y200 (t) dt
0 0
π Z π
= −y1 (t)y20 (t) + y10 (t)y20 (t) dt
0 0
Z π
0 0
= −(y1 (π)y2 (π) − y1 (0)y2 (0)) + y10 (t)y20 (t) dt
0
Z π
= y10 (t)y20 (t) dt.
0
et pour k 6= m, on a :
Zπ
hyk , ym i = cos(kt) cos(mt) dt
0
Zπ
1
= [cos(k + m)t + cos(k − m)t] dt
2
"0 π #
1 sin(k + m)t sin(k − m)t
= + = 0.
2 k+m k−m
0
r
yk (t) 2
Les fonctions ϕk (t) = = cos(kt) (k = 1, 2, · · · ) forment
kyk (t)k π
un espace orthonormal dans C([0, π]).
y(0) = 0, (5.7)
y(π) + y 0 (π) = 0.
√ √
y(t) = c1 sin( λt) + c2 cos( λt).
on a y(0) = 0 on obtient c2 = 0 En utilisant la deuxième condition y(π) +
0
y (π) = 0 pour obtenir
√ √ √
c1 sin( λπ) + c1 λcos( λπ) = 0.
pour c1 6= 0 on a √ √
tg( λπ) = − λ.
√
On pose µ= λπ , alors
µ
tg(µ) = − .
π
28
µ2n µn
λn = 2 et tg(µn ) = − .
π π
µ
Les solutions de tg(µ) = − peuvent être considérées comme les points
π
t
d'intersection graphes des fonctions y = tg(t) et y = . Cette équation a une
π
innité solutions. Les fonctions autonomes respectives sont :
µn √
yn (t) = sin( )t = sin( λn t).
π
Exercice 5.5 Résoudre l'équation d'Euler
t2 y 00 (t) + ty 0 (t) + λ2 y(t) = 0
[3] G. Birkho, G-C. Rota, Ordinary dierential equations. John Wiley &
Sons (1989).
[4] W.E. Boyce and R.C. Diprima, Elementary Dierential Equations with
Boundary Value Problems, John Wiley & Sons (1986).
[5] N. El Hage Hassan, Topologie et espaces normés, Dunod. 2011.
[11] S. Djebali, Problème aux limites associés aux EDO du second ordre,
cours polycopié. E.N.S.-Kouba, 2006.