Fractal de Newton

El fractal de Newton es un conjunto de contornos en el plano complejo que se caracteriza por el método de Newton aplicado a un polinomio fijo $p (z) \in$ $[$ $z$ $]$ o función trascendental . Es el conjunto de Julia de la función meromórfica $z$ $\mapsto$ $z$ $-$ $⁠$ $\mathbb {C}$ $.mw-parser-output .sfrac{white-space:nowrap}.mw-parser-output .sfrac.tion,.mw-parser-output .sfrac .tion{display:inline-block;vertical-align:-0.5em;font-size:85%;text-align:center}.mw-parser-output .sfrac .num{display:block;line-height:1em;margin:0.0em 0.1em;border-bottom:1px solid}.mw-parser-output .sfrac .den{display:block;line-height:1em;margin:0.1em 0.1em}.mw-parser-output .sr-only{border:0;clip:rect(0,0,0,0);clip-path:polygon(0px 0px,0px 0px,0px 0px);height:1px;margin:-1px;overflow:hidden;padding:0;position:absolute;width:1px}p ( z )/p′ ( z )⁠$ que se da por el método de Newton. Cuando no hay ciclos atractivos (de orden mayor que 1), divide el plano complejo en regiones $G k$ , cada una de las cuales está asociada con una raíz $ζ k$ del polinomio, $k = 1, \dots, deg(p)$ . De esta manera, el fractal de Newton es similar al conjunto de Mandelbrot , y como otros fractales exhibe una apariencia intrincada que surge de una descripción simple. Es relevante para el análisis numérico porque muestra que (fuera de la región de convergencia cuadrática ) el método de Newton puede ser muy sensible a su elección del punto de inicio.

Casi todos los puntos del plano complejo están asociados con una de las raíces $deg(p)$ de un polinomio dado de la siguiente manera: el punto se utiliza como valor inicial $z 0$ para la iteración de Newton $z n + 1 := z n - ⁠ p (zn) / p' (z n) ⁠$ , dando una secuencia de puntos $z 1, z 2, \dots,$ Si la secuencia converge a la raíz $ζ k$ , entonces $z 0$ era un elemento de la región $G k$ . Sin embargo, para cada polinomio de grado al menos 2 hay puntos para los cuales la iteración de Newton no converge a ninguna raíz: ejemplos son los límites de las cuencas de atracción de las diversas raíces. Incluso hay polinomios para los cuales conjuntos abiertos de puntos de partida no convergen a ninguna raíz: un ejemplo simple es $z 3 - 2 z + 2$ , donde algunos puntos son atraídos por el ciclo $0, 1, 0, 1\dots$ en lugar de por una raíz.

Un conjunto abierto cuyas iteraciones convergen hacia una raíz o ciclo dado (que no es un punto fijo), es un conjunto de Fatou para la iteración. El conjunto complementario a la unión de todos ellos es el conjunto de Julia. Los conjuntos de Fatou tienen un límite común, es decir, el conjunto de Julia. Por lo tanto, cada punto del conjunto de Julia es un punto de acumulación para cada uno de los conjuntos de Fatou. Es esta propiedad la que causa la estructura fractal del conjunto de Julia (cuando el grado del polinomio es mayor que 2).

Para trazar imágenes del fractal, primero se puede elegir un número específico $d$ de puntos complejos $(ζ 1, \dots, ζ d)$ y calcular los coeficientes $(p 1, \dots, p d)$ del polinomio

p(z)=z^{d}+p_{1}z^{d-1}+\cdots +p_{d-1}z+p_{d}:=(z-\zeta _{ 1})(z-\zeta _{2})\cdots (z-\zeta _{d})

Luego, para una red rectangular

z_{mn}=z_{00}+m\,\Delta x+in\,\Delta y;\cuadrado m=0,\ldots ,M-1;\cuadrado n=0,\ldots ,N-1

de puntos en , se encuentra el índice $k$ $($ $m$ $,$ $n$ $)$ de la raíz correspondiente $ζ$ $k$ $($ $m$ $,$ $n$ $)$ y se utiliza esto para llenar una cuadrícula ráster $M$ $\times$ $N$ asignando a cada punto $($ $m$ $,$ $n$ $)$ un color $f$ $k$ $($ $m$ $,$ $n$ $)$ . Adicionalmente o alternativamente los colores pueden depender de la distancia $D$ $($ $m$ $,$ $n$ $)$ , que se define como el primer valor $D$ tal que $|$ $z$ $D$ $-$ $ζ$ $k$ $($ $m$ $,$ $n$ $)$ $| <$ $ε$ para algún pequeño $ε$ $> 0$ previamente fijado . $\mathbb {C}$

