Méthodes de calcul d'intégrales de contour

En analyse complexe, l'intégration de contour est une technique de calcul d'intégrale le long de chemins sur le plan complexe[1],[2],[3]

L'intégration de contour est fortement liée au calculs de résidus[4], une méthode de calcul utilisée pour évaluer des intégrales curvilignes sur l'axe des réelles, que les outils de la théorie de l'intégration ne permettent pas de calculer par une simple analyse réelle[5]

Les méthodes d'intégration de contour incluent :

Ces méthodes peuvent être combinées pour obtenir les résultats attendus.

Courbes dans le plan complexe[modifier | modifier le code]

Les contours donnent une définition précise des courbes sur lesquelles une intégrale peut être exactement définie. Une courbe du plan complexe est donc une fonction continue renvoyant un intervalle réel fermé sur le plan complexe : z : [a, b] → C.

Cette définition coïncide avec la notion intuitive d'une courbe, mais inclut une paramétrisation par une fonction continue à partir d'un intervalle fermé. Cette définition plus précise permet de réfléchir aux propriétés que doit avoir une courbe pour qu'elle soit utile à l'intégration. Dans les sous-sections suivantes, on réduira l'ensemble de courbes intégrables pour n'inclure que celles qui peuvent être construites à partir d'un nombre fini de courbes continues, auxquelles on peut donner une direction. De plus, on empêchera les "morceaux" de se croiser, et on imposera que chaque morceau ait une dérivée continue finie (non nulle). Ces exigences permettent de ne considérer que les courbes qui peuvent être tracées, comme par un stylo, dans une séquence de traits constants et réguliers, qui ne s'arrêtent que pour commencer un nouveau morceau de la courbe, le tout sans relever le stylo[6].

Courbes lisses orientées[modifier | modifier le code]

Les contours sont définis comme des courbes régulières orientées[6]. On peut ainsi définir précisément un "morceau" d'une courbe lisse, dont on trace un contour.

Une courbe lisse est une fonction z : [a, b] → C avec une dérivée continue qui ne s'annule pas et ne possèdent pas de point double, sauf peut-être aux extrémités (tel que z(a) = z(b)). On parle alors de courbe fermée, et la fonction doit alors être bijective et les dérivées doivent correspondre aux extrémités (z′(a) = z′(b)). Sinon, on parle d'arc lisse ou de courbe ouverte[6].

La paramétrisation d'une courbe donne un ordre naturel des points de la courbe : z(x) est avant z(y) si x < y. On peut ainsi définir une courbe lisse orientée, ce qui est surtout utile pour les courbes indépendantes d'une paramétrisation spécifique. Ceci peut être fait en considérant des classes d'équivalence de courbes lisses avec la même direction. Une courbe lisse orientée est alors un ensemble ordonné de points du plan complexe, image d'une courbe lisse dans leur ordre naturel (selon la paramétrisation). Il faut noter que tous les ordres de points correspondent à l'ordre naturel d'une courbe lisse. En effet, une courbe lisse donnée n'a que deux ordres. De plus, une courbe fermée peut avoir tout point comme extrémité, alors qu'un arc lisse aura deux possibilités.

Contours[modifier | modifier le code]

Les contours sont les classes de courbes sur lesquelles on peut définir l'intégration de contour. Un contour est une courbe orientée qui est formé d'une suite finie de courbes lisses dirigées dont les extrémités sont reliées pour donner une seule direction. Cela nécessite que la suite de courbes γ1, …, γn soit telle que, le point final de γi coïncide avec le point initial de γi+1, i, 1 ≤ i < n. Cela inclut toutes les courbes lisses dirigées. De même, un point singleton du plan complexe peut être considéré comme un contour. Le symbole + est souvent utilisé pour désigner le raccordement des courbes afin de former une nouvelle courbe. Ainsi on peut construire un contour Γ qui est formé de n courbes telles que

Intégrales de contour[modifier | modifier le code]

