es Fractal de Newton

El fractal de Newton es una frontera en el plano complejo delimitada mediante el método de Newton aplicado a un polinomio fijo $p (Z) \in ℂ[Z]$ o a una función trascendente. Es el conjunto de Julia de la función meromorfa $z \mapsto z - p (z) / p' (z)$ , que viene dado 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, al igual que otros fractales, exhibe una apariencia intrincada que surge de una descripción simple. Es relevante en análisis numérico porque muestra que (fuera de la región de orden de convergencia) el método de Newton puede ser muy sensible a la elección del punto de inicio.

Muchos puntos del plano complejo están asociados con una de las raíces $deg(p)$ del polinomio de la siguiente manera: el punto se usa como valor inicial $z 0$ para la iteración de Newton $z n + 1 := z n - p (z n) / p' (z n)$ , produciendo 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 que 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 para el cual las iteraciones convergen hacia una raíz o ciclo dado (que no es un punto fijo), es un conjunto de Julia para la iteración. El conjunto complementario a la unión de todos estos, es el conjunto de Julia. Los conjuntos de Fatou tienen un límite común, a saber, 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 interesantes, 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 retícula rectangular

z_{mn}=z_{00}+m\,\Delta x+in\,\Delta y;\quad m=0,\ldots ,M-1;\quad 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 usa para llenar una cuadrícula de $M \times N$ píxeles, asignando a cada punto $(m, n)$ un color $f k (m, n)$ . Además 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 algunos $ε > 0$ pequeños previamente fijados.

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 elección particular con $a = 1$ corresponde al fractal de Newton.

Los puntos fijos de esta aplicación 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, la aplicación aún exhibe una estructura fractal en el sentido de 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 conjunto de Julia.

Fractal de Newton para tres raíces de grado 3 $p (z) = z 3 - 1$ , coloreado por 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 llegan a la 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 x 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) = sin z$ , coloreado por la raíz alcanzada, sombreado por el número de iteraciones requeridas
$Otro fractal de Newton para p(z) = sin z$

Otro fractal de Newton para $p (z) = sin z$
Fractal de Newton generalizado para $p (z) = z 3 - 1$ , $a = - 1 / 2$ . El color se eligió en función del 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) = z 6 + z 3 - 1$
$p (z) = sin z - 1$
$p (z) = sin z - 1$
$p (z) = cosh z - 1$
$p (z) = cosh z - 1$

$p (z) = z 3 - 1$ , a=1
$p (z) = z 3 - 1$ , a=2
$p (z) = z 4 - 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 7 - 1$ , a=2
$p (z) = z 8 - 1$ , a=1
$p (z) = z 10 - 1$ , a=1

Fractales compuestos de potencias y funciones trigonométricas.

$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 Nova

Artículo principal: Fractal Nova

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 constante $c$ sobre la imagen e inicializa $z 0$ a las coordenadas de píxeles. La variante "Mandelbrot" del fractal Nova inicializa $c$ en las coordenadas de 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$ .

Generación

Para generar el fractal de Newton, es necesario tener una función inicial y 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 Function (float2 z)
{
	return cpow(z, 3) - float2(1, 0); //cpow is an exponential function for complex numbers
}

//3*z^2
float2 Derivative (float2 z)
{
	return 3 * cmul(z, z); //cmul is a function that handles multiplication of complex numbers
}

Ahora es solo cuestión de aplicar el método de Newton usando las funciones dadas.

float2 roots[3] = //Roots (solutions) of the polynomial
{
	float2(1, 0), 
	float2(-.5, sqrt(3)/2), 
	float2(-.5, -sqrt(3)/2)
};
	
color colors[3] =  //Assign a color for each root
{
    red,
    green,
    blue
}

Para cada pixel (x, y) como objetivo, hacer que:
{
	zx = scaled x coordinate of pixel (scaled to lie in the Mandelbrot X scale (-2.5, 1))
    zy = scaled y coordinate of pixel (scaled to lie in the Mandelbrot Y scale (-1, 1))

    float2 z = float2(zx, zy); //z is originally set to the pixel coordinates

	for (int iteration = 0;
	     iteration < maxIteration;
	     iteration++;)
	{
		z -= cdiv(Function(z), Derivative(z)); //cdiv is a function for dividing complex numbers

        float tolerance = 0.000001;
        
		for (int i = 0; i < roots.Length; i++)
		{
		    float2 difference = z - roots[i];
		    
			//If the current iteration is close enough to a root, color the pixel.
			if (abs(difference.x) < tolerance && abs (difference.y) < tolerance)
			{
				return colors[i]; //Return the color corresponding to the root
			}
		}
		
    }
    
    return black; //If no solution is found
}

Véase también

Referencias

↑ Simon Tatham. «Fractals derived from Newton–Raphson».
↑ Damien M. Jones. «class Standard_NovaMandel (Ultra Fractal formula reference)».
↑ Damien M. Jones. «dmj's nova fractals 1995-6». Archivado desde el original el 20 de agosto de 2017. Consultado el 4 de enero de 2022.
↑ Michael Condron. «Relaxed Newton's Method and the Nova Fractal».
↑ Frederik Slijkerman. «Ultra Fractal Manual: Nova (Julia, Mandelbrot)».

Lecturas relacionadas

J. H. Hubbard, D. Schleicher, S. Sutherland: "Cómo encontrar todas las raíces de polinomios complejos por el método de Newton", Inventiones Mathematicae vol. 146 (2001) - con una discusión de la estructura global de los fractales de Newton
Sobre el número de iteraciones para el método de Newton por Dierk Schleicher 21 de julio de 2000
El método de Newton como sistema dinámico por Johannes Rueckert

Enlaces externos

Wikimedia Commons alberga una categoría multimedia sobre Fractal de Newton.

Datos: Q1151001
Multimedia: Newton fractals / Q1151001

[1] Simon Tatham. «Fractals derived from Newton–Raphson».

[2] Damien M. Jones. «class Standard_NovaMandel (Ultra Fractal formula reference)».

[3] Damien M. Jones. «dmj's nova fractals 1995-6». Archivado desde el original el 20 de agosto de 2017. Consultado el 4 de enero de 2022.

[4] Michael Condron. «Relaxed Newton's Method and the Nova Fractal».

[5] Frederik Slijkerman. «Ultra Fractal Manual: Nova (Julia, Mandelbrot)».

[1]

[2]

[3]

[4]

[5]