École polytechnique de l’Université de Nantes

Systèmes électroniques et informatique industrielle

Première année
2006-2007

Analyse fonctionnelle

III. Approximations hilbertiennes

et développements en série

Laurent Guillopé
 TABLE DES MATIÈRES

Liste des figures. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . iii

Prologue. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . iv
1. Espaces de Hilbert. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 1
1.1. Produit scalaire. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 1
1.2. Orthogonalité. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 3
1.3. Forme linéaire. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 4
1.4. Base hilbertienne. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 5
1.5. Exercices. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 6
2. Systèmes orthogonaux. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 9
2.1. Exponentielles de Fourier. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 9
2.2. Polynômes de Legendre. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 13
2.3. Polynômes de Chebyshev. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 18
2.4. Fonctions d’Hermite. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 20
2.5. Système de Haar. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 22
2.6. Exercices. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 25
3. Fonctions de Bessel. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 28
3.1. Modélisation dans un domaine circulaire. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 28
3.2. Les fonctions de Bessel. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 29
3.3. Développements en série de Fourier-Bessel. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 32
3.4. Exercices. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 34
A. La fonction Gamma d’Euler Γ . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 37
Index. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 38
 LISTE DES FIGURES

1 Géométrie euclidienne dans un espace de Hilbert : identité du parallélogramme,

relation de Pythagore, projection d’un point sur un convexe fermé.. . . . . . . . . . . . . . . . . 4

2 Les fonctions cos/sinusoïdales En sur [0, 1] pour n = 0, . . . , 20.. . . . . . . . . . . . . . . . . . . . 10

3 Approximations de la fonction g1 définie sur [0, 1] par P g1 (t) = t si t < 1/2 et g1 (t) =
−1 k+1 −1
t−1 si t ≥ 1/2 par les sommes partielles de Fourier π 1≤2k−1≤N (−1) k sin k2πt
pour N = 2, 6, 10, 14, 18, 22, 50, 100.. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 11
4 Les polynômes de Legendre Ln sur [−1, 1] pour n = 0, . . . , 20.. . . . . . . . . . . . . . . . . . . . . 13
5 Approximations de la fonction f1 définie sur [−1, 1] par f1 (t) = t si |t| < 1/2 et
PN
f1 (t) = t − sgn(t) si |t| ≥ 1/2 par la somme de polynômes de Legendre n=0 (n +
1/2)`n (f1 )Ln pour N = 2, 6, 10, 14, 18, 22, 50, 100.. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 16
6 Les polynômes de Chebyshev Tn sur [−1, 1] pour n = 0, . . . , 20.. . . . . . . . . . . . . . . . . . . 18
7 Approximations de la fonction f1 définie sur [−1, 1] par f1 (t) = t si |t| < 1/2
et f1 (t) = t − sgn(t) si |t| ≥ 1/2 par la somme de polynômes de Chebyshev
PN
n=0 tn (f1 )Tn pour N = 2, 6, 10, 14, 18, 22, 50, 100.. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 20
8 Les fonctions d’Hermite hn pour n = 0, . . . , 19.. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 21
9 Approximations sur R de la fonction f1 à support [−1, 1] et telle que f1 (t) = t si
|t| < 1/2 et f1 (t) = t − sgn(t) si |t| ≥ 1/2 par la somme de polynômes d’Hermite
PN −2
n=0 kHn k hfT , Hn iHn pour N = 2, 6, 10, 14, 22, 50, 100, 200, 400, 700, 1000.. . . . . . 23
10 Les ondelettes de Haar ϕ, ψ , ainsi que ϕ`,m , ψ`,m .. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 24

11 Les fonctions de Bessel Jn , n = 0, 1, 2, 3, 4 sur l’intervalle (0, 25).. . . . . . . . . . . . . . . . . . . 29

12 Les fonctions de Bessel Yn , n = 0, 1, 2, 3, 4 sur l’intervalle (0, 25).. . . . . . . . . . . . . . . . . . . 32
13 Les fonctions de Bessel J0 (x0n t), t ∈ [0, 1] pour n = 0, . . . , 20.. . . . . . . . . . . . . . . . . . . . . . 33
14 Approximations de la fonction f1 définie sur [0, 1] par f1 (t) = t si t < 1/2 et
f1 (t) = t − 1 sinon, par la somme de fonctions de Bessel (46) tronquée à l’ordre N
pour N = 2, 6, 10, 14, 18, 22, 30, 50.. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 34
 PROLOGUE

Le titre général de ces notes insiste sur le point de vue de l’approximation en moindres
carrés d’une fonction, via sa représentation par une série dans un espace de Hilbert
N
X
f = lim hf, en ien
N →∞
n=0
avec
N 2 ∞
X X
f− hf, en ien = |hf, en i|2 →N →∞ 0.
n=0 n=N +1
La famille (en )n∈N , un système de fonctions orthonormées, sera tour à tour celle des expo-
nentielles de Fourier, et leurs compagnons sin et cos, celle de polynômes de divers types
(Legendre, Chebyshev, Hermite), celle des fonctions de Haar (ancêtres des ondelettes) et
celle de fonctions de Bessel.
Le cadre géométrique des séries précédentes est celui de la géométrie euclidienne (celle
du plan ou de l’espace physique) en dimension infinie (le domaine de l’analyse), qui mêle
formules de Pythagore de la géométrie élémentaire et complétude de l’analyse fonctionnelle :
c’est la théorie des espaces de Hilbert dont les résultats de base sont exposés dans le premier
chapitre.
Les deux chapitres suivants se concentrent sur les familles particulières de fonctions don-
nant des bases orthonormées. Leur caractère de catalogue aurait pu imposer le titre général
de Fonctions spéciales au lieu de la référence hilbertienne (qui ne couvre pas les résultats
plus subtils de convergence simple) : pourquoi ces choix, et pas d’autres, tant les fonctions
spéciales sont nombreuses ? Sans risque d’erreur, on peut dire que ce sont les fonctions qui
apparaissent le plus fréquemment dans les sciences de l’ingénieur, après les fonctions ex-
ponentielle, logarithme, sinusoïdales : par ex., les fonctions de Bessel apparaissent lors de
l’étude des fonctions radiales (pour le calcul de leur transformée de Fourier, pour l’expression
du Laplacien en coordonnées polaires et la séparation des variables exposée au début du cha-
pitre 3), les polynômes de Legendre sont aussi intimement liés à l’analyse de fonctions dans
l’espace R3 comme par exemple la formule (20) pour le dipôle ou les harmoniques sphériques
de la sphère S 2 à l’instar des fonctions sinusoïdales sur le cercle S 1 du plan.
L’étude de, et les résultats sur, ces fonctions spéciales sont emmêlés, comme l’exprime par
exemple parfaitement la décomposition (21) attribuée à Rayleigh de l’onde plane ejhk,ri de
vecteur k ∈ R3
r
π Jn+1/2 (kkk krk)
 
jhk,ri
X
n hk, ri
e = (2n + 1)j p Ln , r ∈ R3 .
n≥0
2 kkk krk kkk krk
en série de polynômes de Legendre Ln , avec les amplitudes exprimés en terme de fonctions
de Bessel Jn+1/2 d’ordre demi-entier.
 PROLOGUE v

L’examen des graphes des diverses familles (Fig. 2, 4, 6, 8, 13), et des approximations
des fonctions en dent de scie classique (Fig. 3, 5, 7, 9, 14) montre par ailleurs toute la
similarité des approximations. Si les démonstrations pour les exponentielles sont classiques
(et aisées), les résultats sont repris quasiment mot pour mot pour les familles de polynômes
ou de fonctions de Bessel : ces notes auraient pu être intitulées Développements de Fourier-
Legendre-Chebyshev-Hermite-Bessel ! Aussi simples soient-ils, ces résultats requièrent pour
leur preuve des résultats fins de la théorie moderne, i. e. lebesguienne, de l’intégration (ne
serait-ce que la définition de l’espace L2 (I) associé à un intervalle I de R) : le lecteur
intéressé est renvoyé au cours Linéarité et convergences qui indique au long d’un bref aperçu
de cette théorie quelques résultats utilisés sans vergogne ici.
Ces fonctions classiques ont été introduites et étudiées il y a bien longtemps : Bessel, Che-
byshev, Hermite, Legendre sont des mathématiciens du XIXe siècle (à l’habitude, l’Index
de la version en ligne renvoie pour les mathématiciens cités dans le texte à leurs notices de
l’encyclopédie biographique MacTutor history of mathematics archive de J. J. O’Connor et
E. F. Robertson de l’Université de St Andrews, Écosse). Avec la disponibilité des ordina-
teurs et de leurs bibliothèques de programmes scientifiques (C, maple, matlab, scilab,. . .),
le scientifique ou l’ingénieur du XXIe siècle ont un accès commun et aisé à ces fonctions
classiques.

Nantes, le 7 janvier 2007

Laurent Guillopé

[Link]/˜guillope/seii1-af/
[Link]@[Link]
 CHAPITRE 1

ESPACES DE HILBERT

Un espace de Hilbert est un espace normé complet, non nécessairement de dimension finie,
dont la norme dérive d’un produit scalaire comme la norme euclidienne || ||2 de Rn , ce
qui lui confère une géométrie euclidienne qui permet de généraliser simplement certaines
propriétés de la dimension finie.
L’archétype des espaces de Hilbert (non de dimension finie) est l’espace des suites dénom-
brables (i. e. indexée par N ou, par une bijection convenable, par Z)
X
`2 (N) = {u = (un )n∈N ∈ CN , |un |2 < ∞}
n∈N

avec produit scalaire défini par

X
(1) hu, vi = un v n , u, v ∈ `2 (N)
n∈N

et norme || ||2 par

sX
kuk2 = |un |2 , u ∈ `2 (N).
n∈N

