Ensemble de Mandelbrot
En mathématiques, lensemble de Mandelbrot est une fractale définie comme l'ensemble des points c du plan complexe pour lesquels la suite de nombres complexes définie par récurrence par :

est bornée.

L'ensemble de Mandelbrot a été découvert par Gaston Julia et Pierre Fatou avant la Première Guerre mondiale. Sa définition et son nom actuel sont dus à Adrien Douady, en hommage aux représentations qu'en a réalisées Benoît Mandelbrot dans les années 1980. Cet ensemble permet d'indicer les ensembles de Julia : à chaque point du plan complexe correspond un ensemble de Julia différent. Les points de l'ensemble de Mandelbrot correspondent précisément aux ensembles de Julia connexes, et ceux en dehors correspondent aux ensembles de Julia non connexes. Cet ensemble est donc intimement lié aux ensembles de Julia, ils produisent d'ailleurs des formes similairement complexes.

Les images de l'ensemble de Mandelbrot sont réalisées en parcourant les nombres complexes sur une région carrée du plan complexe et en déterminant pour chacun d'eux si le résultat tend vers l'infini ou pas lorsqu'on y itère une opération mathématique. On considère la partie réelle et imaginaire de chaque nombre complexe comme des coordonnées et chaque pixel est coloré selon la rapidité de divergence, ou si elle ne diverge pas.

Les images de l'ensemble de Mandelbrot exposent une limite élaborée qui révèle progressivement des détails récursifs toujours plus fins en augmentant le grossissement. La limite de l'ensemble est constituée de plus petites versions de la forme principale, donc la propriété fractale de l'autosimilarité s'applique à l'ensemble tout entier (et pas simplement à certaines parties).

L'ensemble de Mandelbrot est devenu populaire hors des mathématiques, comme inspiration artistique et comme exemple de structure complexe venant de l'application de règles simples. C'est l'un des exemples les plus connus de visualisation mathématique.

Historique 

L'ensemble de Mandelbrot tire ses origines de la dynamique complexe, un domaine défriché par les mathématiciens français Pierre Fatou et Gaston Julia au début du .

La première représentation de cet ensemble apparaît en 1978 dans un article de Robert W. Brooks et Peter Matelski.

