Fonctions Speciales Et Polynomes Orthogonaux Cour-75-96

Chapitre 5
Les polynômes orthogonaux
Les polynômes orthogonaux sont très importants en physique. Ce sont les solutions
d’équations qui surviennent très souvent lorsqu’un problème possède une symétrie
sphérique. Les polynômes orthogonaux sont profondément liés à la mécanique
quantique et ils sont présents dans de nombreuses applications, telles que la théorie
électromagnétique, l’hydrodynamique et la conduction thermique. En ingénierie,
les polynômes orthogonaux apparaissent dans de nombreuses applications, par ex-
emple dans la théorie des lignes de transmission, la théorie des circuits électriques,
la physique des réacteurs nucléaires ainsi que la sismologie.
5.1 Polynômes de Legendre

Adrien-Marie Legendre a introduit, en 1784, les polynômes de Legendre, tout en
étudiant l’attraction des sphéroı̈des et des ellipsoı̈des. Ces polynômes sont les solu-
tions d’une équation différentielle ordinaire appelée équation différentielle de Leg-
endre. Cette équation est fréquemment rencontrée en physique et en ingénierie.
En particulier, cela se produit lors de la résolution de l’équation de Laplace en
coordonnées sphériques. L’équation de Legendre est donnée par
d2 y dy
(1 − x2 ) − 2x + l(l + 1)y = 0 (5.1)
dx2 dx
dans le cas générale l ∈ R. L’équation (5.1) est solvable par la méthode de Frobe-
nius (qui a été utilisée pour résoudre les équations de Bessel). Les solutions ont
la forme de séries entières et elles sont convergentes uniquement dans le domaine
−1 < x < 1. Dans le cas où l ∈ N, il est possible d’obtenir des solutions qui sont
régulières aux points x ± 1, et pour lesquelles la série s’arrête au terme de degré l.
Dans le dernièr cas, la solution de l’Eq. (5.1) est appellée polynône de Legendre.
Elle est donnée par
[l/2]
X (2l − 2r)!
Pl (x) = (−1)r xl−2r −1≤ x≤1 (5.2)
r=0
2l r!(l − r)!(l − 2r)!
67
Notons que la solution (5.2) est convergente seulement dans le domaine [−1, 1].
Les points x ± 1 sont des points singuliers réguliers de l’équation différentielle
(5.1).
Démonstration
Avant de chercher la solution de l’Eq. (5.1), cherchons d’abord le domaine de con-
vergence. Pour cela, on utilise le théorème suivant; si on a une équation différentielle
de la forme
d2 y dy
x2 2
− xq(x) + p(x)y = 0 (5.3)
dx dx
et si q(x) et p(x) peuvent être developpées sous les formes
∞
X
q(x) = qn x n ,
n=0
X∞
p(x) = pn x n
n=0
donc, le domaine de convergence de la solution de l’Eq. (5.3) est la même que le

domaine de convergence de q(x) et p(x).
Pour déterminer les fonction p(x) et q(x), on écrit l’Eq. (5.1) sous la forme
d2 y 2x2 dy l(l + 1)x2

x2 − x + y=0
dx2 1 − x2 dx 1 − x2
En comparant avec l’Eq. (5.3), on déduit que
2x2
q(x) = ,
1 − x2
l(l + 1)x2
p(x) =
1 − x2
Le développement en série entière des fonctions q(x) et p(x) est
∞
X
q(x) = 2x 2
x2m ,
m=0
∞
X
p(x) = l(l + 1)x 2
x2m
m=0
donc les fonctions q(x) et p(x) sont convergentes seulement dans l’intervalle −1 <
x < 1. Pour résoudre l’Eq. (5.1), on utilise la méthode de Frobenius qui consiste à
chercher des solutions sous forme de séries entières
∞
X
y(x, s) = ar xr+s
r=0
CHAPITRE 5. Les polynômes orthogonaux 68
Si en insère y, sa dérivée y0 et la dérivée seconde y00 dans l’Eq. (5.1), on trouve

trois équations indiciales
a0 s(s − 1) = 0,
a1 s(s + 1) = 0,
an+2 (s + n + 2)(s + n + 1) − an {(s + n)(s + n + 1) − l(l + 1)} n≥0
Pour a0 et a1 arbitraires et différents de 0, les équations indiciales admettent comme
solutions
s = 0, s=1
et
n(n + 1) − l(l + 1)
an+2 = an
(n + 1)(n + 2)
(n − l)(l + n + 1)
= an
(n + 1)(n + 2)
Pour déduire la forme du terme général an en fonction de a0 (ou a1 ) et n, on con-
sidére quelques cas particuliers
a2 = −a0 l(l + 1)
(1 − l)(l + 2)
a3 = a1
2.3
(l − 2)(l + 3)
a4 = −a2
3.4
l(l − 2)(l + 1)(l + 3)
= a0
1.2.3.4
(l − 3)(l + 4)
a5 = −a3
4.5
(l − 1)(l − 3)(l + 2)(l + 4)
= a1
1.2.3.4.5
A partir des équations précédentes on peut déduire le terme général qui est donné
par
(−1)n a0
a2n = l(l − 2)(l − 4) . . . (l − 2n + 2)(l + 1)(l + 3) . . . (l + 2n − 1)
(2n)!
(−1)n a1
a2n+1 = (l − 1)(l − 3) . . . (l − 2n + 1)(l + 2)(l + 4) . . . (l + 2n)
(2n + 1)!
Donc, la solution générale peut être écrite sous la forme
∞
 
 X l(l − 2)(l − 4) . . . (l − 2n + 2)(l + 1)(l + 3) . . . (l + 2n − 1) 
y(x, 0) = a0  + n 2n 
 
1 (−1) x

(2n)!

 


n=1
∞
 
