Banach et son microscope à trouver des points fixes

Vignette écrite par Christiane Rousseau.
Dans cette vignette, nous allons découvrir, en commençant avec un petit jeu, l’un des plus puissants théorèmes des mathématiques, à savoir le théorème du point fixe de Banach. Ce théorème a de fantastiques applications aussi bien en mathématiques que dans d’autres domaines. Dans la troisième partie de cet article nous verrons une fascinante application en compression d’images.

Mais commençons avec notre jeu, et regardons de plus près le couvercle d’une boîte de Vache Qui Rit.

La boucle d’oreille droite de la Vache est encore une boîte de Vache Qui Rit. A chaque point du couvercle, on peut associer le point correspondant sur la boucle d’oreille droite. On obtient bien sûr une fonction du couvercle dans lui-même, que l’on appelle F. Par exemple, on associe l’extrémité du menton de la Vache à l’extrémité du menton de la petite vache de la boucle d’oreille. On associe le centre de l’œil droit de la vache au centre de l’œil droit de la petite vache de la boucle d’oreille, etc. Et voilà maintenant la question: Y a-t-il un point qui est envoyé sur lui-même par ce procédé? Un tel point, s’il existe, sera appelé un point fixe. Si un point fixe existe, alors alors ce n’est pas un des points dont nous avons parlé précédemment. De plus, si un point fixe existe, alors il doit être dans la boucle d’oreille droite. Mais cette boucle d’oreille est envoyée sur la boucle d’oreille de la petite vache, etc. Visuellement, on devine que ces boucles d’oreilles droites emboîtées semblent converger vers un point, que l’on appelle A, et A est un candidat pour notre solution.

Commençons maintenant avec un point quelconque, par exemple l’extrémité du menton de la vache, et appelons-le P_0. Alors P_0 est envoyé sur P_1=F(P_0), qui est l’extrémité du menton de la petite vache sur la boucle d’oreille droite. Alors P_1 est envoyé sur P_2=F(P_1), qui est l’extrémité du menton de la vache de la boucle d’oreille droite de la petite vache, etc. On remarque trois choses:

  1. On pourrait continuer ce procédé indéfiniment et créer ainsi une suite \{P_n\}, où P_{n+1}=F(P_n).
  2. Visuellement, seul un nombre fini de points semblent distincts, et les autres ne peuvent pas être distingués les uns des autres. Bien sûr, on peut zoomer et voir plus de points. Mais peu importe la puissance du zoom, on ne distinguera qu’un nombre fini de points distincts, et les autres sembleront superposés.
  3. Cette suite semble converger vers le point A défini précédemment.

Si nous avions pris un autre point Q_0 et construit la suite \{Q_n\}, où Q_{n+1}=F(Q_n), nous pourrions voir que la suite \{Q_n\} semble également converger vers le point A. En fait on peut remarquer que cela vient du fait que le singleton \{A\} est l’intersection de toutes les boucles d’oreilles emboîtées, dont le diamètre tend vers 0.

Que nous dira le théorème du point fixe de Banach? Il nous dira que, en effet, la fonction F a un unique point fixe, c’est-à-dire il existe un unique point A du plan tel que F(A)=A. De plus, il établira également que si nous prenons un point P_0 quelconque et que nous construisons la suite \{P_n\}, où P_{n+1}=F(P_n), alors la suite \{P_n\} convergera vers A.

Pourquoi? Est-ce que cela peut arriver pour n’importe quelle fonction F? Bien sûr que non. Par exemple, une translation du plan n’a pas de point fixe. Et la fonction G(x,y)=(x+(x^2-1),y) a deux points fixes, à savoir (\pm 1,0). La fonction F de notre jeu a une propriété spéciale. C’est une contraction. En effet, l’image de la fonction est plus petite que son ensemble de définition. Si deux points P et Q sont écartés d’une certaine distance, alors leurs images F(P) et F(Q) seront plus proches l’un de l’autre que ne le sont P et Q. Cette affirmation a du sens, car dans {\mathbb R}^2 on peut mesurer la distance entre deux points. C’est parce que {\mathbb R}^2 est un espace métrique. C’est le fait que F soit une contraction qui garantit que lorsque l’on construit la suite \{P_n\}, alors peu importe le recul depuis lequel on regarde la suite (c’est-à-dire si on regarde la suite de loin, ou de près, ou à travers un microscope ou un microscope électronique), alors à partir d’un rang N qui dépend du recul choisi, tous les éléments P_n, n>N, de la suite deviennent indistincts. Deux éléments sont distincts si la distance qui les séparent est plus grande qu’un certain seuil. Nous verrons dans la prochaine partie qu’une suite ayant cette propriété est appelée une suite de Cauchy. Dans {\mathbb R}^2, toute suite de Cauchy est une suite convergente. On dit que {\mathbb R}^2 est un espace métrique complet.