Le  mars 1980, au centre Thomas J. Watson de recherche IBM (dans l'État de New York), Benoît Mandelbrot obtient pour la première fois, une visualisation par ordinateur de cet ensemble. Mandelbrot étudie l'espace des paramètres des polynômes quadratiques complexes dans un article publié en 1980.

En 1984, l'étude de l'ensemble de Mandelbrot commence réellement avec les travaux d'Adrien Douady et de John H. Hubbard, qui établissent ses propriétés fondamentales et baptisent l'ensemble en l'honneur de Mandelbrot.

En 1985, les mathématiciens  et Peter Richter popularisent l'ensemble de Mandelbrot par des images de qualité et qui frappent les esprits.

Dans le numéro d'août 1985 du magazine Scientific American, l'ensemble de Mandelbrot est présenté au grand public comme  et présente l'algorithme qui permet de le tracer soi-même. La couverture de ce numéro reprend une image créée par Peitgen.

Les travaux de Douady et Hubbard coïncidaient avec un intérêt considérable pour la dynamique complexe et l'étude de l'ensemble de Mandelbrot a été le centre d'attention de ce domaine depuis. Parmi les mathématiciens qui apportèrent une contribution significative à l'étude de cet ensemble, il faut citer Tan Lei, , Curtis T. McMullen, John Milnor,  et Jean-Christophe Yoccoz.

Propriétés

Définition 
La théorie générale, développée par Pierre Fatou et Gaston Julia au début du , associe à toute fonction (suffisamment régulière) f(z, c) (à arguments et valeurs complexes) les ensembles de Julia Jc, définis (pour un c fixé) comme la frontière de l'ensemble des a complexes tels que la suite définie par z0 = a et zn+1 = f(zn, c) reste bornée (en module) ; pour la fonction particulière  on définit l'ensemble de Mandelbrot M comme l'ensemble des c pour lequel Jc est connexe. Fatou et Julia ont démontré que cette définition est équivalente à celle donnée dans l'introduction, c'est-à-dire que c appartient à M si et seulement si la suite (zn) définie par z0 = 0 et zn+1 = zn2 + c reste bornée (en module) ; restant dans les réels, il s'agit donc des points de coordonnées (a,b) tels que les deux suites (xn) et (yn) définies par la récurrence x0 = y0 = 0 et xn+1 = xn2 – yn2 + a ; yn+1 = 2xnyn + b restent bornées.

Barrière du module égal à 2 
Si la suite des modules des z est strictement supérieure à 2 pour un certain indice alors, cette suite est croissante à partir de cet indice, et elle tend vers l'infini.

Soit α la racine positive de l'équation α = α + |c| (donc α ≥ 1). En posant x = |z| – α, on a :

α + x ≥ (α + x) – |c|,

d'où x ≥ 2αx.

Par conséquent si, pour un certain indice k, |z| > α alors, à partir de cet indice, la suite (x) croît géométriquement.

Ceci a lieu en particulier, pour |c| > 2, dès que k = 1 mais aussi, pour |c| ≤ 2, dès que pour un certain k, |z| > 2.

Pour que la suite (z) soit bornée, il suffit donc qu'elle ne tende pas vers l'infini, et il faut qu'elle reste bornée par 2 (en module). Ceci fournit deux définitions de l'ensemble de Mandelbrot équivalentes à celle donnée en introduction, et prouve de plus que cet ensemble est inclus dans le disque fermé de centre 0 et de rayon 2.

Géométrie élémentaire 

L'ensemble M de Mandelbrot est compact, symétrique par rapport à l'axe réel et contient le disque fermé de centre 0 et de rayon 1/4.

Son intersection avec l'axe réel est le segment [–2, 1/4].

Notons (z(c)) la suite associée à un complexe c et pour tout entier n > 0, M l'ensemble des points c pour lesquels |z(c)| ≤ 2. D'après le paragraphe ci-dessus, M est l'intersection de la suite des M (on peut, accessoirement, en déduire aussi que cette suite est décroissante). Or M est un disque fermé donc un compact du plan, et chaque M est fermé, comme image réciproque de M par la fonction polynomiale c ↦ z(c). Il en résulte que M est compact.

Les fonctions c ↦ z(c) étant des polynômes à coefficients réels, tous les M ainsi que M sont stables par conjugaison complexe, ce qui se traduit géométriquement par la symétrie des M et de M par rapport à l'axe réel.

Si |c| ≤ 1/4 alors, par récurrence, |z(c)| ≤ 1/2 pour tout n donc c appartient à M.

Il reste à élucider la situation si c est un réel tel que 1/4 < |c| ≤ 2.
 si –2 ≤ c ≤ 0 (autrement dit si c + c ≤ – c) alors, par récurrence, c ≤ z(c) ≤ – c pour tout n donc c appartient à M.
 pour tout réel c, z(c) – z(c) ≥ c – 1/4 donc si c > 1/4 alors c n'appartient pas à M.

Remarque. Toutes les racines des polynômes c ↦ z(c) sont dans M. En effet, si z(c) = 0 alors la suite des z(c) est de période r donc bornée.

Son aire est estimée autour de  ± .

Son centre de gravité est estimé à la position réelle –0,286 768 3 ± 0,000 000 1.

La principale structure remarquable est la cardioïde centrale, de centre  et d'équation polaire . On peut aussi citer le cercle de centre  et de rayon .

Connexité 
Douady et Hubbard ont montré, en 1985, que l'ensemble était connexe. Ce résultat n'était pas évident à l'observation des premiers tracés de Mandelbrot, qui faisait apparaître des « îlots » semblant détachés du reste. Pour ce faire, ils ont montré que le complémentaire de l'ensemble de Mandelbrot est conformément isomorphe au complémentaire dans ℂ du disque unité.

On conjecture que l'ensemble de Mandelbrot est localement connexe.

Autosimilarité 

L'ensemble de Mandelbrot est autosimilaire dans le voisinage des . Ces points sont denses sur toute la frontière de l'ensemble. On conjecture qu'il est également autosimilaire, à la limite, autour des points de Feigenbaum (ex. : −1,401 155 ou −0,152 8 + 1,039 7 i).

Des versions réduites de l'ensemble de Mandelbrot apparaissent sur toute sa frontière, jusqu'à des grossissements infinis, avec de légères différences.

L'ensemble de Mandelbrot n'est pas, en général, strictement autosimilaire.

Universalité 
L'ensemble de Mandelbrot M a un caractère universel pour de nombreuses fonctions holomorphes. Des copies de M sont visibles sur les frontières de leurs bassins d'attraction, c'est-à-dire des ensembles des c pour lesquels les itérés f(f(…f(z)))) convergent vers un complexe donné. En voici quelques exemples, avec des fonctions transcendantes :
  (fonction cosinus)
  (fonction de Gauss)

On retrouve également M lors de l'itération d'une famille de fonctions cubiques telles que , par la méthode de Newton. L'ensemble des points  ne convergeant pas vers une racine de ce polynôme prend la forme de M.

Là encore, l'autosimilarité n'est pas stricte.

Dimension de Hausdorff 
La dimension de Hausdorff de la frontière de l'ensemble de Mandelbrot vaut 2. Ce résultat a été démontré en 1990 par Mitsuhiro Shishikura. On ne sait pas si cette frontière a une mesure de Lebesgue (surface) positive.

Lien avec l'équation logistique 

Les paramètres de M, dans l'intervalle réel [–2, 1/4], peuvent être mis en correspondance bijective avec ceux de l'équation logistique : 
la correspondance étant donnée par :

Lien avec les ensembles de Julia 

L'ensemble de Mandelbrot M peut être défini comme l'ensemble des complexes c pour lesquels l'ensemble de Julia correspondant, J, est connexe. Mais, plus précisément, on a (à la limite) identité entre J et M au voisinage de c, lorsque c appartient à la frontière de M. Ainsi, il y a une correspondance proche entre la géométrie de l'ensemble de Mandelbrot à un point donné et la structure de l'ensemble de Julia lui correspondant. Par exemple, un point est dans l'ensemble de Mandelbrot exactement là où l'ensemble de Julia correspondant est connexe.

Ce principe est exploité dans beaucoup de preuves et découvertes profondes sur l'ensemble de Mandelbrot. Par exemple, Mitsuhiro Shishikura a prouvé que la frontière de l'ensemble de Mandelbrot possède une dimension de Hausdorff de 2 et que l'ensemble de Julia correspondant à un point de cette frontière a aussi une dimension 2.

Jean-Christophe Yoccoz a prouvé que l'ensemble de Julia est un espace localement connexe pour certains paramètres, avant de l'établir également pour l'ensemble de Mandelbrot aux points correspondant à ces paramètres.

Bourgeons, antennes et périodicités 

L'ensemble de Mandelbrot fait apparaître nombre de structures en forme de bourgeons entourant une structure principale en forme de cardioïde.
 La cardioïde est l'ensemble des points c qui convergent vers un point fixe (période 1). Ce sont les points de la forme  pour tout  appartenant au disque unité.
 Le bourgeon principal, à gauche de la cardioïde, lui est attaché au point c = –3/4. Il s'agit d'un disque centré en c = –1 et de rayon 1/4. Il s'agit de l'ensemble des points paramètres qui, à la limite, convergent vers un cycle de période 2.
 Les autres bourgeons tangents à la cardioïde sont les points admettant d'autres périodicités. Pour tout nombre rationnel p/q, avec p et q premiers entre eux, il existe un bourgeon presque circulaire tangent à la cardioïde au point . Chacun de ces bourgeons est appelé « bourgeon p/q ». Il s'agit de l'ensemble des points paramètres convergeant vers un cycle de période q.
 Enfin, chaque bourgeon porte lui-même des bourgeons, représentatifs d'une périodicité différente, selon le même schéma. Par exemple, le bourgeon à gauche du grand bourgeon de période 2 a une périodicité de période 4. Le bourgeon immédiatement à sa gauche est de période 8, puis 16, etc., suivant ainsi le motif de doublement de période du diagramme de bifurcation de l'équation logistique.

Les bourgeons sont également surmontés de filaments en forme d'antenne. Le nombre d'antennes est directement lié à la périodicité du bourgeon. Ainsi, compter le nombre d'antennes permet de déterminer la périodicité du bourgeon.

Dessiner l'ensemble 

On a vu plus haut que dès que le module d'un zn est strictement plus grand que 2, la suite diverge vers l'infini, et donc c est en dehors de l'ensemble de Mandelbrot. Cela nous permet d'arrêter le calcul pour les points ayant un module strictement supérieur à 2 et qui sont donc en dehors de l'ensemble de Mandelbrot. Pour les points de l'ensemble de Mandelbrot, le calcul n'arrivera jamais à terme, donc il doit être arrêté après un certain nombre d'itérations déterminé par le programme.
Il en résulte que l'image affichée n'est qu'une approximation du vrai ensemble.

Bien que cela n'ait aucune importance sur le plan mathématique, la plupart des programmes générant des fractales affichent les points en dehors de l'ensemble de Mandelbrot dans différentes couleurs. La couleur attribuée à un point n'appartenant pas à l'ensemble dépend du nombre d'itérations au bout desquelles la suite correspondante est déclarée divergente vers l'infini (par exemple quand le module est strictement supérieur à 2).
Cela donne plusieurs zones concentriques, qui entourent l'ensemble de Mandelbrot. Les plus éloignées sont constituées des points c pour lesquels la suite (zn) tend « plus rapidement » vers l'infini. Ces différentes zones délimitent d'une manière plus ou moins précise l'ensemble de Mandelbrot.

Le nombre maximum d'itérations influe donc fortement sur la représentation de l'ensemble de Mandelbrot. En effet si la forme de base de l'ensemble ne varie pas outre mesure à partir de 100 itérations il n'en est pas de même pour tous les centres et tous les zooms de cet ensemble.

Optimisations

Cardioïde / bourgeon principal 

Un moyen de limiter les calculs est de trouver au préalable les points appartenant à la cardioïde et au bourgeon principal, c'est-à-dire au disque de centre (–1, 0) et de rayon 1/4. Ces points coûtent particulièrement cher en calculs, puisqu'ils appartiennent à l'ensemble de Mandelbrot, et nécessitent donc d'aller jusqu'au maximum d'itérations. Avant de passer un point de coordonnées  dans l’algorithme itératif, il faut donc vérifier :
 si  avec alors le point appartient à la cardioïde, donc à l'ensemble de Mandelbrot ;
 si alors le point appartient au disque de centre (–1, 0) et de rayon 1/4, c'est-à-dire au bourgeon principal, et donc à l'ensemble de Mandelbrot.
Les autres bourgeons n'ont pas de test équivalent, puisqu'ils ne sont pas parfaitement circulaires. Il est néanmoins possible de vérifier si un point appartient à un cercle que l'on sait inscrit dans un de ces bourgeons, permettant à la plupart des points de ces derniers de ne pas lancer l'algorithme itératif.

Algorithme par perturbation 

Il est possible de reformuler l'équation principale, de manière à raisonner sur la différence d'évolution par rapport à un point voisin. Cette méthode permet de ramener ces différences  à une faible précision (vingtaine de chiffres significatifs), dans le domaine des unités de calcul en virgule flottante matérielles de l'ordinateur, qui sont beaucoup plus rapides que les bibliothèques d'arithmétique multiprécision qui ne sont pas autant accélérées matériellement. L'ensemble est ainsi calculé à partir d'un (ou plusieurs) points de référence, calculés de manière absolue.

Si on raisonne en absolu, les nombres mis en œuvre pour le calcul de l'ensemble peuvent atteindre des milliers de chiffres significatifs pour des grands niveaux de zoom, et s'accroissant davantage avec le zoom. Avec cet algorithme, la vitesse de calcul reste relativement constante selon le niveau de zoom, la vitesse de calcul dépendant plus de la complexité de la scène.

Cet algorithme est à la base d'un des programmes de calcul le plus rapide de l'ensemble de Mandelbrot en 2017 : Kalles Fraktaler.

Principe :

Tous les , calculés à partir des  de référence mémorisés, ont une faible précision, permettant d'utiliser les unités de calcul en virgule flottante matérielles de l'ordinateur pour les opérations les plus nombreuses et coûteuses.

Le calcul peut encore être davantage optimisé en exprimant  en utilisant l'approximation :

avec  ne dépendant pas de , mais seulement de  (série du point de référence à partir duquel on calcule les perturbations).

Donc, une fois calculés et mémorisés tous les  pour le point de référence (), on peut rechercher directement pour tout point voisin (en faisant varier ), par recherche dichotomique, le  qui provoque la divergence de la suite (puisqu'on n'a pas besoin du  pour calculer ) et le calcul d'un point se fait avec une complexité en temps  au lieu de ,  étant le nombre maximum d'itérations.

Un zoom commenté 

L'ensemble de Mandelbrot doit beaucoup sa popularité à la variété et la beauté de ses structures et à la profondeur infinie de ses détails, mais aussi à la possibilité de l'explorer soi-même à l'aide des nombreux logiciels aujourd'hui disponibles.

La séquence d'exploration commentée ci-dessous est un zoom profond vers la valeur de c = –0,743643887037151 + 0,13182590420533i, à travers nombre de motifs caractéristiques. Le rapport de grossissement entre la dernière et la première image est d'environ 60 milliards (si la dernière image était en taille réelle, la première ferait une petite dizaine de fois la distance Terre-Lune).
La séquence entière (et sa suite jusqu'à un grossissement d'environ 1030) peut être visionnée sur l'animation ci-contre.

Généralisations et variantes 

L'ensemble de Mandelbrot peut être généralisé pour des puissances d supérieures à 2 pour z ↦ z + c. Ces généralisations sont parfois appelées « multibrot », mais pour certains auteurs (comme McMullen), le terme « ensemble de Mandelbrot » doit également désigner ces généralisations.

Une modification de la méthode de traçage, proposée par Melinda Green en 1993, mène à une variante baptisée « Buddhabrot ». On y voit la densité des points visités par les orbites correspondant aux valeurs de c qui divergent, donc choisies à l'extérieur de l'ensemble de Mandelbrot.

En 3 dimensions, il n'existe pas de structure de corps comparable à celle des nombres complexes, donc pas d'extension « naturelle ». Notons toutefois l'extension de Daniel White, réalisée en 2009, baptisée « Mandelbulb. »

En 4 dimensions, l'extension naturelle au corps des quaternions a été étudiée par John Holbrook en 1987. D'autres résultats théoriques ont été traités pour l'ensemble de Mandelbrot quaternionique  (voir video) defini avec un polynôme quadratique de la forme .

En 2000, Dominic Rochon a utilisé l’anneau commutatif des nombres bicomplexes pour donner une nouvelle version de l’ensemble de Mandelbrot en dimension trois et quatre. De ses travaux émerge le .

Il est aussi possible de généraliser encore, et de considérer plutôt des itérations quelconques .

Références culturelles 
L'ensemble de Mandelbrot est l'un des objets mathématiques les plus populaires et les plus connus. Il est évoqué par divers artistes, dans des livres ou des musiques notamment :
 La musique de Jonathan Coulton Mandelbrot Set est un hommage à l'ensemble de Mandelbrot lui-même et à son « père », Benoît Mandelbrot. Les paroles rappellent entre autres la définition par récurrence de l'ensemble.
 Le second livre de la série  de Piers Anthony, "Fractal Mode", décrit un monde qui est un modèle 3D parfait de l'ensemble.
 Le roman de Arthur C. Clarke "Le Fantôme venu des profondeurs" contient un lac artificiel créé pour répliquer la forme de l'ensemble de Mandelbrot.
 Le chanteur de heavy metal sud-coréen a fait un clip musical Ni pal za ya (ton destin), qui commence avec une vidéo hypnotique de l'ensemble de Mandelbrot.
 L'album Jupiter's Darling du groupe de rock américain Heart présente un ensemble de Mandelbrot sur sa couverture. L'ensemble est tourné de manière que le cardioïde soit en haut et que la figure ressemble à un cœur.

Logiciels générateurs de fractales 
 ChaosPro - (propriétaire - pour Windows)
 FractaNep (gratuit - pour Windows)
 Fractile Plus (gratuit pour iOS)
 Fractint (propriétaire - porté sur de nombreuses plates-formes)
 Gecif
 Gnofract 4D - (puissant et rapide - pour Linux, FreeBSD ou Mac OS X)
 Kalles Fraktaler (gratuit - pour Windows)
 Mandelbulb 3D (gratuit pour de nombreuses plateformes, exemple de réalisation : Video)
 XaoS (Logiciel Libre - Licence GPL, pour Windows, Mac OS X, Linux et systèmes type Unix)

Notes et références

Voir aussi 

 Liste de fractales par dimension de Hausdorff
 Mandelbox
 Mandelbulb
 Fractale de Newton

Mandelbrot (Ensemble de)
Art fractal