L'intégrale de contour d'une fonction complexe f : CC, est une généralisation de l'intégrale pour les fonctions réelles. Pour les fonctions continues du plan complexe, l'intégrale de contour peut être définie de façon analogue à l'intégrale curviligne en définissant d'abord l'intégrale le long d'une courbe lisse orientée en termes d'intégrales sur un paramètre réel. Une définition plus générale peut être donnée en termes de partitions du contour par analogie avec la partition d'un intervalle et l'intégrale de Riemann. Dans les deux cas, l'intégrale sur un contour est définie comme la somme des intégrales sur les courbes lisses orientées qui constituent le contour.

Pour les fonctions continues[modifier | modifier le code]

Pour définir l'intégrale de contour dans ce cas, on dit d'abord considérer l'intégrale, sur une variable réelle, d'une fonction à valeurs complexes. Soit f : RC une fonction à valeurs complexes d'une variable réelle t. Les parties réelles et imaginaires de f sont souvent dénotées u(t) et v(t), respectivement, de sorte que

Alors l'intégrale de f sur l'intervalle [a, b] est donnée par

Soit f : CC une fonction continue sur une courbe lisse orientée γ. Soit z : RC une paramétrisation de γ consistante avec sa direction. Alors l'intégrale le long de γ est notée

et se calcule par[6]

Cette définition est bien définie. Ainsi, le résultat est indépendant de la paramétrisation choisie[6]. Dans le cas où l'intégrale à droite n'existe pas, l'intégrale le long de γ n'est pas définie.

Comme une généralisation de l'intégrale de Riemann[modifier | modifier le code]

La généralisation de l'intégrale de Riemann aux fonctions d'une variable complexe est faite par analogie avec les fonctions d'une variable réelle. La partition d'une courbe lisse orientée γ est définie comme un ensemble de points ordonnés et fini sur γ. L'intégrale sur la courbe est la limite des sommes finies de valeurs de fonction, prise aux points de la partition, où la limite est celle où la distance maximum entre n'importe lequel des deux points successifs sur la partition (dans le plan complexe 2D), connu comme le maillage, tend vers zéro.

Exemples de contours particuliers[modifier | modifier le code]

Arc de Jordan

Un arc de Jordan est une courbe fermée simple du plan complexe qui sépare le plan en deux parties : l'intérieur et l'extérieur. Les contours les plus simples et intuitifs pour les intégrales de contour sont les arcs de Jordan et les lacets. Afin de pouvoir appliquer le théorème des résidus, on construira donc un lacet entourant les pôles de l'intégrande.

Contour de Hankel
Un contour de Hankel, parcouru dans le sens positif.
Une version totalement linéaire du contour de Hankel.

Un contour de Hankel est une courbe ouverte entourant un demi-axe du plan.

La version la plus courante consiste en deux demi-droites parallèles au demi-axe choisi, reliées par un arc de cercle centré à l'origine du demi-axe. Il est parcouru dans le sens horaire.

Ce type de contour est utilisé lorsqu'un logarithme complexe est impliqué, comme la fonction Gamma d'Euler sur le plan complexe.

Contour de Pochhamer
Un contour de Pochhammer tourne dans le sens horaire autour d'un point, puis dans le sens horaire autour d'un autre point, puis dans le sens anti-horaire autour du premier point, et enfin dans le sens anti-horaire autour du second point. La position exacte, la courbure, etc. ne sont pas essentielles dans ce cas, seule la suite de tours importe.

Un contour de Pochhammer, introduit par Camille Jordan[7] et Leo Pochhammer, est utile pour les intégrales de contour autour de deux points. En notant CA et CB les boucles parcourues autour des points A et B respectivement, alors le contour de Pochhammer est ABA−1B−1, où l'exposant −1 désigne un chemin parcouru dans le sens inverse.

Par sa construction, la courbe est homologue à 0 mais pas homotopique à 0. Son indice autour de tout point est de 0 malgré le fait qu'à l'intérieur le plan doublement épointé il ne peut pas être réduit à un seul point.

Un cycle de Pochhammer est homologue à zero : c'est la frontière de la surface verte moins la frontière de la rouge.

Méthodes directes[modifier | modifier le code]

