Voir la source de Fractale de Newton

[[image:Julia-set_N_z3-1.png|thumb|400px|Fractale de Newton {{math|''z''{{3}} -1}} et les 3 bassins d'attraction des racines du polynôme (en couleur).]]
La '''fractale de Newton''' est un [[Frontière (topologie)|ensemble frontière]] défini dans le [[plan complexe]] caractérisé par l’application de la [[méthode de Newton]] à un [[polynôme]] <math>\textstyle{p(z), z\in\mathbb{C}.}</math>

== Définition ==
La fractale de Newton est l’[[ensemble de Julia]] d’une [[fonction méromorphe]] <math>z\mapsto z-\tfrac{p(z)}{p'(z)}</math> qui est donnée par la méthode de Newton. Lorsqu’il n'y a pas de cycles attractifs, il divise le plan complexe en régions {{mvar|G{{ind|k}}}}, chacune d’elles associée à chaque racine de ce polynôme.

La '''fractale de Newton classique''' est ainsi associée au polynôme {{math|''z''{{3}} -1}} et divise le plan en trois régions associées à ses trois racines : <math>{1, -\frac{1}{2}+\frac{\sqrt{3}}{2}\mathrm{i}}</math> et <math>{-\frac{1}{2}-\frac{\sqrt{3}}{2}\mathrm{i}}</math>.

=== Construction ===
De nombreux points du plan complexe sont associés à chaque racine de la manière suivante :

Un point {{math|''z''{{ind|0}}}} du plan complexe est choisi comme point de départ. On applique la méthode itérative de Newton :
:<math>z_{n+1}=z_n-\frac{p(z_n)}{p'(z_n)}.</math>

En particulier, la fractale de Newton classique s'obtient en itérant :
:<math>z_{n+1}=z_n-\frac{z_n^3-1}{3z_n^2}=\frac{2z_n^3+1}{3z_n^2}.</math>

Cette règle mène à une suite de points {{math|''z''{{ind|1}}}}, {{math|''z''{{ind|2}}}} {{etc.}} Si la suite converge vers la racine {{mvar|R{{ind|k}}}} du polynôme, alors {{math|''z''{{ind|0}}}} appartient à la région {{mvar|G{{ind|k}}}}. Cette région est aussi appelée « bassin d'attraction de la racine {{mvar|R{{ind|k}}}} ».

Toutefois, pour tout polynôme de degré égal au moins à 2, il existe des points pour lesquels la suite de Newton ne converge pas, c’est le cas de la frontière des bassins d’attraction de chaque racine.

=== Structure fractale ===
La fractale de Newton présente, à l’instar de toute [[fractale]], une apparence complexe, malgré une description simple, et des auto-similarités visibles à toutes échelles (voir zoom successifs ci-dessous).
<gallery>
Image:Newton1.png|Newton {{math|''z''{{3}} -1}}.
Image:Newton2.png|{{1er|zoom}}.
Image:Newton3.png|{{2e|zoom}}.
</gallery>

Elle suggère également que la méthode de Newton peut être très sensible aux conditions initiales et que deux points initiaux infiniment proches peuvent converger vers des racines différentes.

Elle montre, enfin, que chaque point de la fractale de Newton est un point-frontière multiple, séparant chacun des ''n'' bassins d'attraction. Si deux points infiniment proches convergent vers deux racines distinctes, alors il existe un troisième point, infiniment proche également, qui converge vers la troisième racine. Voir l'article sur les [[lacs de Wada]].

=== Généralisation ===

Une généralisation de l’itération de Newton est :
:<math>z_{n+1}=z_n- a \frac{p(z_n)}{p'(z_n)} </math>
où {{mvar|a}} est un nombre complexe<ref>{{lien web|langue = en |url = http://www.chiark.greenend.org.uk/~sgtatham/newton/ | titre = Fractals derived from Newton-Raphson iteration | auteur = Simon Tatham }}.</ref>. Le cas particulier {{math|''a'' {{=}} 1}} correspond à la fractale de Newton classique.
 
Les points fixes de cette transformation sont stables si {{mvar|a}} appartient au disque centré en 1 de rayon 1. Hors de ce disque les points fixes sont localement instables, toutefois la transformation présente une structure fractale au sens de l’[[ensemble de Julia]]. Si {{mvar|p}} est un polynôme de degré {{mvar|n}}, alors la suite {{mvar|z{{ind|n}}}} est bornée tant que {{mvar|a}} reste dans le disque de rayon {{mvar|n}} centré en {{mvar|n}}.

=== Autres méthodes ===
En analyse numérique, nombre de [[Algorithme de recherche d'un zéro d'une fonction|méthodes de résolution d'équation]] existent.

Les fractales associées partagent des caractéristiques communes avec la fractale de Newton : la frontière triple, des auto-similarités à toutes les échelles, et trois bassins d'attraction non connexes (en couleur). Selon les conditions initiales choisies, la méthode de la sécante crée des zones de non-convergence.

Voir les exemples ci-dessous, appliqués à la fonction polynomiale {{math|''z''{{3}} -1}} :
{| class="wikitable"
! Méthode || Formule || Convergence || Illustration || Remarques
|-
| [[Méthode de la sécante]] || <math>{z_{n+1} = z_n - \frac{z_n-z_{n-1}}{p(z_n)-p(z_{n-1})} p(z_n). }</math> || 1,618  || [[image:Secant fractal.jpg|150px|]] || La [[méthode de la sécante]] permet de s'affranchir du calcul de la dérivée en approximant la dérivée <math>p'(z_n)</math> par <math>\textstyle{\frac{p(z_n)-p(z_{n-1})}{z_n-z_{n-1}}}</math>.
Dans l'illustration on a posé <math>z_{-1}</math> proche de <math>z_0</math>.
|-
| [[Méthode de Newton]] || <math>{z_{n+1} = z_n - \frac{p(z_n)}{p'(z_n)}}</math> || quadratique  || [[image:Newton method fractal.jpg|150px|]] ||
|-
| [[Itération de Householder|Méthode de Householder]] || <math>{z_{n+1} = z_n - \frac{p(z_n)}{p'(z_n)}\times (1+h_n)}</math> avec <math>{h_n = \frac{p(z_n) p''(z_n)}{2 p'(z_n)^2}}</math> || cubique|| [[image:Householder method fractal.jpg|150px|]] || Les méthodes de Householder généralisent les méthodes de Newton et de Halley.
|-
| [[Méthode de Halley]] || <math>{z_{n+1} = z_n - \frac {2 p(z_n) p'(z_n)} {2 {[p'(z_n)]}^2 - p(z_n) p''(z_n)}} </math>  || cubique  || [[image:Halley method fractal.jpg|150px|]] ||
|}

<gallery caption="Fractales de Newton" >
FRACT008.png|Fractale de Newton pour le polynôme {{math|''p''(''z'') {{=}} ''z''{{3}} -1}}, coloré en fonction du nombre d’itérations de convergence.
Newtroot 1 0 0 m1.png|Fractale de Newton pour le polynôme {{math|''p''(''z'') {{=}} ''z''{{3}} -1}},  colorée par racine atteinte.
Newtroot 1 0 0 0 0 m1.png|Fractale de Newton pour le polynôme {{math|''p''(''z'') {{=}} ''z''{{5}} -1}},  colorée par racine atteinte.
Newton-Fractal 2z3-2z+2.png| Fractale de Newton pour {{math|''p''(''z'') {{=}} ''z''{{3}} - 2''z'' + 2}}. Les points en rouge n’atteignent aucune racine.
Newtroot 1 0 m3i m5m2i 3 1.png| Fractale de Newton pour {{math|''p''(''z'') {{=}} ''z''{{5}} - 3i''z''{{3}} -(5+2i)''z''{{2}} + 3''z'' + 1}}, coloré par racine atteinte.
Image:Newton_z3-2z+2.png|Fractale de Newton pour {{math|''p''(''z'') {{=}} ''z''{{3}} - 2''z'' + 2}}. Les points en rouge n’atteignent aucune racine.
Image:Colored Newton Fractal 2.png|Fractale de Newton pour un polynôme du {{7e}} degré, coloré par racine atteinte nuancé en fonction du nombre d’itérations de convergence.
Image:Sin(x)_detail.png|Une autre fractale de Newton pour {{math|sin(''x'')}}.
Image:Newton z6 z3.jmb.jpg|<center><math>p(z) = z^6 + z^3 - 1</math>.</center>
Image:Newton SINUS.jmb.jpg|<center><math>p(z) = \sin(z)- 1</math>.</center>
Image:Newton COSH.jmb.jpg|<center><math>p(z) = \cosh(z)- 1</math>.</center>
</gallery>

== Voir aussi ==
{{Autres projets|commons=Newton fractal|commons titre=Fractale de Newton}}
* [[Fractale]]
* [[Ensemble de Julia]]
* [[Ensemble de Mandelbrot]]
* [[Liste de fractales]]

== Références ==
<references/>
* [http://www.math.sunysb.edu/~scott/ ''{{Langue|en|How to Find All Roots of Complex Polynomials by Newton's Method}}''] par J. H. Hubbard, D. Schleicher, S. Sutherland, ''Inventiones Mathematicae'' vol. 146 (2001) – Avec une discussion sur la structure globale des fractals de Newton.
* [http://citeseer.ist.psu.edu/313081.html ''{{Langue|en|On the Number of Iterations for Newton's Method}}''] par Dierk Schleicher, {{date-|21 juillet 2000}}
* [http://www.math.sunysb.edu/cgi-bin/thesis.pl?thesis06-1 ''{{Langue|en|Newton's Method as a Dynamical System}}''] par Johannes Rueckert
* [http://images.math.cnrs.fr/La-methode-de-Newton-et-son.html ''La méthode de Newton et son fractal''] Une approche didactique, par [[Tan Lei]] (CNRS).

{{Palette|Fractales}}
{{Portail|mathématiques}}

[[Catégorie:Fractale]]
[[Catégorie:Calcul numérique]]