Cours de mathématique

1
 Table des matières

Première partie. Analyse complexe 7

Chapitre 1. Fonctions d’une variable complexe 9
1. Introduction 9
2. Continuité 9
3. Différentiabilité 10
3.1. Conditions de Cauchy-Riemann 10
3.2. Dérivées complexes 10
4. Potentiel et fonction de courant en physique 11
5. La différentiation complexe et l’idée des transformations conformes 12
6. Le théorème de l’application conforme de Riemann 14
7. L’idée des transformations conformes 15
8. L’exemple de la sphère de Riemann 16
Chapitre 2. Intégration 19
1. Intégrales 19
2. Intégrale sur des chemins 20
3. Intégrale de Cauchy 21
4. Conséquences 23
5. Formule intégrale de Cauchy 24
6. Le logarithme 25
7. Réponse linéaire 26
8. Series et développement en série des fonctions analytiques 28
8.1. Domaine de convergence 28
8.2. Convergence uniforme 28
Chapitre 3. Fonctions analytiques 29
1. Fonctions analytiques 29
1.1. Définitions 29
2. Prolongement analytique 29
2.1. Prolongement par la méthode des séries 29
3. Autres théorèmes 31
4. Exemples de fonctions analytiques 35
4.1. Fonctions trigonométriques 35
4.2. Fonction logarithmique 35
4.3. L’idée des coupures et des surfaces de Riemann 36
5. Fonction puissance 37
6. Singularités 37
6.1. Classification des singularités 37
6.2. Théorème du prolongement de Riemann 37
7. Séries de Laurent 38
Chapitre 4. Transformations intégrales 43
1. Séries de Fourier 43
2. Transformée de Fourier 44
3. Transformée de Laplace 45

3
4 TABLE DES MATIÈRES

4. Application : Description hydrodynamique des péhnomènes en matière condensée 46

5. Méthodes du col et de la phase stationnaire 48
5.1. Développement asymptotique 48
5.2. Un exemple 48
5.3. Méthode du col 49
5.4. Méthode de la phase stationnaire 50
Chapitre 5. Résidus 51
1. Calcul des résidus 51
2. Évaluation d’intégrales
R +∞ définies 54
3. Intégrale du type −∞ f (x)dx 54
3.1. Index d’une courbe 54
4. Produits infinis 54
5. Fonction Gamma 55
6. Fonction Zeta de Riemann 56
7. Problèmes 58

Deuxième partie. Fonctions de Green 61

Chapitre 1. Introduction 63
1. Introduction 63
2. L’espace L2 64
Chapitre 2. Distributions 67
1. Solutions élémentaires 69
2. Equations aux dérivées partielles du 1er ordre: Méthode des caractéristiques 69
Chapitre 3. Le problème de Sturm-Liouville 71
1. Existence des fonctions propres (vecteurs propres) 75
1.1. Fonctions de Green 75
1.2. Construction de la fonction de Green 76
1.3. Un problème équivalent 76
2. Fonctions de Green et développement sur une base propre 77
Chapitre 4. Applications des fonctions de Green 79
1. Fonctions de Green pour le problème de Poisson 79
2. Fonctions de Green sur R2 et R3 80
2.1. Si D = R2 80
2.2. Si D = R3 80
3. Fonctions de Green pour des problèmes dépendants du temps 80
Chapitre 5. Polynômes orthogonaux et fonctions de Bessel 81
1. Généralité sur les polynômes orthogonaux 81
2. Polynômes de Legendre 83
3. Polynômes d’Hermite et de Laguerre 84
3.1. Polynômes d’Hermite 84
3.2. Polynômes de Laguerre 84
4. Applications physiques 84
4.1. Équation de Laplace en coordonnées sphériques 84
4.2. Oscillateur harmonique en mécanique quantique 85
5. Fonctions de Bessel 85
5.1. Introduction 85
5.2. Forme intégrale des solutions 86
5.3. La fonction Jν (x) 87
6. Fonction génératrice des fonctions de Bessel 88
 TABLE DES MATIÈRES 5

7. Diffraction de Fraunhofer 89
8. Épilogue 1 90
9. Fonctions de Bessel modifiées 90
10. Épilogue 2 91
Annexe A. Rappels de topologie sur un espace métrique 93
1. Ouverts, fermés, compacts 93
1.1. Espaces métriques 93
1.2. Ouverts, fermés 93
1.3. Adhérence d’un ensemble 93
1.4. Compacts 94
 Première partie

Analyse complexe

7
 CHAPITRE 1

Fonctions d’une variable complexe

Pour cette partie, on pourra consulter les références [?, ?, ?, ?, ?]

1. Introduction
Soit une fonction, pour l’instant supposée uni-valuée, associant à chaque point z du plan
complexe un autre nombre complexe w. Nous notons :
(1) w = f (z)
La notion de fonction étant très générale, nous nous limiterons aux fonctions f (z) données par
la paire de fonctions réelles u(x,y) et v(x,y) avec z = x + iy (i2 = −1), où z est un point d’un
domaine Ω du plan complexe
(2) f (z) = u(x,y) + iv(x,y)
Sans autres hypothèses, on ne peut aller plus en avant. Par contre, si la fonction est f (z)
est continue et si elle est différentiable, cette deuxième propriété impliquant la première,
les conséquences sont importantes et très surprenantes. Voici quelques exemples qui nous
intéresserons par la suite :
(1) Toutes les fonctions différentiables du plan complexe donnent les solutions de l’équation
de Laplace.
(2) Si f (z) est holomorphe sur un domaine Ω, son maximum est sur la frontière de Ω.
(3) Deux fonctions analytiques, c’est-à-dire développables en série entière, qui coïncident
sur un sous-ensemble non discret (segment, voisinage d’un point, etc.) d’un ouvert Ω
coïncident sur Ω. 1 Cette propriété est absolument remarquable. Il suffit de connaître
la fonction dans un voisinage d’un point pour la connaître dans tout le domaine
de définition. On remarque l’analogie avec les polynômes de degré n qui sont eux
parfaitement déterminés par leurs n racines.
Les fonctions holomorphes généralisent les polynômes complexes, mais ce sont des objets
beaucoup plus flexibles : l’inverse d’une fonction holomorphe, son logarithme ou sa racine carrée,
la division de deux fonctions holomorphes etc. sont des fonctions holomorphes alors que ces
opérations ne laissent pas invariant l’ensemble des polynômes à variables complexes.

2. Continuité
Définition 1.
(3) lim f (ζ) existe et lim = f (z)
ζ→z ζ→z

ou de façon équivalente:
Définition 2. Soit z1 ,z2 , . . . ,zn , . . . une suite de points telle que
(4) lim zn = z
n→∞

Alors
(5) lim f (zn ) = f (z)
n→∞

1. Un ouvert Ω et tel que pour chaque point de Ω il existe un disque centré en ce point contenu dans Ω..

9
10 1. FONCTIONS D’UNE VARIABLE COMPLEXE

Il est bien évident que ces définitions impliquent que les fonctions réelles u(x,y) et v(x,y)
soient elles aussi continues 2 .

3. Différentiabilité
Définition 3.
f (ζ) − f (z)
(6) lim existe et est notée f 0 (z)
ζ→z ζ −z
Définition 4. Une fonction qui est définie et différentiable au sens des complexes en tous
points d’une région Ω est dite holomorphe dans cette région. On parle ici de région et cette
notion ne porte que sur les points intérieurs. Dans certaines présentations, on rajoute une
condition : la fonction est holomorphe si sa dérivée est continue. Nous verrons que cela n’a
aucune importance, car la fonction sera nécessairement infiniment dérivable.
3.1. Conditions de Cauchy-Riemann. Les fonctions régulières possèdent une structure
interne surprenante (nous verrons que les fonctions holomorphes sont nécessairement analytiques,
donc C ∞ ). Un point z = x + iy étant donné, il existe deux et seulement deux méthodes pour
approcher z dans (6). Nous pouvons poser
(1) Soit ζ = x +  + y.
(2) Soit ζ = x + i(y + )
et former le quotient dans les deux cas:
f (ζ) − f (z)
(7)
ζ −z
Théorème 1. Soit f (z) = u(x,y) + iv(x,y)
(1) La première méthode donne f 0 (z) = ux + ivx
(2) La deuxième méthode donne f 0 (z) = 1i (uy + ivy )
Mais il faut que le résultat soit unique. D’où le théorème :
Théorème 2. (Cauchy) Si f (z) = u + iv est différentiable en un point, alors les dérivées
partielles obéissent à
∂u ∂v ∂u ∂v
(8) = et =−
∂x ∂y ∂y ∂x
Cette propriété conduit immédiatement au résultat suivant:
Théorème 3. Si f (z) = u(x,y) + iv(x,y) est différentiable, alors u(x,y) et v(x,y) satisfont
tous les deux à l’équation de Laplace
∂ 2ϕ ∂ 2ϕ
(9) + 2 =0
∂x2 ∂y
3.2. Dérivées complexes. Soit f = u + iv une fonction C 1 . Nous pouvons penser f (x,y)
comme une fonction de z et z̄, f (z,z̄). Nous définissons:
 
∂ 1 ∂ ∂
(10) f = −i f
∂z 2 ∂x ∂y
 
∂ 1 ∂ ∂
(11) f = +i f
∂ z̄ 2 ∂x ∂y
Exercice 1. Montrer que la condition ∂f /∂ z̄ = 0 est équivalente à reformuler les conditions
de Cauchy-Riemann.

2. Une fonction dont la kième dérivée est continue est dite de classe C k .
 4. POTENTIEL ET FONCTION DE COURANT EN PHYSIQUE 11

Exercice 2. Vérifiez que chacune des fonctions suivantes est holomorphe dans le domaine
où elle est définie :
z2
(1) f (z) = sin z − z+1 .
3
(2) f (z) = e2z−z − z 2 .
(3) f (z) = zcos z
2 +1 .

(4) f (z) = z (tan z + z).

Exercice 3. Vérifier que les fonctions suivantes ne sont pas holomorphes :
(1) f (z) = |z|4 − |z|2 .
(2) f (z) = z2z̄+1 .
(3) f (z) = z (z̄ 2 − z).
(4) f (z) = z̄(sin z)(cos z̄).
Exercice 4. Vérifier
∂ ∂
(12) z = 1, z=0
∂z ∂ z̄
∂ ∂
(13) z = 0, z̄ = 1
∂ z̄ ∂ z̄
Exercice 5. Vérifier que chacune des fonctions suivantes est harmonique :
(1) f (z) = Re z.
(2) f (z) = x3 − 3xy 2 .
(3) f (z) = |z|2 − 2x2 .
(4) f (z) = ex cos y.
Exercice 6. Vérifier que chacune des fonctions suivantes n’est pas harmonique
(1) f (z)|z|2 .
(2) f (z) = z̄ sin z.
(3) f (z) = f (z) = ez̄ cos z .
Exercice 7. Pour chacune des fonctions réelles et harmoniques suivantes, trouver une
fonction v réelle et harmonique telle que f = u + iv soit holomorphe :
(1) u(z) = ex sin y.
(2) u(z) = 3x2 y − y 3 .
(3) u(z) = e2y sin x cos x.
(4) u(z) = x − y.
Exercice 8. Utilisant ∂/∂x = ∂r/∂x · ∂/∂x + . . ., démontrer que :
   
∂ 1 ∂ sin θ ∂ i ∂ cos θ
(14) = cos θ − − sin θ +
∂z 2 ∂r r ∂θ 2 ∂r r
   
∂ 1 ∂ sin θ ∂ i ∂ cos θ
(15) = cos θ − + sin θ +
∂ z̄ 2 ∂r r ∂θ 2 ∂r r
Exercice 9. Démontrer que les conditions de Cauchy-Riemann en coordonnées polaires
prennent la forme :
(16) r · ur = vθ et r · vr = −uθ

4. Potentiel et fonction de courant en physique

En raison des conditions de Cauchy-Riemann, nous avons la propriété suivante :
∂u ∂v ∂u ∂v
(17) + =0
∂x ∂x ∂y ∂y
12 1. FONCTIONS D’UNE VARIABLE COMPLEXE

qui peut être interprété comme la condition d’orthogonalité des deux vecteurs gradients ∇u·∇v =
0. On en conclut que les isoclines u = const. et v = const. se coupent à angle droit. En
électromagnétisme, u est associé au potentiel électrique et la condition v = cst. donne les lignes
de champ.
Une autre application est donnée par la fonction courant en mécanique des fluides. Nous
considérons un fluide incompressible et non visqueux. Nous sommes intéressés à trouver le
champ de vitesse v(x,y). Cette question n’est pas évidente, car tout dépend des conditions aux
limites, c’est-à-dire que la géométrie du domaine est importante.
Si v est le champ de vitesse d’un fluide irrotationnel, nous pouvons l’écrire comme un
gradient :
vx = ∂x φ
(18)
vy = ∂y φ
où φ définit le potentiel des vitesses. Supposons maintenant le fluide incompressible, alors
div v = 0, et le champ de vitesses peut être écrit comme le rotationel d’une fonction appelée
fonction de courant χ
vx = ∂y χ
(19)
vy = −∂x χ
Les courbes χ = cst sont les lignes de courant, car les particules de fluide se meuvent
tangentiellement aux lignes χ(x,y) = cste. 3 . Les deux représentations, soit par φ, soit par χ,
devant nécessairement être les mêmes, on retrouve les conditions de Cauchy-Riemann pour la
fonction φ(x,y) + iχ(x,y)
(20) ∂x φ = +∂y χ
(21) ∂y φ = −∂x χ
Nous pouvons donc former la fonction complexe
(22) F = φ + iχ
qui est généralement appelée potentiel des vitesses.
En mécanique des fluides, les solution v(x,y) sont données par des conditions aux limites.
La condition la plus naturelle est que la vitesse v(x,y) est parallèle aux frontières du domaine
considéré. Lorsque la géométrie de ce domaine est simple, on peut facilement trouver une
solution. Par exemple, pour une soltuion limitée par un plan, v(x,y) = v0 :
(23) φ(x,y)) = v0 Re(z)
L’outil des transformations conformes que nous allons aborder permet de trouver la solution
dans des géométries beaucoup complexes en procédant par changement de repère. Pour aller de
l’avant, il faut mieux définir ce qu’est une application du plan complexe sur lui-même.

5. La différentiation complexe et l’idée des transformations conformes

Une transformation du plan complexe sur lui-même est dite conforme en un point P si elle
préserve les angles en P et si elle dilate-contracte l’espace de façon uniforme dans toutes les
directions. Une transformation conforme est par essence une transformation locale où le facteur
de dilation et l’angle de rotation sont des fonctions locales du point où la transformation est
effectuée. Nous allons montrer:
Théorème 4. Dans un voisinage infinitésimal d’un point P où f 0 (P ) 6= 0, l’application
z → w = f (z) est la composée d’une rotation d’angle arg |f 0 (P )| et d’une dilatation par un
facteur d’échelle |f 0 (P )|. Les deux application conservant les angles, la propriété sera démontrée.
Lorsque f 0 (P ) = 0, la structure au voisange de P est plus complexe et sera examinée dans le
cadre du prolongement analytique.
3. le vecteur −dχ/dy, dχ/dx est tangent à cette courbe. Prendre un cercle x2 +y 2 = R2 pour s’en convaincre.
 5. LA DIFFéRENTIATION COMPLEXE ET L’IDéE DES TRANSFORMATIONS CONFORMES 13

Fig. 1. Une transformation conforme préserve les angles et revient à une

composition d’une rotation et d’une dilation locales.

Supposons f holomorphe dans un voisinage de P et soient w1,2 deux nombres complexes de

module unité. les deux w1,2 définissent donc des vecteurs unités dans le plan complexe.
Nous définissons les dérivées directionnelles suivant w1,2 par :
f (P + tw1 ) − f (P )
(24) Dw1 f (P ) = lim
t→0 t
f (P + tw2 ) − f (P )
(25) Dw2 f (P ) = lim
t→0 t
On remarque que tout comme f (z), Dw1 et Dw2 sont des nombres complexes. Nous pouvons
identifier les parties réelles et imaginaires aux composantes d’un vecteur. Nous aurons :
(1) |Dw1 f (P )| = |Dw2 f (P )|
(2) Si |f 0 (P )| =
6 0, alors l’angle θ entre w1 et w2 est égal à l’angle entre Dw1 f (P ) et
Dw2 f (P ).

Ces deux propriétés peuvent être démontrées de la manière suivante. Soit f = u + iv. f
définit une application du plan réel sur lui-même (x,y) → (u,v) dont le Jacobien en P est la
matrice
 
ux uy
(26) J(P ) =
vx vy
Cette matrice donne la façon avec laquelle f transforme le plan complexe en lui-même dans un
voisinage local du point P . Dans ce voisinage, l’approximation linéaire de f est satisfaisante et
f peut être représentée par cette matrice. qui peut être ré-écrite en utilisant les conditions de
Cauchy
 
ux uy
(27) J(P ) =
−uy ux
Posons r2 = u2x + vy2 et nous pouvons écrire J(P ) comme le produit d’une matrice orthogonale
de déterminant 1 et d’une dilatation. Posons cos ψ = ux /r sin ψ = −uy /r
   
1 0 cos ψ − sin ψ
(28) J(P ) = r ·
0 1 sin ψ cos ψ
Les matrices orthogonales de déterminant 1 étant associées à des rotations directes dans le sens
horaire, elles préservent les angles tout comme les dilatations.
14 1. FONCTIONS D’UNE VARIABLE COMPLEXE

Note 1. Une application conforme transforme donc un carré élémentaire en un carré

élémentaire. Il s’agit d’un résultat infinitésimal : ce résultat ne peut pas être généralisé à des
pavages grossiers. On verra que tout ouvert simplement connexe du plan complexe peut
être mis en bijection continue avec le disque unité par le théorème de l’application conforme de
Riemann.
Exercice 10. On considère la fonction holomorphe f (z) = z 2 Calculez la dérivée de f
au point P = 1 + i. Calculez le matrice du Jacobien en P comme une matrice 2 × 2. Vérifiez
directement que cette matrice est la composition d’une dilatation et d’une rotation et donner
l’angle de cette rotation.
Exercice 11. Faire le même exercice avec la fonction h(z) = ez au point P = 2 − i.
6. Le théorème de l’application conforme de Riemann
Théorème 5. Soit U un ouvert simplement connexe 4 du plan complexe distinct du plan. Il
existe une unique bijection holomorphe continue entre U et le disque unité D = {z ∈ C : |z| < 1}
(ce qui veut aussi dire que l’application inverse f −1 est aussi continue ). En d’autres mots,
n’importe quel domaine fini de C est topologiquement conforme au disque unité.
Ce théorème a été énoncé pour la première fois par Riemann dans sa thèse. Il semblerait
que la preuve soit beaucoup moins évidente que la preuve de physicien que nous détaillons ici.
Remarquons que la fonction
1 1
(29) − ln z = − [ln|z| + iθ]
2π 2π
a pour partie réelle φ1 (x,y) = −1/2π ln|z|. Elle est solution de l’équation de Laplace avec un
terme source en fonction delta à l’origine w0 = 0
(30) ∆φ = δ(w)
qui donne le potentiel électrostatique en dimension 2 pour une charge électrique unité placée à
l’origine. Remarquons que nous avons volontairement changé de notations en appelant le plan
complexe le plan des w.
Nous pouvons, bien sûr, rajouter une constante pour que φ soit au potentiel 0 sur la frontière
du disque unité. Pour toute fonction φ1 , nous montrerons plus tard que nous pouvons trouver
une fonction conjuguée harmonique φ2 . Supposons que nous puissions trouver la solution du
problème électrique dans l’espace w et appelons Φ(w) la fonction associée. Nous posons
1

(31) φ1 + iφ2 = Φ(w) = − ln zeiα
2π
qui définit une nouvelle variable z. La transformation w → z définie par
(32) z = e−iα e−2πΦ(w)
associe le domaine D (variable w) au disque unité (variable z). En effet, lorsque |z| = 1, Φ(w) = 0,
donc w est sur le bord ∂D et lorsque z → 0, Φ(w) tend bien vers l’infini, donc w tend vers 0,
là où est la charge. La variable α est un angle arbitraire, car nous pouvons toujours faire une
rotation dont l’axe passe par le centre du cercle.
Cette démonstration, en réalité, n’est est pas une. Nous avons supposé que problème électrique
avait toujours une solution, ce qui n’est pas nécessairement le cas pour des frontières suffisamment
irrégulières. Bien qu’il s’agisse d’un résultat « naturel », la démonstration de ce théorème est
loin d’être évidente et nous entrainerait trop loin. Une conséquence immédiate de ce théorème
est que toute question concernant les fonctions holomorphes et analytiques peut être ramenée
soit au disque unité, soit au plan complexe.

4. Simplement connexe : tout parcours fermé peut être réduit à un point sans sortir du domaine. À ne pas
confondre avec connexe dont la définition est : ne peut pas avoir une partition en deux ouverts disjoints. C’est le
cas d’un intervalle.
 Theorem: The interior of any simply connected domain D in C whose bound-
ary consists of more that one point can be mapped conformally one-to-one
and onto the interior DES
7. L’IDéE of theTRANSFORMATIONS
unit circle. It is possibleCONFORMES
to choose an arbitrary 15
interior point w0 of D and map it to the origin, and to take an arbitrary
direction through w0 and make it the direction of the real axis. With these
two choices the mapping is unique.

w z

w0
D f
O

Figure 17.3: The Riemann mapping theorem.

Fig. 2. Le théorème de l’application conforme de Riemann.

This theorem was first stated in Riemann’s PhD thesis in 1851. He re-
garded it as “obvious” for the reason that we will give as a physical “proof.”
Riemann’s argument is not rigorous, however, and it was not until 1912 that
a real proof was obtained by Constantin Carathéodory. A proof that is both
shorter and more in spirit of Riemann’s ideas was given by Leopold Fejér
and Frigyes Riesz in 1922.

Fig. 3. Transformation conforme pour le problème de Dirichlet.

Exemple 1. Méthode des images:

Trouvons l’application qui associe le plan supérieur au disque unité. Il faut donc trouver un
potentiel électrique qui est nul sur l’axe réel. Pour le trouver, considérons un argument de
symétrie. Si une charge ponctuelle de signe + localisée en (0,i) crée un potentiel, nous pouvons
lui associer une charge de signe − localisée en (0, − i). Par symétrie, le potentiel total est
nécessairement nul sur l’axe réel. C’est ce qu’on appelle la méthode des images. Nous avons
donc :
1
φ1 + iφ2 = − [ln(w − i) − ln(w + i)]
(33) 2π
1