Les méthodes directes impliquent le calcul de l'intégrale par des méthodes similaires aux calculs d'intégrales curvilignes en calcul multivarié. Ainsi, on utilise la méthode suivante :

  • paramétrer le contour
    le contour est paramétré par une fonction de variables réelles à valeurs complexes différentiable, ou en morceaux paramétrables
  • substitution de la paramétrisation dans l'intégrande
    substituer la paramétrisation dans l'intégrande permet de revenir à des intégrales de la variable réelle
  • évaluation directe

Exemple[modifier | modifier le code]

Un résultat fondamental en analyse complexe est que l'intégrale de contour de 1/z vaut 2πi, avec un chemin pris sur le cercle unité dans le sens anti-horaire (ou toute courbe de Jordan orientée positivement autour de 0). Dans le cas du cercle unité, il y a une méthode directe de calcul de l'intégrale

Dans le calcul de cette intégrale, le cercle unité est caractérisé par |z| = 1, qu'on paramètre par z(t) = eit, avec t ∈ [0, 2π], ce qui donne dz/dt = ieit, d'où

On retrouve bien la valeur attendue.

Applications de théorèmes intégraux[modifier | modifier le code]

Les théorèmes intégraux sont souvent utilisés les intégrales de contour sur un chemin, en utilisant le fait que l'intégrale d'une variable réelle est calculée en même temps que le calcul sur le chemin dans le plan complexe.

On utilise notamment la formule intégrale de Cauchy ou le théorème des résidus, de la façon suivante :

  • un contour spécifique est choisi :
    le contour est choisi de façon à contenir la partie du plan complexe correspondant à l'intégrale réelle, et entourer les singularités de l'intégrande de sorte que le théorème s'applique
  • application du théorème intégral de Cauchy
    l'intégrale est réduite à une intégration autour d'un petit cercle autour de chaque pole
  • application du théorème intégral de Cauchy ou théorème des résidus
    l'application donne alors une valeur le long du contour
  • division du contour en un contour sur les parties réelle et imaginaire
    l'entièreté du contour peut être divisée en un contour suivant la partie du plan complexe qui décrit l'intégrale réelle telle que choisie auparavant (on la note R), et la partie qui traverse le plan complexe (on la note I). L'intégrale sur le contour entier est la somme des intégrales sur ces deux contours.
  • démonstration que l'intégrale passant par le plan complexe n'a aucun rôle dans la somme
    si l'intégrale sur I peut être réduite à 0, ou si l'intégrale réelle s'avère impropre, alors si on prouve que l'intégrale sur I tend vers 0, l'intégrale sur R va tendre vers l'intégrale sur R + I.
  • conclusion
    si le point suivant est vérifié, alors on peut calculer directement R, l'intégrale réelle.

Exemple 1[modifier | modifier le code]

Illustration du chemin utilisé dans l'exemple.

On considère l'intégrale

Pour le calcul, on passe par la fonction d'une variable complexe

qui présente deux singularités en i et −i. On choisit donc un contour qui contient l'intégrale réelle, par exemple un demi-cercle dont le diamètre est formé par un intervalle sur la droite réelle (allant de a vers a, pour a réel strictement supérieur à 1) conviendrait. On note ce contour C.

Il y a deux façons de procéder, selon le théorème qu'on souhaite utiliser :

Avec la formule intégrale de Cauchy

On note que :

soit

On observe que

Comme la seule singularité contenue dans le contour est celle en i, on peut écrire

ce qui donne la fonction sous une forme qui permet l'application directe de la formule intégrale de Cauchy :

On prend la dérivée première, dans les calculs précédents, car le pole est d'ordre 2. Ainsi, (z − i) est prise au carré, donc on utilise a dérivée première de f(z). Si on avait (z − i) élevé au cube, on aurait utilisé la dérivée seconde et divisé par 2!, etc. Le cas de (z − i) à la puissance 1 correspond à une dérivée d'ordre 0 — soit f(z) elle-même.

On doit montrer que l'intégrale sur le demi-cercle tend vers zéro pour a → ∞, par le lemme d'estimation

avec M un majorant de |f(z)| sur l'arc et L désignant la longueur de l'arc. Ainsi,

Ainsi

Par le théorème des résidus

On considère la série de Laurent de f(z) en i, qui est la seule singularité à considérer. On a

Il est clair par l'observation que le résidu vaut i/4, ainsi, la formule des résidus permet de conclure :