Generalización de los fractales de Newton

Una generalización de la iteración de Newton es

z_{n+1}=z_{n}-a{\frac {p(z_{n})}{p'(z_{n})}}

donde $a$ es cualquier número complejo . ^[1] La opción especial $a = 1$ corresponde al fractal de Newton. Los puntos fijos de este mapa son estables cuando $a$ se encuentra dentro del disco de radio 1 centrado en 1. Cuando $a$ está fuera de este disco, los puntos fijos son localmente inestables, sin embargo el mapa aún exhibe una estructura fractal en el sentido del conjunto de Julia . Si $p$ es un polinomio de grado $d$ , entonces la secuencia $z n$ está acotada siempre que $a$ esté dentro de un disco de radio $d$ centrado en $d$ .

De manera más general, el fractal de Newton es un caso especial de un conjunto de Julia .

Fractal de Newton para tres raíces de grado 3 $p (z) = z 3 - 1$ , coloreado según el número de iteraciones requeridas
Fractal de Newton para tres raíces de grado 3 $p (z) = z 3 - 1$ , coloreado por la raíz alcanzada
Fractal de Newton para $p (z) = z 3 - 2 z + 2$ . Los puntos en las cuencas rojas no alcanzan una raíz.
Fractal de Newton para un polinomio de séptimo orden, coloreado por la raíz alcanzada y sombreado por la tasa de convergencia.
Fractal de Newton para $p (z) = z 8 + 15 z 4 - 16$
Fractal de Newton para $p (z) = z 5 - 3 iz 3 - (5 + 2 i) z 2 + 3 z + 1$ , coloreado por la raíz alcanzada, sombreado por el número de iteraciones requeridas.
Fractal de Newton para $p (z) = sen z$ , coloreado por la raíz alcanzada, sombreado por el número de iteraciones requeridas
$Otro fractal de Newton para p(z) = sen z$
Otro fractal de Newton para $p (z) = sen z$
Fractal de Newton generalizado para $p (z) = z 3 - 1$ , $a = - ⁠ 1 / 2$ El color fue elegido en base al argumento después de 40 iteraciones $.$
Fractal de Newton generalizado para $p (z) = z 2 - 1$ , $a = 1 + i$ .
Fractal de Newton generalizado para $p (z) = z 3 - 1$ , $a = 2$ .
Fractal de Newton generalizado para $p (z) = z 4 + 3 i - 1$ , $a = 2,1$ .
$p (z) = z6 + z3-1 $
$p (z) = sen z - 1$
$p (z) = sen z - 1$
$p (z) = cosh z - 1$
$p (z) = cosh z - 1$
$p (z) = z20-2z + 2 $

Serie: $p (z)= zn - 1$

$p (z) = z 3 - 1$ , a=1
$p (z) = z 3 - 1$ , a=2
$p (z) = z4-1,$ a= $1$
$p (z) = z 4 - 1$ , a=2
$p (z) = z 5 - 1$ , a=1
$p (z) = z 5 - 1$ , a=2
$p (z) = z 6 - 1$ , a=1
$p (z) = z 7 - 1$ , a=1
$p (z) = z 8 - 1$ , a=1
$p (z) = z 10 - 1$ , a=1

Otros fractales donde se multiplican funciones potenciales y trigonométricas. $p (z) = z n *Sin(z) - 1$

$p (z) = z 2 *Sin(Z) - 1$ , a=1
$p (z) = z 2 *Sin(Z) - 1$ , a=1 (zoom)
$p (z) = z 3 *Sin(Z) - 1$ , a=1
$p (z) = z 4 *Sin(Z) - 1$ , a=1
$p (z) = z 4 *Sin(Z) - 1$ , a=1 (zoom)
$p (z) = z 5 *Sin(Z) - 1$ , a=1
$p (z) = z 6 *Sin(Z) - 1$ , a=1
$p (z) = z 6 *Sin(Z) - 1$ , a=1 (zoom)

Fractal nuevo

El fractal Nova inventado a mediados de la década de 1990 por Paul Derbyshire, ^[2]^[3] es una generalización del fractal de Newton con la adición de un valor $c$ en cada paso: ^[4]

z_{n+1}=z_{n}-a{\frac {p(z_{n})}{p'(z_{n})}}+c=G(a,c,z)

La variante "Julia" del fractal Nova mantiene $c$ constante sobre la imagen e inicializa $z 0$ en las coordenadas de los píxeles. La variante "Mandelbrot" del fractal Nova inicializa $c$ en las coordenadas de los píxeles y establece $z 0$ en un punto crítico, donde ^[5]

{\frac {\partial }{\partial z}}G(a,c,z)=0.

Los polinomios de uso común como $p (z) = z 3 - 1$ o $p (z) = (z - 1) 3$ conducen a un punto crítico en $z = 1$ .

Fractal animado "Julia" Nova para $p (z) = z 3 - 1$ con $c$ yendo de −1 a 1, coloreado por la raíz alcanzada.
Fractal animado "Julia" Nova para $p (z) = z 3 - 1$ con $c = ⁠ 1 / 2 ⁠ e iφ$ y $φ$ van de 0 a 2 $π$ , coloreados por la raíz alcanzada.

Implementación

Para implementar el fractal de Newton es necesario tener una función inicial así como su función derivada:

{\begin{aligned}f(z)&=z^{3}-1\\f'(z)&=3z^{2}\end{aligned}}

Las tres raíces de la función son

z=1,\ -{\tfrac {1}{2}}+{\tfrac {\sqrt {3}}{2}}i,\ -{\tfrac {1}{2}}-{ \tfrac {\sqrt {3}}{2}}i

Las funciones definidas anteriormente se pueden traducir en pseudocódigo de la siguiente manera:

//z^3-1 float2 Función ( float2 z ) { return cpow ( z , 3 ) - float2 ( 1 , 0 ); //cpow es una función exponencial para números complejos }         //3*z^2 float2 Derivative ( float2 z ) { return 3 * cmul ( z , z ); //cmul es una función que maneja la multiplicación de números complejos }

Ahora sólo es cuestión de implementar el método de Newton utilizando las funciones dadas.

float2 roots [ 3 ] = //Raíces (soluciones) del polinomio { float2 ( 1 , 0 ), float2 ( -.5 , sqrt ( 3 ) / 2 ), float2 ( -.5 , - sqrt ( 3 ) / 2 ) }; color colors [ 3 ] = //Asigna un color para cada raíz { rojo , verde , azul }              Para cada píxel ( x , y ) en el objetivo , haga : { zx = coordenada x escalada del píxel ( escalada para estar en la escala X de Mandelbrot ( -2,5 , 1 )) zy = coordenada y escalada del píxel ( escalada para estar en la escala Y de Mandelbrot ( -2 , 1 ))                                          float2 z = float2 ( zx , zy ); //z se establece originalmente en las coordenadas de píxeles     para ( int iteración = 0 ; iteración < maxIteración ; iteración ++ ;) { z -= cdiv ( Función ( z ), Derivada ( z )); //cdiv es una función para dividir números complejos             float tolerancia = 0.000001 ; for ( int i = 0 ; i < roots . Length ; i ++ ) { float2 diferencia = z - raíces [ i ]; //Si la iteración actual está lo suficientemente cerca de una raíz, colorea el píxel. if ( abs ( diferencia . x ) < tolerancia && abs ( diferencia . y ) < tolerancia ) { return colores [ i ]; //Devuelve el color correspondiente a la raíz } } } return negro ; //Si no se encuentra ninguna solución }

Véase también

Referencias

Wikimedia Commons tiene medios relacionados con Fractales de Newton .

^ Simon Tatham. "Fractales derivados de Newton-Raphson".
^ Damien M. Jones. "clase Standard_NovaMandel (referencia de fórmula Ultra Fractal)".
^ Damien M. Jones. "Fractales de nova de dmj 1995-6".
^ Michael Condron. "El método relajado de Newton y el fractal Nova".
^ Frederik Slijkerman. "Manual Ultra Fractal: Nova (Julia, Mandelbrot)".

Lectura adicional

JH Hubbard, D. Schleicher, S. Sutherland: Cómo encontrar todas las raíces de polinomios complejos mediante el método de Newton , Inventiones Mathematicae vol. 146 (2001), con un análisis de la estructura global de los fractales de Newton
Sobre el número de iteraciones del método de Newton por Dierk Schleicher 21 de julio de 2000
El método de Newton como sistema dinámico de Johannes Rueckert
El fractal de Newton (del que Newton no sabía nada) por 3Blue1Brown , junto con una demostración interactiva del fractal en su sitio web y el código fuente para la demostración.