= − ln eiα z
2π
d’où
w−i
(34) z = e−iα
w+i
qui marche car si w ∈ R, |z| = 1, et z = 0 si w = i.

7. L’idée des transformations conformes

Il s’agit d’étudier les fonctions holomorphes et bijective h : U → V , où U et V sont
deux ouverts de C. Si il existe une telle application entre U et V , alors ces deux ouverts sont
identiques du point de vue de l’analyse complexe. Sur le plan des applications, ces transformations
permettent de trouver la solution d’équations différentielles aux dérivées partielles.
16 1. FONCTIONS D’UNE VARIABLE COMPLEXE

Exemple 2. Cherchons à trouver la solution du problème de Dirichlet (la solution est donnée
sur les frontières. Pour le problème de Neuman qui obéit au même théorème, c’est la dérivée
normale qui est donnée) :
(35) ∆u = 0 sur Ω et u = f sur ∂Ω
où f est une donnée du problème. Si nous connaissons la transformation h(z) qui associe Ω
au disque unité, nous pouvons résoudre le problème comme suit. Soit u0 (x,y) la solution du
problème sur le disque, alors, comme le bord de Ω a pour image le bord du disque, la solution
cherchée est u(x,y) = u0 ◦ h.

Fig. 4. Effet de z → z + 1/z.

Exercice 12. Considérez l’application:

(36) h : z → w = h(z) = z + 1/z
Cette application envoie le demi-plan supérieur à l’exclusion de la sphère unité sur le demi-plan
supérieur: Pourquoi ? (quelle est l’image du demi-cercle unité ?) En se référant au problème
hydrodynamique vu au début de ce chapitre, donner le potentiel des vitesses et les lignes de
courant.
Rép. si F0 (z) = v0 z, la solution est donnée par
1
(37) F (z) = F0 (h(z)) = v0 (z + )
z
et donc (cf. Figure ??):
1
(38) φ(r,θ) = v0 (r + ) cos θ
r
1
(39) χ(r,θ) = v0 (r − ) sin θ
r

8. L’exemple de la sphère de Riemann

Soit S 2 la sphère de rayon unité dans R3 . Nous considérons l’application S : S 2 \{e3 } → R2
où e3 = (0,0,1). S est définie par
(40) S(x1 ,x2 ,x3 ) = (1 − x3 )−1 (x1 ,x2 )
Cette application associe la sphère (moins le pôle nord) au plan complexe toute entier et est
connue comme la projection stéréographique de S 3 . L’inverse de S est donnée par la formule
1
(41) S −1 (x,y) = (2x,2y,r2 − 1), r2 = x2 + y 2
1 + r2
On montre que S est la seule transformation conforme bijective et continue de S 2 dans C. Par
conséquent, deux courbes ayant une intersection sur S 2 à un angle α ont des images dans le
plan qui se coupent avec un même angle.
 8. L’EXEMPLE DE LA SPHèRE DE RIEMANN 17

Fig. 5. Ligne de courant au voisinage d’un demi-cylindre.

Fig. 6. La projection stéréographique associe la sphère S 2 au plan complexe.