Remarque sur le contour

On aurait pu choisir le demi-cercle dans la partie inférieure du plan complexe, qui entoure la singularité en −i. Pour avoir l'intégrale sur l'axe réel bien orienté, il faudrait parcourir le demi-cercle dans le sens horaire, ce qui change le signe de l'intégrale complète, mais pas le résultat final sur l'intégrale réelle.

Exemple 2 – Distribution de Cauchy[modifier | modifier le code]

L'intégrale

Le contour
Le contour

(qui apparait en théorie des probabilités comme un multiple scalaire de la fonction caractéristique de la loi de Cauchy) résiste aux techniques d'analyse élémentaire. On l'évalue en l'exprimant comme la limite d'intégrales de contour faites sur les demi-cercles dont le diamètre est formé par l'intervalle de la droite réelle entre a et a et sur le demi-cercle centré en 0 et de rayon a, sur le demi-plan complexe des nombres de partie imaginaire positives, parcouru dans le sens anti-horaire. On prendra a plus grand que 1, de façon que l'affixe de i soit à l'intérieur de la courbe. L'intégrale de contour devient alors

Comme eitz est une fonction entière (sans singularité sur le plan complexe), l'intégrande a des singularités aux lieux où le dénominateur z2 + 1 s'annule. La factorisation donne z2 + 1 = (z + i)(z − i), les singularités sont en z = i ou z = −i. Seul un de ces points est dans la zone délimitée par le contour. Le résidu de f(z) en z = i est

Selon le théorème des résidus, on a alors

Le contour C peut être séparé en une partie "droite" et un arc, donc

est donc

Selon le lemme de Jordan, si t > 0 alors

Ainsi, si t > 0 alors

Un argument similaire avec un arc qui tourne autour de −i au lieu de i montre que si t < 0 alors

ainsi, on obtient:

(Si t = 0 alors l'intégrale peut être calculée grâce à l'analyse réelle et vaut π.)

Exemple 3 – intégrales trigonométriques[modifier | modifier le code]

Certaines substitutions peuvent être aux intégrales impliquant des fonctions trigonométriques, ainsi l'intégrande est transformé en une fonction rationnelle d'une variable complexe et on peut appliquer les méthodes vues au-dessus pour calculer l'intégrale. On peut considérer par exemple

On veut faire une substitution de z = eit. On applique la formule d'Euler

et donc

On considère C le cercle unité, ce qui donne par substitution :

Les singularités à considérer sont en Soit C1 un petit cercle autour de et C2 un petit cercle autour de On en déduit :

Exemple 3a – intégrales trigonométriques[modifier | modifier le code]

La méthode décrite ici est utile pour les intégrales de la forme

P et Q sont des polynômes, i.e. une fonction rationnelle de termes trigonométrique est intégrée. On notera que les bornes d'intégration peuvent être π est , comme dans l'exemple, ou les extrémités tout intervalle de longueur .

L'astuce est d'utiliser la substitution z = eitdz = ieit dt et donc

Cette substitution envoie l'intervalle [0, 2π] sur le cercle unité. De plus,

et
de sorte qu'une fonction rational f(z) en z résulte de la substitution, est l'intégrale devient
qui à son tour se calcule en sommant les résidus de f(z)1/iz dans le cercle unité.

On prend par exemple l'intégrale

On doit d'abord voir que

La substitution donne

Les poles de cette fonction sont en 1 ± 2 et −1 ± 2. De ces quatre poles, 1 + 2 et −1 − 2 sont hors du cercle unité (en rouge), mais 1 − 2 et −1 + 2 (en bleu) sont à l'intérieur. Les résidus correspondants sont tous deux égaux à i2/16, ainsi l'intégrale vaut

Exemple 4 – intégrales sur branches[modifier | modifier le code]

On considère l'intégrale réelle

On peut d'abord commencer par poser l'intégrale de contour

Si les résultats utilisés pour le calcul des résidus restent les mêmes, il faut noter que z1/2 = e(Log z)/2, donc z1/2 a un point de branchement, ce qui doit être pris en compte dans le choix du contour C. Normalement le branchement du logarithme est définie comme le demi-axe réel négatif, cependant, ce choix rend le calcul de l'intégrale un peu compliqué, donc on choisira le demi-axe réel positif.