n (l − 1)(l − 3) . . . (l − 2n + 1)(l + 2)(l + 4) . . . (l + 2n) 2n+1 
 X 
+ a1  +
 
x (−1) x

(2n + 1)!

 


n=1
= a0 y1 (x) + a1 y2 (x)
69
Nous soulignons que cette solution est valable dans l’intervalle −1 < x < 1. Afin
d’étendre la solution aux points x = −1 and x = 1, nous utilisons l’observation
suivante: si l est pair, ı.e., l = 2n, on a
a2 , 0, a4 , 0, . . . , a2n , 0,
a2n+2 = a2n+4 = . . . = 0
et si l est impair, ı.e., l = 2n + 1
a3 , 0, a5 , 0, . . . , a2n+1 , 0,
a2n+3 = a2n+5 = . . . = 0
Donc pour l = 2n, le polynôme y2 (x) est nul et y1 (x) devient fini pour tout x à
l’intérieur de l’intervalle [−1, 1] et on a l’inverse pour l = 2n + 1 où le polynôme
y1 (x) est nul tandis que le polynôme y2 (x) devient fini pour tout x à l’intérieur de
l’intervalle [−1, 1]. La solution de l’équation de Legendre dans l’intervalle [−1, 1]
est donc donnée soit par y1 (x) (dans le cas où l est pair) ou par y2 (x) (si l est impair).
y(x) peut être écrit sous la forme générale
(
a0 si l est pair
y(x) = al x + al−2 x + al−4 x + . . . +
l l−2 l−4
a1 si l est impair
[l/2]
X
= al−2r xl−2r (5.4)
r=0
[l/2] est la fonction de plafond (”ceiling function”) et elle définie par

" #  l
l 
 2 si l est pair
=

2 l−1
si l est impair


2
al−2r peut être déduit de la forme générale de an . On a
(n + 2)(n + 1)
an = −an+2
(l − n)(l + n + 1)
En remplaçant n par l − 2 dans l’equation précédente, on trouve
l(l − 1)
al−2 = −al
2(2l − 1)
pour n = l − 4 on obtient
(l − 2)(l − 3)
al−4 = al−2
4(2l − 3)
(l − 1)(l − 2)(l − 3)
= al
2.4(2l − 1)(2l − 3)
A partir des équations de al−2 et al−4 , on peut déduire la forme générale

l(l − 1)(l − 2) . . . (l − 2r + 1)
al−2r = (−1)r
2.4 . . . (2r)(2l − 1)(2l − 3) . . . (2l − 2r + 1)
En remplaçant al−2r dans l’Eq. (5.4), on obtient
[l/2]
X l(l − 1)(l − 2) . . . (l − 2r + 1)
y(x) = al (−1)r xl−2r
r=0
2.4 . . . (2r)(2l − 1)(2l − 3) . . . (2l − 2r + 1)
On peut simplifier l’expression de y(x) en écrivant autrement le numérateur et le

dénominateur
(l − 2r)!
l(l − 1) . . . (l − 2r + 1) = l(l − 1) . . . (l − 2r + 1)
(l − 2r)!
l!
= ,
(l − 2r)!
2.4.6 . . . 2r = (2.1)(2.2)(2.3) . . . (2.r) = 2r r!
et
2l(2l − 1)(2l − 2)(2l − 3) . . . (2l − 2r + 1)(2l − 2r)!
(2l − 1)(2l − 3) . . . (2l − 2r + 1) =
2l(2l − 2) . . . (2l − 2r + 2)(2l − 2r)!
(2l)!
= r
2 l(l − 1) . . . (l − r + 1)(2l − 2r)!
(2l)!(l − r)!
= r
2 l!(2l − 2r)!
donc
[l/2]
X l! 1 2r l!(2l − 2r)! l−2r
y(x) = al (−1)r x
r=0
(l − 2r)! 2r r! (2l)!(l − r)!
[l/2]
X (l!)2 (2l − 2r)! l−2r
= al (−1)r x
r=0
r!(l − 2r)!(2l)!
Choisissant al d’être
(2l)!
al =
2l (l!)2
On obtient le résultat final pour la solution de l’équation de Legendre Eq. (5.1)
[l/2]
X (2l − 2r)!
Pl (x) = (−1)r xl−2r −1≤ x≤1
r=0
2l r!(l − r)!(l − 2r)!
71
Les cinq premiers polynômes sont
P0 (x) = 1
P1 (x) = x
1
P2 (x) = (3x2 − 1)
2
1
P3 (x) = (5x3 − 3x)
2
1
P4 (x) = (35x4 − 30x2 + 3)
8
et ils sont représentés dans la Fig. (5.1).
0.5
Pl (x)
0
P0
P1
-0.5 P2
P3
P4
-1
-1.5 -1 -0.5 0 0.5 1 1.5
x
Figure 5.1: Polynômes de Legendre pour n = 0, 1, 2, 3, 4.
5.1.1 Fonction génératrice
∞
1 X
√ = tl Pl (x) if |t| < 1, |x| ≤ 1 (5.5)
1 − 2tx + t 2
l=0
Démonstration
On a
1 1
√ = [1 − t(2x − t)]− 2
1 − 2tx + t 2
!
1 (−1/2)(−3/2)
= 1 + − [−t(2x − t)] + [−t(2x − t)]2
2 2!
(−1/2)(−3/2) . . . (−(2r − 1)/2)
+ ... + [−t(2x − t)]r + . . .
r!
∞
X 1.3.5 . . . (2r − 1)
= r r!
(−1)r tr (2x − t)r
r=0
2
∞
X (2r)! r
= 2r 2
t (2x − t)r
r=0
2 (r!)
D’aprés la formule du binôme de Newton, on a

r
X r!
(2x − t)r = (2x)r−s (−t) s
s=0
s!(r − s)!
En remplaçant dans l’équation précédente, on trouve

∞ r
1 X (2r)! X r!
√ = 2r 2
(−1) s (2x)r−s tr+s
1 − 2tx + t 2
r=0
2 (r!) s=0
s!(r − s)!
Si on pose n = r + s, on a n qui varie entre 0 et ∞, et puisque s varie entre 0 et r,

donc r varie entre [n/2] et n, alors
∞ X n
1 X (2r)!
√ = tn (−1)n−r n x2r−n
1 − 2tx + t 2
n=0 r=[n/2]
2 r!(n − r)!(2r − n)!
En faisant le changement de variable l = n − r, on obtient finalement
∞ 0
1 X X (2n − 2l)!
√ = tn (−1)l n xn−2l
1 − 2tx + t 2
n=0 l=[n/2]
2 l!(n − l)!(n − 2l)!
X∞
= tn Pn (x)
n=0
5.1.2 Formule de Rodrigues
1 dl 2
Pl (x) = (x − 1)l (5.6)
2l l! dxl
73
Démonstration
D’après la formule du binôme de Newton, on a
l
l X l!
x2 − 1 = (−1)r x2(l−r) (5.7)
r=0
r!(l − r)!
donc
l
1 dl 2 1 X l! r d
l
(x − 1)l
= (−1) x2(l−r)
2l l! dxl 2l l! r=0 r!(l − r)! dxl
mais
dl 2(l−r)
x =0 si 2(l − r) < l, i.e., si r > l/2
dxl
Pl P[l/2]
On peut donc remplacer r=0 par r=0 .
Si r ≤ l/2, on a
dl 2(l−r)
x = {2(l − r)}{2(l − r) − 1}{2(l − r) − 2} . . . {2(l − r) − l + 1}x2(l−r)−l
dxl
(2l − 2r)! l−2r
= x
(l − 2r)!
et on obtient alors
[l/2]
1 dl 2 X (2l − 2r)!
l l
(x − 1)l
= (−1)r l xl−2r
2 l! dx r=0
2 r!(l − r)!(l − 2r)!
= Pl (x)
5.1.3 Représentation intégrale de Laplace
Z π l
1 p
Pl (x) = x+ x2 − 1 cos θ dθ (5.8)
π 0
Démonstration
Considérons la relation (qui peut être prouvé en utilisant le changement de variable
t = tan θ)
π
π
Z
dθ
= √ (5.9)
0 1 + λ cos θ 1 − λ2
√
Posons λ = −u x2 − 1/(1 − ux), En développant les deux côtés de l’Eq. (5.9) en
puissance de u, on obtient
1 1
= √
1 + λ cos θ 1 + u x2 −1
cos θ
1−ux
p −1
= (1 − ux) 1 − u x + x2 − 1 cos θ
∞
X p l
= (1 − ux) ul x + x2 − 1 cos θ
l=0
où on autilisé le théorèrme de binôme de Newton (1 − a)−1 =

P∞ n.
n=0 a
1 1
√ = q
1−λ 2 2 (x2 −1)
1 − u(1−ux) 2
1 − ux
= p
(1 − ux)2 − u2 (x2 − 1)
1 − ux
= √
1 − 2ux + u2
Substitution dans l’Eq. (5.9) donne
Z πX ∞ p l π
ul x + x2 − 1 cos θ dθ = √
0 l=0 1 − 2ux + u2
On insère la formule de fonction génératrice (5.9) dans l’équation précédente, on
obtient
X∞ Z π p l X∞
ul
x + x2 − 1 cos θ dθ = π ul Pl (x)
l=0 0 l=0
Par identification des coefficients de ul on obtient

Z π p l
πPl (x) = x + x2 − 1 cos θ dθ
0
5.1.4 Propriétés des polynômes de Legendre
(a) Pl (1) = 1,
(b) Pl (−1) = (−1)l ,
1
(c) P0l (1) = l(l + 1), (5.10)
2
1
(d) Pl (−1) = (−1)l−1 l(l + 1),
0
2
2l
(e) P2l (0) = (−1)l 2l 2 ,
2 (l!)
( f ) P2l+1 (0) = 0
75
Démonstration
(a) On met x = 1 dans l’Eq. (5.5), on obtient
∞
1 X
= tl Pl (1)
1 − t l=0
Le développement en série de 1/(1 − t) donne

∞
1 X
= tl
1 − t l=0
En identifiant les deux expressions précédentes, il en résulte que Pl (1) = 1, ∀ l.

(b) la preuve est exactement similaire à (a) avec x = −1.
(c) En posant x = 1 dans l’équation de Legendre (5.1), on obtient
−2P0l (1) + l(l + 1)Pl (1) = 0
et puisque Pl (1) = 1, il en résulte que P0l (1) = l(l + 1)/2.

(d) Même preuve que (c) en posant ici x = −1.
Pour montrer (e) et (f), on pose x = 0 dans l’Eq. (5.6). En faisant le développement
en série, on trouve
∞
1 X
√ = tl Pl (0)
1 + t2 l=0
∞
X (2l)!
= (−1)l 2l 2 t2l
l=0
2 (l!)
En identifiant les coefficients de la même puissance de t, on obtient

(2l)!
P2l (0) = (−1)l
22l (l!)2
P2l+1 (0) = 0
5.1.5 Relation d’orthogonalité


1
 0 si l , m
Z 
Pl (x)Pm (x)dx = 

(5.11)
 2 si l = m

−1 2l+1
Démonstration
Pl (x) et Pm (x) satisfont l’équation de Legendre qui peut être écrite sous la forme
( )
d 2 dPl
(1 − x ) + l(l + 1)Pl = 0 (5.12)
dx dx
( )
d dPm
(1 − x2 ) + m(m + 1)Pm = 0 (5.13)
dx dx
En multipliant (5.12) par Pm (x) et (5.13) par Pl (x) et en faisant la soustraction des
deux équations et en intégrant par rapport à x de −1 à 1, on trouve
Z 1" ( ) ( )#
d 2 dPl d 2 dPm
Pm (1 − x ) − Pl (1 − x ) dx
−1 dx dx dx dx
Z 1
+ {l(l + 1) − m(m + 1)} Pl Pm dx = 0
−1
Maintenant, en remplaçant les deux termes de la première intégrale par

( ) ( )
d 2 dPl d 2 dPl dPm dPl
Pm (1 − x ) = Pm (1 − x ) − (1 − x2 )
dx dx dx dx dx dx
( ) ( )
d dP m d dP m dP l dP m
Pl (1 − x2 ) = Pl (1 − x2 ) − (1 − x2 )
dx dx dx dx dx dx
on obtient
" #1 Z 1
dPl dPm
Pm (1 − x2 ) − Pl (1 − x2 ) + (l − m)(l + m + 1) Pl Pm dx = 0
dx dx −1 −1
ce qui montre que si l , m, on devrait avoir

Z 1
Pl (x)Pm (x)dx = 0
−1
R1
Pour montrer que −1
(Pl )2 dx = 2/(2l+1), on utilise la fonction génératrice Eq. (5.6),
∞ 2
1 X 
=  tl Pl (x)
1 − 2tx + t2  l=0 
∞
X ∞
X
= l
t Pl (x) tm Pm (x)
l=0 m=0
En intégrant les deux côtés par rapport à x, on obtient

∞ Z 1 " #1
X 1
t l+m
Pl (x)Pm (x)dx = − ln 1 + t − 2t
2
l,m=0 −1 2t −1
1
= [ln(1 + t) − ln(1 − t)]
t
∞
X t2l
=2
l=0
2l + 1
La derniére ligne résulte du développement en série de ln(1 + t) et ln(1 − t). En

identifiant les coefficients de même puissance de t, il en résulte que si l = m
Z 1
2
[Pl (x)]2 dx =
−1 2l + 1
77
5.1.6 Séries de Legendre

Si f (x) est un polynôme de degré n, alors il peut être écrit sous la forme
n
X
f (x) = cr Pr (x) (5.14)
r=0
avec les coefficients cr donnés par

!Z 1
1
cr = r + f (x)Pr (x)dx (5.15)
2 −1
Démonstration
Si f (x) est un polynôme de degré n, on peut écrire
f (x) = bn xn + bn−1 xn−1 + . . . + b1 x + b0
D’aprés l’Eq. (5.4), le polynôme de Legendre s’écrit sous la forme
Pn (x) = an xn + an−2 xn−2 + . . .
Si on multiplie la dernière expression par bn /an et en le soustrayant de f (x), on

trouve que la différence est un polynôme de degré (n − 1)
f (x) − cn Pn (x) = gn−1 (x)
où cn = bn /an et gn−1 (x) est un polynôme de degré n − 1. En faisant la même chose
pour gn−1 (x), on peut démontrer facilement que gn−1 (x) peut s’écrire sous la forme
gn−1 (x) = cn−1 Pn−1 (x) + gn−2 (x)
donc
f (x) = cn Pn (x) + cn−1 Pn−1 (x) + gn−2 (x)
On fait la même chose pour gn−2 (x), gn−3 (x) et ainsi de suite, on obtient le résultat
désiré
f (x) = cn Pn (x) + cn−1 Pn−1 (x) + cn−2 Pn−2 (x) + . . . + c0 P0 (x)

Xn
= cr Pr (x)
r=0
Les coefficients cn peuvent être calculer de la manière suivante

Z 1 X n Z 1
f (x)Pl (x)dx = cn Pn (x)Pl (x)dx
−1 r=0 −1
2cl
=
2l + 1
donc
!Z 1
1
cl = l + f (x)Pl (x)dx
2 −1
Si f (x) est un polynôme de degré inférieur à l, alors

Z 1
f (x)Pl (x)dx = 0
−1
Si f (x) est de degré n tel que n < l, d’aprés la relation d’orthogonalité on a

!Z 1
1
cl = n + f (x)Pl (x)dx
2 −1
! n Z
1 X 1
= n+ Pr (x)Pl (x)dx
2 r=0 −1
=0
car r est toujours inférieur à l (r < l).
5.1.7 Relations de récurrence
l+1 l
(a) xPl (x) = Pl+1 (x) + Pl−1 (x)
2l + 1 2l + 1
[ 21X
(l+1)]
(b) Pl (x) =
0
(2l − 4r − 1)Pl−2r−1 (x)
r=0
(c) Pl+1 (x) − P0l−1 (x) = (2l + 1)Pl (x)
0
(5.16)
(d) xP0l (x) − P0l−1 (x) = lPl (x)
(e) P0l (x) − xP0l−1 (x) = lPl−1 (x)
Démonstration:
(a) Pl (x) est un polynôme d’ordre l, donc xPl (x) est un polynôme d’ordre l + 1.
D’aprés (5.14), on a
(
c1 P1 (x) si l est pair
xPl (x) = cl+1 Pl+1 (x) + cl−1 Pl−1 (x) + . . . +
c0 P0 (x) si l est impair
où les coefficients cr sont donnés par

!Z 1
1
cr = r + xPl (x)Pr (x)dx
2 −1
!Z 1
1
= r+ Pl (x) {xPr (x)} dx
2 −1
79
xPr (x) est un polynôme d’ordre r + 1, donc cr = 0 si r + 1 < l, ı.e. si r < l − 1. alors
xPl (x) = cl+1 Pl+1 (x) + cl−1 Pl−1 (x) (5.17)
Pour calculer cl+1 et cl−1 , on pose x = 1 dans l’expression (5.17) et dans sa dérivée,
on obtient
Pl (1) = cl+1 Pl+1 (1) + cl−1 Pl−1 (1)

Pl (1) + P0l (1) = cl+1 P0l+1 (1) + cl−1 P0l−1 (1)
En utilisant les propriétés (5.10)(a) et (c), on obtient les équations suivantes
1 = cl+1 + cl−1
" # " #
1 1 1
1 + l(l + 1) = cl+1 (l + 1)(l + 2) + cl−1 (l − 1)l
2 2 2
qui admettent comme solutions

l+1 l
cl+1 = , cl−1 =
2l + 1 2l + 1
On obtient finalement
l+1 l
xPl (x) = Pl+1 (x) + Pl−1 (x)
2l + 1 2l + 1
(b) P0l (x) est un polynôme de degré l − 1, donc on peut l’écrire en série de Legendre
P0l (x) = cl−1 Pl−1 (x) + cl−3 Pl−3 (x) + . . . +cl−2r+1 Pl−2r+1 (x)
(
c P (x) si l est pair
+ ... + 1 1
c0 P0 (x) si l est impair
avec
!Z 1
1
cn = n + P0l (x)Pn (x)dx
2 −1
En intégrant par parties, on trouve

!( Z 1 )
1
cn = n + 1
[Pl (x)Pn (x)]−1 − 0
Pl (x)Pn (x)dx
2 −1
et puisque P0n (x) est un polynôme de degré n−1 < l, donc l’intégrale dans l’équation
précédente est nulle. En utilisant les propriétés (5.10)(a) et (b) des polynômes de
Legendre on obtient
!
1 h i
cn = n + 1 − (−1)n+l
2
n prends les valeurs l − 1, l − 3, l − 5,. . . , donc n + l prends les valeurs 2l − 1, 2l − 3,

2l − 5,. . . . Puisque toutes les valeurs sont impaires, (−1) s+l = −1, donc
cn = 2n + 1
et
cl−2r−1 = 2l − 4r − 1
Il en résulte que
P0l (x) = (2l − 1)Pl−1 (x) + (2l − 5)Pl−5 (x) + . . . + (2l − 4r − 1)P2l−4r−1 (x)
(
3P1 (x) si l est pair
+ ... +
P0 (x) si l est impair
[ 21X
(l−1)]
= (2l − 4r − 1)P2l−4r−1 (x)
r=0
(c) On applique la relation (b) pour P0l+1 (x) et P0l−1 (x), on obtient
P0l+1 (x) = (2l + 1)Pl (x) + (2l − 3)Pl−2 (x) + (2l − 7)Pl−4 (x) + . . .
P0l−1 (x) = (2l − 3)Pl−2 (x) + (2l − 7)Pl−4 (x) + . . .
la soustraction des deux expressions donne
P0l+1 (x) − P0l−1 (x) = (2l + 1)Pl (x)
(d) En multipliant (a) par 2l + 1 et dérivant par rapport à x, on obtient

n o
(l + 1)P0l+1 (x) − (2l + 1) Pl (x) + xP0l (x) + lP0l−1 (x) = 0
On remplace P0l+1 (x) par (2l + 1)Pl (x) − P0l−1 (x) (Eq. (5.16)(c)), et on réarrange les
termes de l’équation on trouve que
xP0l (x) − P0l−1 (x) = lPl (x)
(e) En multipliant (c) par x et en substituant le résultat pour (2l + 1)xPl (x) dans
(5.16)(d), on obtient
(l + 1)Pl+1 (x) + lPl−1 (x) = xP0l+1 (x) − xP0l−1 (x) (5.18)
Si on écrit (d) avec l remplacé par l + 1, on obtient
xP0l+1 (x) − P0l (x) = (l + 1)Pl+1 (x)
Maintenant, en substituant (l + 1)Pl+1 (x) dans l’Eq. (5.18) et en réarrangeant les

termes, on obtient le résultat désiré
P0l (x) − xP0l−1 (x) = lPl−1 (x)

81
5.2 Polynôme associé de Legendre

L’équation associée de Legendre est donnée par
2 m2
" #
2 d y dy
(1 − x ) 2 − 2x + l(l + 1) − y =0 (5.19)
dx dx 1 − x2
L’Eq. (5.19) admet une solution régulière uniquement dans l’intervalle [−1, 1]
avec 0 ≤ m ≤ l. Pour m = 0, l’Eq. (5.19) se réduit à l’équation de Legendre.
La solution de cette équation noté Pm
l (x) est donnée par
dm
l (x) = (1 − x )
2 m/2
Pm Pl (x) (5.20)
dxm
La solution (5.20) est appellée polynôme associé de Legendre. Pour des valeurs
négatives de m, on utilise
(l − m)! m
l (x) = (−1)
P−m m
P (x) (5.21)
(l + m)! l
où encore 0 ≤ m ≤ l.
Démonstration
Pour montrer que (5.20) est une solution de l’Eq. (5.19), on doit montrer que
d2 m2
" #
d m
(1 − x2 ) 2 Pm (x) − 2x P (x) + l(l + 1) − Pm (x) = 0
dx l dx l 1 − x2 l
Puisque Pl (x) est une solution de l’équation de Legendre, on a

d2 d
1 − x2 2
Pl (x) − 2x Pl (x) + l(l + 1)Pl (x) = 0
dx dx
En dérivant l’équation précédente m fois par rapport a x, elle devient
dm
" 2 #
dm
" #
dm
2 d d
m
1−x P l (x) − 2 m
x Pl (x) + l(l + 1) Pl (x) = 0
dx dx2 dx dx dxm
En utilisant la règle de Leibniz pour le m−ième dérivatives

m
dm X m! dn f dn−m g
{ f g} =
dxm n=0
n!(n − m)! dxn dxn−m
l’équation précédente devient

dm+2 d m+1
2 d m(m − 1) d2 m
2 d
1 − x2 Pl (x) + m 1 − x Pl (x) + 1 − x Pl (x)
dxm+2 dx dxm+1 2 dx2 dxm
" m+1
d dm dm
#
d
− 2 x m+1 Pl (x) + m x m Pl (x) + l(l + 1) m Pl (x) = 0
dx dx dx dx
En regroupant les coefficients de dm+2 Pl /dxm+2 , dm+1 Pl /dxm+1 et dm Pl /dxm , on

obtient
dm+2 dm+1
1 − x2 P l (x) − 2x(m + 1) Pl (x) (5.22)
dxm+2 dxm+1
dm
+ [l(l + 1) − m(m + 1) − 2m] Pl (x) = 0 (5.23)
dxm
l (x) = d Pl (x)/dx
qui devient, on dénote ym m m
d2 d
1 − x2 ym (x) − 2x(m + 1) ym (x) + [l(l + 1) − m(m + 1) − 2m] ym
l (x) = 0
dx2 l dx l
Si on remplace maintenant
2 m/2 m

l (x) = 1 − x
Pm yl (x)
dans l’équation précédente, on obtient

d 2 d
2 −m/2 m 2 −m/2 m

1 − x2 1 − x Pl (x) − 2(m + 1)x 1 − x Pl (x)
dx2 dx
−m/2
+ [l(l + 1) − m(m + 1)] 1 − x2 l (x) = 0
Pm (5.24)
mais
d −m/2 dPm (x) −m/2 −(m/2)−1
1 − x2 Pm
l (x) = l
1 − x2 + mx 1 − x2 Pm
l (x)
dx dx
et
d 2 2 −m/2 m

1 − x P l (x) =
dx2
d2 Pml (x)

2 −m/2

2 −(m/2)−1
dPm
l (x)

2 −(m/2)−1 m

1 − x + mx 1 − x + mx 1 − x Pl (x)
dx2 dx
donc, l’Eq. (5.24) devient

2 −(m/2)+1
d2 Pm
l (x)

2 −m/2
dPm l (x)

2 −m/2 m

1−x + 2mx 1 − x + m 1 − x Pl (x)
dx2 dx
−(m/2)−1
+ m(m + 2) 1 − x2 x 2 Pm
l (x)
" m
dPl (x)
#
2 −m/2 2 −(m/2)−1 m

− 2(m + 1)x 1 − x + mx 1 − x Pl (x)
dx
−m/2
+ [l(l + 1) − m(m + 1)] 1 − x2 l (x) = 0
Pm
83
−m/2
On annule le facteur commun 1 − x2 , on obtient
d2 Pm
l (x) dPm
l (x)
1−x 2
2
+ [2mx − 2(m + 1)x]
dx dx
m(m + 1) 2(m + 1)mx2
" #
+ m+ − + l(l + 1) − m(m + 1) Pm
l (x) = 0
1 − x2 1 − x2
qui peut être simplifié à
d2 Pm
l (x) dPm
l (x)
"
m2
#
1−x 2
− 2x + l(l + 1) − Pm (x) = 0
dx2 dx 1 − x2 l
qui est l’équation désiré.

Les polynômes associé de Legendre Pml (x) pour l = 5 et 0 < m < l sont représentés
dans la Fig. (5.2). Les premiers polynômes associés de Legendre, y compris ceux
Plm
1 P50
P51
P52
0 P53
P54
-1 P55
-2
-1 -0.5 0 0.5 1
x
l (x) pour l = 5 et 0 < m < l.

Figure 5.2: Polynômes associés de Legendre Pm
des valeurs négatives de m, sont les suivants

P00 (x) = 1
p
P11 (x) = − 1 − x2
P01 (x) = x
1 1 1p
P−1
1 = − P (x) = 1 − x2
2 1 2
P22 (x) = 3(1 − x2 )
p
P12 (x) = −3x 1 − x2
1
P02 (x) = (3x2 − 1)
2
1 1 1 p
P−1
2 (x) = − P (x) = x 1 − x2
6 2 2
1 2 1
2 (x) =
P−2 P2 (x) = (1 − x2 )
24 8
P3 (x) = −15(1 − x )
3 2 3/2
P23 (x) = 15x(1 − x2 )

3 p
P13 (x) = − (5x2 − 1) 1 − x2
2
1 3
P3 (x) = (5x − 3x)
0
2
1 1 1 2 p
P−1
3 (x) = − P (x) = − (5x − 1) 1 − x2
12 3 8
1 2 1
3 (x) =
P−2 P (x) = x(1 − x2 )
120 3 8
1 1
3 (x) = −
P−3 P33 (x) = (1 − x2 )3/2
720 48
P44 (x) = 105(1 − x2 )2
P34 (x) = −105x(1 − x2 )3/2
15
P24 (x) = (7x2 − 1)(1 − x2 )
2
5 p
P4 (x) = − (7x3 − 3x) 1 − x2
1
2
1
P04 (x) = (35x4 − 30x2 + 3)
8
1 1 p
P4 (x) = − P14 (x) = (7x3 − 3x) 1 − x2
−1
20 8
1 1
4 (x) =
P−2 P24 (x) = (7x2 − 1)(1 − x2 )
360 48
1 3 x
P4 (x) = −
−3
P (x) = (1 − x2 )3/2
5040 4 48
1 1
4 (x) =
P−4 P44 (x) = (1 − x2 )2
40320 384
85
5.2.1 Relation d’orthogonalité
1
2(l + m)!
Z
l (x)Pl0 (x)dx =
Pm δl,l0
m
(5.25)
−1 (2l + 1)(l − m)
Démonstration
Pour démontrer que si l , l0 , l’intégral donne 0, nous suivons les mêmes étapes que
celle pour la relation d’orthogonalité de polynôme de Legendre Eq. (5.11). Pm l (x)
et Pm
l 0 (x) satisfy l’équation de Legendre associée qui peut être écrire sous la form
dPm m2
( ) " #
d 2
(1 − x ) l
+ l(l + 1) − Pm = 0 (5.26)
dx dx 1 − x2 l
dPm0 m2
( ) " #
d
(1 − x2 ) l + l0 (l0 + 1) − Pm0 = 0 (5.27)
dx dx 1 − x2 l
En multipliant Eq. (5.26) par Pm m
l0 (x) et l’Eq. (5.27) par Pl (x), en fiasant la substrac-
tion des deux équations et intégrant par rapport à x de −1 à 1, on trouve
dPm dPm
Z 1" ( ) ( )#
m d 2 l m d 2 l0
Pl0 (1 − x ) − Pl (1 − x ) dx
−1 dx dx dx dx
1 m m
Z
+ l(l + 1) − l (l + 1)
0 0
Pl Pl0 dx = 0

−1
En remplaçant les deux termes du première intégrale par

dPm dPm dPm dPm
( ) ( )
m d d l0
Pl0 2
(1 − x ) l
= m 2
P 0 (1 − x ) l
− (1 − x2 ) l
dx dx dx l dx dx dx
dPm dPm dPm dPm
( ) ( )
m d l0 d l0 l0
Pl 2
(1 − x ) = m
P (1 − x )2
− l
(1 − x2 )
dx dx dx l dx dx dx
il s’ensuit
#1
dPm dPm
" Z 1
l0
m 2
Pl0 (1 − x ) l m 2
− Pl (1 − x ) + (l − l )(l + l + 1)
0 0
l Pl0 dx = 0
Pm m
dx dx −1 −1
donc
Z 1
(l − l )(l + l + 1)
0 0
l (x)Pl0 (x)dx = 0
Pm m
−1
ce qui montre que si l , l0 , on devrait avoir

Z 1
Pml (x)Pl0 (x)dx = 0
m
−1
Montrons maintenant que si l = l0 , on a

2(l + m)!
Z 1h i2
Pml (x) dx =
−1 (2l + 1)(l − m)!
Supposons tout d’abord que m > 0. D’après l’Eq. (5.20), on a

Z 1h Z 1 " m #" m #
2 m d d
i2
Pl (x) dx =
m
1−x Pl (x) Pl (x) dx
−1 −1 dxm dxm
" m−1 ! m
!#1
d 2 m d

= Pl (x) 1 − x Pl (x)
dxm−1 dxm −1
Z 1 " m−1 # " m d m #
d d 2
− m−1 l
P (x) 1−x Pl (x) dx
−1 dx dx dxm
Z 1 " m−1 # " m
#
d d 2 m d

=− m−1 l
P (x) 1−x Pl (x) dx (5.28)
−1 dx dx dxm
où nous avons intégré par parties. A partir de l’Eq. (5.22) avec m est remplacé par
m − 1, on a
dm+1 dm dm−1
1 − x2 P l (x) − 2xm P l (x) + [l(l + 1) − m(m − 1)] Pl (x) = 0
dxm+1 dxm dxm−1
m−1
En multipliant par 1 − x2 , l’équation devient
m dm+1 m
2 m−1 d

1 − x2 P l (x) − 2xm 1 − x Pl (x)
dxm+1 dxm
m−1 dm−1
+ [l(l + 1) − m(m − 1)] 1 − x2 Pl (x) = 0
dxm−1
qui peut être écrit comme
m m−1
" #
d 2 m d

2 m−1 d

1−x P l (x) = −(l + m)(l − m + 1) 1 − x Pl (x)
dx dxm dxm−1
On insère ce résultat dans l’Eq. (5.28), on obtient
Z 1h i2
Pml (x) dx
−1
Z 1 " m−1 # " m−1 #
d
2 m−1 d

= P (x) (l + m)(l − m + 1) 1 − x
m−1 l
Pl (x) dx
−1 dx dxm−1
Z 1 " m−1 #2
2 m−1 d

= (l + m)(l − m + 1) 1−x Pl (x) dx
−1 dxm−1
Z 1h i2
= (l + m)(l − m + 1) Pm−1
l (x) dx
−1
Dans la dernière ligne, nous avons utilisé la définition de polynône de Legendre
associé Eq. (5.20). En appliquant à nouveau ce résultat, on trouve
Z 1h i2 Z 1h i2
Pl (x) dx = (l + m)(l − m + 1)(l + m − 1)(l − m + 2)
m
Pm−2
l (x) dx
−1 −1
Z 1h i2
= (l + m)(l + m − 1)(l − m + 1)(l − m + 2) Pm−2
l (x) dx
−1
87
En répétant le processus m fois, on trouve

Z 1h i2
Pm
l (x) dx
−1
Z 1 h i2
= (l + m)(l + m − 1) . . . (l + 1)(l − m + 1)(l − m + 2) . . . l P0l (x) dx
−1
2
= (l + m)(l + m − 1) . . . (l + 1)l(l − 1) . . . (l − m + 2)(l − m + 1)
2l + 1
(l + m)! 2
=
(l − m)! 2l + 1
qui est le résultat désiré. Supposons maintenant que m < 0, ı.e., m = −n avec
n > 0, donc
Z 1h i2 Z 1h i2
Pl (x) dx =
m
P−n
l (x) dx
−1 −1
Z 1" #2
(l − n)! n
= n
(−1) P P (x) dx
−1 (l + n)! l
" #2 Z 1 h
(l − n)! i2
= P Pnl (x) dx
(l + n)! −1
#2
(l − n)! (l + n)! 2
"
= P
(l + n)! (l − n)! 2l + 1
(l − n)! 2
=
(l + n)! 2l + 1
(l + m)! 2
=
(l − m)! 2l + 1
5.3 Harmoniques sphériques

On appelle harmoniques spheriques, les solutions de l’équation de type
∂ ∂Ψ 1 ∂2 Ψ
!
1
sin θ + 2 + l(l + 1)Ψ = 0 (5.29)
sin θ ∂θ ∂θ sin θ ∂φ2
Les harmoniques sphériques sont particulièrement utiles pour les problémes qui ont
une symétrie de rotation, car elles sont les vecteurs propres de certains opérateurs
liés aux rotations. En physique, elles sont souvent utilisées en acoustique, en
cristallographie, en mechanique quantique, en cosmologie, etc. Cette équation
peut être resolue par la méthode de séparation des variables, ı.e., on cherche les
solutions sous la forme Ψ(θ, φ) = Θ(θ)Φ(φ). En remplaçant dans l’Eq. (5.29), cela
nous donne
Φ Θ d2 Φ
( !)
d dΘ
sin θ + 2 + l(l + 1)ΘΦ = 0
sin θ dθ dθ sin θ dφ2

Fonctions Speciales Et Polynomes Orthogonaux Cour-75-96

Transféré par

Droits d'auteur :

Formats disponibles

Fonctions Speciales Et Polynomes Orthogonaux Cour-75-96

Transféré par

Informations du document

Titre original

Copyright

Formats disponibles

Partager ce document

Partager ou intégrer le document

Options de partage

Avez-vous trouvé ce document utile ?

Ce contenu est-il inapproprié ?

Droits d'auteur :

Formats disponibles

Fonctions Speciales Et Polynomes Orthogonaux Cour-75-96

Transféré par

Droits d'auteur :

Formats disponibles

Chapitre 5

Les polynômes orthogonaux

5.1 Polynômes de Legendre

donc, le domaine de convergence de la solution de l’Eq. (5.3) est la même que le

d2 y 2x2 dy l(l + 1)x2

Si en insère y, sa dérivée y0 et la dérivée seconde y00 dans l’Eq. (5.1), on trouve

et si l est impair, ı.e., l = 2n + 1

[l/2] est la fonction de plafond (”ceiling function”) et elle définie par

al−2r peut être déduit de la forme générale de an . On a

En remplaçant n par l − 2 dans l’equation précédente, on trouve

A partir des équations de al−2 et al−4 , on peut déduire la forme générale

On peut simplifier l’expression de y(x) en écrivant autrement le numérateur et le

2.4.6 . . . 2r = (2.1)(2.2)(2.3) . . . (2.r) = 2r r!

Les cinq premiers polynômes sont

et ils sont représentés dans la Fig. (5.1).

Figure 5.1: Polynômes de Legendre pour n = 0, 1, 2, 3, 4.

5.1.1 Fonction génératrice

D’aprés la formule du binôme de Newton, on a

En remplaçant dans l’équation précédente, on trouve

Si on pose n = r + s, on a n qui varie entre 0 et ∞, et puisque s varie entre 0 et r,

En faisant le changement de variable l = n − r, on obtient finalement

5.1.2 Formule de Rodrigues

5.1.3 Représentation intégrale de Laplace

où on autilisé le théorèrme de binôme de Newton (1 − a)−1 =

Par identification des coefficients de ul on obtient

5.1.4 Propriétés des polynômes de Legendre

Le développement en série de 1/(1 − t) donne

En identifiant les deux expressions précédentes, il en résulte que Pl (1) = 1, ∀ l.

−2P0l (1) + l(l + 1)Pl (1) = 0

et puisque Pl (1) = 1, il en résulte que P0l (1) = l(l + 1)/2.

En identifiant les coefficients de la même puissance de t, on obtient

5.1.5 Relation d’orthogonalité

Maintenant, en remplaçant les deux termes de la première intégrale par

ce qui montre que si l , m, on devrait avoir

En intégrant les deux côtés par rapport à x, on obtient

La derniére ligne résulte du développement en série de ln(1 + t) et ln(1 − t). En

5.1.6 Séries de Legendre

avec les coefficients cr donnés par

f (x) = bn xn + bn−1 xn−1 + . . . + b1 x + b0

D’aprés l’Eq. (5.4), le polynôme de Legendre s’écrit sous la forme

Pn (x) = an xn + an−2 xn−2 + . . .

Si on multiplie la dernière expression par bn /an et en le soustrayant de f (x), on

f (x) − cn Pn (x) = gn−1 (x)

gn−1 (x) = cn−1 Pn−1 (x) + gn−2 (x)

f (x) = cn Pn (x) + cn−1 Pn−1 (x) + gn−2 (x)

f (x) = cn Pn (x) + cn−1 Pn−1 (x) + cn−2 Pn−2 (x) + . . . + c0 P0 (x)

Les coefficients cn peuvent être calculer de la manière suivante

Si f (x) est un polynôme de degré inférieur à l, alors

Si f (x) est de degré n tel que n < l, d’aprés la relation d’orthogonalité on a

car r est toujours inférieur à l (r < l).

5.1.7 Relations de récurrence

où les coefficients cr sont donnés par

xPl (x) = cl+1 Pl+1 (x) + cl−1 Pl−1 (x) (5.17)

Pl (1) = cl+1 Pl+1 (1) + cl−1 Pl−1 (1)

En utilisant les propriétés (5.10)(a) et (c), on obtient les équations suivantes

qui admettent comme solutions

En intégrant par parties, on trouve

n prends les valeurs l − 1, l − 3, l − 5,. . . , donc n + l prends les valeurs 2l − 1, 2l − 3,

la soustraction des deux expressions donne