1. Le théorème du point fixe de Banach

Nous avons maintenant tous les ingrédients pour le cas général, et nous pouvons énoncer le théorème.

Théorème (du point fixe de Banach). Soit \mathcal{K} un espace métrique complet dans lequel la distance entre deux points P et Q est notée d(P,Q). Et soit F:\mathcal{K}\rightarrow \mathcal{K} une contraction, c’est-à-dire il existe c\in ]0,1[ tel que pour tout P,Q\in \mathcal{K}, alors

    \[d(F(P),F(Q))\leq c  \;d(P,Q).\]

Alors F a un unique point fixe, c’est-à-dire il existe un unique A\in \mathcal{K} tel que F(A)=A.

Nous allons définir chaque terme apparaissant dans cet énoncé. Cette partie est plus formelle et peut être sautée si vous préférez vous concentrer sur les applications fascinantes.
Nous savons ce qu’est la distance entre deux points P et Q dans {\mathbb R}^2. Comment la généraliser à un ensemble \mathcal{K}?

Définition. une distance sur un ensemble \mathcal{K} est une fonction d:\mathcal{K}\times \mathcal{K}\rightarrow {\mathbb R} vérifiant

  1. Pour tout P,Q\in \mathcal{K}, d(P,Q)\geq0;
  2. d(P,Q)=0 si et seulement si P=Q;
  3. Pour tout P,Q,R\in \mathcal{K}, d(P,Q)\leq d(P,R)+d(R,Q); (inégalité triangulaire.)

Nous savons que ces propriétés sont vérifiées pour la distance euclidienne usuelle dans {\mathbb R}^2.

Nous rappelons maintenant la définition d’une suite de Cauchy, qui est la formalisation du concept de suite pour laquelle, peu importe quel seuil de précision on choisit, après un nombre fini d’éléments, tous les autres éléments deviennent indistincts. On rappelle également la définition d’une suite convergente.

Définition.

  1. Une suite \{P_n\} d’éléments d’un espace métrique \mathcal{K} est une suite de Cauchy si pour tout \epsilon>0, il existe N\in {\mathbb N} tel que pour tout n,m>N, alors

        \[d(P_n,P_m)<\epsilon.\]

  2. Une suite \{P_n\} d’éléments d’un espace métrique \mathcal{K} converge vers une limite A\in \mathcal{K} si pour tout \epsilon>0, il existe N\in {\mathbb N} tel que pour tout n>N, alors

        \[d(P_n,A)<\epsilon.\]

Définition. Un espace métrique \mathcal{K} est un espace métrique complet si toute suite de Cauchy \{P_n\} d’éléments de \mathcal{K} converge vers un élément A de \mathcal{K}.

Comment démontrer le théorème du point fixe de Banach? L’unicité du point fixe est facile. En effet, supposons que A et B soient deux points fixes. Alors F(A)=A et F(B)=B. De plus, puisque F est une contraction,

    \[d(F(A),F(B))\leq c\;d(A,B),\]

et donc, d(A,B)\leq c\;d(A,B). La seule solution est d(A,B)=0, c’est-à-dire A=B.

Pour l’existence, l’idée de la preuve est également simple: Nous l’avons déjà trouvée dans notre jeu avec la Vache Qui Rit! Prenons un point P_0\in \mathcal{K} quelconque et construisons (comme précédemment) la suite \{P_n\}, où P_{n+1}=F(P_n). Ainsi cette suite est une suite de Cauchy et sa limite est un point fixe. (Bien évidemment, démontrer ces deux assertions requiert du travail, mais nous omettrons les détails techniques. Ce qui importe est le fait que la preuve soit la même dans le cas général d’un espace métrique complet compliqué que dans le cas simple \mathcal{K}={\mathbb R}.)

L’idée de la preuve n’est pas seulement simple et intuitive, mais également très puissante. Elle nous fournit une façon de construire numériquement le point fixe A. Ce qui explique que de nombreuses applications de ce théorème puissent aussi bien être trouvées du côté théorique que du côté appliqué.

2. Les applications en analyse

Une des application théoriques très importantes du théorème du point fixe de Banach est la preuve de l’existence et de l’unicité de solutions d’équations différentielles suffisamment régulières. Dans cette application, l’espace complet \mathcal{K} est un ensemble de fonctions, et la fonction F transforme une fonction en une autre (on dit souvent que F est un opérateur). L’astuce consiste à démontrer que la solution de l’équation différentielle, si elle existe, est un point fixe de l’opérateur F.

Vous avez peut-être étudié des équations différentielles simples et appris des astuces pour trouver des formules pour les solutions. Ces équations différentielles sont l’exception, et pour la plupart des équations différentielles il n’existe pas de formule pour la solution. D’où l’importance d’un théorème affirmant l’existence et l’unicité d’une solution. Vous ne devriez pas être surpris qu’il soit impossible de donner une formule pour les solutions de la plupart des équations différentielles. En effet, considérons l’équation simple

    \[y'= e^{-x^2}.\]

Sa solution est donnée par

    \[y=\int e^{-x^2}dx.\]

Vous vous souvenez peut-être de vos cours de probabilité ou de statistiques qu’il n’existe pas de formule pour la primitive de la fonction e^{-x^2}, et donc pourquoi on utilise des tableaux lorsque l’on étudie la loi normale.

3. Application à la compression d’images

Le meilleur moyen de sauvegarder une image est d’enregistrer la couleur de chaque pixel. Il y a deux problèmes avec cette méthode:

  • Elle requiert une énorme quantité de mémoire.
  • Si on essaye d’agrandir l’image, par exemple en l’utilisant sur une grande affiche, alors les pixels deviendront de grands carrés et nous n’aurons pas assez d’informations pour compléter les détails de ces carrés.

Quel est le principe de la compression d’images? C’est d’encoder moins d’informations que l’image originale, mais en le faisant de manière intelligente, de sorte que l’oeil ne puisse pas voir que l’image observée est détériorée. Internet a augmenté le besoin de bon systèmes de compression d’images. En effet, les images ralentissent de façon significative la navigation sur internet. Ainsi, pour naviguer sur la Toile, il faut des images encodées dans des fichiers aussi petits que possible. Lorsque vous regardez cette image sur votre écran d’ordinateur, vous ne pouvez pas voir qu’elle a été détériorée. Mais si vous essayez de l’agrandir ou de l’imprimer sur une affiche, vous pourriez voir immédiatement que la qualité est mauvaise.

Il existe plusieurs procédés de compression d’images et le plus commun est le JPEG, qui est devenu standard pour les images numériques. L’encodage d’une image en format JPEG est aussi un algorithme mathématique.

Dans cette vignette nous allons nous concentrer sur une autre méthode, qui est restée plus expérimentale. Ce procédé, découvert par Barnsley, a été appelé système de fonctions itérées. L’idée sous-jacente de la méthode est d’approcher une image par des objets géométriques. Pour avoir suffisamment d’objets disponibles, on ne se limitera pas aux objets géométriques usuels comme les droites et les courbes, mais on utilisera également des figures fractales complexes comme la fougère ou le tapis de Sierpiński (voir la figure à gauche).

Nous allons expliquer l’idée du procédé de compression sur le tapis de Sierpiński à gauche. Il semble compliqué à priori. Comment le stocker dans la mémoire d’un ordinateur de façon économique? Le mieux est d’enregistrer un programme qui le reconstruira lorsqu’on en aura besoin.

Et pour construire ce programme on a besoin de comprendre ce qui caractérise cet objet géométrique. Observons notre tapis de Sierpiński: c’est l’union de trois tapis de Sierpiński (c’est-à-dire 3 copies de lui-même), qui sont deux fois plus petits (en largeur et en hauteur). En effet, en partant du tapis de Sierpiński, on peut en construire un second avec le procédé suivant:

  • On réduit le tapis de Sierpiński de moitié en partant du sommet inférieur gauche.
  • On fabrique une seconde copie de ce demi tapis de Sierpiński et on le colle à droite du premier.
  • On fabrique une troisième copie de ce demi tapis de Sierpiński et on le colle au-dessus des deux autres.

La seconde figure que nous avons construite est identique au tapis de Sierpiński initial. Donc, le tapis de Sierpiński est un point fixe de ce procédé.

Écrivons cela en termes mathématiques. Remarquons que la base du tapis de Sierpiński est égale à sa hauteur. On peut donc construire un repère dont l’origine se situe au sommet inférieur gauche du tapis de Sierpiński, et dont les unités des axes sont telles que la base et la hauteur ont toute deux pour longueur 1. On construit également les transformations affines suivantes définies sur {\mathbb R}^2:

    \begin{align*} T_1(x,y)&=\left(\frac{x}2,\frac{y}2\right),\\ T_2(x,y)&=\left(\frac{x}2+\frac12,\frac{y}2\right),\\ T_3(x,y)&=\left(\frac{x}2+\frac14,\frac{y}2+\frac12\right).\end{align*}

Si S est le tapis de Sierpiński, alors on a

    \[S=T_1(S)\cup T_2(S)\cup T_3(S).\]

Y a-t-il d’autres sous-ensembles B du plan qui ont la même propriété, à savoir

(1)   \begin{equation*}B=T_1(B)\cup T_2(B)\cup T_3(B)?\end{equation*}

Nous allons faire une expérience et voir que non! Nous aurons donc caractérisé notre tapis de Sierpiński comme étant le seul sous ensemble B du plan tel que (1) soit satisfaite. Qu’avons nous fait? Nous avons construit une fonction qui associe à un sous ensemble B du plan le sous ensemble T_1(B)\cup T_2(B)\cup T_3(B). Appelons cette fonction W. Elle est définie par

(2)   \begin{equation*}B\mapsto W(B)=T_1(B)\cup T_2(B)\cup T_3(B),\end{equation*}

et nous avons vu que S=W(S), i.e. S est un point fixe de cette fonction.

Nous avons dit que nous allions faire une expérience pour montrer que le tapis de Sierpiński est le seul point fixe de cette fonction. Essayons avec un carré C_0 comme dans la Figure 1(a). Son image est C_1 (Figure 1(b)). On applique le même procédé à C_1 et on obtient C_2, C_3C_5. (Figure 1 (c)-(f)).

Figure 1

On observe trois choses:

  1. Aucun des ensembles C_0, …, C_5 n’est un point fixe pour W.
  2. Nous pourrions continuer indéfiniment ce procédé, et ainsi obtenir une suite infinie d’ensembles \{C_n\}, où C_{n+1}=W(C_n).
  3. La suite \{C_n\} semble converger rapidement vers le tapis de Sierpiński. En effet, notre oeil ne peut distinguer C_{10} de S. Donc, à la place de S sur notre image, le programme de reconstruction peut simplement produire C_{10}. Et si une meilleure résolution était nécessaire, on utiliserait le même programme en lui demandant de s’arrêter à C_{20}, ou C_{30}. Le même minuscule programme peut reconstruire S avec une précision arbitraire!

De plus, on peut refaire cette expérience et voir que cela fonctionne avec n’importe quel ensemble de départ! Un deuxième exemple avec un pentagone est visible dans la Figure 2. Les mêmes remarques 1, 2 et 3 comme ci-dessus peuvent s’appliquer à cet exemple.

Figure 2

Nous avons vu que le théorème du point fixe de Banach s’applique aux contractions des espaces métriques complets. En (2) nous avons défini la fonction W sur les sous-espaces du plan. Pour l’espace métrique, nous allons prendre \mathcal{K} comme étant l’ensemble des sous-ensembles (fermés) bornés du plan. Nous allons introduire une distance sur \mathcal{K} appelée la distance de Hausdorff. La définition de la distance de Hausdorff d_H(B_1,B_2), entre deux sous-ensembles B_1 et B_2 est une formule obscure et compliquée, donc nous expliquerons cette notion d’une façon différente. Commençons par expliquer ce que nous entendons par

    \[d_H(B_1,B_2)\leq \epsilon\]

(même si nous n’avons pas défini ce que d_H(B_1,B_2) signifie!). Cela veut dire que si notre oeil a une précision de \epsilon, alors il ne peut pas distinguer B_1 de B_2. En termes mathématiques, alors d_H(B_1,B_2) \leq \epsilon signifie

(3)   \begin{equation*}\forall P\in B_1\;\exists Q\in B_2 \: d(P,Q)\leq \epsilon \quad\text{et}\quad \forall P'\in B_2\;\exists Q'\in B_1 \: d(P',Q')\leq \epsilon. \end{equation*}

(ici d est la distance euclidienne usuelle sur \mathbb R^2.) Ceci nous permet de donner la définition indirecte suivante.

Définition. La distance de Hausdorff entre deux ensembles fermés et bornés B_1 et B_2 est le minimum de tous les \epsilon\geq0 tels que (3) soit vérifiée.

On peut maintenant admettre que la fonction W est une contraction. En effet, supposons que d_H(B_1,B_2)\leq \epsilon. Alors, on peut prouver que d_H(W(B_1),W(B_2))\leq \frac{\epsilon}{2}. Soit P\in W(B_1). Alors, il existe i\in \{1,2,3\} et P_1\in B_1 tel que P=T_i(P_1). Puisque d_H(B_1,B_2)\leq \epsilon, il existe Q_1\in B_2 de sorte que d(P_1,Q_1)<\epsilon. Soit Q=T_i(Q_1). Alors Q\in W(B_2) et d(P,Q)= d(T_i(P_1),T_i(Q_1))=\frac12 d(P_1,Q_1)\leq\frac{\epsilon}{2}. De même, si on commence avec P'\in W(B_2) alors il existe Q'\in W(B_1) tel que d(P',Q')\leq\frac{\epsilon}{2}.

Ce procédé a été adapté à la compression de vraies images (voir [K] ou [RS]). La méthode produit des images de grande qualité lorsque l’image originale a un caractère fractal. Cependant, le niveau de compression n’est pas aussi flexible et performant qu’avec le format JPEG. De plus, le procédé d’encodage (transformant l’image en un programme censé la reconstruire) est toujours trop complexe pour avoir un intérêt pratique. Cependant, la simplicité de cette idée combinée avec sa puissance reste impressionnante et très séduisante.

4. Une application surprenante: l’algorithme PageRank

Le succès de Google en tant que moteur de recherche vient de son algorithme: l’algorithme PageRank. Dans cet algorithme, on calcule un point fixe d’un opérateur linéaire sur {\mathbb R}^n qui est une contraction, et ce point fixe (qui est un vecteur) engendre le classement des pages. En pratique, le point fixe (qui est un vecteur propre pour la valeur propre 1) est calculé de façon approximative par P_n pour n suffisamment grand. Nous invitons le lecteur intéressé à regarder plus en détail la vignette Comment Google fonctionne.

5. Conclusion

Nous avons vu dans cette vignette comment, en commençant par un simple jeu, on peut découvrir des idées très puissantes qui peuvent conduire à des avancées majeures pour les mathématiques et la technologie. Lorsque nous cherchons une solution unique à un problème, il est désormais normal, dans de nombreux domaines des mathématiques, d’essayer de voir si la solution au problème peut être caractérisée comme l’unique point fixe d’un opérateur construit spécialement à cette fin.

Nous avons vu que l’avantage de cette approche est que le théorème fournit une façon efficace et concrète de construire la solution comme la limite d’une suite, puisque la convergence est rapide.

L’analyse est l’étude des fonctions. Les fonctions sont généralement définies sur les nombres. En analyse à plusieurs variables, on généralise la notion de fonction aux vecteurs, qui sont des éléments de {\mathbb R}^n. Mais pourquoi s’arrêter aux éléments de {\mathbb R}^n? Nous avons aussi vu que les mathématiciens aiment généraliser la notion de fonctions, et s’autorisent à les définir, par exemple, sur des ensembles d’ensembles, des ensembles de fonctions, etc. L’analyse portant sur les ensembles de fonctions est désormais devenue un chapitre important de l’analyse moderne, appelé analyse fonctionnelle, qui est standard dans les études de niveau master.

Nous vous invitons à faire le lien avec des procédés itératifs que vous avez rencontrés. Par exemple avec les procédés itératifs en dimension 1 associés aux suites de Héron pour obtenir des racines carrées. La convergence rapide des suites géométriques peut également être comprise du point de vue expliqué dans cette vignette.

Bibliographie

[B] M. F. Barnsley, Fractals everywhere, San Diego, Academic Press, 1988.

[K] J. Kominek, Advances in fractal compression for
multimedia applications, Multimedia Systems Journal, vol. 5,
n^{o} 4, 1997, 255–270.

[R] C. Rousseau, Point fixe de Banach, Accromath 5, hiver-printemps 2010 (www.accromath.ca).

[R2] C. Rousseau, Comment Google fonctionne: chaînes de Markov et valeurs propres, Klein vignette (www.kleinproject.org).

[RS] C. Rousseau et Y. Saint-Aubin, Mathématiques et technologie, SUMAT Series, Springer-Verlag, 2008 .

Ce post est disponible en: Anglais, Italien, Espagnol, Arabe, Khmère, Portugais - du Brésil

This entry was posted in Vignettes. Bookmark the permalink.

  1. Microscope says:

    Merci pour cet article vraiment complète. Technique mais très instructif.

Leave a Reply

Your email address will not be published. Required fields are marked *