Exercice 13. Considérer les transformations de la figure 7. Vérifiez que les transformations
sont bijectives et que les frontières des domaines sont bien images l’une de l’autre.
Exercice 14. Soit H la partie supérieure du plan, y > 0. Nous désirons trouver la fonction
harmonique sur H telle que
u(x) = c0 si x ∈] − ∞,x1 ]
u(x) = c1 si x ∈]x1 ,x2 [
(42)
...
u(x) = cn si x ∈]xn , + ∞[
(1) Vérifier que la solution du problème est donnée par la partie réelle de
1
(43) cn + {(cn−1 − cn ) ln(z − xn ) + . . . + (c0 − c1 ) ln(z − x1 )}
iπ
(2) Montrer que
1
(44) u(x,y) = cn + {(cn−1 − cn )θn + . . . + (c0 − c1 )θ1 }
π
où on définira les angles θi sur un schéma.
(3) En utilisant ce résultat ainsi que l’un des transformations de la fig. 7 trouver la solution
du problème suivant. Considérer un cercle dont le demi-cercle supérieur est maintenu
au potentiel +1 et le demi-cercle inférieur à φ = 0. Trouver la distribution du potentiel
φ(x,y) à l’intérieur du cercle.
h 2 2i
Rép. φ(x,y) = 1 − π arctan 1−x2y−y
1

Exercice 15. À faire en travaux dirigés:

(1) La fonction f (z) = z 3 /(x2 + y 4 ) a-t-elle une limite quand z → 0?
18 1. FONCTIONS D’UNE VARIABLE COMPLEXE

Fig. 7. Quelques transformations conformes.

(2) Même question pour la fonction f (z) = xy/(x + iy).

(3) En quels points de C les fonctions suivantes sont-elles dérivables?
(a) f (z) = x2 + y 2
(b) f (z) = x2 + ixy
(4) Soit f (z) = u(x) + iv(y). Montrer que f est holomorphe si et seulement si f (z) = az + b
avec a ∈ R et b ∈ C.
(5) Si P et Q, fonctions différentiables de U ouvert de R2 dans R, sont telles que les lignes
P (,y) = constante soient toujours orthogonales aux lignes Q(x,y) = constantes. La
fonction f = P + iQ est-elle forcément holomorphe?
(6) Soit f (z) une fonction holomorphe sur C.
(a) la fonction z → f (z) est-elle holomorphe?
(b) la fonction z → f (z) l’est-elle?
(7) Soit f = P + iQ une fonction holomorphe sur U ouvert connexe. On fait le changement
de variable polaire:
(45) (r,θ) → (x,y) = (r cos θ,r sin θ)
et on désigne pour préciser les notations les fonctions avec un tilde lorsqu’elles sont
considérées comme des fonctions de r et θ
P̃ (r,θ) = P (r cos θ,r sin θ) Q̃(r,θ) = Q(r cos θ,r sin θ)
Montrer que les conditions de Cauchy-Rieman en termes des fonctions P̃ et Q̃ s’écrivent
∂ Q̃ ∂ P̃ ∂ P̃ ∂ Q̃
(46) =r = −r
∂θ ∂r ∂θ ∂r
 CHAPITRE 2

Intégration

1. Intégrales
Définition 5. Soit w = f (z) une fonction continue dans un domaine Ω. Soient z0 et Z
deux points arbitraires de Ω. On définit l’intégrale en se donnant un chemin entièrement contenu
dans Ω. Ce chemin est divisé en n − 1 intervalles z0 ,z1 , . . . ,zn = Z. Sur chaque division [zν−1 ,zν ],
choisir un point ζv et former la somme
ν=n
X
(47) Jn = (zν − zν−1 )f (ζν )
ν=1

La limite
(48) lim Jn = J
n→∞

existe et est indépendante de la division z0 ,z1 , . . . ,zn à condition que la longueur des intervalles
tende vers zéro lorsque n tend vers l’infini. Cependant, la valeur de la limite J dépend du chemin.
Le problème est maintenant de calculer cette intégrale. Nous nous limitons au cas suivant :
Supposons que le chemin soient donné par x(t), y(t) avec t parcourant l’intervalle [α,β]. Si les
deux fonctions ont des dérivées x0 (t), y 0 (t) continues, l’arc est dit rectifiable. Nous pouvons
subdiviser l’intervalle [α,β] comme
α = t0 < t1 < t2 . . . < tn = β
et poser la notation
(49) zν = z(tν ) pour ν = 0, 1, 2, . . . , n
(50) ζν = z(τν ) pour τν ∈ [tν−1 ,tν ], ν = 1, 2, . . . , n
Allégeons l’écriture en posant
(51) u(x(t),y(t)) = ū(t)
(52) v(x(t),y(t)) = v̄(t)
Les sommes Jn deviennent
ν=n
X ν=n
X
(53) Jn = (zν − zν−1 )f (ζν ) = [xν − xν−1 + i(yν − yν−1 )] [ū(τν ) + iv̄(τν )]
ν=1 ν=1

qui, après développement des crochets, se décompose en 4 sommes de type Σ. Lorsque le

refinement n tend vers l’infini, chacune des sommes approche une limite connue. Par exemple
( ν=n ) Z
X β
(54) lim [xν − xν−1 ] ū(τν ) = ū(t)x0 (t)dt
n→∞ α
ν=1

D’où
Théorème 6. Lorsque l’arc Γ est rectifiable avec t ∈ [α,β]
Z Z β
(55) J = f (z)dz = (ū + iv̄)(x0 (t) + iy 0 (t))dt
Γ α

19
20 2. INTéGRATION

Exemple 3.
1
f (z) = , Γ : z(t) = cos t + i sin t, t ∈ [0,2π]
z
Z Z 2π Z 2π
dz − sin t + i cos t
(56) = dt = i dt = 2πi
Γ z 0 sin t + i cos t 0

résultat qui sera constamment utilisé par la suite.

Exemple 4.
f (z) = Re(z); z0 = 0, Z = 1 + i
où l’intégrale est évaluée sur deux chemins :
(1) Suivant la ligne droite z(t) = (1 + i)t. Réponse : J = 12 (1 + i).
(2) Suivant le parcours de 0 à 1 en ligne droite puis de 1 à 1 + i, toujours en ligne droite.
Réponse J = 12 + i.
Conclusion : Deux chemins différents peuvent donc donner deux résultats différents.

2. Intégrale sur des chemins

Il est souvent utile d’associer une fonction continue γ d’un intervalle [a,b] dans R2 ' C à
un chemin. Lorsque γ(a) = γ(b), le chemin est fermé. Pour faciliter l’écriture, nous écrirons
γ(t) = (γ1 (t),γ2 (t)) ou γ(t) = γ1 (t) + iγ2 (t).
Prenons comme exemple le cercle de rayon R. Nous avons x = R cos θ, y = R sin θ Le
changement de notation t = θ est ici évident. On a x0 (t) = −R sin t, y 0 (t) = R cos t d’où
γ(t) = R cos t + i sin t et γ 0 (t) = −R sin t + iR cos t.
Définition 6. Une courbe est de classe C 1 , c’est-à-dire continûment différentiable, si γ1 (t)
et si γ2 (t) le sont toutes les deux. On écrira:
dγ dγ1 dγ2
(57) = +i
dt dt dt
Théorème 7. Si f est de classe C 1 sur un ouvert U contenant la courbe γ, on a la formule
classique d’analyse
Z b  
∂f dγ1 ∂f dγ2
(58) f (γ(b)) − f (γ(a)) = dt +
a ∂x dt ∂y dt
Si la fonction f est holomorphe, cette égalité peut-être réécrite en utilisant les conditions de
Cauchy-Riemann
Z b
∂f dγ
(59) f (γ(b)) − f (γ(a)) = dt (γ(t)) ·
a ∂z dt
et cette intégrale est généralement notée :
I b
∂f
(60) dz (z)
a ∂z
Exercice 16. Pour chacun des cas suivants, calculer l’intégrale complexe sur le chemin γ(t)
(1) f (z) = z z̄ 2 − cos z , γ(t) = cos(2t) + i sin(2t) , 0 ≤ t ≤ π/2.
(2) f (z) = z̄ 2 − sin z , γ(t) = t + it2 , 0 ≤ t ≤ 1.
(3) f (z) = z 3 + z/(z + 1) , γ(t) = et + ie2t , 1 ≤ t ≤ 2.
 3. INTéGRALE DE CAUCHY 21

Fig. 1.
3. Intégrale de Cauchy
Dans l’un des exemples précédents, nous avons montré que l’intégrale d’une fonction suivant
un parcours dépendait a priori du parcours. Pour une fonction continue, on ne peut rien
démontrer de plus. Par contre, si la fonction est holomorphe dans un domaine où l’arc peut être
déformé, l’intégrale entre deux points ne dépend pas du parcours.
Théorème 8. Soit w = f (z) une fonction régulière ou holomorphe dans un domaine Ω
supposé simplement connexe. Alors l’intégrale
Z Z
f (z)dz
z0
est indépendante du chemin d’intégration.
Corollaire 1. Prenons 2 chemins Γ1,2 . Alors
Z Z
(61) =
Γ1 Γ2
Soit Z Z
(62) − =0
Γ1 Γ2

Remarque 1. On remarque que la condition de régularité sur la fonction w = f (z) est

cruciale, car, sans cette condition, les exemples précédents ont montré que l’intégrale dépend du
chemin choisi.
D’où l’énoncé équivalent du théorème de Cauchy :
Théorème 9. Théorème intégral de Cauchy: Soit w = f (z) une fonction régulière dans un
domaine Ω supposé simplement connexe. Alors :
Z
(63) f (z)dz = 0
C
où C est un chemin arbitraire, pas nécessairement simple, mais fermé de Ω.
La démonstration est faite en trois étapes : le chemin est un triangle, le chemin est une ligne
polygonale et, enfin, le chemin est arbitraire.

Partie 1: Cas du triangle. Divisons le triangle en quatre sous-triangles congruents T I , T II , T III , T IV .

L’intégrale revient à faire la somme sur tous les côtés en sommant éventuellement les parcours
22 2. INTéGRATION

effectués en sens contraire (et qui donnent une contribution nulle une fois sommés).
Z Z Z Z Z
(64) = + + +
T TI T II T III T IV
Parmi les quatre triangles, il y en a un dont l’intégrale est supérieure en valeur absolue aux trois
autres. Nous le notons T1
Z  Z 
   
(65)   ≤ 4 
   
T T1
Le sous-triangle T1 peut à son tour est divisé en quatre sous-triangles avec
Z  Z 
   
(66)   ≤ 4 
   
T1 T2
Soit au bout de n sous-divisions Z  Z 
   
(67)  ≤4 
n 
   
T Tn
Nous savons qu’il existe un seul point z0 commun à tous les triangles Ti , i = 0,1, . . . ∞.
Soit alors  > 0. Il existe un voisinage Dδ de z0 tel que
(68) ∀z ∈ Dδ , f (z) = f (z0 ) + (z − z0 )f 0 (z0 ) + η(z − z0 )
(69) avec |η(z − z0 )| ≤ 
Une fois donné , nous pouvons choisir n tel que Ti ⊂ Bδ , i ≥ n.
D’où : Z Z Z Z Z
0 0
(70) f (z)dz = f (0)dz − z0 f (z)dz + zf (z)dz + η(z − z0 )dz
Tn Tn Tn Tn Tn
Z Z
(71) f (z)dz = 0 + 0 + 0 + η(z − z0 )dz
Tn Tn
Prenant la norme, sn étant le périmètre du triangle Tn
Z  Z
 
(72)  f (z)dz  ≤  |z − z0 | dz

Tn Tn
 2
(73) ≤ s
2 n
À chaque itération, le périmètre sn est divisé par 2. D’où
Z 
  1 n s0
(74)   ≤ 4
  2 4n
T
qui peut être aussi petit que désiré. D’où le théorème.

Partie 2: Cas de la ligne polygonale. Il suffit alors de décomposer le chemin dur la ligne en
une somme finie sur des triangles composant la ligne polygonale.

Partie 3: Cas du chemin arbitraire. On prend la limite d’un chemin polygonal où le nombre de
segments tend vers l’infini.
Il existe une autre démonstration utilisant le théorème de Green (que nous verrons dans un
chapitre ultérieur);
Théorème 10. Soit un domaine Ω fini de R2 de frontière (non singulière) δΩ. On suppose
que f (x,y) et g(x,y) sont de classe C 1 dans Ω. Alors
Z Z   Z
∂g ∂f
(75) − dx dy = (f dx + gdy)
Ω ∂x ∂y δΩ
 1
4. CONSéQUENCES 23

82 CHAPTER 3. 3THE CAUCHY THEORY

2

z
z

Figure 3.16: A triply connected domain.

[Link]
2. Le3.12:contour
A curve peut
that can be déformé
être deformed to
ena un
point.
point.

We shall see in Section 4.1.1 that the Cauchy integral formula gives an
easy proof that a holomorphic function is infinitely differentiable. Thus the
Cauchy-Goursat theorem is swept under the rug: holomorphic functions are
1
as smooth as can be, and we can differentiate them at will.
3
2
3.1.5 More General Versions of the Cauchy Formula
A more general version of the Cauchy formula—the one that is typically used
in practice—is this:

THEOREM 1 Let U ⊆ C be a domain. Let γ : [0, 1] → U be a simple,

closed curve that can be continuously deformed to a point inside U . See
Figure 3.12. If f is holomorphic on U and z lies in the region interior to γ,
then !
1 contours
f(ζ) into one.
Fig. 3. . 3Figure 3.17: peuvent
contours Turning
f(z) =three
être considérés
2πi γ ζ − z
dζ . comme 1 seul contour.

Nous appliquons ce théorème

The proof à l’intégrale
is nearly identical surthat
to the one le chemin δΩ :
we have presented above in the
Z
special case. We omit the details.
Z
(76) f (z)dz = f (x,y)(dx + idy)
Example 34 Let δΩ ∈ C : −2 < x < 2, 0 < y < 3} \ D(−1 +
U = {z = x + iy δΩ
Posant g = if , il(7/4)i,
vient1/10).
: Let γ(t) = cos t + i(3/2 + sin t). Then the curve γ lies in U .
Z Z to a point inside
Z certainly be deformed
The curve γ can  U . Thus if f is any
∂f ∂f
(77) f (z)dz = i − dx dy
δΩ Ω ∂x ∂y
qui est nulle, car f (z) étant holomorphe, elle satisfait aux conditions de Cauchy.

4. Conséquences
Théorème 11. Si l’anneau défini par les chemins C1 et C2 est tel que la fonction f (z) y
soit régulière et si C1 et C2 sont orientés dans le même sens alors
Z Z
(78) f (z)dz = f (z)dz
C1 C2

même si la fonction f (z) n’est pas régulière dans l’intérieur de la région définie par C2 .

Théorème 12. Si f (z) est régulière dans une région simplement connexe 1 Ω et si f (z)
possède une primitive F (z), alors
Z z1
(79) f (z)dz = F (z1 ) − F (z0 )
z0
à condition que z0,1 et que le chemin d’intégration appartiennent à Ω.
1. Une région est dite simplement connexe si tout chemin peut être rétrécis à un point (en tirant sur le
lacet !) sans rencontrer d’obstacle, c’est-à-dire de trou.
24 2. INTéGRATION

Exemple 1. D’après ce théorème, la fonction

Z z
dζ
1 ζ
est une fonction régulière de z dans chaque domaine simplement connexe incluant le point +1
mais pas le point 0. Le domaine Ω est donc le demi-plan Re(z) > 0.
Exercice 17. (1) Soit f (z) = z 2 − z et γ(t) = (cos t, sin t), 0 ≤ t ≤ 2π. Vérifiez le
théorème de Cauchy en calculant l’intégrale.
(2) Le théorème de Cauchy ne marche pas pour la fonction f (z) = cot z sur l’anneau
{z ∈ C : 1 < z < 2}. Calculer l’intégrale et vérifiez que l’intégrale n’est pas zéro. Quelle
hypothèse supposée dans le théorème de Cauchy n’est pas vérifiée dans ce cas?

5. Formule intégrale de Cauchy

Le théorème suivant est remarquable. Si un chemin délimite une région où une fonction
est régulière, alors la valeur de cette fonction en tous points de cette région est parfaitement
déterminée par ses valeurs sur le chemin, c’est-à-dire par ses valeurs sur la frontière.
Théorème 13. Soit f (z) une fonction régulière dans un domaine Ω. Pour chaque chemin
orienté C et pour chaque point z ∈ Ω, nous avons :
Z
1 f (ζ)
(80) f (z) = dζ
2πi C ζ − z
si le chemin C et son intérieur sont tous les deux contenus dans Ω.
Nous avons :
Z Z Z
1 f (ζ) 1 f (z) 1 f (ζ) − f (z)
(81) dζ = dζ + dζ
2πi C ζ −z 2πi C ζ −z 2πi C ζ −z
La première intégrale a déjà été calculée :
Z Z
1 f (z) f (z) 1
(82) dζ = dζ = f (z)
2πi C ζ − z 2πi C ζ − z
La deuxième est équivalente à l’intégrale calculée sur un cercle Cz de centre z. Pour chaque
 > 0, nous pouvons choisir le rayon du cercle ρ tel que :
|f (z) − f (ζ)| < 
en raison de la continuité de f (z). Nous avons alors :
Z  Z 
1  f (ζ) − f (z)  1  f (ζ) − f (z) 
(83) dζ  = dζ 
2π  C ζ −z 2π  Cr ζ −z
1 
(84) ≤ 2πρ
2π ρ
qui peut donc être aussi petit que désiré. D’où le théorème.
Théorème 14. Avec les mêmes hypothèses :
Z
(n) n! f (ζ)
(85) f (z) = dζ
2πi C (ζ − z)n+1
Commençons par examiner l’intégrale
Z
1 φ(ζ)
f (z) = dζ
2πi Γ ζ −z
où φ(z) est définie sur le chemin Γ.
 6. LE LOGARITHME 25

Considérons une suite {zn }n≥0 de Ω approchant z mais pour l’instant arbitraire. Nous avons
Z
f (zn ) − f (z) 1 φ(ζ)
(86) = dζ
zn − z 2πi Γ (ζ − z)(ζ − zn )
D’où :
Z Z Z
f (zn ) − f (z) 1 φ(ζ) 1 φ(ζ) 1 φ(ζ)
An =
(87) − 2
dζ = dζ − dζ
zn − z 2πi Γ (ζ − z) 2πi Γ (ζ − z)(ζ − zn ) 2πi Γ (ζ − z)2
Z
zn − z φ(ζ)
(88) = dζ
2πi Γ (ζ − z)2 (ζ − zn )
Soit M une borne supérieure de |φ(z)|. Nous avons :
Z
|zn − z| 1
(89) An ≤ M 2
dζ
2π Γ (ζ − z) (ζ − zn )

où l’intégrale est majorée car, pour n suffisamment grand, aucun des points n’est sur le contour.
La suite {zn }n≥0 tendant vers z, An tend vers zéro.
Comme les dérivées d’ordre supérieur peuvent être traitées de la même manière, le théorème
est démontré.
2
Remarque 2. Nous venons de démontrer qu’un fonction holomorphe est infiniment dérivable,
ce qui est évidemment faux pour des fonctions de variable réelle.

6. Le logarithme
Théorème 15. Soit Ω un région simplement connexe et supposons 0 ∈ / Ω. Alors il existe une
fonction analytique f (z) unique à 2πi près telle que exp[f (z)] = z. Nous écrirons f (z) = ln z et
le choix de la constante arbitraire qui multiplie 2πi est appelé coupure
L’exsitence de ectte fonction nous est assurée en prenant f 0 (z) = 1/z. Choisissons un point
z0 arbitraire. La fonction ln est donc définie par l’intégrale
Z z
(90) ln z = f 0 (z) dz
z0

où l’integration ne dépend du chemin si 0 ∈ / Ω. On ajuste la constante en se fixant une valeur

de f (z0 ). Reste à démontrer que exp[f (z)] = z si z 6= 0. Posons g(z) = exp[f (z)]/z. La fonction
g(z) est analytique dans tout domaine n’incluant pas l’origine.
Nous avons

(91) g 0 (z) = 0

D’où le résultat.
Nous aurons besoin dans la suite de la définition suivante:
Définition 7. Soit c une singularité d’une fonction réelle. On appelle Partie Principale de
l’intégrale de f (x) sur [a,b] avec c ∈ [a,b] la limite
Z c− Z b 
(92) lim f (x)dx + f (x)dx
→0 a c+
Rb
On remarque que l’intégrale a
dx f (x) n’existe pas forcément. Quand elle existe, les deux
définitions coïncident.
 26 2. INTéGRATION

Im(ω’)

ω
Re(ω’)
Fig. 4. Contour utilisé dans (100).
Figure 1: Contour C for Eq. (12 – Kramers-Kronig relation
7. Réponse linéaire
Considérons une quantité mesurable x(t). La réponse de x(t) à un champ F (t) est linéaire :
where C could be any contour confined to a region Z +∞ where χ(ω) is analytic, and not containing
he point z(93)
= ω. Let C denote the contour x(t) =shownχ(t t0 )F (t0 )Because χ vanishes as |z| goes to
in−fig.1.
−∞
nfinity, the large arc0
of C contributes nothing as it recedes to infinity. The remaining part
où χ(t − t ) est une suceptibilité. Nous avons supposé le système invariant sous translation dans
f the contour can: laberéponse
le temps separated intot the
à l’instant ne peutstraight
dépendrepart
quealong the real taxis,
de la différence − t0 oùwhich becomes a
t0 est l’instant
rinciple-value integral
où le champ F (tas
0 theappliqué.
) est small arc
Si F shrinks, and àanunintegral
(t) est assimilé overnous
pic de Dirac, theavons
smallx(t)arc= which
χ(t). is
iφ 0 0
arameterized by z =
Le principe de ωcausalité
+ ϵe .exigeThus − t )(12)
χ(tEq. = 0, tbecomes
< t : La réponse est provoquée par la force et ne
la précède pas.
! ∞ ! 0
la′ réponse du sysème
′ χ(ω )
Nous considérons maintenant iϵeiφ dφ
à un champ alternatif
iφ
F (t) = F0 e−iωt .
D’où : 0=P dω ′ + lim iφ
χ(ω + ϵe ). (13)
−∞ ω −ω ϵ→0 π ϵe
(94) x(t) = χ̃(ω)F0 e−iωt
This becomes the general relation
où la susceptbilité χ̃(ω) a une partie réelle et une partie imaginaire :
! Z ∞ ′
P ∞ ′ χ(ω ) iωt
(95) χ(ω) == χ1 + iχ2dω
χ̃(ω) = .
dt χ(t)e (14)
iπ −∞ ω0 ′ − ω
Le principe de causalité entraîne que l’intégrale est prise sur [0, + ∞[.
inally, separating
Commeinto real and
un champ imaginary parts
est nécessairement une and using
quantité theseule
réelle, results from réelle
la partie Eqs. de
(8,9)
(94)that
1 is even compte
in ω and
: χ2 is odd, this becomes
(96) (ω)F!0 ∞
x(t) = χ12P cos ωt +
ωχ′χ2 (ω)F
(ω ′0)sin ωt
′ 2
χ1en
où le premier terme oscille (ω) = avec F0 cos
phase dωωt, ′2alors que (15)
π 0 ω − ω 2 le deuxième est en quadrature.
Nous montrons que χ1 et χ2 ne sont pas indépendants en raison du principe de causalité.
! ∞
Posons z = x + iy avec x = ω et 2ωP
définissons χ1 (ωle′ )plan complexe supérieur :
(95)′ dans
χ2 (ω) = − Z ∞ dω ′2 2
. (16)
(97) π
χ̃(x + iy) = 0 ω −
−yt ixt ω
dt χ(t)e e , ∀y > 0
0
car χ(t) est nullle pour t <= 0. Cette fonction est clairement holomorphe dans le demi-plan
3 Damped
supérieur.
Harmonic Oscillator
The HarmonicExercice [Link]
oscillator directementexample.
the canonical que les conditions de Cauchyresults
The quantum s’appliquent
are essentially the
(98) dχ
ame as the classical results, so consider the 1 /dx = dχ
Newtonian 2 /dy equation
(99) dχ2 /dx = −dχ1 /dy
mẍ + mẋ/τ + mω02 x = F (t) (17)
where 1/τ is a phenomenological damping rate, and mω02 is the spring constant K. If F (t)
 7. RéPONSE LINéAIRE 27

Le théorème intégral de Cauchy donne que l’intégrale sur le contour C de la figure 4 est nul :
I
χ(z)
(100)
c z −ω

ξ(z) tendant vers 0 sur le grand cercle quand R → ∞, nous avons

Z ω− Z +∞  Z 0 iφ
χ(ω 0 ) 0 χ(ω 0 ) 0 ie dφ
(101) lim 0
dω + 0
dω + lim iφ
χ(ω + eiφ )
→0 −∞ ω − ω ω+ ω − ω →0 π e
qui devient :
Z −∞
1 χ(ω 0 )
(102) χ(ω) = P 0
dω 0
iπ −∞ ω −ω
où P est la partie principale définie par le crochet dans (101). Séparant parties réelles et
imaginaires, et utilisant le fait que χ1 (ω) est paire et χ2 (ω) est impaire, on a les relations de
Kramers-Kronig :
Z +∞
1 ω 0 χ2 (ω 0 )
(103) χ1 (ω) = P dω 0 02
π −∞ ω − ω2
Z +∞
ω χ1 (ω 0 )
(104) χ2 (ω) = − P dω 0 02
π −∞ ω − ω2
Démontrons cette propriété: Par définition de la partie principale et en utilisant le fait que χ̃1 (ω)
est paire et en séprant partie réelle et imaginaire dans l’équation (102) :
1
(105) χ̃1 (ω) = [χ̃1 (ω) + χ̃1 (−ω)]
2
soit
Z −∞ Z −∞ Z −∞
1 χ2 (ω 0 ) 0 1 χ2 (ω 0 ) 0 1 χ2 (ω 0 ) 0
(106) P 0
dω = P 0
dω + P dω
π −∞ ω −ω 2π −∞ ω −ω 2π −∞ ω0 + ω
Mais
1 1 2ω 0
(107) + =
ω0 − ω ω0 + ω ω 02 − ω 2
D’où le résultat.
Exemple 5. Nous considérons l’oscillateur hamrmonique amorti de fréquence propre ω0
(108) mẍ + mẋ/τ + mω02 x = F (t)
Lorsque F (t) n’a qu’une seule composante de Fourier F (t) = F0 e−iωt , alors :
(109) x(t) = χ(ω)F0 e−iωt
avec
−1/m −1/m
(110) χ(ω) = 2
=
ω 2 + iω/τ − ω0 (ω − ω1 )(ω − ω2 )
où
q
(111) ω1,2 = −i/2τ ± ω̃ avec ω̃ = ω02 − 1/4τ 2
Nous pouvons calculer χ(t) en faisant la transformée de Fourier inverse
Z +∞
−1 e−iωt
(112) χ(t) =
2πm −∞ (ω − ω1 )(ω − ω2 )
L’intégrand possède deux pôles, tous eux dans le demin plan complexe inférieur. Nous verrons
que l’intégrale est nécessairement nulle pour t < 0.
28 2. INTéGRATION

8. Series et développement en série des fonctions analytiques

8.1. Domaine de convergence. Considérons une séquence infinie f0 (z), f1 (z), f2 (z), . . . , fn (z), . . .,
où z est un point appartenant aux domaines de définition de toutes les fonctions fn (z). La série
X
fn (z)
n≥0

converge éventuellement. Soit Ω le domaine de convergence : on sait que Ω est l’intérieur d’un
certain cercle dont le centre est un point z0 . Nous avons le théorème :
P
Théorème 16. (Hadamard) Soit la suite des nombres an , n > 0. La série n an (z − z0 )n
converge absolument 2 si |z − z0 | < r et diverge si |z − z0 | > r où r est donné par
1
(113) r = lim p
n→∞ n |a |
n

On remarque que rien n’est dit sur ce qui se passe sur le cercle de convergence lui-même. Tout
est cas d’espèce.
Souvent plus commode, on a la règle de d’Alembert
Théorème 17. (D’Alembert) Le rayon du disque de convergence est déterminé par la limite
 
 an 
(114) r = lim  
n→∞ an+1 

8.2. Convergence uniforme. Il existe Pen fait une convergence plus forte que la convergence
simple d’une série. Supposons que la série an z n converge dans un domaine z ∈ Ω, alors pour
chaque  et pour chaque point z1 de Ω, il existe un n1 = n1 (,z1 ) tel que
|fn+1 (z1 ) + fn+2 (z1 ) + . . . fn+p (z1 )| < , ∀n > n1 , ∀p > 0
On remarque que n1 dépend non seulement de , mais aussi de z1 . D’où la définition suivante :
P
Définition 8. La série an z n converge uniformément si n1 ne dépend que de  et non
de z1 dans tout le domaine Ω.
Cette convergence uniforme nous assure des théorèmes suivants
Théorème 18. Toute série converge de façon uniforme dans un cercle concentrique à son
domaine de convergence et de rayon strictement plus petit.
P
Théorème 19. Considérons la série fn (z) où toutes les fonctions sont régulières dans la
même région simplement connectée Ω.
P
– La série F (z) = fn (z) représente une fonction continue dans Ω.
– Pour chaque chemin Γ
XZ Z
fn (z)dz = F (z)dz
n Γ Γ

– La série est dérivable terme à terme et est égale à F (p) (z).

2. C’est-à-dire si X n
an |z − z0 |
n
converge.
 CHAPITRE 3

Fonctions analytiques

1. Fonctions analytiques
1.1. Définitions. Nous avons (enfin) la définition suivante :
Définition 9. La fonction f (z) est dite développable en série entière en z0 si et si seulement
il existe une série entière de rayon de convergence r > 0 tel que
X
f (z) = an (z − z0 )n
n≥0

Une fonction développable en série entière dans un voisinage est dite analytique.
Théorème 20. Si f (z) est une fonction holomorphe sur un ouvert Ω. Alors f (z) est
nécessairement analytique en tout point z0 de Ω avec
Z
1 f (ζ)
(115) an =
2πi C(z0 ,r) (ζ − z0 )n+1
En particulier
Théorème 21. Si f (z) est une fonction holomorphe sur un ouvert Ω, elle est infiniment
différentiable en tout point de Ω avec
(116) ∀z ∈ Ω, f (n) (z) = an n!
Il y a donc équivalence entre les principes d’holomorphie et d’analycité, une fonction
analytique étant clairement holomorphe. Ce résultat explique une terminologie fluctuante.
Suivant les sources, on parle de façon équivalente de fonctions holomorphes, analytiques,
différentiables au sens des complexes, régulières ou, enfin, conformes. Enfin, on appelle fonction
entière une fonction analytique sur tout le plan complexe.
Exemple 6. Il est évident que tout polynôme est localement la somme de son développement
de Taylor.
Exemple 7. f (z) = 1/z est analytique sur C/ 0, car pour z0 6= 0, Z = z − z0 et
 
1 X (−1)n Z 
n  <1
(117) f (z) = = n+1 (z − z0 ) si  z0 
z0 + Z n≥0
z0

2. Prolongement analytique
2.1. Prolongement par la méthode des séries.
Théorème 22. Soient deux fonctions f1 (z) et f2 (z) régulières dans une région Ω. Nous
supposons que f1 (z) et f2 (z) coïncident dans une région de z0 de Ω. Cette région peut aussi être
un chemin incluant z0 . Alors ces deux fonctions coïncident sur Ω.
Pour démontrer ce théorème, nous aurons besoin du lemme
Lemme 1. On se donne deux séries :
X X
an (z − z0 )n et bn (z − z0 )n
n≥0 n≥0

ayant un rayon de convergence. Si les deux séries coïncident dans un voisinage de z0 sur un
nombre infini de points avec comme point limite z0 , alors les deux séries sont identiques.

29
30 3. FONCTIONS ANALYTIQUES

Fig. 1. Illustration du prolongement anaytique par la méthode des disques de

convergence. Chaque disque a un et un seul domaine d’intersection avec le précédent.
Le théorème sur l’identité des fonctions analytiques nous assure qu’on peut définir
de proche en proche une seule fonction analytique dans le plan complexe. Il peut
arriver que les valeurs trouvées dans le domaine commun au dernier et premier
cercle ne coïncident pas. C’est le cas en particulier si la fonction analytique définie
par le processus est muli-valuée (voir exemple de ln z. La contradiction avec le
théorème d’unicité n’est qu’apparente. La fonction analytique a été définie de
proche en proche sur des cercles qui s’enroulent autour de l’origine. Le plan
complexe en en réalité feuilleté pour éviter les coupures (voir surfaces de Riemann).

Supposons que ce lemme soit vrai: Soit C0 le plus grand cercle de centre z0 contenu dans Ω.
Les deux fonctions étant analytiques, elles sont développables en série et ces séries convergent
dans C0 . D’après le lemme précédent, les deux fonctions sont identiques sur C0 .
Donnons-nous maintenant un point arbitraire ζ de C0 . Les deux points z0 et ζ peuvent être
connectés par un chemin entièrement contenu dans Ω. Divisons ce chemin de façon arbitraire
z0 , z1 , z2 , . . . zm = ζ. Pour chaque zν , nous traçons le cercle Cν entièrement contenu dans Ω. La
division est choisie de telle sorte que chaque cercle contient le centre du suivant.
Les deux fonctions coïncidant dans C0 , elles coïncident en z1 et dans un voisinage de z1 .
Elles sont développables en série en z1 et leur séries sont identiques sur C1 . Les deux fonctions
sont donc identiques sur C1 . Itérant le processus, nous arrivons à zm−1 . Nous avons donc
f1 (ζ) = f2 (ζ).
2
Reste à démontrer le lemme :
(1) On fait z = z0 , d’où a0 = b0 .
(2) Supposons que nous ayons démontré la même propriété jusqu’au rang m, ai = bi , 0 ≤
i ≤ m. Dans le voisinage de z0 , nous avons donc
am+1 + am+2 (z − z0 ) + . . . = bm+1 + bm+2 (z − z0 ) + . . .
Les séries de gauche et de droite représentant des fonctions continues, on peut passer à
la limite z = z0 . D’où am+1 = bm+1
2
Ce théorème peut être interprété comme une extension aux fonctions analytiques du théorème
des polynômes : deux polynômes coïncidant en n + 1 points sont identiques. Ici, il faut que les
fonctions coïncident en un nombre infini de points avec un point d’accumulation.
Remarque 3. Ce théorème n’a de sens que pour les fonctions analytiques et il est en règle
général faux pour des fonctions infiniment dérivables. Il suffit de prendre f (x) = 0, si x ≤ 0 et
2
f (x) = e−1/x si x > 0 qui ne coïncide pas avec la fonction nulle.
 3. AUTRES THéORèMES 31

Définition 10. Soit Ω un domaine et V un ouvert de U . Soit f (z) une fonction définie sur
V. On appelle prolongement analytique de f (z) sur U toute fonction définie sur Ω et coïncidant
avec f (z) sur V.
Conclusion : pour démontrer que deux fonctions analytiques sont identiques, il suffit de le
démontrer sur l’axe réel ou sur un segment [0,1].
Remarque 4. feuillets de Riemann : L’énoncé de ce théorème cache en fait une très
grande subtilité quant au domaine de définition de la fonction analytique. Considérons un
prolongement analytique obtenu par la méthode des séries en définissant de proche en proche
des disques de convergence (voir Fig. 3). Chaque disque de convergence Di , 1 ≤ i ≤ n, a une
intersection non-nulle avec le précédent. Il peut arriver que la valeur de la fonction analytique
définie dans le dernier disque Dn ne coïncide pas avec celle qui a été définie dans le premier
D1 . En apparence, nous pouvons donc avoir deux valeurs pour un même point, ce qui contredit
notre conclusion sur l’unicité du prolongement unique. Dans ce cas, il faut assimiler le plan
complexe à une surface construite de proche en proche suivant un processus géométrique inventé
par Riemann.
Imaginons que nous empilions une série de feuillets 1 ≤ i ≤ n. Chaque feuillet est coupé
suivant une ligne que nous pouvons assimiler à l’axe des réels négatifs. En collant par les bords de
la coupure, x < 0, deux feuillets successifs, nous construisons une surface qui s’enroule comme
une spirale autour de l’origine (voir l’exemple du stationnement). Projetons cette surface spirale
sur le plan et nous avons défini une fonction multi-valuée sur plan. L’apparente discontinuité de
la fonction analytique est donc uniquement due au fait que nous avons oublié le feuilletage de la
surface de Riemann. Ce sera le cas de la fonction logarithme étudiée plus bas.

3. Autres théorèmes
Théorème 23. Si Ω est un domaine ouvert et connexe (c’est-à-dire fait d’une seul tenant).
Si f est holomorphe dans D et si |f (z)| atteint son maximum dans Ω, et non sur la frontière
∂Ω, alors f (z) est une fonction constante. En particulier, une fonction holomorphe atteint sa
borne supérieure quelque part sur le bord δΩ si Ω est fini
(118) sup |f (z)| = sup |f (z)|
z∈Ω̄ z∈δΩ

Supposons qu’il existe z0 ∈ Ω où |f (z)| atteind son maximum. Définissons

(119) O = {ζ ∈ Ω : f (ζ) = f (z0 )}
et, évidemment, z0 ∈ O. La fonction f étant continue, O est fermé (toute suite de O converge
dans 0) . Considérons ζ0 ∈ O.
Rappelons la formule integrale de Cauchy(Cf. Théorème 5). Nous avons :
Z 2π
1
(120) f (z0 ) = f (z0 + ρeiθ )dθ
2π 0
car f est supposée holomorphe dans le disque de rayon ρ centré en z0 .
Exercice 19. Démontrer cette formule à l’aide du théorème 5.
Notons :
ZZ Z 2π Z ρ
(121) f (z) dx dy = f (z0 + seiθ )s ds dθ = πρ2 f (z0 )
Dρ (z0 ) 0 0

qui revient à (c’est une autre forme du théorème de la moyenne):

ZZ
1
(122) f (ζ0 ) = f (z) dxdy
Aρ
32 3. FONCTIONS ANALYTIQUES

Comme |f (ζ0 )| ≥ |f (z)|, nous avons f (ζ0 ) = f (z), ∀z ∈ Dρ (ζ0 ). Tout point de O possède donc
un voisinage dans O. O est donc ouvert. Mais il est aussi fermé. La seule possibilité est O = Ω
et f est donc constante.
Définition 11. Une fonction f (z) est dite entière si elle est définie et holomorphe sur C.
Par exemple, f (z) = sin z, ez etc. sont des fonctions entières. La fonction f (z) = 1/z n’est pas
entière, car elle est singulière en z = 0.
Théorème 24. (Liouville) Si f est une fonction holomorphe sur C et si il existe une
constante M telle que |f (z)| ≤ M, ∀z ∈ C, alors f (z) est une constante.
Introduisons la fonction auxiliaire g(z) = (f (z) − f (0)/z. Même si le dénominateur tend vers
0 quand z → 0, g est holomorphe sur C, car f est C ∞ . On a évidemment g(0) = f 0 (0).
Posons |z| = R. Nous avons :
|f (z)| + |f (0)| 2M
(123) |g(z)| ≤ =
|z| R
Soit w un point du plan et R un nombre réel choisi tel que |w| ≤ R.
Utilisant la formule intégrale de Cauchy :
Z
1 g(z)
(124) g(w) = dz
2πi |z|=R z − w
Mais
Z
1 g(z) 1 sup|g(z)|
(125) | dz| ≤ 2πR
2π |z|=R z − w 2π inf|z − w|
Mais |g(z)| ≤ 2M/R et, pour z sur le cercle, la distance minimale entre z et w est obtenue
quand les deux vecteurs sont colinéaires. D’où inf|z − w| = R − |w|, cf. Fig. 2. D’où
2M
(126) |g(w)| ≤
R − |w|
Mais R peut être aussi grand que voulu. Faisant R → ∞, nous avons g(w) = 0. Revenant à la
définition de g, nous avons f (z) = f (0), ∀z ∈ C. Donc, f est bien constante.

Remarque 5. Ce théorème ne s’applique pas aux fonctions sur la droite réelle: heureusement !
La fonction sin(x) est bornée sur R et elle n’est pas constante.
Nous avons une autre démonstration : Pour z ∈ C, nous avons pour tout R ∈ [0,∞[
Z
0 1 f (z)
(127) f (z) = dζ
2πi ∂DR (z) (ζ − z)2
Nous pouvons paramétriser ∂DR (z) par ζ(t) = z + Reit , 0 ≤ t ≤ 2π. L’intégrale précédente
devient
Z 2π
0 1 f (z + Reit )
f (z) = 2 2it
iReit dt
2πi 0 R e
(128) Z 2π
1
= f (z + Reit )e−it dt
2πR 0
Donc si : ∀z ∈ C, |f (z)| ≤ M :
M
(129) |f 0 (z)| ≤
R
Faisant R → ∞, nous avons f 0 (z) = 0 quelque soit z et donc f est constante.

Ce théorème a pour conséquence :
 3. AUTRES THéORèMES 33

z
Z

Fig. 2. La distance minimale entre z sur le cercle et w est obtenue quand z est
colinéaire avec w.

Théorème 25. (dit théorème fondamental de l’algèbre) Si p(z) est un polynôme de degré
n > 0, alors p(z) s’annule quelque part dans C.
Supposons que p(z) ne s’annule pas. Alors g(z) = 1/p(z) est une fonction entière. Nous
avons aussi |p(z)| → ∞ quand |z| → ∞. Donc |g(z)| → 0 quand |z| → ∞ et g(z) est bornée.
Donc g(z) est constante. D’où la contradiction.
Définition 12. Lorsque qu’une fonction de classe C 2 obéit à l’équation de Laplace
∂ 2f ∂ 2f
(130) + 2 =0
∂x2 ∂y
on dit que f (x,y) est harmonique. On sait que si f = u + iv est holomorphe, alors les parties
réelles et imaginaires u et v sont harmoniques.
Nous avons :
Théorème 26. Si u(x,y) est une fonction réelle et harmonique de classe C 2 dans un
domaine Ω, alors il existe une fonction v(x,y) de classe C 2 telle que la fonction f = u + iv soit
holomorphe.
Exercice 20. Nous démontrons ce théorème.
(1) Nous démontrons d’abord que si Ω est simplement connexe, alors il existe une fonction
analytique f sur Ω telle que u = Re f .
(a) Considérer la fonction g(z) = ∂u/∂x − i∂u/∂y Démontrer que g(z) est analytique.
(b) Nous désirons démontrer le théorème suivant :
Théorème 27. Soit g une fonction analytique sur un domaine simplement connexe
Ω. Alors il existe une fonction f analytique sur Ω telle que f 0 (z) = g(z). La fonction
f est unique à une constante additive près.

Choisissons un point z0 arbitraire dans Ω et définissons

Z z
(131) f (z) = g(z) dz
z0

où le chemin d’intégration est quelconque.

(i) Démontrer que pour tous points z et z0
Z z
f (z) − f (z0 ) 1
(132) − g(z0 ) = (g(w) − g(z0 )) dw
z − z0 z − z0 z0
34 3. FONCTIONS ANALYTIQUES

(ii) la fonction f étant analytique, le chemin d’intégration est arbitraire. En

choisissant comme chemin d’intégration la droite allant de z0 à z, et en
utilisant que g est nécessairement continue, démontrer que
f (z) − f (z0 )
(133) lim = g(z)
z→z0 z − z0
(iii) Attention. La preuve suppose le domaine Ω simplement connexe. À titre
d’exemple, démontrer qu’il n’existe pas de fonction analytique définie sur C\0
telle que f 0 z) = 1/z. La fonction ln z est effectivement une primitive, mais
elle analytique sur C auquel on a retranché la droite Re x ≤ 0.
(2) Utiliser ce résultat pour démontrer le théorème.
Soit Ω un ouvert connexe de C. Si f, Ω → C est holomorphe, alors le théorème intégral de
Cauchy s’applique: Finalement f ∈ C ∞ (Ω) et f 0 (z) est holomorphe Nous avons la réciproque :
Théorème 28. (Morera). Si g(z) est continue dans un ouvert D et si
Z
(134) g(z)dz = 0
∂Γ
pour tout rectangle contenu dans Ω, alors g(z) est une fonction holomorphe.
Soit z0 ∈ Ω donné. Définissons
Z z
(135) f (z) = g(ζ)dζ
z0
où le chemin d’intégration est une suite de droites soit verticales soit horizontales. Nous pouvons
approcher z0 suivnat deux chemins et, d’après les hypothèses, le résultat sera le même.
Si nous approchons z par un chemin qui se termine par une droite horizontale
∂f
(136) = g(z)
∂x
Par contre, si le chemin se termine par uen droite verticale
∂f
(137) = ig(z)
∂y
Donc f est c1 et satisfait aux conditions de Cauchy-Riemann. Par conséquent, f (z) est
holomorphe. Comme g(z) = f 0 (z), g est aussi holomorphe. 
Théorème 29. Soit f (z) la continuation analytique de f (x) où f (x) est une fonction réelle
définie sur l’axe réel ou sur un segment, disons [a,b]. On a (b∗ désigne le conjugué complexe
comme le symbole b̄). Nous avons le principe de réflexion de Schwarz
(138) f (z ∗ ) = [f (z)]∗
Supposons f = u + iv. Alors g = f¯ = u − iv. Donc si f satisfait aux conditions de Cauchy-
Riemann, g aussi. Donc, g est analytique dans le demi-plan inférieur. Reste à démontrer que g
est analytique sur le segment [a,b] !
Nous allons le faire en utilisant le théorème de Morera (cf. Théorème 28): à faire !
Théorème 30. Soit f une fonction holomorphe sur un ouvert Ω. Alors l’ensemble des
points où f s’annule ne peut pas être quelconque. Ce théorème nous rappelle encore une fois
que les fonctions holomorphes se comportent comme des polynômes. Plus précisément : Soit
Z0 = z ∈ Ω : f (z) = 0. Si il existe z0 ∈ Ω et {zj }j≥1 ⊂ Ω\z0 tel que zj → z0 , alors f = 0 sur Ω.
Théorème 31. (Jordan) Soit CR une partie ou le tout d’un arc circulaire de contour Reiθ
où 0 < θ ≤ π. Supposons lim|z|→∞ f (z) = 0 . Soit λ un nombre positif. Considérons l’intégrale :
Z
(139) I(R) = f (z)eiλz dz
CR
 4. EXEMPLES DE FONCTIONS ANALYTIQUES 35

Nous avons : limR→∞ |I(R)| = 0 qui sera constamment utilisé dans les application du théorème
des résidus.
Exercice 21. Nous allons démontrer ce théorème.
(1) Choisir R suffisamment grand pour que |f (z)| < .
(2) Utilisez
2
(140) θ ≤ sin θ θ ∈ [0,π/2]
π
pour majorer cette intégrale (faire un dessin pour démontrer cette propriété).
(3) En déduire le résultat.

4. Exemples de fonctions analytiques

4.1. Fonctions trigonométriques. Les fonctions trigonométriques et hyperboliques sont
définies par leurs séries
X 2n
n z
(141) cos(z) = (−1)
n≥0
(2n)!
X z 2n+1
(142) sin(z) = (−1)n
n≥0
(2n + 1)!
X z 2n
(143) cosh(z) =
n≥0
(2n)!
X z 2n+1
(144) sinh(z) =
n≥0
(2n + 1)!
et la plupart des égalités trigonométriques s’étendent à C par prolongement analytique. Par
exemple sin2 x + cos2 x = 1 sur la droite réelle, la fonction constante 1 étant analytique, cette
égalité est vraie partout.
4.2. Fonction logarithmique. Nous avons déjà démontré que la fonction
Z z
dζ
(145) f (z) =
1 ζ
est une fonction holomorphe pour Re z > 0. Comme le logarithme naturel peut être défini pour
x ∈ R+ \0 par
Z x
dζ
(146) ln x =
1 ζ
la fonction f (z) est l’unique prolongement analytique de ln x dans le domaine complexe. Quel
est alors le domaine de définition de f (z) et quelles sont ses valeurs?
Comme f (z) est toujours défini si le chemin d’intégration évite l’origine, la fonction f (z)
est holomorphe dans C\0. Cependant cette fonction n’est pas uni-valuée, car prenant comme
chemin d’intégration soit le demi-cercle supérieur ou soit le demi-cercle inférieur, on obtient πi
ou −πi. D’après le théorème de Cauchy, ces deux valeurs seront identiques pour n’importe quel
chemin d’intégration contenu ou dans le demi-plan supérieur ou dans le demi-plan inférieur.
Choisissant un chemin d’intégration partant de +1 et se terminant à −1, après avoir fait m
tours dans le sens positif autour de l’origine, nous avons :
(147) ln(−1) = i + 2mπi
et
(148) ln(−1) = i − 2mπi
si le chemin tourne dans le sens négatif.
36 3. FONCTIONS ANALYTIQUES

(a) (b) (c)

Fig. 3. (a) Tracé de la partie imaginaire et réelle, (b), du ln. En (c), la surface
de Riemann du ln.

En conséquence, suivant le chemin, nous avons une infinité de valeur de ln(−1). La fonction
est multi-valuée. Chaque valeur de ln z pouvant être obtenue à partir de l’autre par l’addition
d’un multiple de 2πi. Chacune de ces déterminations constitue une branche de la fonction ln et
chaque branche est bien holomorphe dans n’importe quelle région qui ne contient pas l’origine.
Théorème 32. On considère le plan complexe auquel on a retranché la partie négative de
l’axe réel. En d’autres mots, on considère C\{x + iy|x ≤ 0 et y = 0}. La fonction
(149) ln z = ln|z|+i arg z − π < z < +π
est appelée la branche principale du logarithme. Cette fonction est analytique sur C\{x + iy|x ≤
0 et y = 0} et
d 1
(150) ln z =
dz z
4.3. L’idée des coupures et des surfaces de Riemann. Lorsque la fonction f (z) est
multi-valuée, il existe deux méthodes pour l’étudier :
(1) Introduire des coupures dans le plan complexe.
(2) Remplacer le plan complexe par une surface. À partir de fonction f originalement
définie dans le plan complexe, on définit une nouvelle fonction f˜ sur une surface où la
fonction devient uni-valuée. Cette surface est généralement compliquée dans un espace
de dimension 4.
Commençons par tracer un graphe. Considérons la fonction z = f (w) = ew , w ∈ C. Cette
fonction est périodique, car
f (w + 2kπi) = f (w)
et f (w) ne peut donc être surjective. Par contre, si le domaine de variation est restreint au
domaine
(151) −∞ < Re w < ∞ et − π < Im w ≤ π
la fonction est clairement surjective et on peut définir son inverse, le ln ! D’où l’idée de la coupure
qui est mieux visualisée en traçant les deux surfaces donnant les parties imaginaires et réelles
du ln, c.f. Fig. 3, sous forme paramétrée :
(1) (Re f (w), Im f (w), Im w).
(2) (Re f (w), Im f (w), Re w).
Prenant Im w = 0, la deuxième surface revient à tracer la courbe (ew ,w) qui n’est autre que le
graphe du logarithme donné sous forme paramétrée. Les deux surfaces correspondent donc aux
parties imaginaires et réelles du ln dans le plan complexe. La partie réelle, cas (b), est clairement
 6. SINGULARITéS 37

univaluée et donne un comportement singulier au voisinage de l’origine comme d’habitude.

Quant à elle, la partie imaginaire, cas (a), montre clairement qu’il n’existe pas de prolongement
analytique continu de ln dans le plan complexe. L’axe R− donne une discontinuité apparaissant
comme une coupure. Celle-ci peut être éliminée en collant les différents feuillets correspondants
à −π + 2kπ < Im w ≤ π + 2kpi, k ∈ Z par leur bord. On obtient ainsi une surface de Riemann
qui est l’union de différent feuillets ou “clones” de C où la fonction ln est continue lorsqu’on se
déplace sur la surface.
Exercice
√ 22. √ (1) √ Calculer le logarithme complexe pour les les cas suivants : 3 − 3i, −
3 + i, − 2 − 2i
(2) Expliquer pourquoi le logarithme complexe de 0 n’est pas bien défini.
(3) Calculer : (1 + i)1−i , i1−i , (−3)4−i .

5. Fonction puissance
Soient a et b deux nombres complexes. Nous avons ab = eb ln a et cette fonction est multivaluée
car ln est multivaluée. Nous avons :

Théorème 33. (1) Quelque soit le choix de la branche choisie pour la fonction ln, la
z
fonction z → a est entière, c-à-d. analytique dans le plan complexe, et a pour dérivée
z → ln a az .
(2) Fixons la branche de la fonction ln, par exemple, la branche principale. Alors la fonction
z → z b est analytique dans le domaine de la branche du ln qui a été choisie et sa dérivée
est z → bz b−1 .
Exercice 23. Donner la dérivée des fonctions suivantes en précisant les régions où ces
fonctions sont analytiques:
(1) ez .
(2) sin ez
(3) ez /(z 2 + 3).
√
(4) ez + 1.
(5) cos z̄.
(6) 1/(ez − 1).
(7) ln(ez + 1).

6. Singularités
6.1. Classification des singularités. Il y a trois possibilités pour une fonction f au
voisinage d’un point P :
(1) Si la singularité est effaçable : |f (z)| est bornée dans un disque ayant P pour centre.
(2) Si la fonction f a un pôle : limz→P |f (z)| = +∞.
(3) Si la singularité est essentielle : Aucune de ces possibilités.
6.2. Théorème du prolongement de Riemann.
Théorème 34. Soit f : D(P,r)\P → C un fonction holomorphe et bornée. Alors nous
pouvons effacer la singularité. Celle-ci est dite effaçable :
(1) limz→P f (z) existe.
(2) la fonction fˆ(z) définie par
(152) fˆ(z) = f (z) si z 6= P et lim f (z)sinon
z→P

est holomorphe. Autremement dit, nous pouvons prolonger de façon analytique f en P .

38 3. FONCTIONS ANALYTIQUES

Exercice 24. Discuter les singularités en 0 dans les cas suivants :

(1) z 2 /(1 − cos z).
(2) sin z/z.
(3) (sec z − 1)/ sin2 z
(4) ln(1 + z)/z 2 .
(5) e1/z .
Exercice 25. On constate qu’une solution numérique dans une simulation de mécanique
des fluides possède une signularité effaçable. Quelle est votre interprétation physique? Même
question si cette singularité est essentielle.
Théorème 35. (Casorati-Weierstrass) Soit f : D(P,r0 )\{P } → C une fonction holomorphe
et P une singularité essentielle de f , alors f (D(P,r)\{P } est dense dans C quel que soit
0 < r < r0 .
La preuve est une application du théorème du prolongement de Riemann. Supposons que
f (D(r,P ))\{P } ne soit pas dense dans C. Il exsite donc un disque de centre Q et de rayon s,
D(Q,s), qui n’est pas couvert par f . Considérons la fonction
1
(153) g(z) =
f (z) − Q
Le dénominateur ne peut pas s’annuler. Donc, g(z) est bornée dans un voisinage de P . Par
conséquent, g(z) peut être prolongée analytiquement en P et, comme limz→P f (z) 6= ∞, g(z)
est prolongeable de façon analytique. Donc f (z) = 1/g(z) + Q est prolongeable, ce qui contredit
nos hypothèses.
.
7. Séries de Laurent
Les séries de Laurent généralisent le développement en série puissance pour des fonctions
holomorphes sur un anneau et non sur un disque. Par exemple, le développement de Taylor
2
ne s’applique pas à des fonctions comme 1/z ou 1/zez , car ces fonctions son singulières en
z = 0, mais nous verrons que nous pourrons faire leur développement de Laurent -. Pour de
telles fonctions, il existe une autre façon de faire des développments qui utilisent des puissances
négatives de z.
Soit l’anneau de rayon intérieur r0 et de rayon extérieur r1 :
(154) A = {z ∈ C : r0 < |z − z0 | < r1 }
Soient γj les cercles parcourues dans le sens antihoraire, |z − z0 | = rj de telle sorte que
(155) ∂A = γ1 − γ0
1
Si f ∈ C (Ā), le théorème de Cauchy donne
Z Z
1 f (ζ) 1 f (ζ)
(156) f (z) = −
2πi γ1 ζ − z 2πi γ0 ζ − z
Le point z étant donné, nous avons sur chaque cercle :
 n
1 1 X z − z0
(157) = , ζ ∈ γ1
ζ −z ζ − z0 n≥0 ζ − z0
 m
1 1 X z − z0
(158) = − , ζ ∈ γ0
ζ −z z − z0 m≥0 ζ − z0
qui peuvent être utilisées dans (156) pour donner
n=+∞
X
(159) f (z) = an (z − z0 )n , z ∈ ∂A
n=−∞
 7. SéRIES DE LAURENT 39

avec
Z
1 f (ζ)
(160) an = dζ , n ≥ 0
2πi γ1 (ζ − z0 )n+1
Z
1
(161) an = f (ζ)(ζ − z0 )m dζ , n = −m − 1 < 0
2πi γ0
qui donne le développement de Laurent. Déformant les contours :
Z
1
(162) an = f (ζ)(ζ − z0 )n+1 dζ , n ∈ Z
2πi γ
où γ est n’importe quel cercle parcouru dans le sens anti-horaire et entourant z0 .
Un cas intéressant est r0 = 0.
Proposition 1. Supposons f (z) holomorphe sur le disque DR (z0 )\z0 . Si :
(1) Si il existe un n fini et < 0 tel que am = 0, ∀m < n, alors f (z) a un pôle en z0 .
(2) Dans le cas contraire, la singularité est essentielle.
Théorème 36. Le j-ième coefficient de la série de Laurent d’une fonction f ayant un pôle
d’ordre k en P est
1 dk+j  k

(163) aj = (z − P ) f (z) |z=P
(k + j)! dz k+j
40 3. FONCTIONS ANALYTIQUES

Exercice 26. À faire en classe:

(1) Déterminer la série de Laurent de z/(z 2 + 1) au voisinage de z = i.
(2) Déterminer l’ordre des pôles dans chacun des cas suivants :
(a) cos z/z 2 .
(b) (ez − 1)/z 2
(3) Quelle est la série de Laurent de f (z) = 1/(z − 1) dans les domaines z : |z| < 1 et
z : |z| > 1)?
(4) Même question pour g(z) = 1/(z − 2) dans les domaines z : |z| < 2 et z : |z| > 2.
(5) Considérer la série de Laurent:
+∞
X
z+1/z
(164) e = an z n
n=−∞

(a) Quel est le domaine de convergence de cette série?

(b) Démontrer:
Z
1 π 2 cos t
(165) an = e cos nt dt
π 0
(c) Démontrer:
X 1
(166) ak = a−k =
m≥0
m!(m + k)!

(6) Considérer maintenant la fonction yn (x) définie par:

n=+∞
X
x
(z−1/z)
(167) e 2 = yn (x)z n
n=−∞

Nous verrons plus tard que les fonctions yn (x) sont des fonctions de Bessel Jn (x) et
que la fonction du membre de gauche est appelée fonction génératrice.
(a) Exprimer yn (x) comme une intégrale de Cauchy.
(b) Faire la substitution t = 2u/x et démontrer :
 x n X (−1)k  x 2k
(168) yn (x) =
2 k≥0 (n + k)!k! 2

(c) Choisir le contour t = eiφ et démontrer

Z
1 π
(169) yn (x) = cos(x sin φ − nφ) dφ
π 0
(d) La méthode du col est en général valable pour les intégrales dy type:
Z
(170) I(λ) = eλf (z) dz
C
où λ est grand et positif.
(i) Tracer le graphe de la fonction g(t) = tx e−t .
(ii) Utiliser la méthode du col pour l’intégrale
Z ∞
(171) Γ(x + 1) = g(t) dt x  1
0
et montrer
√
(172) Γ(x + 1) ' 2πxxx e−x
 7. SéRIES DE LAURENT 41

(e) La méthode de la phase stationnaire est en général valable pour les intégrales du
type:
Z
(173) I(λ) = eiλf (z) dz
C
où λ est grand et positif et où f (z) est un nombre réel le long du contour.
(i) Nous considérons l’intégrale d’Airy:
Z
1 ∞ 1
(174) cos( s3 + zs)ds
π 0 3
(ii) Faire le changement de variable s = z 1/2 t et x = z 3/2 . Quelle est alors
l’intégrale à évaluer?
(iii) Montrer que le chemin d’intégration peut être déformé en n’importe quel
chemin tel que:
– Le chemin part dans le secteur 2π/3 < t < π.
– Le chemin se termine dans le secteur 0 < t < π/3.
(iv) Quels sont les points col?
(v) Posant t = x + iy, montrer que le chemin de plus grand gradient est obtenu
pour
(175) x(x2 − 3y 2 + 3) = 0
(vi) On trace sur la figure la cubique solution de cette équation. Les flèches
indiquent les directions où Re h[t] décroît. Quelle branche choisissez-vous
pour faire l’intégration?
(vii) Pour information, le résultat final est
Γ(1/2) 1  3/2

(176) Ai(x) ' exp −2/3x
4π x1/4

Fig. 4.
 CHAPITRE 4

Transformations intégrales

1. Séries de Fourier
Nous nous concentrons sur l’intervalle [0,2π[ en identifiant les point 0 et 2π. Sur le plan
géométrique, notre analyse portera donc sur un cercle. Si f (θ) est une fonction intégrable sur
[0,2π[, nous définissons :
Z 2π
˜ 1
(177) fn = f (t)e−int dt
2π 0
On appelle série de Fourier de f l’expression :
+∞
X
(178) f (t) ∼ f˜n eint
n=−∞

où le signe ∼ signifie converge dans un sens qu’il reste à définir.

Théorème 37. Si f est intégrable sur [0,2π[ et si f est différentiable en t, alors
n=+N
X
(179) f (t) = lim f˜n eint
N →∞
n=−N

Si la dérivée de f est continue, la convergence est absolue et uniforme.

Exemple 8. Nous désirons calculer les coefficients de Fourier de la fonction f (t) = e2i sin t .
Nous avons
1
(180) 2i sin t = eit − it
e
D’où l’idée de considérer la fonction de variable complexe :
(181) ez−1/z = ez · e−1/z
Chaque fonction étant développable en série, nous pouvons développer :
" #" #
X z k X (−1)l z −l
z −1/z
(182) F (z) = e · e =
k≥0
k! l≥0
l!
soit :
n=+∞
" ∞
#
X X (−1)m
(183) F (z) = zn
n=−∞ m=n
m!(m − n)!
et nous avons finalement les coefficients de Fourier de la fonction f :
X∞
˜ (−1)m
(184) fn =
m=n
m!(m − n)!
Exemple 9. Nous considérons ici l’équation de Laplace en 2d sur un disque de rayon R
(185) ∆u = 0
avec les condition aux limites :
(186) u |∂D = f

43
44 4. TRANSFORMATIONS INTéGRALES

où la fonction f (θ) est une donnée du problème. C’est ce qu’on appelle une condition de Dirichlet.
Pour information, la condition de Neumann imposerait de trouver une solution ayant un flux
nul, c’est-à-dire une dérivée nulle, à la frontière du disque.
En 2d le laplacien est (formule qu’il faut connaître) :
∂2 1 ∂ 1 ∂2
(187) ∆= + +
∂r2 r ∂r r2 ∂θ2
et nous cherchons une solution de (185) sous la forme : w(r,θ) = u(r) · v(θ). Nous avons :
1 1
(188) u00 (r)v(θ) + u0 (r)v(θ) + 2 u(r)v 00 (θ) = 0
r r
D’où :
v 00 (θ) r2 u00 (r) + ru0 (r)
(189) − =
v(θ) u(r)
Le membre de gauche ne dépend que de la variable θ. Le membre de droite, quant à lui, ne
dépend que de la variable r. Ces deux variables étant indépendantes, les deux membres sont
constants. Soit λ cette constante :
(190) v 00 (θ) + λv(θ) = 0
(191) r2 u00 (r) + ru0 (r) − λu(r) = 0
On demande que v(θ) soit continue et périodique : On a donc λ = n2 où n est un entier positif.
(192) vn (θ) = an cos nθ + bn sin nθ
où qn et bn sont deux constantes d’intégration. Supposons n > 0. Dans l’équation pour u(r), on
fait le changement de variable r = ez qui transforme cette équation en une équation linéaire à
coefficients constants :
(193) u00 (z) − n2 u(z) = 0
avec la solution u(z) = e±nz lorsque n 6= 0. D’où la solution générale : un (r) = An rn +Bn r−n pour
n > 0. Si n = 0, il faut revenir à l’équation originale, dont la solution est u0 (r) = A0 + B0 ln r.
On a donc la soltution générale
X
(194) w(r,θ) = A0 + B0 ln r + un (r)vn (θ)
n>0
Dans le cas qui nous intéresse, les solutions sont régulières à l’origine. Cela élimine le ln et les
puissances n < 0.
X
(195) w(r,θ) = A0 + rn (An cos nθ + Bn sin nθ)
n>0

Enfin, la fonction f (θ) est supposée régulière. Elle est donc développable en série de Fourier :
1 X
(196) f (θ) = a0 + (an cos nθ + bn sin nθ)
2 n>0

Ce qui permet d’identifier les coefficients entre (195) et (196).

2. Transformée de Fourier
Soit S l’espace vectoriel des fonctions infiniment différentiables f : R → C telles que
limx→±∞ |x|n f (x) → 0, ∀n ∈ N. Nous définissons la transformée de Fourier par
Z +∞
1
(197) f˜(y) = √ e−ixy f (x) dx
2π −∞
lorsque l’intégrale converge pour y ∈ C. Lorsque f (x) est intégrable, sa transformée de Fourier
f˜(y) est continue.
Exemple 10. Si f (x) = e−x alors f˜(y) = e−y .
2 2
 3. TRANSFORMéE DE LAPLACE 45

Exercice 27. Démontrer le résultat précédent en dérivant sous le signe somme donnant
une équation pour f˜0 (y)/f˜(y).
Théorème 38. Nous avons la transformée inverse :
Z +∞
1
(198) f (x) = √ eixy f˜(y) dy
2π −∞
Théorème 39. Nous avons la formule de Poisson :
n=+∞
X n=+∞
1 X ˜
(199) f (x + 2πn) = √ f (n)einx
n=−∞
2π n=−∞
P
Définissons h(x) = +∞n=−∞ f (x + k). La fonction h(x) est clairement périodique de période
1. Elle est donc développable en séries de Fourier. Ses coefficients de Fourier sont an = f˜(n).
D’où le résultat.
2 /4t
Exemple 11. Si f (x) = e−x , nous avons l’identité de Jacobi :
l=+∞
X k=+∞
−πl2 τ 1 X −πk2 /τ
(200) e =√ e
l=−∞
τ k=−∞

3. Transformée de Laplace
Soit f : R+ → C une fonction integrable sur [0,R], ∀R > 0 telle que
Z +∞
(201) |f (t)| e−at dt < ∞, ∀a > A
0

où A est une nombre réel quelconque.

Définition 13. On définit la transformée de Laplace de f par
Z +∞
(202) Lf (s) = f (t)e−st dt, Re s > A
0

Cette intégrale est absolument convergente quelque soit s dans le demi-plan HA = {z ∈ C :

Re z > A} et elle définit une application continue Lf : HA → C.
Théorème 40. Lf (s) est une fonction holomorphe.
Considérons l’intégrale de Lf (s) sur un contour fermé γ. Interchangeant les deux intégrales,
nous avons :
Z Z ∞Z
(203) Lf (s) ds = f (t)e−st dsdt = 0
γ 0 γ

Nous pouvons alors appliquer le théorème de Morera et en conclure le théorème.

La transformation de Laplace permet de trouver la solution de nombreuses équations
différentielles. Nous nous intéressons ici à trouver sa transformée inverse.
Théorème 41. En fait, si B > A, nous avons :
√
(204) Lf (B + iζ) = 2πg̃(ζ), ζ∈R
où g(x) est définie par
(205) g(x) = f (x)e−Bx , x ≥ 0 et g(x) = 0 si x < 0
et où g̃ désigne la transformée de Fiurier de g.
46 4. TRANSFORMATIONS INTéGRALES

Fig. 1. Chemin pour la transformée de Laplace inverse.

Appliquant la transformée de Fourier inverse pour g̃, nous avons :

Z
eBt +∞
f (t) = Lf (B + iζ)eiζt dζ
2π −∞
(206) Z
1
= Lf (s)esζ dζ
2πi γ
où γ est la ligne verticale γ(ζ) = {B + iζ, − ∞ < ζ < +∞} laissant tous les pôles de Lf (s)
sur sa gauche.
On verra que les transformées inverses de Fourier et de Laplace se calculent par la méthode
des résidus.

4. Application : Description hydrodynamique des péhnomènes en matière

condensée
Let us consider a magnetic system with spin variable σ(x). The total magnetization is
X
(207) M= σ(x − xi (t))
i
where the spin on the particle i moves with the particle i. The spin variable is general and can
thought as a density ρ(x,t). Typically, we have the equation of conservation
(208) ∂t M(x,t) = −∇ · j(x,t)
where j(x,t) is the magnetic current. j is simply the magnetic moment times the velocity:
X
(209) j= σ i vi
i

The equation of conservation (208) is equivalent to the condition:

Z
d
(210) dx M(x,t) = 0
dt
Eq. (208) does not give a complete description. We must know the relations between the
current and the magnetization. The traditional way to solve the problem is to postulate a
constitutive equation
(211) j = −D∇M(x,t)
where, in general, D is a phenomenological transport coefficient. Here, it is a simple diffusion
coefficient as for the classical diffusion problem for particles. Using the equation of conservation,
we must solve:
(212) ∂t M − D∇2 M = 0
 4. APPLICATION : DESCRIPTION HYDRODYNAMIQUE DES PéHNOMèNES EN MATIèRE CONDENSéE 47

To solve this equation, we assume some initial magnetization modulation and Fourier transform
the initial condition:
Z
(213) M(k,0) = dkeik·x M̄(k,0)

Without loss of generality, we take:

(214) M̄(k = 0,0) = 0

since, if M (x,t) solve (212), M (x,t) − Cste. is also a solution.

Fourier transform in space the last equation

Z
(215) M̄(k,t) = dxe−ik·x M(x,t)

where we use the “bar” symbol for the spatial Fourier transform. We also Laplace transform the
last equation in time.

Z ∞
(216) M̃(k,z) = dt M̄(k,t)eizt
0

where z is a C number in the half plane Im z ≥ 0.

Because
Z Z
 +∞  
(217) dze ∂t M̄(k,t) = eizt M̄(k,t) t=0 −
izt
dt∂t eizt M̄(k,t)

We get

i
(218) M̃(k,z) = M(k,0)
z + iDk2

Since M(k,0) is given by the initial conditions, the last equation tell us how the initial condition
spreads with time through the system.
From (218), the solution M̃(k,z) is an analytic function in the C-plane but at the the point
z = iDk2 . The pole for z = −iDk2 is characteristic of an hydrodynamic mode. Taking the
inverse Laplace transform, we get:

2
(219) M̄(k,t) = e−Dk t M̄(k,t = 0)

This equation is characteristic of an hydrodynamic mode. It is exponentially damped with a

time scale 1/Dk2 tending to infinity in the long wavelength limit k → 0.
Eq. (219) is special. The system has a priori many channels to relax a perturbation. We
have, however, a conservation law. This conservation law tells us that for a fluctuation to relax,
matter or spin must be transported in space via diffusion. In this way, the perturbation spreads
and relaxes.
48 4. TRANSFORMATIONS INTéGRALES

Fig. 2. An arbitrary initial condition can be decomposed into Fourier modes.

Each mode will evolve independently from the other ones and will spread by
diffusion.

Assuming the initial condition M(x,0) = M0 δ(x), we get the classical result (in 3 dimensions)
M0  
(220) < M(x,t) >= 3/2
exp −x2 /6Dt
(4πDt)

5. Méthodes du col et de la phase stationnaire

5.1. Développement asymptotique. Les développements asymptotiques permettent
d’obtenir des expressions approchées. Par exemple, la formule de Stirling donne le développement
asymptotique de (n − 1)! = Γ(n)
√
(221) Γ(x) ∼ 2πe−x xx−1/2
Définition 14. Une fonction φ(z) est o(1/z n ) si et seulement si limz→∞ z n φ(z) = 0. Si
z n φ(z) est simplement borné, φ(z) est O(z n ). On a : O(z n ) ⇒ o(z n+1 ), mais pas l’inverse.
P
Considérons maintemant la série S = l≥0 al /z l et les sommes partielles
X al
(222) Sn =
0≤l≤n
zl
Si Sn est bien défini pour z 6= 0, nous n’exigeons pas que la série converge. En fait, dans la plus
part des cas, elle diverge.
Définition 15. On dit que f admet S pour développement asymptotique, i.e. f ∼ S, si
(223) f − Sn = o(1/z n )
Pour deux fonctions arbitraires, nous écrivons f (z) ∼ g(z) lorsque nous avons :
 a1 a2 
(224) f (z) ∼ g(z) 1 + + 2 + ...
z z
5.2. Un exemple. Il est souvent important de pouvoir approximer des intégrales (même
celles que l’on sait calculer de façon exacte !). Ces méthodes sont connues sous le nom de méthode
du point de selle, ou méthode du col.
À titre d’exemple, considérons la fonction Γ(x) définie par
Z ∞
(225) Γ(x + 1) = tx e−t dt
0
et cherchons une approximation pour x → +∞.
 5. MéTHODES DU COL ET DE LA PHASE STATIONNAIRE 49

20

15

10

5 10 15 20 25

Fig. 3. Graphe de tx e−t pour x = 5.

Le graphe de la fonction tx e−t est porté en Fig. 3 et fait clairement apparaître un maximum.
La contribution dominante à l’intégrale (225) viendra donc du voisinage de ce maximum, car
c’est dans cette partie que l’aire par “unité de courbe” est la plus concentrée.
Pour trouver le maximum à x donné, on calcule t qui annule la dérivée d(tx e−t )/dt = 0 et
on trouve t = x. Pour faire l’approximation qui nous intéresse, on pose
(226) ef (t) = tx e−t = e−t+x ln t
avec f 0 (t = x) = 0 et f 00 (t = x) = −x/t2 .
Développant f (t) 1 au deuxième ordre :
Z +∞  
1 2
Γ(x + 1) ≈ exp x ln x − x − (t − x)
0 2x
(227) Z +∞  
x ln x−x 1 2
≈e exp (t − x)
−∞ 2x
où, pour x suffisant grand, nous pouvons étendre l’intégrale jusqu’à −∞. L’intégrale sur t est
maintenant une gaussienne 2 . Cela donne :
√
(228) Γ(x + 1) ≈ 2πxxx e−x x → ∞
qui est le premier terme de la formule de Stirling
√  n n  1 1

(229) n! = 2πn 1+ + + ...
e 12n 288n2
5.3. Méthode du col. Soit γ =] − ∞, + ∞[→ C une courbe (γ peut aussi être définie sur
un intervalle). Soit h(ζ) une fonction analytique sur γ. On suppose que :
R
(1) f (z) = γ ezh(ζ) dζ converge absolument.
(2) h0 (ζ0 ) = 0 et h00 (ζ0 ) 6= 0.
(3) Im(zh(ζ)) constant pour ζ sur γ dans un voisinage de ζ0 .
1. Il faut développer la fonction apparaissant dans l’argument de l’exponentielle : le développement de la
fonction exp [f (t)] n’apporterait rien.
R +∞ 2
2. Pour calculer −∞ e−x dx, on fait :
Z +∞ 2 Z +∞ Z +∞ Z +∞ Z +∞
−x2 −x2 −y 2 2 2
e dx = e dx e dy = e−x −y dxdy
−∞ −∞ −∞ −∞ −∞
Passer en coordonnées polaires donne :
Z 2π Z +∞
2
dθ dr re−r = π
0 0
D’où le résultat : Z +∞ √
2
e−x dx = π
−∞
50 4. TRANSFORMATIONS INTéGRALES

(4) Re(zh(ζ)) possède un maximum absolu sur γ en entier.

Alors :
√
2πezh(z)
(230) f (z) ∼ p
−h00 (ζ0 )z
lorsque z → ∞.
Sous les mêmes hypothèses, nous pouvons considérer :
Z
(231) f (z) = ezh(ζ) g(ζ) dζ
γ

où g(ζ) est continue. Alors :

√
2πezh(z) g(ζ0 )
(232) f (z) ∼ p
−h00 (ζ0 )z
lorsque z → ∞.
5.4. Méthode de la phase stationnaire. Nous avons la même approximation pour des
intégrales complexes dans le cas où l’argument tend vers l’infini :
Théorème 42.
Z b
(233) f (z) = eizh(t) g(t) dt
a
√
eizh(ζ0 ) 2π −iπ/4
(234) f (z) ∼ √ p e g(ζ0 )
z −h00 (ζ0 )
lorsque z → ∞. Si h(z) possède un minimum avec h00 (ζ0 ) > 0, cette formule devient
√
eizh(ζ0 ) 2π iπ/4
(235) f (z) ∼ √ p e g(ζ0 )
z h00 (ζ0 )
 CHAPITRE 5

Résidus

1. Calcul des résidus

Donnons-nous f supposée holomorphe dans un domaine ouvert Ω à l’exception des points
pj ∈ Ω supposés en nombre fini. La fonction est alors dite méromorphe can holomorphe partout
sauf en des points isolés où elle a des résidus Au voisinage de chaque pôle sur un disque Dj , f
est développable en série de Laurent
+∞
X
(236) f (z) = an (z − pj )n
n=−∞

Définition 16. Le coefficient a−1 est appelé résidu et noté Respj (f ).

Par définition
Z
1
(237) Respj (f ) = f (z)dz
2πi ∂j D

Le théorème de Cauchy donne donc une formule très simple pour évaluer l’intégrale de f sur le
bord ∂Ω
Z XZ X
(238) f (z)dz = f (z)dz = 2πi Respj (f )
∂Ω j ∂Dj j

qui permet de calculer de très nombreuses intégrales.

La justification de cette propriété est la suivante. Soit γ(θ) un lacet encerclant P . Le théorème
intégrale de Cauchy permet de déformer le parcours en un cercle de rayon r qui peut être
paramétré par µ(θ) = P + reiθ . Ce cercle est à l’intérieur du disque de convergenece de la série
de Laurent de f . On rappelle que dz = ireiθ dθ. Le coefficient −1 de la série de Laurent donne
bien
Z 2π
1 a−1
(239) dz = a−1
2πi 0 z − P
alors que pour n 6= −1
Z Z 2π
n n+1
(240) (z − P ) dz = r i ei(n+1)θ dθ = 0
γ 0

Théorème 43. On appelle ordre du pôle, le dernier entier k tel que (z − P )k f (z) est borné
en P . Autrement dit,
 
 ln |f (z)| 
(241) Ordre du pôle = lim  
z→P ln |z − P | 

Le résidu d’une fonction f ayant un pôle d’ordre k en P est calculé à partir de la formule
 k−1
1 ∂  
(242) Resf (P ) = z − P )k f (z) z=P
(k − 1)! ∂z
Exemple 12. Calculons
Z +∞
dx
(243)
−∞ 1 + x2

51
52 5. RéSIDUS

qui peut être calculée directement (= π!). Les pôles sont simples, car les zéros son simples, et ils
sont localisés en z = ±i. Si les pôles sont simples, le calcul des résidus est immédiat
(244) Resp (f ) = lim (z − p)f (z)
z→p

D’où
1
(245) Resi (f ) =
2i
Soit γR le chemin allant de −R à +R > 1 sur l’axe réel et complété par le demi arc-de-cercle de
rayon R dans le demi-plan supérieur. Ce chemin encercle le pôle +i et par conséquent
Z
dz
(246) 2
= 2πi Resi (f ) = π
γR 1 + z

D’un autre côté, sur le demi-cercle βR de rayon R, nous avons en majorant l’intégrale
Z
dz
(247) lim =0
R→∞ β 1 + z 2
R
car
Z  Z 
 dz   +π Reiθ 
(248)  = dθ
 2  2 2iθ
βR 1 + z 0 1+R e
 
d’où le résultat, car 1 + R2 e2iθ  ≥ |1 − R2 |.
Exercice 28. Montrer que
Z +∞
dx π
(249) = √
−∞ 1 + x4 2
Exemple 13. Nous désirons calculer :
Z +∞
sin x
(250) dx
−∞ x
Cette intégrale converge en 0, car sin x ≈ x, x → 0. Nous considérons la fonction intermédiaire :
eiz
(251) g(z) =
z
Le pôle de cette fonction est en x = 0 et son résidu en 0 est 1. Pour appliquer le théorème des
résidus, il faut choisir un contour. Mais le contour d’intégration doit être tel que ce pôle n’y
réside pas ! Nous formons donc le contour µR = ∪i=1,4 µiR (t) comme sur la Fig. 1 avec
(252) µ1R (t) = t, − R ≤ t ≤ −1/R
(253) µ2R (t) = eit /R, π ≤ t ≤ 2π
(254) µ3R (t) = t, 1/R ≤ t ≤ R
(255) µ4R (t) = Reit , 0 ≤ t ≤ π
et donc :
I XI
(256) g(z) dz = g(z) dz = 2πi
j=1,4 µjR

Nous avons
"I I # Z Z
+∞ +∞
eix eiz
(257) lim g(z) dz + g(z) dz = dx dz =
R→∞ µ1R µ3R −∞ −∞ x z
√ √
Pour le contour µ4R , nous coupons le parcours en deux suivant que Im z > R (A) ou Im z < R
(B).
√
e− R
(1) Sur (A), g(z) ≤ R
.
 1. THE
5.5. APPLICATIONS TO CALCUL DES RéSIDUSOF INTEGRALS
CALCULATION 151 53

curve R

4
R

3
- 1/R 1/R R
-R 1 R
R 2
R

Figure 5.15: The curve µR in Section 5.5.4.

Fig. 1. Contour utilisé pour calculer l’intégrale (250).
Clearly
! 4 !
" √ √
R (5.48) R
g(z) dz = 0≤
g(z) dz
(2) Sur (B), l’arc de cercle parcourt les domaines
j
. θ ≤ arcsin
R
et π − arcsin R
≤
µ µR
θ≤π . j=1

Comme On the one hand, for R > 0,

! Z Z
g(z) dz = 2πiResg (0) · Indµ (0) = 2πi · 1 · 1 = 2πi . (5.49)
(258) µ g(z) dz = g(z) dz
µ4R A∪B
On the other hand,
alors ! ! # ∞
eix
Z g(z) dz + Z g(z) dz → Z dx as R → Z∞ . (5.50)
 1  3 −∞ x
 µR
 µR
(259)  g(z) dz  ≤ |g(z)| dz = |g(z)| dz + |g(z)| dz

Furthermore,
4
µR  µ4R A B
$ $ $ $
avec $! $ $! $ $! $
$ $ $ $ $ $
$ $ $ $ $ $
$Z g(z) dz $ ≤ $ g(z) dz $ + $ g(z)
√ dz $ (5.51)
$ µ4R $ $ $ $ − R $ √
$ µ4
R √ $ $ R √e
µ4 $ − R
(260) |g(z)| dz ≤ ( Longueur de AIm)y ×
Im y < R ≥ R ≤ πe
A R
Z ≡ A+√ B. (5.52)
R
(261) |g(z)| dz ≤ 2 × arcsin
B R

de telle sorte que les deux intégrales tendent vers 0 quand R → +∞. Enfin
I Z 2π it  
eie /R i it
g(z) dz = e dt
µ2R π eit /R R
(262) Z 2π
it /R
=i eie dt → πi, quand R → +∞
π

Nous avons donc

Z +∞
eix
(263) dx + πi = 2πi
−∞ x

D’où en prenant la partie imaginaire

Z +∞
sin x
(264) dx = π
−∞ x
54 5. RéSIDUS

P P
Ind = +1
Ind = +2
P
P

Ind = -3 Ind = -1

Fig. 2. Différents index d’une courbe encerclant un point P .

2. Évaluation d’intégrales définies

R +∞
3. Intégrale du type −∞ f (x)dx
Théorème 44. Si |f (z)| ≤ M/|z|α où α > 1
Z +∞ X
f (x)dx = 2πi (residus dans le plan supérieur)
(265) −∞
X
= −2πi (residus dans le plan supérieur)

3.1. Index d’une courbe. Considérons un arc fermé γ : [a,b] → C supposé de classe C 1
encerclant un point P (pas nécessairement une fois). On suppose que P n’appartient pas à γ.
La quantité
I
1
(266) Indγ (P ) = dz
γ z −P

ne prend que des valeurs entières et mesure le nombre de fois que l’arc entoure le point P , Cf.
Fig. 2. Elle donne aussi le sens. Si l’arc γ entoure k-fois P dans le sens anti-horaire, son index
vaut +k. Si le sens est anti-horaire, son index est −k. Ce nombre est l’index d’une courbe ou
son “winding number” (indice d’un point autour d’un cycle).

4. Produits infinis
Considérons le produit :
k=N
Y
(267) (1 + ak )
k=0

où pour des raisons techniques aucun des ak n’est supposé égal à −1. Que peut-on dire sur la
limite n → ∞?
Théorème 45. Si :
X
(268) |ak | < ∞
k

alors
Y Y
(269) (1 + |ak |) et (1 + ak ) convergent
k k
 5. FONCTION GAMMA 55

10

-4 -2 2 4

-5

-10

Fig. 3. . Fonction Γ(x).

D’où le résultat :

Théorème 46. Soit f : Ω → C une fonction holomorphe telle que

X
(270) |fk (z)|k converge ∀z ∈ Ω
k

alors
Y
(271) g(z) = (1 + fk (z))
k

est holomorphe sur Ω.

5. Fonction Gamma
La fonction Γ(z) intervient dans la transformée de Laplace de f (t) = tz−1 en raison de
définition :
Z ∞
(272) Γ(z) = e−t tz−1 dt, Re z > 0
0

Intégrant par parties

Z ∞
Γ(z + 1) = e−t tz−1 dt
0
Z ∞
(273) d −t z−1
=− e t dt
0 dt
= zΓ(z)

La définition de la fonction Γ(z) montre que Γ(1) = 1 et nous en déduisons

(274) Γ(k) = (k − 1)Γ(k − 1) = (k − 1)!

Si la définition de la fonction Gamma n’est valable que pour Re z > 0, l’identité (273) montre
que Γ(z) est aussi méromorphe dans le dommaine {z ∈ C : z > −1} avec un pôle simple
en −1. Itérant la même identité, on démontre que Γ(z) est en fait définie dans tout le plan
complexe sauf aux entiers négatifs où elle a un pôle simple.

Proposition 2. Pour z ∈ C\Z :

πz
(275) Γ(z)Γ(1 − z) =
sin πz
56 5. RéSIDUS

Il suffit d’établir cette identité pour 0 < Re z < 1 pour qu’elle soit vraie dans le plan
complexe. Dans ce cas
Z ∞Z ∞
Γ(z)Γ(1 − z) = e−(s+t) s−z tz−1 dsdt
Z0 ∞ Z0 ∞
(276) = e−u v z−1 (1 + v)−1 dudv
Z0 ∞ 0
= (1 + v)−1 v z−1 dv
0

après le changement de variable u = s + t et v = t/s. La dernière intégrale est égale à

Z +∞
(277) (1 + ex )−1 exz dx
−∞

qui est une fonction holomorphe dans le domaine 0 < Re z < 1. Pour démontrer la dernière
égalité, il suffit de se concentrer sur la ligne z = 1/2 + iξ, ξ ∈ R
Z +∞ 
x −1 ixξ π
(278) 2 cosh e dx =
−∞ 2 cosh(πξ)
Or
π π
(279) =
sin π( 12 + iξ) cosh πξ
Ce qui achève la démonstration.

6. Fonction Zeta de Riemann

La fonction zeta de Riemann est définie par :
X
(280) ζ(s) = n−s , Re s > 1
n≥1

En particulier
π2 π4
(281) ζ(2) = , ζ(4) =
6 90
Cette fonction intervient dans de nombreaux calculs de physique statistique, voir en particulier
la condensation de Bose-Einstein. Elle intervient aussi en théorie des nombres dû au théorème
suivant :
Théorème 47. Soit {pj : j ≥ 1} = {2,3,5,7,11, . . .} l’ensemble des nombres premiers. Alors,
pour Re s > 1
∞
Y
−1
(282) ζ(s) = 1 − p−s
j
j=1

De façon surprenante, la démonstration est très simple. Nous pouvons développer le membre
de droite de (282)
X X X
(283) 1+ psj + (pj1 pj2 )s + (pj1 pj2 pj3 )s + . . .
j j1 <j2 j1 <j2 <j3

qui est équivalente à la définition (281), car chaque nombre peut être décomposé de façon unique
en produit de facteurs premiers.
Théorème 48. La fonction de Riemann définit une fonction méromorphe dans C avec un
pôle simple en 1.
 6. FONCTION ZETA DE RIEMANN 57

Définissons la fonction
X 2 πt
(284) g(t) = e−n
n≥1

Cette fonction est reliée à la fonction ζ. Nous avons en effet :

Z ∞ X Z ∞
−2s −s
s−1
g(t)t dt = n π e−t ts−1
(285) 0 n≥1 0

= ζ(2s)π −s Γ(s)
Rappelons l’identité de Jacobi, cf. (200), qui implique
1 1
(286) g(t) = − + t−1/2 + t−1/2 g(1/t)
2 2
Nous avons :
Z ∞ Z 1 Z ∞
s−1 s−1
(287) g(t)t dt = g(t)t dt + g(t)ts−1 dt
0 0 1
58 5. RéSIDUS

7. Problèmes

(1) Appliquer le théorème intégral de Cauchy (en le justifiant) pour calculer (γ est un
cercle entourant l’origine)
Z Z
cos z sin z
(288) dz et 2
dz
γ z γ z

(2) On rappelle qu’une coupure dans le plan complexe est une « barrière » qui empêche un
chemin quelconque d’être déformé lorsque les fonctions sont multivaluées.
(a) Démontrer que f (z) = ln(z 2 ) est analytique sur A = {z|z 6= 0 et arg z 6= ±π/2}.
√
(b) On cherche à déterminer la région du plan complexe où la √ fonction z → √ ez + 1
est analytique. On choisira la coupure de la fonction w → w telle que w soit
analytique dans le domaine C\{x + iy|y = 0, x ≤ 0}. Quelle est sa dérivée?
√
(3) Démontrer que la √ branche de la fonction w → w peut être choisie de telle sorte que
2
la fonction z → z − 1 soit analytique dans la région schématisée dans la figure.

(4) Démontrer que :

(a)
Z +∞
dx π
(289) 4
=√
−∞ 1+x 2
(b)
Z +∞
eixζ
(290) dx = πe−|ζ|
−∞ 1 + x2
(5) Nous désirons calculer l’intégrale de Fourier
Z +∞
eixζ
(291) A= dx
−∞ 2 cosh(x/2)

(a) Quels sont les pôles de l’intégrant? Pourquoi ne peut-on pas appliquer le théorème
des résidus directement?
 7. PROBLèMES 59

(b) Pour calculer cette intégrale, nous allons la comparer avec celle prise sur le chemin
γ(x) = x − 2πi. Démontrer que :
 
2πζ eiπζ
(292) A + e A = 2πiRes−iπ
2 cosh(z/2)
En déduire A.
(6) Nous désirons calculer :
Z ∞
xα
(293) B= dx α ∈ [0,1]
0 1 + x2
(a) On pose z α = rα eiαθ 0 ≤ α ≤ 2π. Quelle est la coupure introduite dans le plan
complexe par cette définition?
(b) Soit γR le chemin défini par 1) 0 < r < R, θ = 0 2) r = R, 0 ≤ θ ≤ 2π, 3)
R > r > 0, θ = 2π. Démontrer
Z
zα
(294) lim dz = B − e2πiα B
R→∞ γ 1 + z 2
R

(c) Calculer l’intégrale par le théorème des résidus. En déduire B.

(7) Soit z ∈ C donné. On considère la fonction (appelée fonction génératrice)
(295) f (ζ) = ez(ζ−1/ζ)/2
Démontrer que f (ζ) peut être écrit comme une série
n=+∞
X
(296) f (ζ) = Jn (z)ζ n
n=−∞

où la fonction de Bessel Jn (z) est définie par

Z
1
(297) Jn (z) = ζ −n−1 ez(ζ−1/ζ)/2 dζ
2πi γ
γ étant un cercle de rayon quelconque centré sur l’origine.
(8) En se plaçant sur le cercle unité, déduire de l’expression précédente :
Z
1 +π
(298) Jn (z) = cos(nθ − z sin θ) dθ
π 0
En déduire que Jn (z) est bornée quelque soit z ∈ R.
(9) Dans l’intégrale (297), on fait P le changement
R de variable ζ → 2ζ/z. En déduire en
justifiant l’inversion des signes et
Z
1 X (−1)k  z n+2k
∞
eζ
(299) Jn (z) = n+k+1
dζ
2πi k=0 k! 2 γ ζ

Utilisant le théorème des résidus, démontrer le développement suivant

X∞
(−1)k 1  z n+2k
(300) Jn (z) = n≥0
k=0
k! (n + k)! 2
et
∞
X (−1)(k−n) 1  z 2k−n
(301) Jn (z) = n≤0
k=0
k! (k − n)! 2
(10) En déduire une relation entre Jn (z) et J−n (z).
(11) Démontrer que Jn (z) satisfait à l’équation différentielle :
 
d2 φ 1 dφ n2
(302) + + 1 − 2 φ(z) = 0
dz 2 z dz z
60 5. RéSIDUS

(12) Utilisant le theorème de la phase stationnaire, démontrer le comportement aysmptotique

r 
2 π π
(303) Jn (z) ∼z→∞ cos z − n −
πz 2 4
 Deuxième partie

Fonctions de Green

61
 CHAPITRE 1

Introduction

1. Introduction
Nous sommes intéressés à trouver la solution générale des équations aux dérivées partielles
linéaires. À titre d’exemple, considérons l’équation de la chaleur en une dimension (D est le
coefficient de diffusion) :
∂u ∂ 2u
(304) = D 2 + f (x,t), 0 < x < L
∂t ∂x
avec les conditions :
(305) u(0,t) = u(L,t) = 0 (conditions aux frontières)
(306) u(x,0) = u0 (x) (condition initiale)
où f (x,t) et u0 (x) sont des données du problème. La méthode générale est celle de la séparation
des variables. On cherche u(x,t) sous la forme :
X πx
(307) u(x,t) = an (t) sin(n )
n≥1
L
avec :
Z
P πx 2 L πx
(308) f (x,t) = n≥1 fn (t) sin(n L )et fn (t) = f (x,t) sin(n ) dx
L 0 L
Z L
P πx 2 πx
(309) u0 (x) = n≥1 gn sin(n L ) et gn = u0 (x) sin(n ) dx
L 0 L
d’où
n2 π 2
(310) a0n (t) + an (t) = fn0 (t) et an (0) = gn
L2
La méthode de la variation de la constante permet d’intégrer cette dernière équation :
 Z t 
−n2 tπ 2 /L2 +n2 π 2 s/L2
(311) an (t) = e gn + ds e fn (s)
0
La solution générale est donc :
(312)  Z 
Z t Z L
2 X −n2 π2 /L2 t 2 L πy +n2 π 2 s/L2 πy
u(x,t) = e u0 (y) sin(n ) dy + + ds e f (y,t) sin(n ) dy
L n≥1 L 0 L 0 0 L
P R
qui peut être mis sous la forme en intervertissant les signes et :
Z L Z tZ L
(313) u(x,t) = G(x,t; y,0)u0 (y) dy + G(x,t; y,s)f (y,s) dy ds
0 0 0
où la fonction de Green est définie par :
2 X −n2 π2 /L2 (t−s) nπy nπx
(314) G(x,t; y,s) = e sin( ) sin( )
L n≥1 L L
Par conséquent :
– Si nous connaissons G(. . .), nous pouvons trouver la solution quelques soient u0 et f .
– G(. . .) est une “fonction” de deux paires (x,t) et (y,s).

63
64 1. INTRODUCTION

– G(. . .) s’exprime en fonction des fonctions propres du Laplacien avec les mêmes
conditions aux limites.

2. L’espace L2
Pour toutes fonctions à valeurs complexes et continues sur [a,b], nous pouvons définir le
produit scalaire
Z b
(315) ∀f,g ∈ C([a,b]) (f,g) = f (x)g(x)dx
a
ainsi que la norme
p
(316) kf k = (f,f )
Mais cet espace C([a,b]) n’est pas nécessairement un espace approprié pour la « suite des
opérations », car il n’est pas fermé : une suite de fonctions continues peut ne pas converger
dans C([a,b]), alors qu’elle converge dans un sens que nous définirons plus bas. Là encore, il y
a d’autres difficultés. Si deux focntions f, g sont intégrables au sens de Riemann, leur produit
f × g ne l’est pas nécessairement. Cette propriété n’est vraie que si l’intégrale est prise dans
un sens plus large (intégrale de Lebesgue). Comme les fonctions qui nous intéresserons seront
toutes intégrables au sens de Riemann, nous oublierons ce problème 1 .
Définition 17. On appelle L2 (a,b) l’espace vectoriel des fonctions telles que
Z b
(317) |f (x)|2 dx < ∞
a
Comme
(318) kαf + βgk ≤ kαf k + kβgk = |α|kf k + |β| kgk ∀α, β ∈ C, ∀f, g ∈ L2 (a,b)
αf + βg est dans L2 (a,b) si f et g sont toutes les deux dans L2 (a,b).
Remarque 6. Dans L2 (a,b), l’égalité kf k = 0 n’implique pas f (x) = 0 sur tous les points
de l’intervalle [a,b]. La fonction nulle sauf en un nombre fini de points a une norme nulle. La
fonction nulle sur L2 (a,b) définit donc une classe d’équivalence dont la fonction strictement
nulle est la seule représentante qui soit continue. De même, nous dirons que f et g g sont égales
dans L2 (a,b) si kf − gk = 0. Cela ne veut pas dire que f et g sont partout identiques, mais
qu’elles sont simplement égales presque partout.
Définition 18. Une suite fn de fonction dans L2 (a,b) converge si il existe une fonction f
de L2 (a,b) telle
(319) lim kfn − f k = 0
n→∞

Définition 19. Un suite fn de fonctions de L2 (a,b) est une suite de Cauchy si et seulement
si
(320) ∀ > 0 ∃N ∈ N tel que ∀m,n > N kfn − fm k < 
Théorème 49. L2 (a,b) est complet. Autrement dit, toute suite de Cauchy converge dans
L2 (a,b).
Exercice 29. (1) Considérons la suite de fonctions {fn (x)}n∈N définie par :
(321) fn (x) = xn x ∈ [0,1]
(a) Étudier la convergence presque partout sur [0,1].
(b) Démontrer que xn −→ L2 0.
1. «La plus grande jouissance dans l’existence, c’est de vivre dangereusement ...», F. Nietzsche.
 2. L’ESPACE L2 65

(2) Considérer la suite de fonctions



0, x = 0
(322) fn (x) = n, 0 < x < 1/n

0,1/n < x ≤ 1

(a) Déterminer si cette suite converge presque partout sur [0,1].

(b) Démontrer que cette suite ne converge pas dans L2 .
Conclusion : il n’y a pas équivalence entre les deux types de convergence.
Un théorème utile est le suivant (test de Weierstrass)
Théorème 50. Soit fn une séquence de fonctions sur un intervalle I. On suppose qu’il
existe une séquence de nombre réels positifs tels que
(323) |fn (x)| < Mn
Si
X
(324) Mn
converge, alors
X
(325) fn (x)
converge de façon uniforme ate absolument sur I.
 CHAPITRE 2

Distributions

Définition 20. Une distribution est une fonctionnelle continue sur l’espace des fonctions
test noté C0∞ (R). Les fonctions test sont des fonctions infiniment différentiables sur R et nulles
en-dehors d’un intervalle [−a,a], où a n’est pas précisé. Autrement dit, une distribution associe
un nombre réel à toutes fonctions test.
Exemple 14. Un exemple de fonction de test est :
2 −x2 )
(326) φ(x) = e−1/(a si − a < x < a et 0 sinon
Pour une distribution f , on note généralement
Z
(327) f (ϕ) =< f,φ >= f (x)φ(x) dx

bien que l’intégrale ne soit pas une “vraie” intégrale, car f n’est pas une fonction ! f est en fait
l’abstraction mathématique d’une densité.
D’autre part, cette application doit être linéaire
(328) f (aφ + bψ) = af (φ) + bf (ψ)
pour tous nombres réels (a,b).
Exemple 15. Un premier exemple de distribution est donné par n’importe quelle fonction
f (x) qui peut être intégrée. L’opération
Z
(329) f (φ) = f (x)φ(x) dx

est effectivement une opération linéaire. Les distributions incluent donc toutes les fonctions
intégrables, d’où le nom de fonctions généralisées.
Exemple 16. La distribution Delta de Dirac notée δ(z) ou δ(z − x) (distribution centrée en
z = x) n’est pas une fonction. Son action sur l’espace des fonctions est
Z
δ(φ) = φ(0) = δ(z)φ(z) dz
(330) Z
δx (φ) = φ(x) = δ(z − x)φ(z) dz

-1.5 -1.0 -0.5 0.5 1.0 1.5

√ 2
Fig. 1. Graphes de G(j,x) = j/ πe−jx pour j = 2,3,5,7.

67
68 2. DISTRIBUTIONS

Il est possible de généraliser aux distributions les opérations que l’on fait d’habitude sur des
fonctions :
(1) Différentiation
Si f est une fonction usuelle, une intégration par parties montre que
Z Z
(331) < f ,φ >= f (x)φ(x)dx = − f (x)φ0 (x) dx = − < f 0 ,φ >
0 0

Pour une distribution, l’accroissement relatif f 0 (x) = ∆f /∆x n’a aucun sens, mais on
peut définir une dérivée par
(332) < f 0 ,φ >= − < f,φ0 >
et
(333) < f (n) ,φ >= (−1)n < f,φ(n) >
À titre d’exemple, considérons la distribution de Heavyside définie par
(334) H(x) = 1 si x > 0 et 0 sinon
H(x) est une fonction usuelle, mais elle n’est pas dérivable au sens usuel. Mais sa
dérivée au sens des distributions est connue
Z +∞
0 0
(335) < H ,φ >=< H, − φ >= − φ0 (x) dx = φ(0)
0
car φ est nulle en-dehors de l’intervalle. On a donc
d
(336) H(x) = δ(x)
dx
(2) Multiplication par des fonctions suffisamment régulières et ayant une
inverse :
Soit g : R → R une fonction différentiable et bijective, avec g 0 (x) > 0. Nous avons
en faisant le changement de variable y = g(x)
Z Z
dy
(337) f ◦ g,φ) = f (g(x))φ(x)dx = f (y)φ(g −1 (y)) 0 −1
g (g (y))
Exercice 30. Démontrer : δ(cx) = 1c δ(x).
Exercice 31. Que devient la fonction δ(x − x0 ) = δ(x − x0 )δ(y − y 0 ) en coordonnées
polaires? Rép. 1/r0 δ(r − r0 )δ(θ − θ0 ).
(3) Convergence des distributions :
Soit {fj }j≥0 une suite de distributions. Nous dirons que fj tend vers f si et seulement
si pour toutes fonctions test
(338) lim (fj ,φ) = (f,φ)
j→∞

Ce type de convergence est appelée convergence au sens faible.

R
Exercice 32. Soit ψ(x) une fonction suffisamment régulière telle ψ(x) dx = 1.
On définit la suite ψj (x) = jψ(jx).
R
(a) Démontrer ψj (x)dx = 1.
(b) Démontrer limj→∞ ψj (x) = δ(x) au sens faible.
(4) Transformée de Fourier des distributions On définit la transformée de Fourier
des distributions régulières à condition de prendre des précautions.
En effet, la défénition naturelle de la transformée de Fourier d’une distribution T
devrait être
(339) ∀φ, (FT,φ) = (T,F(φ))
 2. EQUATIONS AUX DéRIVéES PARTIELLES DU 1ER ORDRE: MéTHODE DES CARACTéRISTIQUES 69

Mais il y a un problème. Il existe en effet des fonctions à support borné, c’est-à-dire

nulle en dehors d’un intervalle, dont la TF n’est pas à support borné. D’où la nécéssité
de de réduire les contraintes imposées aux fonctions test.
On considère donc un nouvel ensemble de fonctions test : l’ensemble S des fonctions
indéfiniment dérivables à décroissance rapide à l’infini.
Définition 21. Un fonction est dite à décroissance rapide si et seulement si
(340) ∀i,k ∈ N2 lim xk φ(i) (x) = 0
|x|→∞
2
Ainsi e−x fait partie de cet ensemble, alors que 1/(1 + x2 ) n’en fait pas partie. On
appelle distribution tempérée toute application linéaire définie sur S et à valeur dans C.
1. Solutions élémentaires
Commençons par un exemple. On cherche à déterminer la solution de l’équation suivante:
(341) D∆u − bu = δ(x)
où D et b sont des constantes. On traite le problème en d = 2. Cette équation est une
équation de diffusion. L’espèce u diffuse dans le plan avec un coefficient de diffusion D. Elle
est irréversiblement détruite avec un taux b. Le delta de Dirac agit comme un terme source:
l’espèce u est créée avec un taux 1 en x = 0.
On résoud l’équation (341) en passant dans l’espace de Fourier.
(342) −Dk2 ũ(k) − bũ = 1
d’où
1
(343) ũ(k) = −
Dk2 +b
La transformée de Fourier inverse est alors:
Z +∞ Z 2π
1 eix·k
(344) u(x) = − k dk dθ
(2π)2 −∞ 0 Dk2 + b
La definition du produit scalaire donne x · k = |x||k| cos θ.
Définissons les fonctions
Z
1 π z cos θ
(345) In (z) = e cos(nθ) n ∈ N
π 0
(attention, n ∈ N). Ce sont des fonctions de Bessel modifiées.
2. Equations aux dérivées partielles du 1er ordre: Méthode des caractéristiques
 CHAPITRE 3

Le problème de Sturm-Liouville

Dans cette section, nous allons étudier un cadre général qui permet de montrer que certaines
fonctions sont orthogonales entre elles. Nous allons obtenir ces fonctions comme les solutions
d’une équation différentielle qui appartient à une classe très générale[?].
Définition 22. L’équation de Sturm-Liouville est l’équation différentielle ordinaire
 
d dy
(346) r(x) + (q(x) + λp(x))y = 0 où x ∈ [a,b]
dx dx
Ici p(x), r(x), q(x) sont des fonctions à valeurs réelles sur [a,b] et où y(x) et la constante λ ∈ C
sont toutes les deux à à déterminer. L’équation de Schrödinger est une équation de type
Sturm-Liouville, où la fonction q(x) joue le rôle de potentiel (avec r(x) = p(x) = 1).
Exemple 17. (1) L’équation de Bessel : x2 y 00 (x) + cy 0 (x) + (x2 − ν 2 )y(x) = 0
h i
d −x2 0 2
(2) L’équation d’Hermite : dx e y (x) + λe−x y = 0
hp i
d λ
(3) L’équation de Chebyshev : dx 1 − x2 )y 0 (x) + √1−x 2y = 0

Exemple 18. Conidérons l’équation différentielle

(347) y 00 + q(x)y = λy(x) x ∈ [a,b]
avec les conditions aux limites
(348) u(a) cos α + u0 (a) sin α = 0
(349) u(b) cos β + u0 (b) sin β = 0
Cette équation est du type (346) et elle peut être mise sous la forme
(350) Ly = λy
où L est un opérateur. Le paramètre λ suggère que le problème soit apparenté à une sorte de
probème à valeurs propres. Ce stype de problème est un problème spectral de type Sturm-Liouville
où λ joue le rôle des niveaux d’énergie dans le l’équation de Schrödinger.
La raison principale qui nous conduit à considérer le problème des valeurs propres est le
résulat suivant :
Théorème 51. Spit A une matrice n × n dans une espace vectoriel E de dimension fini n
muni d’un produit scalaire. Supposons
(351) ∀(x,y) ∈ E 2 (x,Ay) = (Ax,y)
Alors
(1) E possède une base orthonormale de vecteurs propres associés aux valeurs propres
λ1 , λ2 , . . . ,λn .
(2) A peut être décomposé comme
X
(352) A= λ i Pi
i

où les Pi sont des projections orthogonales.

71
72 3. LE PROBLèME DE STURM-LIOUVILLE

Procédant par analogies, nous chercherons à :

(1) Déterminer les vecteurs propres correspondants à l’opérateur L et décrire le développement
de y sur cette base.
(2) Détermier comment L peut être écrit dans cette base.
Bien que l’équation de Sturm-Liouville ait une forme spéciale, beaucoup d’équations d’ordre
2 se ramènent à cette équations. En particulier :
Théorème 52. Supposons que la fonction α1 (x) ne s’annule jamais. La définition
Z x 
α2 (t)
(353) r(x) = exp dt
0 α1 (t)
α3 (x)
(354) q(x) = r(x)
α1 x)
r(x)
(355) p(x) =
α1 (x)
montre que l’équation de Sturm-Liouville est équivalente à l’équation :
d2 y dy
(356) α1 (x) 2
+ α2 (x) + (α3 (x) + λ)y(x) = 0
dx dx
Définition 23. On appelle problème de Sturm-Liouville (SLP) le problème qui consiste à
déterminer la solution de l’équation :
 
d dy
(357) r(x) + (q(x) + λp(x))y = 0 où x ∈ [a,b]
dx dx
avec les conditions aux limites
(358) k1 y(a) + k2 y 0 (a) = 0
(359) l1 y(b) + l2 y 0 (b) = 0
où (k1 ,k2 ) 6= (0,0) et (l1 ,l2 ) 6= (0,0) sont donnés et où λ est un paramètre quelconque qui est
déterminé par la condition que le problème ait une solution autre que la solution nulle. C’est
donc la solution du problème qui détermine λ. On appelle λ la valeur propre et la solution
associée y(x,λ) la fonction propre associée.
Exercice 33. Trouver la solution du problème de Sturm-Liouville pour les problèmes
suivants:
(1) y 00 + λy = 0, y(0) = y 0 (π) = 0
(2) y 00 + λy = 0, y(0) − y 0 (π) = 0
Définition 24. Il y a 3 types de problèmes de SLP :
(1) Régulier si p(x) > 0 et r(x) > 0 sur [a,b].
(2) Singulier si p(x) > 0 et r(x) ≥ 0 sur ]a,b[ avec p(a) = p(b) = 0.
(3) Périodique y(a) = y(b) = 0.
Dans la suite, nous nous intéresserons principalement aux problèmes réguliers ou périodiques.
Remarque 7. L’une des propriétés remarquables est qu’une matrice symétrique possède des
valeurs propres réelles et que ses vecteurs propres forment une base orthonormée. Nous verrons
que le problème de Sturm-Liouville a des propriétés similaires et que l’on peut définir l’adjoint
d’un opérateur comme on définit l’adjoint d’une matrice et que ceux-ci possèdent une base propre
orthogonale.
Remarque 8. Le spectre des valeurs propres de l’équation de Schrödinger sur l’axe réel est
toujours continu. Il existe des valeurs de λ pour lesquelles les solutions sont à la fois non bornées
et ne décroissent pas à l’infini.
 3. LE PROBLèME DE STURM-LIOUVILLE 73

Anderson () a démontré qu’un potentiel aléatoire stationnaire pour les milieux désordonnés
peut avoir un spectre dense de fonctions propres associées à des fonctions ayant une décroissance
exponentielle. Ce phénomène est appelé localisation d’Anderson.
Comme exemple, considérons l’équation de Schrödinger sous la forme :
" #
X
(360) u00 (x) + Qn f (x − n) u(x) = λu(x)
n∈Z

où f (x) est un potentiel donné et où les Qn sont des variables aléatoires indépendantes distribuées
suivant la même loi de probabilité. Le spectre de cette équation est dense et dénombrable dans
λ ∈ [0,∞[ avec des solutions exponentiellement décroissantes.
Théorème 53. Pour un problème régulier :
(1) Les valeurs propres sont réelles.
(2) Les fonctions propres correspondantes à des valeurs propres distinctes sont orthogonales
avec le poids p(x) :
Z b
dx p(x)ui (x)uj (x) = 0 si i 6= j
a
(3) Les valeurs propres sont simples : l’espace associé à une valeur propre est de dimension
1.
(4) L’opérateur associé est auto-adjoint.
(5) Le spectre de l’opérateur −L défini plus lin est borné inférieurement.
(6) L’ensemble des vecteurs propre constitue une bas complète sur laquelle toute fonction
suffisament régulière et ayant les bonnes conditions aux limites peut être développée.
Point 1 :
Nous avons en notant ȳ le conjugué complexe de y :
 
d dy
(361) ȳ r + ȳ(q + λp)y = 0 + Conditions aux limites
dx dx
 
d dȳ
(362) y r + y(q + λ̄p)ȳ = 0 + Conditions aux limites
dx dx
d’où
0
(363) [r (y 0 ȳ − y ȳ 0 )] + (λ − λ̄)p|y|2 = 0
Soit en intégrant
Z b
0 0 b
(364) [r (y ȳ − y ȳ )]a = −(λ − λ̄) dx p(x)|y(x)|2
a
En raison des conditions aux limites, le membre de gauche est nul. La fonction p(x) est supposée
strictement positive et l’intégrale est nécessairement non nulle si la fonction y(x) n’est pas la
fonction nulle. D’où λ = λ̄. 
Point 2 :
Soient ui,j deux fonctions propres associées à deux valeurs propres λi 6= λj . En procédant comme
au point 1 :
Z b
 0 0

b
(365) r ui uj − ui uj a = −(λi − λj ) dx p(x)ui (x)uj (x) = 0
a
Comme auparavant, le membre de gauche est nul en raison des conditions aux limites. D’où le
résultat. 
Point 3 :
Ce point est le plus délicat. Rappelons le théorème d’algèbre :
74 3. LE PROBLèME DE STURM-LIOUVILLE

Théorème 54. Soit A une matrice n × n définie dans un espace vectoriel V avec un produit
scalaire (x,y). Si (Ax,y) = (x,Ay) alors
(1) V possède une base orthonormale (v1 ,v2 , . . . ,vn ) de vecteurs propres correspondants à
des valeurs propres distinctes (λ1 ,λ2 , . . . ,λn ).
P
(2) La matrice A peut être décomposée comme i=1,n λi Pi où les Pi sont les projecteurs
sur la base (v1 ,v2 , . . . ,vn ).
Pour définir l’équivalent d’une matrice symétrique, il faut définir ce qu’est un opérateur
auto-adjoint.
Considérons l’opérateur L défini par :
d2 u du
(366) Lu = a0 (x) 2
+ a1 (x) + a2 (x)u(x)
dx dx
sur l’espace des fonction définies sur [x0 ,x1 ] avec les conditions aux limites u(x0 ) = u(x1 ) = 0.
On peut prendre d’autres conditions aux limites mais celles-ci simplifieront l’écriture. Nous
désirons trouver la solution de
(367) Lu = f (x)
sous la forme
Z x1
(368) u(x) = G(x; z)f (z) dz = (Gx ,f )
x0

où Gx (z) = G(x; z) est une fonction de x quand z est donné. On a

(369) u(x) = (Gx ,f ) = (Gx ,Lu)
Définition 25. On définit l’adjoint de L, l’opérateur L? tel que
(370) ∀(u,v), (v,Lu) = (L? v,u) + BT
où BT désigne l’ensemble des termes qui ne dépendent que des extrémités de l’intervalle [x0 ,x1 ].
Exemple 19. Ici
Z x1
(v,Lu) = dx v [a0 u00 + a1 u0 + a2 u]
(371) x0
x
= (a0 v 00 + (2a00 − a1 )v 0 + (a2 + a000 − a01 )v,u) + [a0 (vu0 − v 0 u) + (a1 − a00 )uv]x10
Le dernier terme est un terme qui ne dépend que des points x0 et x1 . Il est nul en raison des
conditions aux limites. Le premier terme, quant à lui, donne la définition de l’adjoint de L
(372) L? v = a0 v 00 + (2a00 − a1 )v 0 + (a2 + a000 − a01 )v
Définition 26. Un opérateur est auto-adjoint si L = L? . En particulier, si a00 = a1 , il est
de type Sturm-Liouville (et nécessairement auto-adjoint).
Théorème 55. On a :
(373) L? Gx = δ(x − z)
au sens des distributions. En effet :
(374) u(x) = (Gx ,Lu) = (L? Gx ,u) + BT
Mais BT = 0 si Gx (x0 ) = Gx (x1 ) = 0. Donc
Z x1
(375) u(x) = Gx (z)u(z) dz
x0

D’où la propriété.
 1. EXISTENCE DES FONCTIONS PROPRES (VECTEURS PROPRES) 75

Remarque 9. Un opérateur de type Sturm-Liouville peut être mis sous la forme

(376) Lu = (p(x)u0 )0 + q(x)u

avec les conditions aux limites

(377) α1 u(a) + α2 u0 (a) = 0

(378) β1 u(b) + β2 u0 (b) = 0

Théorème 56. Les valeurs propres de l’opérateur −L possèdent une borne inférieure.

Nous démontrons le théorème lorsque α2 = β2 = 0 (dans le cas contraire, la démonstration

est moins évidente). Nous avons :
Z b
 
(−Lu,u) = −(pu0 )0 ū − r|u|2 dx
a
(379) Z b
 02 
= p|u | − r|u|2 dx + p(a)u0 (a)u(a) − p(b)u0 (b)u(b)
a

Mais u(a) = u(b) = 0 en raison des conditions aux frontières. Si u est un vecteur propre associé
à la valeur propre λ alors :
Z b Z b
2 0 2
(380) λkuk = p(x)|u | dx − r(x)|u(x)|2 ≥k u k2 max {|r(x)| : a ≤ x ≤ b}
a a

Ce qui achève la démonstration. 

1. Existence des fonctions propres (vecteurs propres)

1.1. Fonctions de Green. Considérons l’opérateur auto-adjoint :
d2 d
(381) L=p 2
+ p0 +r
dx dx
avec les conditions aux limites

(382) α1 u(a) + α2 u0 (a) = 0, |α1 | + |α2 | > 0

(383) β1 u(b) + β2 u0 (b) = 0, |β1 | + |β2 | > 0

Définition 27. Une fonction de Green pour cette opérateur est une fonction :
(1) De [a,b] × [a,b] → R.
(2) G(x,ξ) satisfait les conditions aux limites (382) pour les deux variables x et ξ. G(x,ξ) =
G(ξ,x).
(3) G(x,ξ) est continue sur [a,b] × [a,b] est est de classe C 2 sur [a,b] × [a,b] sauf sur la
droite x = ξ. En dehors de cette droite, G satisfait à Lx G(x,ξ) = 0.
(4) Enfin, la derivée (par rapport à x ) de G(x,ξ) subit un saut au voisinage de x = ξ, voir
Fig. 1
∂G ∂G 1
(384) (x = ξ + ,ξ) − (x = ξ − ,ξ) =
∂x ∂x p(ξ)

Exercice 34. Démontrer que la dernière propriété revient à exiger :

(385) Lx G(x,ξ) = δ(x − ξ)

76 3. LE PROBLèME DE STURM-LIOUVILLE

Fig. 1. Saut de la dérivée des fonctions de Green.

1.2. Construction de la fonction de Green. Nous considérons à nouveau les conditions

aux limites (382). L’équation différntielle étant uen équation du deuxième ordre, nous savons
que si nous pouvons trouver deux solutions indépendantes v1 et v2 telles que :
(386) v1 (a) = α2 v10 (a) = −α1
(387) v2 (b) = β2 v20 (b) = −β1
alors ce couple est unique. Nous supposerons que ces deux solutions sont linéairement indépendantes,
c’est-à-dire que le Wronskien W (v1 ,v2 )(x) 1

(389) W (v1 ,v2 ) = v1 (x)v20 (x) − v10 (x)v2 (x)

ne s’annule jamais. Définissons G(x,ξ) comme :
(
c−1 v1 (ξ)v2 (x) a ≤ ξ ≤ x ≤ b
(390) G(x,ξ) =
c−1 v1 (x)v2 (ξ) a ≤ x ≤ ξ ≤ b
où
(391) c = p(x) [v1 (x)v20 (x) − v10 (x)v2 (x)] = p(x)W (v1 ,v2 )
est une constante non-nulle. c est une constante, car
(392) [p(v1 v20 − v10 v2 )] = v1 Lv2 − v2 Lv1 = 0
G(x,ξ) défini par (390) possède toutes les propriétés des fonctions de Green. Nous avons en
particulier
∂G ∂G 1
(393) (x,x + ) − (x,x − ) = [v1 (x)v20 (x + ) − v10 (x − )v2 (x)]
∂ξ ∂ξ c
tend vers 1/p(x) quand  → 0.
1.3. Un problème équivalent.
Définition 28. Pour toutes fonctions f de C([a,b]), on construit l’opérateur T agissant sur
f:
Z b
(394) (T f )(x) = G(x,ξ)f (ξ) dξ
a
1. Soient f,g deux fonctions. Le déterminant
 
 f (x) g(x) 
(388) W (f,g)(x) =  0 
f (x) g 0 (x)
est appelé le Wronskien de f et g.
Si y1 (x) et y2 (x) sont deux solutions différentes d’une équation homogène sur un intervalle, alors leur
wronskien est strictement non-nul sur cet intervalle.
 2. FONCTIONS DE GREEN ET DéVELOPPEMENT SUR UNE BASE PROPRE 77

Théorème 57. L’opérateur T joue le rôle de “l’inverse” de L

(395) L(T f )(x) = u(x)
(396) T (Lu)(x) = u(x)
Théorème 58. Si u est une vecteur propre du problème de Sturm-Liouville avec la valeur
propre λ, c’est-à-dire
(397) Lu = λu
Alors u est valeur propre de T avec la valeur propre 1/λ et λ 6= 0.

2. Fonctions de Green et développement sur une base propre

Considérons à nouveau le problèmede Surm-Liouville
 
d dy
(398) r(x) + (q(x) + λp(x))y = f (x) où x ∈ [a,b]
dx dx
avec les conditions aux limites:
(399) k1 y(a) + k2 y 0 (a) = 0
(400) l1 y(b) + l2 y 0 (b) = 0
où (k1 ,k2 ) 6= (0,0) et (l1 ,l2 ) 6= (0,0) Ce problème possède une séquence ordonnée de valeurs
propres:
λ0 < λ1 < λ2 < λ3 . . .
avec des fonctions propres
(
Lφj = λj φj
(401)
φj satisfaisant aux conditions aux limites
qui peuvent orthonormales au sens du produit scalaire:
(402) (φi ,φj ) = δi,j
Comme cette base est complète, toute fonction peut être décomposée comme (faire l’analogie
avec des coefficients de Fourier)
∞
X
(403) g(x) = cj φj (x) avec cj = (φj ,g)
j=0
2
Cette série converge au sens de L avec
Z b X
(404) lim |g(x) − cj φj (x)|2 dx = 0
N →∞ a j=0,N

Théorème 59. Nous pouvons exprimer la fonction de Green G(x; z) = Gx (z) comme une
série:
∞
X
(405) G(x; z) = Gx (z) = cj (x)φj (z)
j=0

où
1
(406) cj (x) = φj (x)
λj
Remarque 10. (1) La formule
X 1
(407) G(x; z) = φj (x)φj (z)
j
λj
est clairement symétrique en x et z comme attendu.
78 3. LE PROBLèME DE STURM-LIOUVILLE

(2) Si l’une des valeurs propres λj est nulle, cette expression n’a pas de sens. Il n’y a pas
de fonction de Green au sens usuel dans ce cas.
(3) La solution de (398) est écrite comme
Z bX∞ X∞
1 1
(408) u(x) = (Gx ,f ) = φj (x)φj (z)f (z) dz = φj (x)(φj ,f )
a j=0 λj λ
j=0 j

et si λj = 0 alors nous exigeons la condition (φj ,f ) = 0.

La démonstration du théorème est simple. Nous avons :
X∞ ∞
X
(409) δx = LGx = cj (x)(Lφj )(z) = cj (x)λj φj (z)
j=0 j=0

d’où le résultat.

 CHAPITRE 4

Applications des fonctions de Green

1. Fonctions de Green pour le problème de Poisson

Il s’agit d’un problème elliptique. Trouver la solution du potentiel u(x) défini par une
distribution de charges f (x) dans une région D avec une frontière ∂D où u(x) est donné par
g(x).
(
∆u(x) = f (x) x ∈ D
(410)
u = g(x) x ∈ ∂D
Nous allons trouver la soltuion de ce problème en utilsant les focntions de Green.
Cherchons à trouver la solution élémentaire du problème suivant:
(
∆Gx (y) = δx (y) y∈D
(411)
Gx (y) = 0 y ∈ ∂D
‘où δx (y) est la fonction delta bi ou tridimensionnelle
δ(y − x) = δ(y1 − x1 )δ(y2 − x2 )δ(y3 − x3 )
Pour la suite, il sera utile d’employer une notation. Nous définissons la dérivée normale de u
à la frontière ∂D de D à partir de la normale n̂ orientée de l’intérieur vers l’extérieur.
Baptisons :
∂u
(412) = n̂ · ∇u
∂n
Nous aurons besoin du résultat suivant:
Lemme 2. Lemme de Green: Si u1 et u2 sont deux fois continûment différentiables
Z Z  
∂u2 ∂u1
(413) (u1 ∆u2 − u2 ∆u1 ) dV = u1 − u2 dS
D ∂D ∂n ∂n
Preuve : Nous pouvons soustraire les deux égalités :
(414) ∇ · (u1 ∇u2 ) = ∇u1 · ∇u2 + u1 ∆u2
(415) ∇ · (u2 ∇u1 ) = ∇u2 · ∇u1 + u2 ∆u1
et utiliser le théorème de la divergence pour arriver à
Z Z Z  
∂u2 ∂u1
(416) (u1 ∆u2 − u2 ∆u1 ) dV = n̂ · (u1 · ∇u2 − u2 ∆u1 ) = u1 − u2 dS
D ∂D ∂D ∂n ∂n

Prenons maintenant u1 = u comme la solution du problème de Poisson (410) et u2 al fonction
de Green. Nous avons :
Z Z
∂Gx
(417) (u(y)δx (y) − Gx (y)f (y)) dy = g(y) dS(y)
D ∂D ∂x
Soit enfin
Z Z
∂G(x; y)
(418) u(x) = G(x; y)f (y) dy + g(y) dS(y)
D ∂D ∂n(y)

79
80 4. APPLICATIONS DES FONCTIONS DE GREEN

2. Fonctions de Green sur R2 et R3

En règle générale, il n’y a aucun espoir de trouver une formule explicite des fonctions de
Green pour un domaine arbitraire. Lorsque D = R2 ou R3 , les choses sont simples.
2.1. Si D = R2 . Si y 6= x
 
∂2 1 ∂ 1
(419) ∆G = 0 = 2
+ h(r) = (rh0 (r))0
∂r r ∂r r
D’où h0 = c1 /r soit h = c1 ln r + c2 . La constante c2 est une simple constanteadditive, donc ous
pouvons pose c2 = 0. La constante c1 , quant à elle, est évaluée en appliquant le théorème de la
divergence
Z Z
1= δx (y) dy = ∇ · ∇Gx (y) =
|y−x|≤ |y−x|
(420) Z Z
y−x y−x
= ∇Gx (y) · n̂dS(y) = c1 2
· dS(y) = 2πc1
|y−x|= |y−x|= |y − x| |y − x|
D’où c1 = 1/2π et
1
(421) GR2 (x; y) = ln|y − x|
2π
2.2. Si D = R3 . Nous avons
1
(422) GR3 (x; y) = −
4π|y − x|
Exercice 35. Utiliser
 
∂2 2 ∂ 1
(423) ∆G = 0 = 2
+ h(r) = 2 (r2 h0 (r))0
∂r r ∂r r
et démontrer cette égalité.
3. Fonctions de Green pour des problèmes dépendants du temps
Nous nous intéressons au problème de la chaleur qui décrit l’évolution de la température
u(x,t). Le problème est de trouver la solution de l’équation au dérivées partielles (dérivée
première par rapport au temps, dérivée seconde par rapport aux variables d’espace) :
∂u
(424) − D∆u = f (x,t) x ∈ V, t > 0
∂t
avec les conditions aux limites et les conditions initiales (2 conditions)
(425) u(x,t) = g(x,t), t > 0 x ∈ ∂V
(426) u(x,0) = u0 (x) x ∈ V
Nous désirons représenter la solution à l’aide des fonctions de Green. L’opérateur L n’est pas
auto-adjoint dans ce cas.
 CHAPITRE 5

Polynômes orthogonaux et fonctions de Bessel

Nous considérons d’abord 3 exemples de problème de type Sturm-Liouville irréguliers

où les vecteurs propres sont des polynômes. Chaque exemple revient à choisir des coefficients
particuliers pour le problème
(427) (pu0 )0 + ru + λρu = 0
où toutes les paires de fonctions propres doivent satisfaire à l’équation :
(428) p(u0 v − uv 0 )|ba = 0
Supposant qu’il existe une base propre orthogonale {φn }n∈N , nous pouvons développer n’importe
qu’elle fonction f ∈ L2 [a,b] sur cette base
X∞
(f,φn )
(429) f (x) = φn (x)
n=0
kφn k2
qui est l’analogue d’un développment en série de Fourier.

1. Généralité sur les polynômes orthogonaux

Soit ρ(x) une fonction inégrable de signe constant (ρ(x) peut être discontinue). Rappelons
que nous pouvons définir un produit scalaire
Z b
(430) (P,Q) = P (x)Q(x)ρ(x) dx
a
La famille de polynôme {Pn }n≥0 de degré n constitue une famille de polynômes orthogonaux
si et seulement si
(431) (Pn ,Pm ) = δn,m
Théorème 60. Il est possible de construire une suite de polynôme orthogonaux Pn (x) sur
un intervale [a,b] relatovement à une densité arbitraire ρ(x), définie sur cet intervalle.
La preuve de ce théorème s’appuie sur l’algorithme de Gram-Schmidt qui donne constructivement
l’existence d’une base orthonormée. Soit {vn }n∈N un ensemble de vecteurs indépendants. La
base othogonale est obtenue en retranchant les projections successives :
(432) u1 = v1
(433) u2 = v2 − Proju1 (v2 )
(434) u3 = v3 − Proju1 (v3 ) − Proju2 (v3 )
(435) ...
(436)
où le projecteur sur le vecteur u s’obtient en prenant le produit scalaire avec u.
(u,v)
(437) Proju (v) = u
(u,u)
Théorème 61. Nous avons
Z b
(438) ρ(x)xk Pn (x) dx = 0
a

81
82 5. POLYNôMES ORTHOGONAUX ET FONCTIONS DE BESSEL

pour 0 ≤ k < n.
Théorème 62. Si {Pn (x)}n∈N est une famille de polynômes orthogonaux, il existe trois
nombres An , Bn , Cn tels que :
(439) xPn (x) = An Pn+1 (x) + Bn Pn (x) + Cn Pn−1 (x)
Le polynôme xPn (x) est un polynôme de degré n + [Link] peut donc le développer sur la base
des polynômes Pn (x) comme suit :
k=n+1
X
(440) xPn = ck Pk (x)
k=0

Or (Pk ,xPn ) = (Pn ,xPk ) = 0 si

(441) k > n + 1 et n > k + 1
Ce qui veut dire que seuls les termes k < n − 1 et k > n + 1 sont non nuls dans (440). 
Théorème 63. Les zéros des polynômes Pn (x) sont réels, simples et contenus dans l’intervalle
ouvert ]a,b[
Comme pour n > 0,
Z b
(442) ρ(x)Pn (x) dx = 0
a
Pn (x) doit changer de signe au moins une fois sur ]a,b[. Soient xi , 1 ≤ i . . . ≤ p la suite des
points où Pn (x) s’annule. Nous avons 1 ≤ p ≤ n. Formons le polynôme
Y
(443) Q(x) = (x − xi )
i=1,...,p

Le produit P (x)Q(x) garde un signe constant sur l’intervalle et donc

Z b
(444) ρ(x)P (x)Q(x) dx 6= 0
a
si p > n, ce n’est pas possible, car Q(x) est un polynôme de degré strictement inférieur à n et
qui est donc nécessaierement orthogonal à P (x) (par définition de l’algorithme de Schmidt).
Donc p = n. 
Définition 29. Soit une fonction G(x,t) développable en série entière de t dans un certain
intervalle
X
(445) G(x,t) = φn (x)tn
n≥0

G(x,t) est appelée la fonction génératrice des φn (x).

Rappelons enfin la définition :
Définition 30. Un problem de type SL est dit régulier si et seulement si :
(1) [a,b] est un intervalle fini.
(2) p, p0 , q sont continues sur [a,b].
(3) p(x) garde un signe constant.
En règle pratique, les problèmes SL deviennent des problèmes singuliers dans les cas suivants :
(1) p(x) > 0 pour x ∈]a,b[ mais p(a) = 0 et/ou p(b) = 0.
(2) L’intervalle est infini.
Les points a ou b sont alors singuliers si a = −∞ ou p(a) = 0 (respectivement pour b).Lorsque
p(a) = p(b) = 0, toutes les conditions aux limites disparaissent. Les conclusions du chapitre
précédent restent valables (bien que nous ne l’ayons pas démontré). Lorsque l’intervalle est infini,
il faut invoquer des conditions suplémentaires.
 2. POLYNôMES DE LEGENDRE 83

2. Polynômes de Legendre
Le polynômes de Legendre ont été introduits dans le cadre du potentiel newtonien. Un
corps situé au point A(0,0,1) crée au point P de coordonnées (r, θ) un potentiel inversement
proportionnel à la distance kAP k:
1 1
(446) =
kAP k (1 − 2r cos θ + r2 )1/2
qu’il est intéressant de développer en puisance de r. Nous aurons :
X ∞
1 1
(447) = 1/2
= Pl (cos θ)rl
kAP k 2
(1 − 2r cos θ + r ) l=0

Exercice 36. Les polynômes de Legendre dont solutions de l’équation différentielle

(448) (1 − x2 )u00 − 2xu0 + λu = 0

qui est donc un SL singulier, x ∈ [−1, + 1].

(1) En supposant
∞
X
(449) u(x) = ck x k
k=0

Démontrer que
k(k + 1) − λ
(450) ck+2 = ck
(k + 1)(k + 2)
En déduire que si λ = n(n + 1), n ∈ N, alors l’une des solutions est un polynôme (on
choisira les constantes c0 et c1 ) .
(2) L’autre solution est la fonction de Legendre Qn (x). Démontrer que
 
1 1+x
(451) Q0 (x) = ln
2 1−x

(3) On considère la fonction génératrice

(1 − 2xt + t2 )−1/2

Cette fonction génératrice est une fonction analytique de la variable t ∈ [−1,1]. Nous
pouvons poser
X
(452) f (t,x) = an (x)tn
n≥0

(4) Démontrer que cette fonction satifait à l’équation différentielle

∂f
(453) (1 − 2xt + t2 ) = (x − t)f
∂t
En déduire l’équation de récurrence

(454) (n + 1)an+1 (x) + nan−1 (x) = (2n + 1)xan (x)

qui est une relation de récurrence pour des polynômes orthogonaux.

(5) Faire x = 1 et x = −1 dans (36). En déduire Pn (1) = 1 et Pn (−1) = (−1)n .
84 5. POLYNôMES ORTHOGONAUX ET FONCTIONS DE BESSEL

3. Polynômes d’Hermite et de Laguerre

3.1. Polynômes d’Hermite. Nous considérons la fonction génératrice
2
X 1
(455) e2xt−t = Hn (x)tn
n≥0
n!

Exercice 37. (1) Démontrer que

dn −x2
2
(456) Hn (x) = (−1)n ex e
dxn
qui est manifestement un polynôme de degré n.
(2) En intégrant par parties, démontrer que (noter la définition du produit scalaire):
Z +∞
2
(457) (Hn ,Hm ) = e−x Hn (x)Hm (x) dx = 0
−∞

Les polynômes Hn (x) sont solutions de l’équation :

(458) u00 − 2xu0 + 2nu = 0
3.2. Polynômes de Laguerre.
4. Applications physiques
4.1. Équation de Laplace en coordonnées sphériques. On cherche à résoudre l’équation
de Laplace
∂ 2u ∂ 2u ∂ 2u
(459) + + =0
∂x2 ∂y 2 ∂z 2
et on se place en coordonnées polaires
   
∂ 2 ∂u 1 ∂ 2u 1 ∂ ∂u
(460) r + + sin θ =0
∂r ∂r sin2 θ ∂ϕ2 sin θ ∂θ ∂θ
Supposant le système invariant par rotation autour de l’axe z, la solution ne peut dépendre de
ϕ. Autrement dit,
   
∂ 2 ∂u 1 ∂ ∂u
(461) r + sin θ =0
∂r ∂r sin θ ∂θ ∂θ
Comme d’habitude, nous utilisons la mtéhode de séparation des variables et nous cherchons une
soltuion sous la forme
(462) u(r,θ) = v(r)w(θ)
Après divisoion par vw, nous avons :
   
1 d 2 dv 1 d dw
(463) r =− sin θ
v(r) dr dr w sin θ dθ dθ
Le membre de gauche ne dépend que de r alors que le membre de droite ne dépend que de θ. Il
existe donc une constante λ telle que
(464) r2 v 00 /2rv 0 − λv = 0
1
(465) (w0 sin θ)0 + λw = 0
sin θ
Posant ξ = cos θ, la deuxième équation peut être mise sous la forme
 
d

2 dw
(466) 1−ξ + λw = 0
dξ dξ
qui est une équation de Legendre. Lorsque λ = n(n + 1), la solution est donnée par les polynômes
de Legendre Pn (cos θ).
 5. FONCTIONS DE BESSEL 85

Exercice 38. Montrer que la solution pour r est alors

(467) vn (r) = an rn + bn r−(n+1)
où an et bn sont deux constantes.
4.2. Oscillateur harmonique en mécanique quantique. En mécanique quantique,
l’état d’une particule est représentée par une fonction d’onde
(468) Ψ(x,t) = ψ(x)e−iEt/~
où la fonction ψ(x) satisfait à l’équation
2m
(469) ψ 00 + 2 [E − V (x)] ψ(x) = 0
~
qui est de la forme Sturm-Liouville. Ce problème est singulier, car l’intervalle est nécessairement
infini, x ∈] − ∞, + ∞[. Suppsons le potentiel harmonique V (x) = ~2 x2 /(2m). Nous sommes
alors ramenés à étudier :
(470) ψ 00 + (λ − x2 )ψ = 0
2 /2
Le changement de variable u(x) = ex ψ(x) ramène cette équation à:
(471) u00 − 2xu0 + (λ − 1)u = 0.
qui est une équation de type Hermite si λn = 2n + 1, n ∈ N. Les fonctions propres sont donc
2 /2
(472) ψn (x) = e−x Hn (x)

5. Fonctions de Bessel
5.1. Introduction. Les fonctions de Bessel se rencontrent aussi fréquemment que les
fonctions sinus ou cosinus. Elles interviennent chaque fois que le problème possède une symétrie
cylindrique. Citons pêle-mêle :
(1) La physique des membranes vibrantes et donc de tous les instruments acoustiques.
(2) La physique de la diffraction par un trou circulaire en optique (tache d’Airy, cf. Fig 1).
(3) La propagation des ondes dans les guides d’onde et dans les fibres optiques.
(4) Tout phénomène où le Laplacien bi-dimensionnel intervient (équation de Schrödinger,
physique de la diffusion etc.)
L’équation différentielle de Bessel est du type :
 
00 1 0 ν2
(473) y + y + 1− 2 y =0
x x
et on désigne habituellement la solution de (473) par Zν (x).
Exemple 20. On considère l’équation de Helmholtz en symétrie cylindrique:
 
2 1 d 1 dψ
(474) ∆ψ + k ψ = + k2ψ = 0
r dr r dr
Cette équation se ramène à la précédente, car
1
(475) ψ 00 + ψ 0 + k 2 ψ = 0
r
et on cherchera la solution sous la forme ψ(r) = y(kr) où y est solution de (473).
De nombreuses équations se ramènent à cette équation par changement de variable et de
fonction. Par exemple, l’équation différentielle
 
00 1 − 2α 0 β−1 2 α2 − ν 2 β 2
(476) y + y + (Kβx ) + y=0
x x2
86 5. POLYNôMES ORTHOGONAUX ET FONCTIONS DE BESSEL

Fig. 1. Tache d’Airy par un trou circulaire

a pour solution
(477) y(x) = xα Zν (Kxβ )
En particulier :
(478) x(y + y 00 ) + (2ν + 1)y 0 = 0
qui a pour solution
(479) y(x) = x−ν Zν (x)

5.2. Forme intégrale des solutions. Chercchons une solution de (477) en la représentant
par une intégrale de Laplace très générale
Z
(480) v(x) = exz f (z) dz
l

où le chemin L et la fonction f sont tous les deux à déterminer. On a :

Z
0
(481) v (x) = exz zf (z) dz
l

D’où
Z
00
(482) x(v + v ) = exz x(1 + z 2 )f (z) dz
L

d’où en intégrant par parties :

Z
00 xz 2 d  
(483) x(v + v ) = |e (1 + z )f (z)|L − exz (1 + z 2 )f (z) dz
L dz
Posons pour alléger l’écriture
(484) g(z) = (1 + z 2 )f (z)
Pour que v soit solution de (477), il suffit que :
Z  
xz 2 xz d z
(485) |e (1 + z )f (z)|L + e g(z) − (2ν + 1) 2 g(z) dz = 0
L dz z +1
qui est évidemment vérifié si les deux termes sont égaux à zéro
(486) |exz (1 + z 2 )f (z)|L = 0
d z
(487) g(z) − (2ν + 1) 2 g(z) = 0
dz z +1
 5. FONCTIONS DE BESSEL 87

1.0

0.8

0.6

0.4

0.2

5 10 15 20

-0.2

-0.4

Fig. 2. Graphe des fonctions J0 (x), J1/2 (x) et J1 (x). La seule fonction qui ne
s’annule pas à l’origine est J0 (x). Lorsque x est grand, on retrouve bien que


2 1/2


Jν (x) ≈ πx cos x − π2 (ν + 1/2) .

La première condition ne dépend que des points limite du parcours L. Celui-ci sera choisi une
fois la fonction f déterminée par la seconde condition qui s’écrit aussi
dg 1 2z
(488) = (ν + ) 2
dz 2 z +1
d’où :
(489) g(z) = Cte · (z 2 + 1)ν+1/2 dz
soit
Z
ν
(490) Zν (x) = Cte · x exz (1 + z 2 )ν−1/2 dz
L

Le parcours L sera déterminé dan la section suivante. L’équation différentielle étant du deuxième
ordre, il existe aussi une autre famille de solutions. Celle-ci sera explicitée plus loin.

5.3. La fonction Jν (x). La fonction (z 2 + 1)ν+1/2 s’annule pour z = ±i à condition que

ν > −1/2. Les deux points limites de L sont donc ±i. Il est donc simple de choisir pour L la
droite qui interpolle entre ces deux points. Faisons les changement de variable z = it et en
multipliant par un constante appropriée, nous définissons les fonction de Bessel Jν (x) comme :
1  x ν Z +1
(491) Jν (x) = 1/2 dt eixt (1 − t2 )ν−1/2 dt
π (ν − 1/2)! 2 −1

à condition que ν > −1/2.

Exercice 39. Utilisant :
 
1 1 √
(492) ! = (− )! = π
2 2
Démontrer
Z 2π Z π
1 ix cos θ 1
(493) J0 (x) = e dθ = cos(x sin(θ)) dθ
2π 0 π 0
88 5. POLYNôMES ORTHOGONAUX ET FONCTIONS DE BESSEL

Pour s’affranchir de cette condition, constatons que x−ν Jν (x) est une fonction paire de x,
analytique et donc développable en série entière. Écrivons pour l’intégrale:
Z +1 Z +1
ixt 2 ν−1/2
dt e (1 − t ) dt = dt cos(xt)(1 − t2 )ν−1/2 dt
−1 0
(494) ∞
X Z +1
x2kk
=2 (−1) dt tk (1 − t2 )ν−1/2
k=0
(2k)! 0

qui devient après le changement de variable t = sin ϕ

X∞ Z +π
x2k
(495) 2 dϕ (sin ϕ)2k (cos ϕ)2ν
k=0
(2k)! 0

Cette intégrale est tabulée. Tout calcul fait, on trouve:

X∞
1  x ν+2k
k
(496) Jν (x) = (−1)
k=0
k!(ν + k)! 2
qui est définie quelque soit ν, car la série converge quelque soit ν. Lorsque ν n’est pas entier, on
constate aussi que
1  x ν
(497) Jν (x) ≈x→0
ν! 2
Comme l’une des deux fonction Jν (x) ou J−ν (x) diverge lorsque x → 0 lorsque n n’est pas
entier, les deux solutions sont linéairement indépendantes et nous obtenons la solution générale
comme une combinaison linéaire de ces deux fonctions (voir épilogue quand ν est un entier).
Théorème 64. Lorsque n est entier :
(498) Jn (x) = (−1)n Jn (x)

6. Fonction génératrice des fonctions de Bessel

Comme les polynômes orthogonaux, les fonction de Bessel ont une fonction génératrice.
Théorème 65.
n=+∞
X
(499) g(x,t) = e(x/2)(t−1/t) = Jn (x)tn
n=−∞

En effet, nous pouvons faire le produit

" #" #
X  x r tr X  x s t−s
(500) ext/2 e−x/(2t) = (−1)s
r≥0
2 r! s≥0
2 s!
Posons n = r − s et sommons sur les n et les s au lieu de sommer sur les r et les s: n varie donc
de −∞ à +∞. Nous avons
n=+∞
X X  x 2s  x n tn
n=+∞
X
(501) ext/2 e−x/(2t) = (−1)s = tn Jn (x)
n=−∞ s≥0
2 2 (n + s)!s! n=−∞


Comme d’habitude, la fonction génératrice permet d’obtenir les relations de récurrence pour
les Jn (x). On a
 
∂g 1 1
= x 1 + 2 e(x/2)(t−1/t)
∂t 2 t
(502) n=+∞
X
= nJn (x)tn−1
n=−∞
 7. DIFFRACTION DE FRAUNHOFER 89

0.25

0.20

0.15

0.10

0.05

-10 -5 5 10

Fig. 3. Fonction (J1 (x)/x)2 utilisé pour calculer l’intensité difractée par une
ouverture circulaire. La fonction est maximale en 0 et est rapidement négligeable.

En substituant la série pour e(x/2)(t−1/t) et en identifiant les puissances de t (on rappelle que la
représentation en série est unique), nous démontrons
2n
(503) Jn−1 (x) + Jn+1 (x) = Jn (x)
x
Exercice 40. Dériver la fonction génératrice g(x,t) par rapport à x. Démontrer :
(504) Jn−1 (x) − Jn+1 (x) = 2Jn0 (x)
avec en particulier:
(505) J00 (x) = −J1 (x)

7. Diffraction de Fraunhofer
Dans la théorie de la diffraction par une ouverture circulaire, nous avons à évaluer l’intégrale
Z a Z 2π
(506) Φ∝ rdr eibr cos θ dθ
0 0
où Φ est l’amplitude de l’onde diffractée. L’angle θ est l’ange azimuthal dasnnle plan de
l’ouverture circulaire de rayon a et l’angle α est l’angle défini par le point sur l’écran avec la
normale à l’ouverture. Le pramètre b est défini par:
2π
(507) b= sin α
λ
où λ est la longueur d’onde incidente.
Nous avons donc:
Z a
(508) Φ ∝ 2π rJ0 (br)
0
soit en utilisant (refeq:bessel0),
 
2πab λa 2πa
(509) Φ ' 2 J1 (ab) = J1 sin α
b sin α λ
D’où l’intensité :
  2
2 λa 2πa
(510) Φ ' J1 sin α
sin α λ
qui est tracée en Fig. 3.
90 5. POLYNôMES ORTHOGONAUX ET FONCTIONS DE BESSEL

2.5

2.0

1.5

1.0

0.5

0.0 0.5 1.0 1.5 2.0 2.5 3.0

Fig. 4. . Graphe des fonctions de Bessel Iν (x), ν = 0,1 et Kν (x). Les fonctions
Kν (x) décroissent exponentiellement et divergent en x = 0 y compris pour ν = 0.

8. Épilogue 1
Nous avons démontré que les fonctions de Bessel Jν (x) donnaient une famille de fonctions
solutions de (473). Lorsque ν n’est pas une entier, les fonctions Jν (x) et J−ν (x) sont indépendantes.
Que se passe-t-il lorsque ν est entier?
Théorème 66. Les fonctions de Bessel modifiées définies par
(
(Jν (x) cos(νx) − J−ν (x)) / sin(νπ) ν 6= 0, 1, 2 . . .
(511) Yν (x) =
limν→n Yν n = 0, 1, 2 . . .
donnent la deuxième famile de solutions de cette équation différentielle. On remarque bien
que ν = n entier donne la forme indéterminée 0/0 dans la première équation.

9. Fonctions de Bessel modifiées

L’équation de Helmholz
(512) ∆ψ + k 2 ψ = 0
a pour solution une combinaison linéaire des fonctions de Bessel Jν et des fonctions de Neuman
Yν avvec les conditions aux limites appropriées.
Nous désirons maintenant connaître les solutions d’une équation de type diffusion
(513) ∆ψ − k 2 ψ = 0
où nous supposons, là encore, une symétrie de révolution. L’analogue de l’équation (473) est ici :
 
00 1 0 ν2
(514) y + y − 1+ 2 y =0
x x
Pour passer de (512) à (513), il suffit de prendre k → ik (i2 = 1). On pose donc :
(515) Iν (x) = i−ν Jν (ix)
où iν est un facteur de normalisation. On pose généralement
(516) Iν (x) = e−νiπ/2 Jν (eiπ/2 x)
 10. ÉPILOGUE 2 91

L’autre famille de solutions est obtenue comme les solutions de Neuman. On pose :
π I−ν (x) − Iν (x)
(517) Kν (x) =
2 sin νπ
La famille des Kν (x) (fonctions de Bessel modifiées) permet d’obtenir les fonctions de Green du
problème de la diffusion. En particulier, la fonction de Green du problème
(518) ∆ψ − ψ = 0
est −1/2πK0 (x).
Exercice 41. En utilisant l’identité
Z ∞
(519) rdr K0 (r) = 1
0
démontrer que la fonction de Green est −1/2πK0 (r).
Enfin, les fonctions Kν (x) décroissent exponentiellement pour x → ∞
r
π −x
(520) Kν (x) 'x→∞ e (1 + . . .)
2x
Cette décroissance exponetielle est souvent utlisée pour définir des longueurs caractéristiques en
physique statistique.
10. Épilogue 2
 ANNEXE A

Rappels de topologie sur un espace métrique

1. Ouverts, fermés, compacts

1.1. Espaces métriques. Soit E un ensemble. On dit que l’application d : E × E → R
est une distance si elle vérifie les trois propriétés suivantes :
(1) Séparation : d(x,y) = 0 ⇒ x = y.
(2) Symétrie : d(x,y) = d(y,x).
(3) Inégalité triangulaire : ∀x,y,z ∈ E, d(x,z) ≤ d(x,y) + d(y,z).
et on baptise le couple (E,d) espace métrique. R et C sont des espaces métriques munis de la
distance habituelle d(x,y) = |x − y|.
1.2. Ouverts, fermés.
Définition 31. Pour tout x0 ∈ E et quelque soit r > 0, on appelle boule ouverte de centre
x0 et de rayon r l’ensemble:
(521) B(x0 ; r) = {x ∈ E, d(x,x0 ) < r}
La boule fermée est par définition:
(522) B̄(x0 ; r) = {x ∈ E, d(x,x0 ) ≤ r}
Définition 32. (1) Une partie U de E est un ouvert de E si pour chaque x de U on
peut trouver  > 0 tel que B(x,) ⊂ U .
(2) Une partie U de E est un fermé si son complémentaire dans E est un ouvert de E.
Théorème 67. Une partie F de E est un fermé si pour toute suite (xn )n d’éléments de F
convergeant vers un élément x de E, alors x ∈ F .
Remarque 11. Pour les espaces “ressemblants” à Rn , toutes les boules ouvertes sont des
ouverts et toutes les boules fermées sont des fermés.
Théorème 68. Soit I un sous-ensemble de N. Soient (Ui )i∈I une famille éventuellement
infinie d’ouverts et (Fi )i∈I une famille de fermés.
(1) ∪i∈I Ui est un ouvert.
(2) Si I est fini, ∩i∈I Ui est aussi un ouvert (penser à ] − n1 , n1 [).
(3) ∩i∈I Fi est un fermé.
(4) Si I est fini, ∪i∈I Fi est un fermé (penser à ∪n>0 [0,1 − n1 ]) = [0,1[ qui n’est pas un
fermé).
1.3. Adhérence d’un ensemble.
Définition 33. Soit P une partie de E. On dit que x est adhérent à P si
(523) ∀ > 0, B(x,) ∩ P 6= ∅
On appelle adhérence de P l’ensemble noté P̄ des points adhérents à P .
Théorème 69. (1) Tout élément de P̄ est nécessairement la limite d’une suite de
d’éléments de P .
(2) Une partie est fermée si elle est sa propre adhérence.
(3) L’adhérence de P est l’intersection de tous les fermés contenant P .

93
94 A. RAPPELS DE TOPOLOGIE SUR UN ESPACE MéTRIQUE

1.4. Compacts.
Définition 34. Si K est une partie de R ou C. K est un compact si et seulement si K
est un fermé borné. Dans un compact K, toute suite d’éléménts de K possède une sous-suite
convergente.