Si (ek )k∈N est la famille de suites où ek est la suite de `2 (N) dont le seul coefficient non nul
2
est le k -ème valant 1, on peut écrire tout
P u = (un )n∈N de2 ` (N) P avec un nombre fini de
coordonnées un non nulles comme u = k∈N uk ek , avec kuk2 = k∈N |uk |2 = k∈N hu, ek i2 .
P
De telles sommes gardent un sens comme série pour des vecteurs quelconques u ∈ `2 (N),
justifiant le qualificatif de base orthonormée pour la famille (en )n∈N . C’est ce que développe
ce chapitre.

1.1. Produit scalaire

Définition 1.1. Un produit scalaire sur l’espace vectoriel E complexe(1) est la donnée d’une
application de E 2 dans C dont l’image de (u, v) ∈ E 2 dans C est notée hu, vi vérifiant les
propriétés
– pour v fixé, l’application u → hu, vi est linéaire,
– hu, vi = hv, ui pour tout u, v ∈ E ,
– hu, ui > 0 pour tout vecteur u non nul.
. Exemple 1.1. L’espace C([0, T ]) est muni du produit scalaire
1 T
Z
hf, gi =
T 0
f (t)g(t)dt, f, g ∈ C([0, T ]). /
(1)
Si l’espace E est réel, le produit scalaire est une fonction à valeurs réelles. Pour un espace de Hilbert
complexe, on précise souvent produit scalaire hermitien.
2 CHAPITRE 1. ESPACES DE HILBERT

p E espace vectoriel muni du produit scalaire h , i. La fonction || ||2

Proposition 1.1. Soit
définie par kuk2 = hu, ui, u ∈ E est une norme sur E , la norme dérivant du produit
scalaire h , i.
Démonstration. — Le polynôme du second degré en λ ∈ R
Tu,v (λ) = kλu + vk22 = hλu + v, λu + vi = λ2 hu, ui + λ(hu, vi + hv, ui) + hv, vi
= kuk22 λ2 + 2<e hu, viλ + kvk22
est toujours positif, ainsi son discriminant réduit
∆ = (<e hu, vi)2 − kuk22 kvk22
est négatif, soit
(2) |<e hu, vi)| ≤ kuk2 kvk2 ,
inégalité connue sous le nom de Cauchy-Schwarz(-Bunyakovsky). Par suite
ku + vk22 = kuk22 + 2<e hu, vi + kvk22 ≤ (kuk2 + kvk2 )2 ,
i. e. l’inégalité triangulaire exigée comme une des propriétés de la fonction norme, les autres
propriétés (homogénéité, positivité et caractérisation du vecteur nul comme seul vecteur de
norme nulle) étant vérifiées aisément.
4 Remarque 1.1. Pour x et y vecteurs quelconques de E , on a l’égalité kx + yk2 + kx −
yk2 = 2kxk2 + 2kyk2 , dite du parallélogramme. 5
Définition 1.2. Un espace de Hilbert est un espace normé (E, || ||) dont la norme || ||
dérive d’un produit scalaire sur E et qui est complet relativement à cette norme.
4 Remarque 1.2. Un espace normé (E, || ||) est dit complet P si toute suite de Cauchy y
est convergente.
P C’est équivalent au fait que toute série k∈N vk absolument convergente
(i. e. k∈N kvk k) est convergente 5
. Exemples 1.2.
1. L’espace `2 (N) est un espace de Hilbert.
Complétude de `2 (N) . — Soit (uk )k∈N une suite de Cauchy dans `2 (N) . On a uk = (ukn )n≥0 où ukn est un
complexe. Pour tout p, q , on a |upn − uqn | ≤ kup − uq k2 , ainsi, à n fixé la suite de complexes (ukn )k∈N est de
Cauchy : soit u∞n sa limite et u_inf ty la suite u∞ = (u∞k )n∈N . Soit ε > 0 et K tel que si k, ` > K on a
kuk − u` k2 ≤ ε . On a donc pour tout N ,
N
X
|upn − uqn |2 ≤ kuk − u` k2 ≤ ε2 .
n=0
PN 2
En passant à la limite lorsque ` → ∞ , on obtient n=0 |ukn − u∞n | ≤ ε2 , puis en faisant N → ∞ ,
P∞ 2 2
k=0 |ukn − u∞n | ≤ ε . On en déduit que la suite u∞ est dans `2 (N) , puis que kuk − u∞ k ≤ ε . On vient
donc de montrer que u∞ = limk→∞ uk , d’où la complétude annoncée.

2. L’espace `f (N), muni du produit scalaire (1) (où la somme est toujours finie) n’est pas
complet : en effet la suite (xk )k≥1 telle que xkn = 1/(n + 1) si n ≤ k et zéro sinon est
une suite de Cauchy, qui n’est pas convergente dans `f (N).
3. L’espace (C([0, 1]), || ||2 ) de l’exemple 1.1 n’est pas complet. La suite (gn )n∈N définie
par 
0,
 si 0 ≤ t ≤ 1/2,
gn (t) = n(t − 1/2), si 1/2 ≤ t ≤ 1/2 + 1/n,

1, si 1/2 + 1/n ≤ t < 1,
est de Cauchy sans être convergente dans cet espace normé.
 1.2. ORTHOGONALITÉ 3

RT
4. L’espace L2 (0, T ) des fonctions f sur (0, T ) de caré intégrable, i. e. 0 |f (t)|2 dt < ∞,
muni du produit scalaire
1 T
Z
hg1 , g2 iT = g1 (t)g2 (t)dt, g1 , g2 ∈ L2 (0, T )
T 0
est un espace de Hilbert, dont la norme d’un vecteur f , une fonction de carré intégrable,
sera notée kf kL2 (0,T ) ou simplement kf k2 . Cet Hilbert contient C([0, T ]) comme sous-
espace dense. Son sous-espace C([0,
e T ]) des fonctions telles que f (0) = f (T ) et dont le
prolongement T -périodique à R est continu, l’est aussi : il contient le sous-espace des
polynômes trigonométriques comme sous-espace dense. /

1.2. Orthogonalité
Définition 1.3. Deux vecteurs u et v de l’espace E sont dits orthogonaux si hu, vi = 0.
Une famille (ei )i∈I de vecteurs de E est dite orthogonale si tous les vecteurs sont deux à
deux orthogonaux, elle est dite orthonormée si elle est orthogonale, avec tous ses vecteurs
unitaires.
Soit A une partie de E . Son orthogonal, noté A⊥ , est la partie définie par
A⊥ = {u ∈ E, hu, ai = 0, a ∈ A}.
4 Remarque 1.3. Si u et v sont orthogonaux, alors on a la relation de Pythagore
ku + vk2 = kuk2 + kvk2 .
Cette égalité, ainsi que l’égalité du parallélogramme de la Rem. 1.1, ne sont nullement
valables pour une norme en général, comme on le vérifie en particulier pour les normes
|| ||1 et || ||∞ telles que
n
X
kzk1 = |zi |, kzk∞ = sup |zi |, z ∈ Cn .
i=1,...,n
i=1

On vérifie simplement que l’orthogonal A⊥ d’une partie A est un sous-espace vectoriel de

E. 5
Introduisons deux notions nécessaires au théorème de projection illustré dans sa version
plane par la Fig. 1
Définition 1.4. La partie F de l’espace vectoriel normé (E, || ||) est dite fermée si toute
suite convergente d’éléments de E a sa limite dans E .
. Exemples 1.3.
1. Ni Q, ni son complémentaire R \ Q, ne sont fermés dans R.
2. En dimension finie, un sous-espace vectoriel de E est fermé, alors qu’en dimension
infinie rien ne peut être dit a priori : le sous-espace `f (N) de `2 (N) n’est pas fermé,
alors que le sous-espace `N (N) des suites de `2 (N) dont les termes sont tous nuls à
partir du rang N (compris) en est un.
3. La boule (ouverte) B|| || (v, r) = {kx − vk < r|} n’est pas fermée : son complémentaire
E \ B|| || (v, r) est fermé, de même que la boule (fermée) B || || (v, r) = {|w − vk ≤
r, w ∈ E}. /
Définition 1.5. Une partie C de l’espace vectoriel E est dite convexe si pour tous vecteurs
u, v de E , le segment [u, v] = {u + λ(v − u) : λ ∈ [0, 1]} est contenu dans C .
. Exemple 1.4. Une boule (ouverte ou fermée) est convexe, de même que tout sous-espace
linéaire. /
4 CHAPITRE 1. ESPACES DE HILBERT

C
x−y v
u+v w u
x+y
y v
x u

Figure 1 . Géométrie euclidienne dans un espace de Hilbert : identité du parallé-

logramme, relation de Pythagore, projection d’un point sur un convexe fermé.

Théorème ∗ 1.1. Soit C une partie fermée convexe d’un espace de Hilbert E . Alors, pour
tout u ∈ E , il existe un unique vecteur v de C tel que ku − vk = inf w∈C ku − wk. Pour un
tel v , on a <e hu − v, w − vi ≤ 0 pour tout w ∈ C . Le vecteur v s’appelle la projection de
u sur la partie C .
C’est le corollaire suivant qui est particulièrement important.
Corollaire 1.1. Soit F un sous-espace vectoriel fermé de E . Alors, pour tout vecteur u ∈ E ,
il existe un unique vecteur v dans E , dit projection orthogonale de u sur E , tel que
(3) ku − vk = inf ku − wk.
w∈F

Le vecteur v ⊥ = u − v est dans l’orthogonal F ⊥ et on a la décomposition en somme directe

⊥ ⊥
E = F⊥ ⊕F , i. e. tout vecteur u de E s’écrit comme somme u = v + v avec v ∈ F et
⊥ ⊥
v ∈ F , et ceci de manière unique.
Démonstration. — Un sous-espace vectoriel étant convexe, il suffit d’appliquer le théorème
précédent pour avoir l’existence et unicité de la projection v . Pour w ∈ F et ε avec ε4 = 1,
le vecteur wε = v + εw est dans F et les quatre inégalités <e hu − v, wε − vi ≤ 0 donnent
l’égalité hu − v, wi = 0, ce qui exprime l’orthogonalité de u − v et F .
4 Remarque 1.4. Le corollaire précédent donne la solution géométrique de problèmes de
minimisation. Par ex., si f est un élément L2 ([0, 1]), chercher à minimiser le défaut kf − P k
de f à être un polynôme de degré au plus n revient à prendre la projection orthogonale
πn (f ) de f sur le sous-espace Pn des polynômes de degré au plus n :
kf − πn (f )k2 = inf kf − Pn k2 = inf kf − a0 − a1 t − . . . − an tn k2 . 5
P ∈Pn a0 ,...,an

1.3. Forme linéaire

Si v est un vecteur d’un espace E muni d’un produit scalaire, alors l’application u ∈ E →
hu, vi ∈ C est linéaire continue, vu que |hu, vi| ≤ kvk kuk. La réciproque est vraie dans un
espace de Hilbert, au sens où toute forme ` (i. e. une application linéaire ` : E → C) est
représentée de cette manière.
Théorème 1.2 (Fischer-Riesz). Soit E un espace de Hilbert. Pour toute forme linéaire
continue ` sur E , il existe un unique vecteur v de E tel que `(u) = hu, vi, u ∈ E .

Démonstration. — L’unicité d’un v résulte du fait que si hw, ui = 0 pour tout u ∈ E , alors le vecteur w est nul. Reste
donc à prouver l’existence. Si la forme ` est nulle, on peut prendre v = 0 . Sinon, le sous-espace K = {v ∈ E, `(v) = 0}
n’est pas E tout entier. Soit u0 un vecteur hors de K , qu’on peut supposer être dans K ⊥ quitte à soustraire à u0
son projeté orthogonal sur K . Alors, tout vecteur u s’écrit comme somme
„ «
hu, u0 i hu, u0 i
u= u− u0 + u0 ,
hu0 , u0 i hu0 , u0 i
 1.4. BASE HILBERTIENNE 5

où le premier terme est dans K . Ainsi, son évaluation par ` est celle du dernier terme, soit

hu, u0 i `(u0 )
`(u) = `(u0 ) = hu, u0 i.
hu0 , u0 i hu0 , u0 i

1.4. Base hilbertienne

La définition suivante généralise les repères orthonormés du plan ou de l’espace

Définition 1.6. Soit E un espace de Hilbert. La famille orthonormée

P (en )n∈N est une base
orthonormée de E si tout vecteur v est somme d’une série n∈N xn en , i. e.
N
X
lim v− xn en = 0.
N →∞
n=0

Les complexes (xn )n∈N sont appelées les coordonnées du vecteur v dans la base (en )n∈N et
on écrit
X ∞
X
v= xn en ou xn en .
n∈N n=0

. Exemple 1.5. Dans l’espace `2 (N), la famille (en )n∈N , où en est la suite de `2 (N) dont
le seul coefficient non nul est le n-ème valant 1, est une base orthonormée. /
On a le théorème d’existence

Théorème ∗ 1.3. Soit E un espace de Hilbert séparable non de dimension finie. Alors E
admet une base orthonormée.

Démonstration. — Le principe de la construction est simple. On définit une suite de vecteurs

unitaires par récurrence. On prend tout d’abord un premier vecteur unitaire, soit e1 . Puis,
e1 , . . . , en étant construit, on considère un vecteur en+1 de l’orthogonal Vn⊥ du sous-espace
Vn = {x1 e1 + . . . + xn en , xi ∈ C, i = 1, . . . , n} engendré par les e1 , . . . , en . Si une telle
construction ne peut aller au delà du rang N , c’est que VN = E et E est de dimension finie,
ce qui est impossible. Sinon elle se poursuit à l’infini. Un argument convenable permet de
conclure (on utilise ici l’hypothèse de séparabilité de l’espace).

. Exemple 1.6. La famille des exponentielles √ (e−2jπnt/T√)n∈Z est une famille orthonormée
2
de L ([0, T ]), ainsi que la famille (1, ( 2 cos(2πnt/T ), 2 sin(2πnt/T ))n∈N∗ ) déduite par
combinaisons linéaires finies de la première. On montre que ce sont des bases orthonormées
de L2 ([0, T ]) : le sous-espace qu’elles engendrent y est dense.
Ainsi, si cn (g), an (g), bn (g) sont les coefficients de Fourier de g
1 T
Z
cn (g) = f (t)e−2jπnt/T dt, n ∈ Z,
T 0
1 T
Z √
an (g) = f (t) 2 sin(2πnt/T )dt, n ∈ N,
T 0
1 T
Z √
bn (g) = f (t) 2 cos(2πnt/T )dt, n ∈ N∗ ,
T 0
on a la décomposition suivant ces bases
X Xh √ √ i
(4) g(t) = cn (g)e2jπnt/T = c0 (g) + an (g) 2 sin(2πnt/T ) + bn (g) 2 cos(2πnt/T ) .
n∈Z n∈N∗
6 CHAPITRE 1. ESPACES DE HILBERT

et la relation de Parseval (5), dite parfois de Bessel-Parseval

/
X X
kgk22 = |cn (g)|2 = |a0 (g)|2 + |an (g)|2 + |bn (g)|2 .

n∈Z n∈N∗

Proposition 1.2. Soit E un espace de Hilbert muni d’une base orthonormée (en )n∈N et v
un vecteur de E . Alors, les coordonnées de v dans la base (en )n∈N sont données par les
produits scalaires hv, en i, soit
X
v= hv, en ien ,
n∈N

et on a l’égalité de Parseval
X
(5) kvk2 = |hv, en i|2 .
n∈N

Démonstration. — Soit v un vecteur de E de coordonnées (xn )n∈N Par continuité de la

forme linéaire v → hv, em i, on a
XN
< v, em >= lim h xn en , em i = xm .
N →∞
n=0

PN 2 PN 2
L’application de Pythagore donne n=0 hv, en ien = n=0 |hv, en i| . D’après l’inégalité
PN
≤ v− N
P
triangulaire kvk − n=0 hv, en ien n=0 hv, en ien , on a donc

N
X
kvk − hv, en ien →0
n=0

lorsque N → ∞. L’égalité de Parseval s’en déduit.

1.5. Exercices

1. Démontrer l’identité d’Apollonius, valable pour tout x, y, z vecteurs de l’espace

de Hilbert E
2
1 x+y
kz − xk + kz − yk = kx − yk2 + 2 z −
2 2
.
2 2

2. Soit E , espace vectoriel sur C, muni d’un produit scalaire. Montrer que
kx + yk2 − kx − yk2 + jkx + jyk2 − jkx − jyk2 = 4hx, yi, x, y ∈ E.

3. Soit P un polynôme P (t) = an tn + · · · + a1 t + a0 avec a2n + . . . + a21 + a20 = 1.

Montrer que
Z 1
|P (t)|dt ≤ π/2
0

On pourra utiliser que 1 + t2 + t4 + . . . + t2n ≤ (1 − t2 )−1 , t ∈ (0, 1).

 1.5. EXERCICES 7

4. Soit N un entier non nul. Montrer que

√
Z 2π
1
1 + cos N t dt ≤ 1,
2π 0
puis, pour deux entiers N et M distincts,
1
Z 2π √
ejM t + ejN t dt ≤ 2,
2π 0
et plus généralement, si N1 , . . . , Nk sont des entiers tous distincts deux à deux,
1
Z 2π √
ejN1 t + . . . + ejNk t dt ≤ k.
2π 0

5. Soit P la partie de `2 (N) définie par P = {(xn ), xn ∈ [0, ∞), n ∈ N}. Montrer
que P est convexe.

6. Comment définir un produit scalaire sur `2 (Z) “naturel”. qui en fasse un espace
de Hilbert ?

7. Soit (en )n=1,...,N une famille orthonormée d’un espace de Hilbert H et VN le sous-
espace vectoriel engendré par cette famille, i. e. VN = {x1 e1 + . . . + xN eN , xi ∈
C}. Montrer que la projection orthogonale de u ∈ H sur VN est le vecteur
v= N
P
n=1 hu, ek iek .

8. Soit, pour n ∈ Z, Vn le sous-espace de L2 (R) constitué des fonctions constantes

sur chaque intervalle (k2−n , (k + 1)2−n ). Soit ϕ la fonction de V0 à support dans
[0, 1] et y valant 1.
2
√ Pour a > 0 et f ∈ L (R), on note par Ua f la fonction définie par Ua f (t) =
af (at), t ∈ R,
(a) Dessiner le graphe de ϕ et de la fonction χ = U2 ϕ − τ1/2 U2 ϕ.
(b) Montrer que Vi ⊂ Vi+1 pour tout entier i. Montrer que U2 induit un
isomorphisme de Vi sur Vi+1 .
(c) Montrer que la famille des translatées (τn ϕ)n∈Z est une base orthonormée
de V0 .
(d) Montrer que f ∈ V1 est orthogonal à V0 si et seulement si f est une
combinaison linéaire (éventuellement infinie) de χ et de ses translatées (τk χ)k∈Z .
(e) Soit pour t ∈ R la fonction ϕt définie par ϕt (u) = ϕ(u − [t]), u ∈ R où [t]
désigne la partie entière de t. Soit f dans V0 . Montrer que, pour t ∈ R \ Z
f (t) = hf, ϕt i.

9. Soient les systèmes Cπ = (cos nt)n≥0 , Sπ = (sin nt)n≥1 et Sπ/2 = (sin(2n +

1)t)n≥0 .
(a) Montrer que les systèmes Cπ , Sπ et Sπ/2 sont orthogonaux dans L2 (0, π),
L (0, π) et L2 (0, π/2) resp. Sont-ils orthonormés ?
2

(b) Soit f une fonction de L2 (0, π). En considérant le prolongement pair f+

(resp. impair f− ) de f à (−π, π), i. e. la fonction f± telle que
f± (t) = ±f± (−t) = f (t), t ∈ (0, π), f (0) = 0,
montrer le développement de f en cosinus
1 π X2Z π
Z
f (t) = f (τ )dτ + f (τ ) cos nτ dτ cos nt,
π 0 n>0
π 0
8 CHAPITRE 1. ESPACES DE HILBERT

resp. le développement en fonction sinus

X2Z π
f (t) = f (τ ) sin nτ dτ sin nt.
n>0
π 0

(c) Soit f une fonction définie sur (0, π/2) et g son prolongement pair à (0, π),
i. e. la fonction g telle que
g(π/2 + t) = g(π/2 − t), t ∈ (0, π/2).
En utilisant le développement en sinus de g , montrer que
X 4 Z π/2
f (t) = f (τ ) sin(2n + 1)τ dτ sin(2n + 1)t, t ∈ (0, π/2).
n>0
π 0

10. Montrer que

∞
2 4 X cos 2nt
−
| sin t| = , t ∈ (−π, π).
π π n=1 4n2 − 1
En déduire une valeur de ∞ 2 −2
P
n=1 (4n − 1) .
11. Soit a, T réels avec 0 < a < T . Calculer dans L2 (−T, T ) les coefficients de
Fourier de la fonction χ[−a,a] , fonction caractéristique de l’intervalle [−a, a]. En
utilisant l’identité de Parseval, en déduire l’égalité
X∞
−1
α =1+2 [sinc(πnα)]2 ,
n=1

où α ∈ (0, 1).
12. Soit H2 l’espace des séries entières f (z) = ∞
P n
P∞ 2
n=0 an z telle que n=0 |an | /(1+n)
soit finie. On le munit du produit scalaire
Z
hf, gi = f (z)g(z)dxdy.
|z|<1
n
(a) Montrer que (z )n∈N est une famille orthogonale de H2 . Quelle est la
norme de z n ?
(b) Soit, pour w complexe avec |w| < 1, la fonction Sw définie sur {|z| < 1}
par Sw (z) = π1 (1 − wz)−2 . Montrer que Sw appartient à H2 , puis que, pour w
avec |w| < 1,
f (w) = hf, Sw i.
 CHAPITRE 2

SYSTÈMES ORTHOGONAUX

Après avoir rappelé comment les exponentielles à fréquences entières constituent un sys-
tème orthogonal dans L2 (0, 1) et quelques résultats de convergence pour les séries de Fourier
d’une fonction, nous introduisons les polynômes de Legendre, de Chebyshev(1) et d’Hermite.
Même si la définition de ces familles de polynômes n’est pas familière, les systèmes orthogo-
naux qu’ils constituent permettent de faire des analyse/synthèse en tout point similaires à
celles de Fourier : les tracés, et quelques démonstrations, témoignent de ces similarités pro-
fondes. Ce chapitre se termine par la présentation du système de Haar introduit au début du
XXe siècle, ancêtre des bases d’ondelettes dont l’efficacité a été plébiscitée dans de nombreux
domaines rapidement après leur introduction à la fin du siècle dernier.

2.1. Exponentielles de Fourier

La famille (en )n∈Z des fonctions exponentielles et la famille (En )n≥0 construite à partir
de leurs parties réelles et imaginaires sont à la base de la théorie de Fourier. Rappelons leur
définition et quelques propriétés importantes.

Théorème/Définition 2.1. Les fonctions en sont définies par

en (t) = e2jπnt , t ∈ [0, 1], n ∈ Z,
les fonctions En par

1,

√
si n = 0
En (t) = 2 cos πnt, si n est pair non nul
√
 2 sin π(n + 1)t, si n est impair
La fonctions en (resp. En ) vérifie l’équation différentielle
y 00 + 4π 2 ω(n)2 y = 0,
avec ω(n) = |n| ([(n + 1)/2] resp.).
La famille (en )n∈Z et la famille (En )n≥0 sont des bases orthonormées de l’espace de Hilbert
L2 (0, 1).

(1)
La translittération de Qebyxv est variable. Ainsi, dans la base de références bibliographiques Zen-
tralblatt Math, 8 000 recensions contiennent référence à Qebyxv, avec 35 translittérations différentes :
Cebisev Cebysev Chebichev Chebisheff Chebishev Chebyschev Chebyshev
Tchebichef Tchébichef Tchebichev Tchebicheff Tchébicheff Tchebischef Tchébischeff
Tchebycev Tchebychef Tchébychef Tchebycheff Tchébycheff Tchebychev Tchébychev
Tchébychew Tchebyschef Tchébyschef Tchebyscheff Tchebyschev Tchebyshef Tchebysheff
Tchebyshev Tchebytchev Tschebicheff Tschebychev Tschebyscheff Tschebyschew Tschebyshev
10 CHAPITRE 2. SYSTÈMES ORTHOGONAUX

Figure 2 . Les fonctions cos/sinusoïdales En sur [0, 1] pour n = 0, . . . , 20.

Les théorèmes de représentation d’une fonction f définie sur [0, 1] comme somme de sa
série de Fourier abondent : le type de convergence dépend fortement des propriétés de ré-
gularité de la fonction. Les résultats suivants rappellent les énoncés pour des fonctions L2
(théorie hilbertienne), indéfiniment dérivable, continue (le théorème de Dirichlet) et généra-
lisée (théorie des distributions de Schwartz). Les coefficients de Fourier sont définis pour une
fonction f intégrable
Z 1
cn (f ) = e−2jπnt f (t)dt,
0

Z 1 Z 1f (t)dt,
√ si n = 0
Cn (f ) = En (t)f (t)dt = 2 cos(πnt)f (t)dt, si n est pair non nul
0 √2 sin(π(n + 1)t)f (t)dt, si n est impair
0 

Théorème 2.1. Soit f une fonction de carré intégrable sur (0, 1). Alors, pour la convergence
dans L2 (0, 1),

X ∞
X
(6) f= cn (f )en = Cn (f )En .
n∈Z n=0
 2.1. EXPONENTIELLES DE FOURIER 11

et on a la relation de Bessel-Parseval
Z 1 X ∞
X
|f (t)|2 dt = |cn (f )|2 = |Cn (f )|2 .
−1 n∈Z n=0

Figure 3 . Approximations de la fonction g1 définie sur [0, 1] par g1 (t) = t

< 1/2 et g1 (t) = t − 1 si t ≥ 1/2 par les sommes partielles de Fourier
si t P
π −1 1≤2k−1≤N (−1)k+1 k −1 sin k2πt pour N = 2, 6, 10, 14, 18, 22, 50, 100.

. Exemple 2.1. Soit g1 la fonction définie sur [0, 1] telle que g1 (t) = t si t < 1/2 et
g1 (t) = t − 1 si t ≥ 1/2. On a alors, au sens L2 dans l’espace de Hilbert L2 (0, 1),
1 X (−1)n+1 2jπnt 1 X (−1)n+1
(7) g1 (t) = e = sin(2πnt)
2jπ n∈Z∗ n π n∈N∗ n

et la relation de Bessel-Parseval (5) s’écrit

Z 1
1 X 1
= g12 (t)dt =
12 0 n∈Z∗
(2πn)2

qui donne la somme des inverses des carrés des entiers(2)

X 1 π2
(8) 2
= .
n∈N ∗
n 6

La fonction P définie sur [0, 1] par P (t) = t2 , t ∈ (0, 1), a comme développement de Fourier
∞ ∞
2 1 X cos(2πnt) X sin(2πnt)
t = + − , t ∈ (0, 1),
3 n=1 (πn)2 n=1
πn

(2)
Trouver la valeur p de cette série a été énoncé comme le Problème de Bâle, problème proposé par Pietro
Mengoli en 1644. L. Euler, dont on célèbre le trois-centième anniversaire de la naissance en 2007, le résolut
en 1735 en affirmant, à raison mais sans Q∞en voir la justification nécessaire, que le sinus cardinal sinc(x)
était égal au produit infini p∞ (x) = k=1 (1 − (x/(πk))2 ), d’où la valeur p de la série en égalant les
coefficients du développement de Taylor à l’ordre 2 en x = 0 : restait à démontrer que sur sinc /p∞ ,
fonction holomorphe sur C sans zéro (donc de la forme eg(x) ) était la fonction constante égale à 1. Euler
donna une autre démonstration en 1743 basée sur le développement en série de arsin2 (dont la dérivée
est 2(1 − x2 )−1/2 arcsin), preuve considérée comme exacte. Si d’autres preuves basées sur le calcul intégral
utilisent les mathématiques du XVIe siècle, celle dérivée ici des séries de Fourier est résolument moderne.
 12 CHAPITRE 2. SYSTÈMES ORTHOGONAUX

et l’identité de Parseval donne, en tenant compte de la somme (8),

X 1 π4
n4
= .
90
/
n∈N∗

Théorème 2.2. Soit h fonction C ∞ sur [0, 1] avec h(k) (0) = h(k) (1), k ≥ 0. Alors les coef-
ficients de Fourier (cn (h))n∈Z , (Cn (h))n≥0 de h sont à décroissance rapide et
X ∞
X
h(t) = cn (h)en (t) = Cn (h)En (t), t ∈ [0, 1]
n∈Z n=0

avec convergence absolue des séries, ainsi que pour leurs dérivées terme à terme.

Démonstration. — En intégrant par parties, on a, pour tout entier k ,

cn (h) = (jn)−k cn (h(k) ), n 6= 0,
ce qui assure la décroissance rapide des coefficients de Fourier cn (h) et des Cn (h) (obtenus
par combinaison linéaire simple des cn (h)).

Théorème ∗ 2.3 (Dirichlet). Soit f continue par morceaux sur (0, 1), ayant des limites à
droite f (t+ ) et à gauche f (t− ) en tout point t ∈ (0, 1) et de carré intégrable sur [0, 1]. Alors
N N
X X f (t− ) + f (t+ )
lim cn (f )en (t) = lim Cn (f )En (t) = , t ∈ (0, 1).
N →∞
n=−N
N →∞
n=0
2

. Exemple 2.2. Le développement (7) pris aux valeurs t = 0, 1/2 et 1 donne bien (g1 (t−)+
g1 (t+))/2. Pour t = 1/4, on obtient
∞ ∞
1 1 X sin(π(2k + 1)/2) 1 X (−1)k )
= =
4 π k=0 2k + 1 π k=0 2k + 1
d’où
1 1 1 1 π
1− + − + − ... = .
3 5 7 9 4
/
4 Remarque 2.1. Pour la base des exponentielles (e2jπn· )n∈Z , la première égalité dans (4)
signifie la convergence
q
X
g(t) = lim cn (g)e2jπnt
p,q→+∞
n=−p

dans L2 (0, 1), i. e.

q
X
lim cn (g)e2jπnt − g(t) = 0.
p,q→+∞
n=−p 2
Les sommes sont asymétriques, au contraire du théorème de Dirichlet (convergence simple
vers la demi-somme des limites à droite et à gauche) où il est essentiel d’avoir une sommation
entre −N et N . 5

Théorème ∗ 2.4. Soit U une distribution tempérée à support dans (0, 1). Alors
X X
U= cn (U )en = Cn (U )En
n∈Z n≥0

où la convergence a lieu au sens distribution.

 2.2. POLYNÔMES DE LEGENDRE 13

2.2. Polynômes de Legendre

Le théorème de Weierstrass(3) énonce que l’espace des fonctions polynômes est dense dans
l’espace des fonctions continues sur [−1, 1] pour la convergence uniforme, cet espace étant
dense dans l’espace L2 (−1, 1) pour la convergence en moyenne quadratique. Ainsi, l’espace
des fonctions polynômes est dense dans l’espace de Hilbert L2 (−1, 1) : toute famille ortho-
normée de polynômes (pn )n≥0 , avec pn de degré n, est une base hilbertienne de√L2 (−1, 1).
Une telle famille est obtenue simplement, aux signes près, en prenant p0 = 1/ 2, puis pn
de degré n, unitaire et orthogonal au sous-espace des polynômes de degré n − 1. La famille
(Ln )n≥0 des polynômes de Legendre est une famille de polynômes orthogonaux : Ln , n ≥ 0
n’est pas unitaire dans L2 (−1, 1), mais normalisé avec Ln (1) = 1.

Figure 4 . Les polynômes de Legendre Ln sur [−1, 1] pour n = 0, . . . , 20.

Théorème/Définition 2.2. La famille des polynômes de Legendre (Ln )n≥0 est caractérisée
par l’une des trois présentations suivantes
1. Les Ln sont donnés par la fonction génératrice
∞
1 X
(9) L(x, t) = √ = Ln (t)xn
1 − 2xt + x 2
n=0
où x est suffisamment petit.
2. Le polynôme Ln , n ≥ 0 est défini par la formule de Rodrigues
1 dn 2 n


(10) Ln (t) = n (t − 1)
2 n! dtn
3. La suite (Ln )n≥0 vérifie la relation de récurrence
(11) (n + 1)Ln+1 (t) = (2n + 1)tLn (t) − nLn−1 (t), n≥1
avec L0 (t) = 1 et L1 (t) = t.
La suite (Ln )n≥0 est orthogonale dans L2 (−1, 1), avec
1 1
Z
1
(12) Ln (t)2 dt = , n ≥ 0.
2 −1 2n + 1
Le polynôme Ln admet la représentation intégrale
1 πh
Z √ in
(13) Ln (t) = t + t2 − 1 cos ϕ dϕ
π 0
(3)
Ce théorème, ainsi que l’assertion suivante, est démontré dans le fascicule Linéarité et convergences,
disponible sur le oueb à des côtés.
14 CHAPITRE 2. SYSTÈMES ORTHOGONAUX

et vérifie l’équation différentielle

(14) (1 − t2 )y 00 − 2ty 0 + n(n + 1)y = 0.

1
Démonstration. — Pour t fixé, la fonction x → √ est holomorphe au voisinage de x = 0 , ainsi, d’après
1−2xt+x2
Cauchy,
∞
1 X
√ = cn (t)xn
1 − 2xt + x2 n=0
avec Z
1 1 dx
cn (t) = √
2jπ1 − 2xt + x2 xn+1
γ

où γ est un petit cercle de centre x = 0 parcouru dans le sens positif. Introduisons u = u(x) vérifiant
p
(15) 1 − ux = 1 − 2xt + x2 .
Ainsi √
1− 1 − 2xt + x2 2
u= = √ (t − x/2)
x 1 + 1 − 2xt + x2
parcourt un contour fermé γ e autour de u = t lorsque x parcourt γ . Développant au carré de (15), on obtient
x = (2(u − t))/(u2 − 1) , et en prenant la différentielle de (15)
dx du
√ = .
x 1 − 2xt + x2 u−x
On en déduit
(u2 − 1)n
Z
1
(16) cn (t) = du
2jπ γ
e 2n (u− t)n+1
soit, par la formule de Cauchy
dn (u2 − 1)n
» –
1
cn (t) = = Ln (t).
n!2n dun u=t
p √
En prenant dans (16) comme contour γ e le cercle de centre t et de rayon |t2 − 1| , soit u(ϕ) = t + t2 − 1ejϕ , ϕ ∈
[−π, π] on obtient
Z π » 2 √ –n
1 t + 2t t2 − 1ejϕ + (t2 − 1)e2jϕ − 1
Ln (t) = √ dϕ
2π −π 2 t2 − 1ejϕ
1 πh
Z p in
= t + t2 − 1 cos ϕ dϕ
π 0
La dérivée ∂x L vérifie
(1 − 2tx + x2 )∂x L + (x − t)L = 0,
d’où en substituant la série et sa dérivée (obtenue par dérivation terme à terme), on a
X X
(1 − 2tx + x2 ) nLn (t)xn−1 + (x − t) Ln (t)xn = 0
n≥0 n≥0

soit en prenant les termes de degré n en x ,

(n + 1)Ln+1 (t) = (2n + 1)tLn (t) − nLn−1 (t), n ≥ 1,
i. e. la relation de récurrence (11).
En développant la formule de Rodrigues, on obtient
[2]n
X (−1)k (2n − k)!
Ln (t) = xn−2k
2n k!(n − k)!(n − 2k)!
k=0
n
ainsi Ln (t) ≤ Ln (ja)/j pour t complexe de module au plus a . Ainsi, x étant fixé, il y a convergence uniforme de
la série L(x, t) sur {t ∈ R, |t| ≤ a} : s’agissant de fonctions holomorphes, on peut dériver terme à terme. On peut
dériver par rapport à t la série donnant la fonction génératrice L . La relation (1 − 2tx + x2 )∂t L = xL donne en
dérivant terme à terme (ce qui est légitime)
(17) L0n+1 − 2tL0n + L0n−1 − Ln = 0
Éliminant entre cette équation et celle obtenue en dérivant (11) le terme L‘n−1 , puis L0n+1 on obtient
(18) L0n+1 − tL0n = (n + 1)Ln ,
(19) tL0n − tL0n−1 = nLn .
En remplaçant n par n − 1 dans (18) et éliminant L0n−1 avec (19), on obtient
(1 − t2 )L0n = nLn−1 − ntLn ,
 2.2. POLYNÔMES DE LEGENDRE 15

qui dérivé, et après élimination de L0n−1 en utilisant (19), donne

[(1 − t2 ]Ln ]0 + n(n + 1)Ln = 0.
C’est l’équation (14).
L’orthogonalité des polynômes de Legendre en résulte : la différence des équations différentielles (14) pour les
indices m et n est
[(1 − x2 )(L0m Ln − L0n Lm )]0 + (m − n)(m + n + 1)Lm Ln = 0
qui intégrée donne
Z 1
(m − n) Lm (t)Ln (t)dt = 0
−1
R1
d’où l’orthogonalité de la famille (Ln ) . Pour calculer −1 L2n (t)dt , on élimine le terme xLn Ln+1 des équations
obtenues en multipliant (11) par Ln+1 et cette même équation (11) au rang n + 1 par Ln :
(n + 1)(2n + 3)L2n+1 + n(2n + 3)Ln−1 Ln+1 − (n + 2)(2n + 1)Ln Ln+2 − (n + 1)(2n + 1)L2n = 0
d’où par intégration, et en utilisant l’orthogonalité démontrée précédemment,
Z 1 Z 1
(n + 1)(2n + 3) L2n+1 (t)dt = (n + 1)(2n + 1) L2n (t)dt
−1 −1

soit Z 1 Z 1
2n + 1
L2n+1 (t)dt = L2n (t)dt
−1 2n + 3 −1

et la formule (12)

4 Remarque 2.2. La fonction génératrice (9) permet l’expression du dipôle de l’électro-

statique ou de la gravitation newtonienne en terme de série de polynômes de Legendre. Soit
→
− r 2 ∈ R3 . En prenant dans la formule (9) x = k−
r 1, −
→ →r 1 k/k−→r 2 k et t = cos(− →
r\ →
−
1 , r 2 ) le

cosinus de l’angle (−→r\ →

− →
− →
−
1 , r 2 ) et en ayant supposé k r 1 k < k r 2 k, on obtient

1 1
→
− →
− =q
k r 1 − r 2k
k−→
r 1 k2 − 2k−
→r 1 k k−
→r 2 k cos(−
→r\ →
− → 2
−
1, r 2) + k r 2k
1 1
= − → q
k r 2k
1 − 2k−
→r 1 k k−
→r 2 k−1 cos(−
→
r\ →
− → −
− → −1 2
1 , r 2 ) + (k r 1 kk r 2 k )

soit, grâce à (9) l’expression du dipôle

 →
k− Ln (cos(−
→ →
−
n
1 X r 1k r\1 , r 2 ))
(20) = .
k−
→
r 1−−→
r 2 k n≥0 k−
→
r 2k k−
→r 2k 5

L’expression des polynômes de Legendre à partir de la fonction génératrice L(x, t) donne

par ailleurs quelques valeurs particulières
Corollaire 2.1. Le polynôme de Legendre Ln est (im)pair si n est (im)pair, et vérifie
1 · 3 · · · (2n − 1)
Ln (1) = 1, Ln (−1) = (−1)n , L2n (0) = (−1)n , L2n−1 (0) = 0.
2 · 4 · · · 2n
Définition 2.1. Les coefficients de Legendre(4) (`n (f ))n≥0 d’une fonction f intégrable sur
(−1, 1) sont définis par
Z 1
`n (f ) = f (t)Ln (t)dt, n ≥ 0.
−1

(4)
On prendra grade aux normalisations : dans l’espace de Hilbert L2 (−1, 1) avec le produit scalaire
−1
R1 √
2 −1
f (t)g(t)dt , les coordonnées de f dans la base orthonormée ( 2n + 1Ln )n≥0 sont les coefficients
p
n + 1/2`n (f ).
16 CHAPITRE 2. SYSTÈMES ORTHOGONAUX

De manière plus générale, si U est une distribution à support dans (−1, 1), ses coefficients
de Legendre sont définis par
`n (U ) = hU, Ln i, n ≥ 0.
√
La suite ( 2n + 1Ln ) étant une base orthonormée de L2 (−1, 1), on a le développement
en série de polynômes de Legendre suivant la convergence quadratique pour toute fonction
de carré intégrable sur [-1,1] :
Théorème 2.5. Soit f une fonction de carré intégrable sur [−1, 1]. Alors, dans L2 (−1, 1),
∞
X
f= (n + 1/2)`n (f )Ln ,
n=0

et on a la relation de Bessel-Parseval
Z 1 ∞
X
2
|f (t)| dt = (n + 1/2)|`n (f )|2 .
−1 n=0

Figure 5 . Approximations de la fonction f1 définie sur [−1, 1] par f1 (t) = t si

|t| < 1/2 et f1 (t) = t − sgn(t) si |t| ≥ 1/2 par la somme de polynômes de Legendre
PN
n=0 (n + 1/2)`n (f1 )Ln pour N = 2, 6, 10, 14, 18, 22, 50, 100.

Comme pour les séries de Fourier, on a des résultats de convergence ponctuelles si la

fonction est plus régulière que simplement de carré intégrable.
Théorème
R1 2.6. Soit h fonction C ∞ sur [−1, 1]. Alors les coefficients de Legendre `n (h) =
−1
h(u)Ln (u)du de h sont à décroissance rapide et
∞
X
f (t) = (n + 1/2)`n (h)Ln (t), t ∈ [−1, 1]
n=0

avec convergence absolue de la série, ainsi que pour ses dérivées terme à terme.
Démonstration. — L’équation différentielle (14) s’écrit ((1 − t2 )L0n )0 + n(n + 1)Ln = 0, ainsi
en intégrant par partie
Z 1 Z 1
1 2 0 0 1
`n (h) = − Ln (t)[(1 − t )h (t)] dt = Ln (t)DLk (h)dt
n(n + 1) −1 (n(n + 1))k −1
où on a introduit l’opérateur différentiel DL défini pour une fonction f deux fois dérivable
sur [−1, 1] par
DL (f )(t) = −[(1 − t2 )f 0 (t)]0 , t ∈ [−1, 1].
La convergence de la série en résulte, en observant que Ln est borné sur (−1, 1) d’après la
formule de représentation intégrale (13) et des majorations analogues pour ses dérivées.
Le théorème de Dirichlet (cf. Théorème 2.3) pour les séries de Fourier a son analogue pour
les séries de Legendre.
 2.2. POLYNÔMES DE LEGENDRE 17

Théorème ∗ 2.7. Soit f continue par morceaux sur (−1, 1), ayant des limites à droite et à
gauche en tout point. Alors
N
X f (t− ) + f (t+ )
lim (n + 1/2)`n (f )Ln (t) = , t ∈ (−1, 1).
N →∞
n=0
2

On a même un théorème de développement de Legendre de distributions.

Théorème ∗ 2.8. Soit U une distribution tempérée à support dans (−1, 1). Alors
X
U= (n + 1/2)hU, Ln iLn ,
n≥0

avec convergence au sens distribution.

. Exemples 2.3.
1. En utilisant la formule de Rodrigues (10) et n intégrations par parties, on obtient pour
u ∈ R,
Z 1
dn (−ju)n 1 jtu 2
Z
jtu jtu 1 2 n
e (t − 1)n dt


he , Ln i = e n n! dtn
(t − 1) dt = n
−1 2 2 n! −1
n Z 1
(ju)
= n cos(tu)(1 − t2 )n dt.
2 n! −1
Vu la représentation intégrale (41) ci-après, les coefficients de Legendre de l’exponen-
tielle ejtu s’expriment en terme de fonction de Bessel
√ Jn+1/2 (u)
(21) hejtu , Ln i = jn 2π √ ,
u
soit la décomposition de Rayleigh
r
jut
X
n π Jn+1/2 (u)
e = (2n + 1)j √ Ln (t), t ∈ (−1, 1).
n≥0
2 u

2. Soit, pour α ∈ (−1, 1), la fonction caractéristique fα de l’intervalle [α, 1]. La somme
de (17) et (18) donnent
L0n+1 − L0n−1 = (2n + 1)Ln , n ≥ 1,
ainsi
Z 1 Z 1  0
Ln+1 (t) − L0n−1 (t) dt = Ln−1 (α) − Ln+1 (α),

(2n + 1) Ln (t)dt = n ≥ 1,
α α
et
∞
1 1X
fα = (1 − α) + [Ln−1 (α) − Ln+1 (α)]Ln .
2 2 n=1
avec convergences comme p décrites dans les théorèmes 2.5 et 2.7.
3. Pour la fonction f (t) = (1 − t)/2, on multiplie la définition (9) de la fonction généra-
trice par f , puis on intègre les deux membres. L’intégration terme à terme du membre
de droite est justifiée par la majoration |L
√ n (t)| ≤ 1, t ∈ [−1, 1] obtenue à partir de (13)
et on utilise le changement de variable 1 − t = (1 − x)(2x)−1/2 sh v pour le membre
de gauche :
√  X n Z 1r
(1 − x)2

1 1+ x n 1−t
1+x− √ log √ = x Ln (t)dt, |x| < 1.
2x 2 x 1− x n=0 −1 2
18 CHAPITRE 2. SYSTÈMES ORTHOGONAUX

Le membre de gauche admet le développement

∞
4 X xn
−4
3 n=1
(4n2 − 1)(2n + 3)
On en déduit
r ∞
1−t 2 Ln (t)
/
X
= L0 (t) − 2 .
2 3 n=1
(2n − 1)(2n + 3)

Les graphes de la figure 4 ont été tracés en utilisant la relation de récurrence (11) : la
commande
with(orthopoly) ; for i to 6 P(i,t) od ;
adressée à maple TM donne les expressions des polynômes de Legendre L0 , . . . , L6 .
1 1 1
L0 (t) = 1, L1 (t) = t, L2 (t) = (3t2 − 1), L3 (t) = (5t3 − 3t), L4 (t) = (35t4 − 30t2 + 3),
2 2 8
1 1
L5 (t) = (63t5 − 70t3 + 15t), L6 (t) = (231t6 − 315t4 + 105t2 − 5).
8 16
Les polynômes de Legendre sont bien sûr aussi disponibles sur matlab TM .

2.3. Polynômes de Chebyshev

Dans l’espace de Hilbert T1 défini par
 Z 1 
2 2 −1/2 2 dt
T1 = L ([−1, 1], (1 − t ) dt) = f : (−1, 1) → C, |f (t)| √ < +∞
−1 1 − t2
avec le produit scalaire
1 1
Z
dt
hf, gi = f (t)g(t) √ ,
π −1 1 − t2
le sous-espace des fonctions polynomiales est dense. Comme pour L2 (−1, 1) et les polynômes
de Legendre, toute famille orthonormée de polynômes (pn )n≥0 , avec pn de degré n, est une
base hilbertienne de l’espace T1 . La famille (Tn )n≥0 des polynômes de Chebyshev est une
famille de polynômes orthogonaux : les polynômes Tn sont normalisés avec Tn (1) = 1, n ≥ 0.

Figure 6 . Les polynômes de Chebyshev Tn sur [−1, 1] pour n = 0, . . . , 20.

Théorème/Définition ∗ 2.3. La famille des polynômes de Chebyshev (Tn )n≥0 est caracté-
risée par l’une des trois présentations suivantes
 2.3. POLYNÔMES DE CHEBYSHEV 19

1. Les Tn sont donnés par la fonction génératrice

∞
1 − xt X
(22) T (x, t) = = Tn (t)xn
1 − 2xt + x2 n=0

où x est suffisamment petit.

2. Le polynôme Tn , n ≥ 0 est défini par la formule de Rodrigues
(−1)n 22n−1 (n − 1)! √ dn
1 − t2 n (1 − t2 )n−1/2


Tn (t) =
(2n)! dt
3. La suite (Tn )n≥0 vérifie la relation de récurrence
(23) Tn+1 (t) = 2tTn (t) − Tn−1 (t), n≥1
avec T0 (t) = 1 et T1 (t) = t.
√
La suite T0 , ( 2Tn )n>0 est orthonormée dans l’espace de Hilbert T1 .
Le polynôme Tn vérifie la relation
(24) Tn (t) = cos(n Arcos(t)) t ∈ [−1, 1],
a comme développement et factorisation
[n/2]   n   
X n n−2m 2 m n−1
Y (2k − 1)π
Tn (t) = t (t − 1) = 2 t − cos
m=0
2m k=1
2n

et satisfait l’équation différentielle

(25) (1 − t2 )y 00 + ty 0 + n2 y = 0.

L’expression des polynômes de Chebyshev à partir de la fonction génératrice donne

quelques valeurs particulières

Corollaire 2.2. Le polynôme de Chebyshev Tn est (im)pair si n est (im)pair, et vérifie

Tn (1) = 1, Tn (−1) = (−1)n , T2n (0) = 1, T2n−1 (0) = 0.

On définit comme pour les polynômes de Legendre la suite (tn (f ))n≥0 des coefficients de
Chebyshev d’une fonction f intégrable sur (−1, 1) relativement à la mesure (1 − t2 )−1/2 ou
d’une distribution U à support dans (−1, 1)
1 1
Z
dt
tn (f ) = f (t)Tn (t) √ , tn (U ) = hU, π −1 Tn (t)(1 − t2 )−1/2 i.
π −1 1−t 2

Les théorèmes 2.5, 2.6, 2.7, 2.8 ont leur équivalent pour les coefficients de Chebyshev. Ils
concernent la convergence de la série suivant diverses topologies, et son égalité avec la fonction
(voire la distribution) f avec comme suite de coefficients de Chebyshev (tn (f ))n≥0
∞
X
f (t) = t0 (f ) + 2 tn (f )Tn (t), t ∈ (−1, 1).
n=1

Il est inutile de les reproduire ici.

4 Remarque 2.3. En effectuant le changement de variable t = cos(πu), u ∈ (0, 1), on in-

troduit un isomorphisme isométrique entre l’espace T1 et l’espace L2 (0, 1), le développement
20 CHAPITRE 2. SYSTÈMES ORTHOGONAUX

Figure 7 . Approximations de la fonction f1 définie sur [−1, 1] par f1 (t) = t si

|t| < 1/2 et f1 (t) = t−sgn(t) si |t| ≥ 1/2 par la somme de polynômes de Chebyshev
PN
n=0 tn (f1 )Tn pour N = 2, 6, 10, 14, 18, 22, 50, 100.

en série de polynômes de Chebyshev se ramenant en un développement en série de cosinus

(cf Exo. 1.5.8)
1 1 X1Z 1
Z
dτ dτ
f (t) − f (τ ) √ =2 f (τ )Tn (τ ) √ Tn (t)
π −1 1 − τ2 n≥1
π −1 1 − τ 2

XZ 1
=2 f (cos πν)Tn (cos πν)dν Tn (cos πu)
n≥1 0

XZ 1
=2 f (cos πν) cos(nπν)dν cos(nπu)
n≥1 0
Z 1
= f (cos πu) − f (cos πν)dν
0
Cela éclaire la simplicité de la forme (24) des polynômes de Chebyshev : la théorie analogue
développée dans l’espace Tα = L2 ([−1, 1], (1 − t2 )−α/2 dt), α < 2 introduit des systèmes de
polynômes orthogonaux de Gegenbauer (Tn (α))n≥0 , cas particuliers des familles polynomiales
dites Jacobi, sans expression aussi simple que pour les polynômes de Chebyshev.

2.4. Fonctions d’Hermite

L’espace de Hilbert
 Z 
2 −t2 2 −t2
L (R, e ) = f : R → C, |f (t)| e dt < +∞
R
contient les fonctions polynômes.
Proposition 2.1. L’espace des fonctions polynômes est dense dans l’espace de Hilbert
2
L2 (R, e−t ).
2
Démonstration. — Soit f ∈ L2 (R, e−t ) orthogonale à tous les polynômes tn , n ≤ 0. Vu
Z 2 Z Z
−ts −t2 −2t<e s −t2 2
e f (t)e dt ≤ e e dt |f (t)|2 e−t dt,
R R R
la transformée de Laplace
Z
−t2 2
L(f (t)e )(p) = e−tp f (t)e−t dt
R
est définie holomorphe sur C. Ses dérivées en p = 0
2
dn L(f (t)e−t )
Z
2

n
(0) = (−t)n f (t)e−t dt
dp R
 2.4. FONCTIONS D’HERMITE 21

Figure 8 . Les fonctions d’Hermite hn pour n = 0, . . . , 19.

22 CHAPITRE 2. SYSTÈMES ORTHOGONAUX

2
sont toutes nulles, L(f (t)e−t ) est nulle et donc aussi f , d’après l’injectivité de la transfor-
mation de Laplace L.
Ce qui a été énoncé ci-dessus pour les polynômes de Legendre et l’espace de Hilbert
2
L2 ([−1, 1]) se développe de manière analogue dans l’espace de Hilbert L2 (R, e−t ), comme
l’indique le théorème suivant
Théorème/Définition ∗ 2.4. La famille des polynômes de Hermite (Hn )n≥0 est caractéri-
sée par l’une des trois présentations suivantes
1. Les Hn sont donnés par la fonction génératrice
∞
−2tx−x2
X
(26) H(x, t) = e = Hn (t)xn .
n=0
2. Le polynôme Hn , n ≥ 0 est défini par la formule de Rodrigues
n
2 d 2
Hn (t) = (−1)n et n e−t .
dt
3. La suite (Hn )n≥0 vérifie la relation de récurrence
(27) Hn+1 (t) = 2tHn (t) − 2nHn−1 (t), n ≥ 1,
avec H0 (t) = 1 et H1 (t) = 2t.
2
La suite (Hn )n≥1 est orthogonale dans L2 (R, e−t dt) avec
√
Z
2
(28) Hn2 (t)e−t dt = 2n n! π.
R
e n = Hn / 2n n!√π est une base
p
La famille (H
e n )n≥0 des polynômes d’Hermite unitaires H
2
orthonormée de L2 (R, e−t dt).
Le polynôme Hn vérifie l’équation différentielle
(29) Hn00 − 2xHn0 + 2nHn = 0
Définition 2.2. Les fonctions d’Hermite, définies suivant,
e n (t)e−t2 /2 , t ∈ R, n ≥ 0.
hn (t) = H
constituent une base orthonormée de l’espace de Hilbert L2 (R).
Les polynômes de Hermite sont disponibles dans la bibliothèque de fonctions spéciales de
maple TM , mais apparemment pas dans celle de matlab TM .

2.5. Système de Haar

Pour a < b, soit χ[a,b] la fonction caractéristique de l’intervalle [a, b] définie sur R par
(
1, si a ≤ t ≤ b,
χ[a,b] (t) =
0, sinon.
L’ensemble des fonctions χ[a,b] , a, b ∈ R engendre linéairement l’espace des fonctions en
escalier à support compact, qui est dense dans l’espace de Hilbert L2 (R).
Il en est donc de même du sous-espace V∞ engendré par les fonctions χ[a,b] avec les extré-
mités a, b dans l’ensemble des rationnels dyadiques D∞ = {k2−` , k ∈ Z, ` ∈ N}. Ainsi toute
famille orthonormée engendrant linéairement le sous-espace V∞ est une base hilbertienne du
Hilbert L2 (R).
Si, pour N ∈ Z, VN désigne le sous-espace des fonctions en escalier à discontinuités dans
l’ensemble DN = {k2−N , k ∈ Z} des dyadiques de niveau N , la famille (VN )N ∈Z vérifie les
propriétés suivantes, dont la preuve est rapide,
 2.5. SYSTÈME DE HAAR 23

Figure 9 . Approximations sur R de la fonction f1 à support [−1, 1] et telle que

f1 (t) = t si |t| < 1/2 et f1 (t) = t − sgn(t) si |t| ≥ 1/2 par la somme de polynômes
PN −2
d’Hermite n=0 kHn k hfT , Hn iHn pour N = 2, 6, 10, 14, 22, 50, 100, 200, 400,
700, 1000.

1. VN ⊂ VN +1 pour N ∈ Z,
2. ∪N ∈Z VN = V∞ , ∩N ∈Z VN = {0},
3. : τ1 (f )(t) = f (t − 1), t ∈ R.
L’espace V0 est stable par la translation entière τ1 √
4. Le changement d’échelle h définie par h(f )(t) = 2f (2t) applique isométriquement
VN sur VN +1 .
 24 CHAPITRE 2. SYSTÈMES ORTHOGONAUX

1 1

0 1 0 1 m2−` (m + 1)2−` m2−` (m + 1)2−`

−1

Figure 10 . Les ondelettes de Haar ϕ, ψ , ainsi que ϕ`,m , ψ`,m .

√
Théorème 2.9. Soient ϕ = χ[0,1] , ψ = 2[hϕ − hτ1 ϕ] et ϕ`,m = h` τm ϕ, ψ`,m = h` τm ψ pour
m, ` ∈ Z.
Alors la famille (ϕ`,m )m∈Z est une base hilbertienne de V` et la famille (ψ`,m )`,m∈Z , dit
système de Haar, est une base orthonormée de L2 (R).
Démonstration. — La famille (ϕ0,m )m∈Z est une base de V0 ; par changement d’échelle h` il
en est de même pour (ϕ`,m )m∈Z et V` .
Soit une fonction f de l’orthogonal W0 de V0 dans V1 . Alors, il existe des complexes
fm , m ∈ Z tels que X
f= fm ϕ1,m .
m∈Z
L’orthogonalité avec ϕ0,M donne fM + fM +1 = 0, d’où
X X √
f= f2M [ϕ1,2M − ϕ1,2M +1 ] = f2M 2ψ0,M .
m∈Z m∈Z

On en déduit que (ψ0,M )M ∈Z est une base orthonormée de W0 . Comme précédemment, par
changement d’échelle, (ψ`,m )m∈Z est une base orthonormée de W` , orthogonal de V` dans
V`+1 . On conclut en remarquant que L2 = ⊕`∈Z W` .
Ainsi, si pour f ∈ L2 (R) les scalaires f`,m désignent les coefficients de f dans le système
de Haar, à savoir, (ψ`,m )`,m∈Z ,
Z Z
2`/2 ψ 2` t − m f (t)dt,


f`,m = ψ`,m (t)f (t)dt =
R R
2
on a, pour la convergence L ,
X X
f`,m 2`/2 ψ 2` t − m


f (t) = f`,m ψ`,m (t) =
`,m∈Z `,m∈Z

Si f est à support compact, l’approximation fL , L ∈ N définie par

L−1 X
X
fL (t) = f`,m ψ`,m (t)
`=−∞ `∈Z

est une somme finie, avec discontinuités dans l’ensemble DL des dyadiques de niveau L.
C’est un élément de VL , en fait la projection de f sur ce sous-espace vu que
L−1 X
X
hfL , ϕL,n i = hf, ψ`,m ihψ`,m , ϕL,n i = hf, ϕL,n i.
`=−∞ m∈Z

Proposition 2.2. Soit f fonction continue sur R à support compact. Alors la suite (fM )M ≥0
définie par
L−1 X
X
fL (t) = f`,m ψ`,m (t)
`=−∞ `∈Z
converge vers f uniformément.
 2.6. EXERCICES 25

Démonstration. — Soit ε > 0. La fonction f , continue et à support compact, est uniformé-

ment continue. Ainsi, il existe ` tel que si |t − t0 | ≤ 2−` , alors |f (t) − f (t0 )| ≤ ε. Ainsi, pour
t dans [m2−` , (m + 1)2−` ], on a
Z (m+1)2−`
`/2
f` (t) = 2 f (u)du2`/2 ϕ(2` t − m)
m2−`

puisque tous les autres termes de la série sont nuls dans cet intervalle, et par suite
Z (m+1)2−`
`
|f` (t) − f (t)| = 2 [f (u) − f (t)] ≤ ε,
m2−`

ce qui achève la preuve.

La proposition précédente approche une fonction continue par des fonctions en escalier :
c’est un des inconvénients du système de Haar. La théorie des ondelettes introduit des couples
de fonctions φ, ψ à support compact, dérivables autant de fois que désiré et qui donnent lieu
à des analyses du même type que le système de Haar : le système (h` τm ψ)`∈Z,m∈Z est une
base orthonormée de L2 (R) et les espaces V` , dont (h` τm ψ)m∈Z est une base orthonormée,
vérifient les propriétés du théorème 2.9.

2.6. Exercices

R1
1. Soit f une fonction définie sur (0, 1) telle que 0 |f (t)|2 dt < ∞. En considé-
rant les développements de Fourier d’une fonction fp (resp. fi ) paire (resp. im-
paire) sur (−1, 1) qui coïncide avec f sur (0, 1), montrer les développements dans
L2 (0, 1) en série de cosinus et sinus
Z 1 XZ 1
f (t) = f (u)du + 2 f (u) cos(nπu)du cos(nπt),
0 n>0 0

XZ 1
f (t) = 2 f (u) sin(nπu)du sin(nπt).
n>0 0

Formuler sous des hypothèses appropriées des résultats de convergence simple

pour les représentations précédentes et appliquer les à la fonction f0 définie par
(
t, si 0 ≤ t ≤ 1/2,
f0 (t) =
t − 1/2, si 1/2 < t ≤ 1.

2. Soient L l’application définie sur l’ensemble C ∞ ([−1, 1]) des fonctions indéfini-
ment dérivables sur [−1, 1] par
Lf (t) = (1 − t2 )f 00 (t) − 2tf 0 (t), t ∈ [−1, 1], f ∈ C ∞ ([−1, 1])
Montrer que
hLf, gi = hf, Lgi f, g ∈ C ∞ ([−1, 1])
R1
où hp, qi = −1 p(t)q(t)dt pour p, q ∈ C ∞ ([−1, 1]).
En déduire que pour fi , i = 1, 2 tels qu’il existe un scalaire λi avec Lfi =
λi fi , i = 1, 2 et λ1 6= λ2 , alors hf1 , f2 i = 0.
26 CHAPITRE 2. SYSTÈMES ORTHOGONAUX

3. Soit Q0 , Q1 définis par sur l’intervalle (−1, 1) par

   
1 1+t 1 1+t
Q0 (t) = log , Q1 (t) = −1 + t log .
2 1−t 2 1−t
Montrer que
(1 − t2 )Q000 (t) − 2tQ00 (t) = 0, (1 − t2 )Q001 (t) − 2tQ01 (t) = −2Q1 .

4. En utilisant la fonction génératrice définissant la suite de polynômes de Legendre

(Ln )n≥0
X
(1 − 2xt + x2 )−1/2 = Ln (t)xn
n≥0
R1
calculer −1
Ln (t)dt et montrer que (n + 1)Ln (t) = L0n+1 (t) − tL0n (t).

5. (a) Montrer que

cos(nt) = 2 cos t cos((n − 1)t) − cos((n − 2)t).
En déduire l’existence, pour tout entier n ≥ 0, d’un unique polynôme Tn tel que
Tn (cos(t)) = cos(nt), t ∈ R.
On pourra, après avoir donner explicitement T0 , T1 , T2 , montrer l’existence de Tn
par une récurrence. Quel est le degré de Tn ?
(b) En calculant la partie réelle de ejnt = (cos t + j sin t)n , montrer que pour
n≥0
[n/2] 
X n
Tn (X) = (X 2 − 1)k X n−2k .
k=0
2k
(c) Montrer que, pour tout entier n ≥ 0, il existe un polynôme Qn tel que
sin t Qn (cos(t)) = sin(nt), t ∈ R.
On pourra dériver la relation cos(nt) = Tn (cos(t)). Donner explicitement
Q0 , Q1 , Q2 . Quel est le degré de Qn ?
Montrer que pour n ≥ 2
Qn (X) = 2XQn−1 (X) − Qn−2 (X).
(d) Montrer que
(1 − X 2 )Qn (X)2 + Tn (X)2 = 1,
Tn0 (X) = nQn (X)
(1 − X 2 )Q0n (X) = XQn (X) − nTn (X).

6. Soit Tn le polynôme de Chebyshev d’ordre n tel que Tn (x) = cos(n arccos x)

pour x ∈ [−1, 1]. Montrer que
1  √ n  √ n 
Tn (x) = x + j 1 − x2 + x − j 1 − x2 , |x| ≤ 1,
2
et en déduire
Tn (x) = ch(n arcch x), x ≥ 1.
 2.6. EXERCICES 27

7. La suite de polynômes de Hermite (Hn )n≥0 est caractérisée par la fonction géné-
ratrice
2
X Hn (t)xn
e2xt−x = ,
n≥0
n!
En dérivant par rapport à t l’égalité précédente, montrer que
Hn0 (t) = 2nHn−1 (t), n ≥ 1.
En dérivant par rapport à x la définition précédente par fonction génératrice,
montrer que
Hn+1 (t) = 2tHn (t) − 2nHn−1 (t), n ≥ 1.
Montrer
  n
t1 + t2 X n!
Hn √ = Hk (t1 )Hn−k (t2 ).
2 k=0
2−n/2 k!(n − k)!
2
En utilisant la formule de Taylor pour la fonction e−x , montrer que
2 dn 2
Hn (t) = et (−1)n n (e−t )
dt
En déduire que
(
si n 6= m,
Z
2 0,
Hn (t)Hm (t)e−t dt = n
√
R 2 n! π, sinon.

8. Soit (Ln )n≥0 (resp. (Hn )n≥0 ) la suite des polynômes de Legendre (resp. Hermite)
caractérisées par leur fonction génératrice
X
(1 − 2xt + x2 )−1/2 = Ln (t)xn ,
n≥0

2
X Hn (t)xn
e2xt−x = .
n≥0
n!
Montrer que
√
Z
2
τ n e−τ Hn (tτ )dτ = πn!Ln (t), t ∈ R.
R
 CHAPITRE 3

FONCTIONS DE BESSEL

3.1. Modélisation dans un domaine circulaire

Divers phénomènes physiques portant sur le disque D = {(x, y) ∈ R2 , x2 +y 2 ≤ 1}, comme
le déplacement vertical u d’une membrane circulaire, la propagation de la chaleur u dans
un disque, sont modélisés par une équation aux dérivées partielles portant sur la fonction
u(x, y, t), (x, y) ∈ D, t ∈ R décrivant l’évolution temporelle de l’observable u :
(30) ∆u = A∂t2 u + B∂t u + Cu
avec une condition sur le bord ∂D = {(x, y) ∈ R2 , x2 + y 2 = 1}, par exemple la condtion
dite de Dirichlet
(31) u(x, y, t) = 0, (x, y) ∈ ∂D.
Dans l’équation (30), A, B sont des constantes et le laplacien ∆ est l’opérateur différentiel
∆ = ∂x2 + ∂y2
en coordonnées cartésiennes, soit
∆ = r−1 ∂r (r∂r ) + r−2 ∂ϕ2
en coordonnées polaires (r, ϕ) avec (x = r cos ϕ, y = r sin ϕ).
Dans le cas de la vibration de la membrane circulaire, on a A = 1 et B = 0, d’où
l’équation des ondes ∆u = ∂t2 u, pour l’évolution de la température on a l’équation de la
chaleur ∆u = ∂t u
Les solutions de l’équation (30) du type R(r)Φ(ϕ)T (t) vérifient
 
d dR d2 Φ d2 T dT
r A + B + CT
dr dr dϕ 2
dt 2 dt
+ =
rR Φr2 T
Les variables r, ϕ, T étant indépendantes, les membres de l’équation précédente sont
constants, soit −K , d’où les deux équations
d2 T dT
(32) A 2
+B + CT = −KT
dt dt 2
d dR dΦ
r
dr dr dϕ2
+ = −K
rR Φr2
où la seconde  
d dR d2 Φ
 
 dr r
dr  = − dϕ
 2
r2  + K
 rR  Φ
 3.2. LES FONCTIONS DE BESSEL 29

se sépare en
   
1 d 2 dR
(33) r r + KR = LR
r dr dr
d2 Φ
(34) + LΦ = 0
dϕ2
Les équations (34) et (32), au contraire de (33), sont résolubles par des fonctions classiques.
En particulier, (34) a des solutions du type Φ(ϕ) = Φ0 einϕ avec n ∈ N, soit L = −n2 . Alors,
l’équation (33) prend alors la forme
d2 R −1 dR 2 −2


+ r + K − n r R=0
dr2 dr
√
soit pour la fonction y(r) = R(r/ K) l’équation différentielle
y 00 + r−1 y 0 + 1 − n2 r−2 y = 0

dite équation de Bessel d’ordre n. √ √

On obtient donc cos(n(ϕ − ϕ0 ))Jn ( Kr) K(t − t0 )) comme solution particulière de
√ cos(−K(t−t
0)
l’équation des ondes, cos(n(ϕ − ϕ0 ))Jn ( Kr)e pour celle de la chaleur. La théorie
générale établit que la solution du problème (30)-(31) peut s’exprimer suivant une série à
termes du type précédent, ainsi pour l’équation de la chaleur A = C = 0, B = 1 dans (30)
une solution générale de la forme
X √
u(r, ϕ, t) = un,p cos(n(ϕ − ϕn,p ))Jn ( K n,p r)e−Kn,p (t−t0 )
n,p

où Kn,p est un zéro de la fonction de Bessel Jn (cf. Section 3.3 ci-dessous) et un,p , ϕn,p des
coefficients appropriés.

3.2. Les fonctions de Bessel

J0
J1
J2 J
3 J
4

Figure 11 . Les fonctions de Bessel Jn , n = 0, 1, 2, 3, 4 sur l’intervalle (0, 25).

Pour n entier, la fonction de Bessel Jn est définie par une série et est solution d’une
équation différentielle (comme la fonction classique et )
Proposition/Définition 3.1. Soit n entier. La série
∞  n+2k
X (−1)k t
(35) Jn (t) =
k=0
k!(n + k)! 2
est absolument convergente pour t ∈ C et est solution de l’équation différentielle
(36) t2 y 00 (t) + ty 0 (t) + (t2 − n2 )y(t) = 0, t ∈ R.
30 CHAPITRE 3. FONCTIONS DE BESSEL

La fonction Jn est appelée fonction de Bessel de première espèce d’ordre n.

Démonstration. — La convergence de la série résulte de la convergence vers 0 du quotient

(|t|/2)2 /((k +1)(n+k +1)) de deux termes consécutifs lorsque k → ∞. La série est dérivable
terme à terme, ce qui permet de vérifier que Jn est solution de (36).
La proposition suivante indique que Jn s’obtient aisément en terme de J0 et J1 (ou J00 )

Proposition 3.1. Les fonctions de Bessel Jn , n ≥ 1 vérifient les relations de récurrence

2n
Jn−1 (t) + Jn+1 (t) = Jn (t),
t
Jn−1 (t) − Jn+1 (t) = 2Jn0 (t).
De plus, J00 = −J1 .

Démonstration. — On vérifie, pour n ≥ 1

d n d −n
[t Jn (t)] = tn Jn−1 (t), [t Jn (t)] = −t−n Jn+1 (t)
dt dt
soit
n n
Jn0 (t) +Jn = Jn−1 (t), Jn0 (t) − Jn = −Jn+1 (t),
t t
qui donnent les relations de la proposition.
Les fonctions de Bessel ont des propriétés pour certaines analogues à celles des polynômes
classiques (Legendre : Th. 2.2, Chebyshev : Th. 2.3, Hermite : Th. 2.4)

Théorème ∗ 3.1. La famille des fonctions de Bessel (Jn )n≥0 est caractérisée par l’une des
deux présentations suivantes
1. Les Jn sont donnés par la fonction génératrice
∞
t
(x−x−1 )
X
(37) J(x, t) = e 2 = J0 (t) + Jn (t)[xn + (−x)−n ], t > 0.
n=1

2. La fonction Jn est donnée par la série entière

∞  n+2k
X (−1)k t
(38) Jn (t) =
k=0
k!(n + k)! 2

De plus, la suite (Jn )n≥0 vérifie la relation de récurrence

tJn+1 (t) = 2nJn (t) − tJn−1 (t)
La fonction Jn a comme équivalents aux bornes de l’intervalle (0, +∞)
 1/2
tn
 
2 1 π
(39) Jn (t) ∼t→0 n , Jn (t) ∼t→+∞ cos t − nπ − .
2 n! πt 2 4
4 Remarque 3.1. En posant, pour n ≥ 0, J−n (t) = (−1)n Jn (t), t ∈ R et en prenant la
variable θ telle que x = ejθ , la relation (37) s’écrit comme un développement en série de
Fourier X
ejt sin θ = Jn (t)ejnθ , θ ∈ R.
n∈Z
Ainsi Z 2π Z π
1 jt sin θ −jnθ 1
Jn (t) = e e dθ = cos(t sin θ − nθ)dθ.
2π 0 π 0
 3.2. LES FONCTIONS DE BESSEL 31

√
Vu que Jn vérifie l’équation différentielle (36), la fonction Un définie par Un (t) = Jn (t) t,
t > 0 vérifie l’équation différentielle
 
00 1
u +u+ − n t−3/2 u = 0.
2
4
√
Ainsi, le terme de droite dans l’asymptotique (39) lorsque t → ∞ vérifie, au facteur t près,
l’équation différentielle précédente où le terme (1/4 − n2 )t−3/2 u a été négligé. 5
En remplaçant dans (35) les factorielles d’entiers `! par les valeurs Γ(` + 1) de la fonction
eulérienne Γ (cf. (47) dans l’appendice), on définit la fonction Jν pour un indice ν non
entier sur R+
∞  ν+2k
X (−1)k t
(40) Jν (t) = , t > 0,
k=0
k!Γ(ν + k + 1)! 2
fonction qui conserve certaines des propriétés des fonctions de Bessel Jn d’ordre n entier,
notamment l’équation différentielle (36) et les relations de récurrence de la Prop. 3.1.
En particulier, pour n = 1/2, U1/2 satisfait à l’équation différentielle u00 + u = 0. En fait,
la définition (35) de J1/2 donne
∞  1/2+2k  1/2 X ∞  1/2
X (−1)k t 2t (−1)k t2k 2
J1/2 (t) = = = sin t,
k=0
Γ(k + 1)Γ(1/2 + k + 1) 2 π k=0
Γ(2k + 2) πt
où on a utilisé la formule de duplication de la fonction Γ (cf. (50) de l’appendice A). L’asymp-
totique pour t → ∞ de (39) est donc une égalité pour J1/2 .
Plus généralement, les séries (40) pour ν demi-entier définissent des fonctions élémentaires
dérivées du sinus cardinal sinc
r  k+1
2 k+1/2 1d
Jk+1/2 (t) = t − (cos t), t > 0.
π t dt
p
soit, si jk est définie par jk (t) = π/(2t)Jk+1/2 (t), t > 0,
 
sin t cos t 3 1 3 cos t
j0 (t) = sinc t, j1 (t) = 2 − , j2 = 3 − sin t − .
t t t t t2
Terminons par une représentation intégrale de Jn .
Proposition 3.2.
π
(t/2)n
Z
(41) Jn (t) = √ cos(t cos θ) sin2n θdθ
πΓ(n + 1/2) 0

Démonstration. — En utilisant l’expression (51) de la fonction beta de l’appendice, on a la formule de

Z 1
1 1
= u2k (1 − u2 )n−1/2 du
Γ(k + n + 1) Γ(k + 1/2)Γ(n + 1/2) −1
qu’on va utiliser dans la représentation de Jn en série
∞ Z 1
X (−1)k (t/2)n+2k 1
Jn (t) = u2k (1 − u2 )n−1/2 du
Γ(k + 1) Γ(k + 1/2)Γ(n + 1/2) −1
k=0
Z 1 ∞
(t/2)n X (−1)k (tu)2k
= (1 − u2 )n−1/2 2k
du
Γ(n + 1/2) −1 2 Γ(k + 1)Γ(k + 1/2)
k=0
Z 1
(t/2)n
= (1 − u2 )n−1/2 cos tudu
Γ(1/2)Γ(n + 1/2) −1
La formule de duplication de la fonction Γ a été utilisée dans la dernière égalité. Cela donne la formule de la
proposition.
32 CHAPITRE 3. FONCTIONS DE BESSEL

Y0
Y1 Y
2 Y
3 Y4

Figure 12 . Les fonctions de Bessel Yn , n = 0, 1, 2, 3, 4 sur l’intervalle (0, 25).

L’équation différentielle (36) a deux solutions indépendantes sur R+ ; la fonction Yn définie

par
Jν (t) cos νπ − J−ν (t)
(42) Yn (t) = lim , t > 0,
ν→n sin νπ
est indépendante de Jn : au contraire de Jn , elle ne se prolonge pas à R, du fait de sa
singularité en t = 0
(43)
n−1  2k−n ∞  n+2k
t 1 X (n − k − 1)! t 1 X (−1)k [ψ(k + 1) + ψ(n + k + 1)] t
Yn (t) = 2Jn (t) log − − ,
2 π k=0 k! 2 π k=0 k!(n + k)! 2

où ψ est la dérivée logarithmique de Γ (cf. l’appendice A).

Les fonctions de Bessel Jn , n ∈ N sont présentes dans maple TM , matlab TM ainsi que dans
la bibliothèque math du compilateur C (utilisé pour le tracé des figures 11 et 12).
#include <math.h>
double jn (int n, double x);
double yn (int n, double x);

3.3. Développements en série de Fourier-Bessel

L’ensemble des zéros des fonctions de Bessel Jν a été bien étudié.

Théorème ∗ 3.2. Soit ν réel positif. La fonction Jν a une infinité de zéros sur R+ , tous
simples. Si ν est un entier n, Jn a pour seuls zéros ses zéros réels, t = 0 étant d’ordre n.

Notons par (xνi )i≥1 la suite croissante des zéros non nuls de Jν .
Pour α réel, la fonction Jν,α définie par Jν,α (t) = Jν (αt), t > 0 est solution de l’équation
différentielle
J00ν,α + t−1 J0ν,α + (α2 − ν 2 t−2 )Jν,α = 0
 3.3. DÉVELOPPEMENTS EN SÉRIE DE FOURIER-BESSEL 33

Soient α+ , α− des réels distincts. En multipliant par α∓ l’équation vérifiée par Jν,α± et
faisant la différence, on obtient
Z 1
1
2 2
(α+ − α− ) tJν (α+ t)Jν (α− t)dt = [t (Jν (α+ t)α− Jν0 (α− t) − α+ Jν0 (α+ t)Jν (α− t))]0
0
soit
1
Jν (α+ )α− Jν0 (α− ) − α+ Jν0 (α+ )Jν (α− )
Z
tJν (α+ t)Jν (α− t)dt = 2 2
0 α+ − α−
Ainsi, si i 6= j ,
Z 1
(44) tJν (xνi t)Jν (xνj t)dt = 0
0
et, en passant à la limite α+ , α− → xνi
Z 1
1 1
(45) tJν2 (xνi t)dt = Jν0 (xνi )2 = Jν+1 (xνi )2
0 2 2
Les identités (44) sont des relations d’orthogonalité : pour un ν positif, la suite de fonc-

Figure 13 . Les fonctions de Bessel J0 (x0n t), t ∈ [0, 1] pour n = 0, . . . , 20.

tions (Jνi )i≥1 avec Jνi (r) = Jν (xνi r), r ∈ (0, 1), est une famille orthogonale dans l’espace
L2 ((0, 1), rdr), dont les normes sont précisées par l’identité (45). On montre que la famille
orthonormée obtenue par normalisation c’est une base hilbertienne, qui donne des dévelop-
pements en série de Fourier-Bessel pour les fonctions définies sur [0, 1], séries au sens L2 et
avec des propriétés de convergence simple sous des hypothèses analogues à celle du théorème
de Dirichlet 2.3 pour les séries de Fourier.
∗
Théorème
R1 3.3. Soit ν ≥ −1/2. Soit f une fonction définie sur (0, 1) telle que
0
r|f (r)| dr soit finie. Alors, dans l’espace L2 ((0, 1), rdr), on a le développement en
2

série
∞ Z 1
X 2
f (t) = uf (u)Jν (xν` u)du Jν (xν` t)
J 2 (xν` ) 0
`=1 ν+1
R1√
Si la fonction f est C 1 par morceaux sur (0, 1) avec l’intégrale 0 t|f (t)|dt finie, alors
∞ Z 1
X 2 f (t+) + f (t−)
(46) 2
uf (u)Jν (xν` u)du Jν (xν` t) =
J (xν` ) 0
`=1 ν+1
2
34 CHAPITRE 3. FONCTIONS DE BESSEL

Figure 14 . Approximations de la fonction f1 définie sur [0, 1] par f1 (t) = t si

t < 1/2 et f1 (t) = t − 1 sinon, par la somme de fonctions de Bessel (46) tronquée à
l’ordre N pour N = 2, 6, 10, 14, 18, 22, 30, 50.

4 Remarque 3.2. Pour ν = 1/2, les séries de Fourier-Bessel de la fonction f √du théorème
précédentp se ramènent à des séries sinusoïdales de Fourier pour la fonction tf , vu que
J1/2 (t) = π/2t sin t, soit
∞ Z
√ X 1 √
tf (t) = 2 uf (u) sin(π`u)du sin(π`t).
`=1 0

3.4. Exercices

1. Si f est à dérivée seconde continue sur [0, 1], on note Lf l’application Lf (t) =
f 00 + t−1 f 0 (t). Montrer que, si f, g sont deux fois dérivables avec f (1) = g(1) = 0
hLf, gi1 = hf, Lgi1
R1
où hp, qi1 = 0 p(t)q(t)tdt.
En déduire que pour λi , i = 1, 2 zéros distincts de la fonction de Bessel J0 et
fi les fonctions définies par fi (t) = J0 (λi t), alors hf1 , f2 i1 = 0.

d p
2. En partant de la définition de Jp comme série entière, montrer que dt
(t Jp (t)) =
tp Jp−1 (t)

3. En dérivant par rapport à x, puis par rapport à t la série génératrice

∞
t
(x−x−1 )
X
e 2 = J0 (t) + Jn (t)[xn + (−x)−n ]
n=1

montrer que
t
nJn (t) = (Jn−1 (t) + Jn+1 (t)) ,
2
1
Jn0 (t) = (Jn−1 (t) − Jn+1 (t)) .
2
 3.4. EXERCICES 35

4. Pour n entier, soit Jn définie par

∞  2k+n
X (−1)k t
Jn (t) = , t ∈ R.
k=0
k!(k + n)! 2
(a) Montrer que, pour n ≥ 1 et t ∈ R,
2n
Jn−1 (t) + Jn+1 (t) = Jn (t),
t
d n
[t Jn (t)] = tn Jn−1 (t).
dt
(b) Soit Γ définie par
Z ∞
dt
Γ(s) = e−t ts .
0 t
Montrer que Γ(s) est bien définie pour s > 0. En intégrant par parties, montrer
que
Γ(s + 1) = sΓ(s), s > 0.
En déduire Γ(n + 1) = n! pour n entier. R
√ 2 √
Montrer que Γ(1/2) = π (on utilisera R e−u du = π ).
(c) Montrer que, pour n entier positif et t réel
∞  2k+n
X (−1)k t
Jn (t) =
k=0
Γ(k + 1)Γ(k + n + 1) 2
(d) Pour t, ν réels positifs, on pose
∞  2k+ν
X (−1)k t
Jν (t) =
k=0
Γ(k + 1)Γ(k + ν + 1) 2
Montrer que r
2
J1/2 (t) = sin(t), t > 0.
πt
5. Soit n un entier relatif et Jn la fonction définie par
1 π
Z
Jn (t) = cos(nu − t sin u)du.
π 0
(a) Montrer que Jn est définie sur R, paire si n est pair et impaire si n est
impair.
Montrer que J−n = (−1)n Jn .
(b) Montrer que pour tout entier p ≥ 0,
1 π
Z
J2p (t) = cos(2pu) cos(t sin u)du,
π 0
1 π
Z
J2p+1 (t) = sin((2p + 1)u) sin(t sin u)du.
π 0
(c) Montrer que Jn est indéfiniment dérivable et que
1 π 1 π
Z Z
0
Jn (t) = sin u sin(nu − t sin u)du = cos u (n − t cos u) cos(nu − t sin u)du,
π 0 π 0
1 π 2
Z
00
Jn (t) = − sin u cos(nu − t sin u)du.
π 0
En déduire que Jn est solution de l’équation différentielle
t2 y 00 (t) + ty 0 (t) + (t2 − n2 )y(t) = 0.
36 CHAPITRE 3. FONCTIONS DE BESSEL

(d) Soit ν réel. Soit y solution sur R de

(EDν ) t2 y 00 (t) + ty 0 (t) + (t2 − ν 2 )y(t) = 0
Montrer que z définie par z(t) = tν y(t), t > 0 est solution de l’équation différen-
tielle
t2 z 00 (t) + (2ν − 1)z 0 (t) + t2 z(t) = 0.
En déduire que y1/2 définie par y1/2 (t) = t−1/2 sin t, t > 0 est solution de l’équa-
tion différentielle (ED 1/2 ) sur R+ .
 APPENDICE A

LA FONCTION GAMMA D’EULER Γ

La fonction Γ (1) d’Euler est définie par

Z ∞
Γ(z) = e−t tz−1 dt, <e z > 0.
0
Elle vérifie
√
(47) Γ(z + 1) = zΓ(z), <e z > 0, Γ(n + 1) = n! si n ∈ N∗ , Γ(1/2) = π,
la factorisation de Weierstrass
∞ 
1 γz
Y z  −z/n
(48) = ze 1+ e ,
Γ(z) n=1
n
la formule des compléments
π
(49) Γ(z)Γ(1 − z) = ,
sin πz
la formule de duplication
22z−1 Γ(z)Γ(z + 1/2)
(50) Γ(2z) = √ ,
π
et la relation avec la fonction beta
Z 1
Γ(x)Γ(y)
(51) B(x, y) = tx−1 (1 − t)y−1 dt = .
0 Γ(x + y)
Sa dérivée logarithmique ψ = Γ0 /Γ a comme valeurs aux entiers
n
X 1
(52) ψ(1) = γ, ψ(n + 1) = −γ + , n ≥ 1.
k=1
k
Dans les formules (48) et (52), γ désigne la constante d’Euler-Mascheroni
X∞
γ= [1/k − log(1 + 1/k)].
k=1

L’intégrale exponentielle se ramène à l’évaluation de Γ en t = 1/2 :

Z ∞ Z ∞ −u
√
Z
−t2 −t2 e
e dt = 2 e dt = 2 √ du = Γ(1/2) = π.
R 0 0 u

(1)
R1
Dans une lettre à Goldbach datée du 13 octobre 1729, Euler introduit la fonction z → 0
(− log u)z−1 du ,
fonction qui fut notée Γ par Legendre en 1809.
 INDEX

base orthonormée, 5 A. Haar, 24

F. Bessel, 6, 28 C. Hermite, 22
boule, 3 D. Hilbert, 1
V. Bunyakovsky, 2 A.-L. Legendre, 37
A.-L. Cauchy, 2, 14 A.-M. Legendre, 13
P. Chebyshev, 18 L. Mascheroni, 37
J. Dirichlet, 16 P. Mengoli, 11
espace de Hilbert, 2 ondelettes, 25
L. Euler, 11, 37 M.-A. Parseval, 6
famille partie
orthogonale, 3 convexe, 3
orthonormée, 3 fermée, 3
E. Fischer, 4 orthogonale, 3
fonction S. Poisson, 31
de Bessel, 17, 29 polynôme
beta, 31, 37 de Chebyshev, 18
de carré intégrable, 3 de Jacobi, 20
exponentielle, 5, 9 de Legendre, 13
de Haar, 22 d’Hermite, 22
d’Hermite, 22 produit scalaire, 1
forme linéaire, 4 projection orthogonale, 4
formule Pythagore, 3
des compléments, 37 J. Rayleigh, iv
de duplication, 31, 37 F. Riesz, 4
de Poisson, 31 L. Schwartz, 10
de Rodrigues, 13 H. Schwarz, 2
J. Fourier, 9 vecteurs orthogonaux, 3
L. Gegenbauer, 20 K. Weierstrass, 37
C. Goldbach, 37