On utilise alors le contour du "trou de serrure", formé d'un petit cercle ouvert de rayon ε, s'étendant sur un segment parallèle au demi-axe jusqu'à un autre cercle ouvert.

Notons que z = −2 et z = −4 sont dans le grand cercle. Ce sont les deux poles restants, dérivables en factorisant le dénominateur de l'intégrande. Le point de branchement en z = 0 est évité en le contournant.



Soit γ le petit cercle de rayon ε, Γ le plus grand, de rayon R, alors

On peut voir que les intégrales sur Γ et γ tendent toutes deux vers 0 pour ε → 0 et R → ∞, par une estimation donnée au-dessus, ce qui laisse deux termes. Comme z1/2 = e(Log z)/2, sur le contour hors du branchement, on a gagné en argument le long de γ. (Par l'identité d'Euler, e représente le vecteur unité, donc a un log égal à π, ce qui correspond à l'argument de z. Le coefficient de 1/2 force à utiliser .) Ainsi

Ainsi :

Avec le théorème des résidus ou la formule intégrale de Cauchy (d'abord en la méthode des fractions partielles pour dériver une somme de deux intégrales de contour simples) on obtient

Exemple 5 – le carré du logarithme[modifier | modifier le code]

Cette section traite des intégrales semblables à l'exemple suivant :

Pour le calcul des intégrales, on utilise la fonction

et la branche du logarithme correspond à −π < arg(z) ≤ π.

On va calculer l'intégrale de f(z) le long du contour "trou de serrure" décrit en image. Comme il s'avère que cette intégrale est un multiple de l'intégrale initiale et par le théorème des résidus, on a

Soit R le rayon du grand cercle, et r le rayon du petit. On note le segment supérieur M, et N le segment inférieur. Comme auparavant, on prend la limite en R → ∞ et r → 0. Les contributions sur les deux cercles tendent à s'annuler. Par exemple, on peut majorer, par le lemme d'estimation avec :

Afin de calculer les contributions de M et N on pose z = −x + iε sur M et z = −x − iε sur N, avec 0 < x < +∞:

ce qui donne

Exemple 6 – logarithmes et le résidu à l'infini[modifier | modifier le code]

On veut évaluer

ce qui implique de considérer

On construit f(z) telle qu'il y ait un branchement sur [0, 3], en rouge dans le schéma. Pour ce faire, on choisit deux branchements du logarithme, ce qui donne

et

Le branchement de z3/4 est donc ]−∞, 0] et le branchement de (3 − z)1/4 est ]−∞, 3]. Il est simple de voir que le branchement du produit des deux, i.e. f(z), est [0, 3], car f(z) est en fait continue sur ]−∞, 0[. Ceci apparait pour z = −r < 0 et quand on approche le branchement par au-dessus, f(z) a la valeur

En l'approchant par en-dessous, f(z) a la valeur

Mais

donc il y a continuité à travers le branchement. C'est illustré dans le schéma, où les deux cercles noirs orientés sont notés avec la valeur correspondante de l'argument du logarithme dans z3/4 et (3 − z)1/4.

On utilise le contour en vert dans le schéma. Pour faire cela, on doit calculer la valeur de f(z) le long des segments juste au-dessus et en-dessous du branchement.

Soit z = r (à la limite, i.e. quand les deux cercles verts se réduisent vers un point), où 0 ≤ r ≤ 3. Le long du segment d'au-dessus, on trouve que f(z) a la valeur

et le long du segment d'en-dessous,

Il suit que l'intégrale de f(z)/5 − z le long du segment vaut −iI à la limite, et le long du segment d'en-dessous, I.

Si on peut montrer que les intégrale sur les deux cercles s'annulent à la limite, alors on a également la valeur de I, par le théorème des résidus. On note ρ le rayon des cercles verts, avec ρ < 0,001 et ρ → 0, pour appliquer le lemme d'estimation. Pour le cercle de gauche CL on obtient

De même, sur le cercle à droite CR, on a

Avec le théorème des résidus, on a

où le signe moins vient que le chemin autour des résidus est parcouru dans le sens horaire. En utilisant la branche du logarithme défini plus tôt, on trouve

Le pole est montré en bleu sur le schéma. La valeur se simplifie en

On utilise la formule suivante pour les résidus à l'infini :

Par substitution, on trouve

et
en utilisant l'égalité −1 = e sur la deuxième branche du logarithme. Ensuite, on applique le développement binomial, qui donne

On en conclut

Enfin, il suit que la valeur de I est de

dont on tire

Évaluation par le théorème des résidus[modifier | modifier le code]

Par le théorème des résidus, on peut évaluer des intégrales sur un contour fermé. On donne ici des exemples d'évaluation d'intégrales de contour par ce résultat.

Un premier exemple étudié sera

On rappelle que le théorème dit que

avec Res désignant le résidu de f(z).

Ici, f(z) a un unique pole en 0. De là, on peut déterminer le résidu de f(z) qui vaut 1/2.

Le théorème des résidus permet de conclure :

Intégrales de contour multivariées[modifier | modifier le code]

Pour calculer des intégrales de contour multivariées (i.e. des intégrales de surface, des intégrales de volume complexes ou des intégrales d'ordre élevé), on doit utiliser le théorème de la divergence. On supposera d'abord qu'on peut permuter avec Div, qui apparaissent dans la divergence du champ de vecteurs noté F. Ce théorème dit :

De plus, on doit aussi évaluer ∇ • F, qui est une notation alternative de div(F). La divergence de toute dimension peut être décrite comme

Exemple 1[modifier | modifier le code]

Soit le champ de vecteurs F = sin(2x) + sin(2y) + sin(2z), dans le pavé

L'intégrale de contour double correspondante sera définie comme:

On évalue maintenant ∇ • F, par l'intégrale triple correspondante :

Exemple 2[modifier | modifier le code]

On considère le champ de vecteurs F = u4 + x5 + y6 + z-3, et n désigne la quatrième dimension. On borne ce champ de vecteurs par le pavé :

On utilise donc le théorème de la divergence. Un élément de volume est donc dV = dx dy dz du.

On peut ainsi évaluer une intégrale de contour de la quatrième dimension.

Représentation intégrale[modifier | modifier le code]

Une représentation intégrale de fonction est une expression de la fonction impliquée dans une intégrale de contour. Plusieurs représentations intégrales sont connus pour de nombreuses fonctions spéciales. Les représentations intégrales peuvent être importantes pour des raisons théoriques, e.g. obtenir des prolongements analytiques ou des équations fonctionnelles, ou parfois des évaluations numériques.

Le contour de Hankel

Par exemple, la définition originale de la fonction zêta de Riemann ζ(s) par une série de Dirichlet,

est valide seulement pour Re(s) > 1. Mais

où l'intégration est faite sur le contour de Hankel H, et est valide pour tout complexe s différent de 1.

Voir aussi[modifier | modifier le code]

Références[modifier | modifier le code]

  1. John Stalker, Complex Analysis: Fundamentals of the Classical Theory of Functions, Springer, (ISBN 0-8176-4038-X, lire en ligne), p. 77
  2. Joseph Bak et Donald J. Newman, Complex Analysis, Springer, , 130–156 p. (ISBN 0-387-94756-6), « Chapters 11 & 12 »
  3. Steven George Krantz, Handbook of Complex Variables, Springer, (ISBN 0-8176-4011-8), « Chapter 2 »
  4. Dragoslav S. Mitrinović et Jovan D. Kečkić, The Cauchy Method of Residues: Theory and Applications, Springer, (ISBN 90-277-1623-4, lire en ligne), « Chapter 2 »
  5. Dragoslav S. Mitrinović et Jovan D. Kečkić, The Cauchy Method of Residues: Theory and Applications, (ISBN 90-277-1623-4, lire en ligne), « Chapter 5 »
  6. a b c d et e Edward B. Saff et Arthur David Snider, Fundamentals of Complex Analysis with Applications to Engineering, Science, and Mathematics, , 3rd éd. (ISBN 0-1390-7874-6, lire en ligne), « Chapter 4 »
  7. Jordan (1887), pp. 243–244

Liens externes[modifier | modifier le code]