Articulo de referencia

Julia estableció

Un mapa del comportamiento de los puntos iniciales z {\displaystyle z} en el plano complejo, mediante la iteración del mapa cuadrático F ( z ) = z 2 + do {\displaystyle f(z)=z^{...

Un mapa del comportamiento de los puntos inicialesz{\displaystyle z}en el plano complejo, mediante la iteración del mapa cuadráticoF(z)=z2+do{\displaystyle f(z)=z^{2}+c}donde c es aproximadamente -0,5125 + 0,5213i. Las regiones verde y azul son el conjunto de Fatou, y el conjunto de Julia consta de los puntos amarillos centrales.
Amplíe la vista de un conjunto de Julia en el plano z de valores complejos con la función polinómica de segundo grado de valores complejos.pag(z)=z2+do{\displaystyle p(z)=z^{2}+c}y los parámetros c re = c im = -0,5251993
Cortes tridimensionales a través del conjunto de Julia (cuatridimensional) de una función sobre los cuaterniones.

En dinámica compleja , el conjunto de Julia y el conjunto de Fatou son dos conjuntos complementarios (Julia "encajes" y Fatou "polvo") definidos a partir de una función . De manera informal, el conjunto de Fatou de la función consta de valores que se comportan de forma similar ante iteraciones repetidas de la función, mientras que el conjunto de Julia consta de valores tales que una perturbación arbitrariamente pequeña puede provocar cambios drásticos en la secuencia de valores iterados de la función. Por lo tanto, el comportamiento de la función en el conjunto de Fatou es "regular", mientras que en el conjunto de Julia su comportamiento es " caótico ".

El conjunto de Julia de una función f se denota comúnmenteJ(F),{\displaystyle \operatorname {J} (f),}y el conjunto de Fatou se denotaF(F).{\displaystyle \operatorname {F} (f).}[ a ] ​​Estos conjuntos reciben su nombre de los matemáticos francesesGaston Julia [ 1 ] yPierre Fatou [ 2 ] cuyo trabajo inició el estudio de ladinámica complejaa principios del siglo XX.

Definición formal

DejarF(z){\displaystyle f(z)}ser una función meromorfa no constante de la esfera de Riemann sobre sí misma. Tales funcionesF(z){\displaystyle f(z)}son precisamente las funciones racionales complejas no constantes , es decir,F(z)=pag(z)/q(z){\displaystyle f(z)=p(z)/q(z)}dóndepag(z){\displaystyle p(z)}yq(z){\displaystyle q(z)}son polinomios complejos . Supongamos que p y q no tienen raíces comunes y que al menos uno tiene grado mayor que 1. Entonces existe un número finito de conjuntos abiertos.F1,...,Fr{\displaystyle F_{1},...,F_{r}}que quedan invariables porF(z){\displaystyle f(z)}y son tales que:

  1. La unión de los conjuntosFi{\displaystyle F_{i}}es denso en el plano y
  2. F(z){\displaystyle f(z)}se comporta de manera regular e igual en cada uno de los conjuntosFi{\displaystyle F_{i}}.

La última afirmación significa que los extremos de las secuencias de iteraciones generadas por los puntos deFi{\displaystyle F_{i}}son o bien exactamente el mismo conjunto, que entonces es un ciclo finito, o bien son ciclos finitos de conjuntos de forma circular o anular que se encuentran concéntricamente. En el primer caso el ciclo es atractor , en el segundo caso es neutral .

Estos conjuntosFi{\displaystyle F_{i}}son los dominios Fatou deF(z){\displaystyle f(z)}y su unión es el conjunto FatouF(F){\displaystyle \operatorname {F} (f)}deF(z){\displaystyle f(z)}Cada uno de los dominios de Fatou contiene al menos un punto crítico deF(z){\displaystyle f(z)}, es decir, un punto (finito) z que satisfaceF(z)=0{\displaystyle f'(z)=0}, oF(z)={\displaystyle f(z)=\infty }si el grado del numeradorpag(z){\displaystyle p(z)}es al menos dos veces mayor que el grado del denominadorq(z){\displaystyle q(z)}, o siF(z)=1/gramo(z)+do{\displaystyle f(z)=1/g(z)+c}para algún c y una función racionalgramo(z){\displaystyle g(z)}satisfaciendo esta condición.

El complemento deF(F){\displaystyle \operatorname {F} (f)}¿Es el conjunto de Julia?J(F){\displaystyle \operatorname {J} (f)}deF(z){\displaystyle f(z)}Si todos los puntos críticos son preperiódicos, es decir, no son periódicos pero finalmente caen en un ciclo periódico, entoncesJ(F){\displaystyle \operatorname {J} (f)}es toda la esfera. De lo contrario,J(F){\displaystyle \operatorname {J} (f)}es un conjunto no denso en ninguna parte (carece de puntos interiores) y un conjunto no numerable (de la misma cardinalidad que los números reales). ComoF(F){\displaystyle \operatorname {F} (f)},J(F){\displaystyle \operatorname {J} (f)}queda invariante porF(z){\displaystyle f(z)}y en este conjunto la iteración es repulsiva, lo que significa que|F(z)F(w)|>|zw|{\displaystyle |f(z)-f(w)|>|zw|}para todo w en un vecindario de z (dentro deJ(F){\displaystyle \operatorname {J} (f)}). Esto significa queF(z){\displaystyle f(z)}El comportamiento es caótico en el conjunto de Julia. Si bien existen puntos en el conjunto de Julia cuya secuencia de iteraciones es finita, solo hay un número contable de dichos puntos (y constituyen una parte infinitesimal del conjunto de Julia). Las secuencias generadas por puntos fuera de este conjunto se comportan de forma caótica, un fenómeno denominado caos determinista .

Se han realizado numerosas investigaciones sobre el conjunto de Fatou y el conjunto de Julia de funciones racionales iteradas , conocidas como mapas racionales. Por ejemplo, se sabe que el conjunto de Fatou de un mapa racional tiene 0, 1, 2 o un número infinito de componentes . [ 3 ] Cada componente del conjunto de Fatou de un mapa racional puede clasificarse en una de cuatro clases diferentes . [ 4 ]

Descripciones equivalentes del conjunto Julia

  • J(F){\displaystyle \operatorname {J} (f)}es el conjunto cerrado más pequeño que contiene al menos tres puntos y que es completamente invariante bajo f .
  • J(F){\displaystyle \operatorname {J} (f)}es el cierre del conjunto de puntos periódicos repulsivos .
  • En todos los casos, salvo como máximo en dos puntos.zincógnita,{\displaystyle \;z\in X\;,}El conjunto de Julia es el conjunto de puntos límite de la órbita completa hacia atrás.norteFnorte(z).{\displaystyle \bigcup _{n}f^{-n}(z).}(Esto sugiere un algoritmo sencillo para representar gráficamente conjuntos de Julia; véase más abajo).
  • Si f es una función entera , entoncesJ(F){\displaystyle \operatorname {J} (f)}es el límite del conjunto de puntos que convergen al infinito bajo iteración.
  • Si f es un polinomio, entoncesJ(F){\displaystyle \operatorname {J} (f)}es el límite del conjunto de Julia lleno ; es decir, aquellos puntos cuyas órbitas bajo iteraciones de f permanecen acotadas.

Propiedades del conjunto Julia y del conjunto Fatou

El conjunto de Julia y el conjunto de Fatou de f son ambos completamente invariantes bajo iteraciones de la función holomorfa f : [ 5 ]

F1(J(F))=F(J(F))=J(F),{\displaystyle f^{-1}(\operatorname {J} (f))=f(\operatorname {J} (f))=\operatorname {J} (f),}
F1(F(F))=F(F(F))=F(F).{\displaystyle f^{-1}(\operatorname {F} (f))=f(\operatorname {F} (f))=\operatorname {F} (f).}

Ejemplos

ParaF(z)=z2{\displaystyle f(z)=z^{2}}El conjunto de Julia es el círculo unitario y en este la iteración viene dada por la duplicación de ángulos (una operación que es caótica en los puntos cuyo argumento no es una fracción racional de2π{\displaystyle 2\pi }). Hay dos dominios de Fatou: el interior y el exterior del círculo, con iteración hacia 0 e ∞, respectivamente.

Paragramo(z)=z22{\displaystyle g(z)=z^{2}-2}El conjunto de Julia es el segmento de línea entre −2 y 2. Hay un dominio de Fatou : los puntos que no están en el segmento de línea iteran hacia ∞. (Aparte de un desplazamiento y escalado del dominio, esta iteración es equivalente aincógnita4(incógnita12)2{\displaystyle x\to 4(x-{\tfrac {1}{2}})^{2}}en el intervalo unitario, que se usa comúnmente como ejemplo de sistema caótico.)

Las funciones f y g son de la formaz2+do{\displaystyle z^{2}+c}donde c es un número complejo. En esta iteración, el conjunto de Julia no es, en general, una curva simple, sino un fractal, y para algunos valores de c puede adoptar formas sorprendentes. Véanse las imágenes a continuación.

Conjunto de Julia (en blanco) para la función racional asociada al método de Newton para f  : zz 3 −1. Coloreado del conjunto de Fatou en tonos rojo, verde y azul según los tres atractores (las tres raíces de f ).

Para algunas funciones f ( z ) podemos decir de antemano que el conjunto de Julia es un fractal y no una curva simple. Esto se debe al siguiente resultado sobre las iteraciones de una función racional:

Teorema : Cada uno de los dominios de Fatou tiene la misma frontera, que en consecuencia es el conjunto de Julia.

Esto significa que cada punto del conjunto de Julia es un punto de acumulación para cada uno de los dominios de Fatou. Por lo tanto, si hay más de dos dominios de Fatou, cada punto del conjunto de Julia debe tener puntos de más de dos conjuntos abiertos diferentes infinitamente cerrados, y esto significa que el conjunto de Julia no puede ser una curva simple. Este fenómeno ocurre, por ejemplo, cuando f ( z ) es la iteración de Newton para resolver la ecuación.PAG(z):=znorte1=0 : norte>2{\displaystyle \;P(z):=z^{n}-1=0~:~n>2\;}:

F(z)=zPAG(z)PAG(z)=1+(norte1)znortenorteznorte1 .{\displaystyle f(z)=z-{\frac {P(z)}{P'(z)}}={\frac {\;1+(n-1)z^{n}\;}{nz^{n-1}}}~.}

La imagen de la derecha muestra el caso n = 3.

Polinomios cuadráticos

Un sistema dinámico complejo muy popular viene dado por la familia de polinomios cuadráticos complejos , un caso especial de aplicaciones racionales . Dichos polinomios cuadráticos pueden expresarse como

Fdo(z)=z2+do ,{\displaystyle f_{c}(z)=z^{2}+c~,}

donde c es un parámetro complejo. Fijar algunosR>0{\displaystyle R>0}lo suficientemente grande como paraR2R|do|.{\displaystyle R^{2}-R\geq |c|.}(Por ejemplo, si c está en el conjunto de Mandelbrot , entonces|do|2,{\displaystyle |c|\leq 2,}así que podemos simplemente dejarR=2 .{\displaystyle R=2~.}) Entonces, el conjunto de Julia lleno para este sistema es el subconjunto del plano complejo dado por

K(Fdo)={zdo:nortenorte,|Fdonorte(z)|R} ,{\displaystyle K(f_{c})=\left\{z\in \mathbb {C} :\forall n\in \mathbb {N} ,|f_{c}^{n}(z)|\leq R\right\}~,}

dóndeFdonorte(z){\displaystyle f_{c}^{n}(z)}es la enésima iteración deFdo(z).{\displaystyle f_{c}(z).}El conjunto JuliaJ(Fdo){\displaystyle J(f_{c})}de esta función es el límite deK(Fdo){\displaystyle K(f_{c})}.

El plano de parámetros de los polinomios cuadráticos  – es decir, el plano de los posibles valores de c  – da lugar al famoso conjunto de Mandelbrot . De hecho, el conjunto de Mandelbrot se define como el conjunto de todos los c tales queJ(Fdo){\displaystyle J(f_{c})}está conectado . Para parámetros fuera del conjunto de Mandelbrot, el conjunto de Julia es un espacio de Cantor : en este caso, a veces se le denomina polvo de Fatou .

En muchos casos, el conjunto de Julia de c se parece al conjunto de Mandelbrot en entornos suficientemente pequeños de c . Esto es cierto, en particular, para los llamados parámetros de Misiurewicz , es decir, parámetros c para los cuales el punto crítico es preperiódico. Por ejemplo:

  • En c = i , la punta más corta del antepié, el conjunto Julia parece un rayo ramificado.
  • En c = −2, la punta de la larga cola puntiaguda, el conjunto de Julia es un segmento de línea recta.

En otras palabras, los conjuntos de JuliaJ(Fdo){\displaystyle J(f_{c})}son localmente similares alrededor de los puntos de Misiurewicz . [ 6 ]

Generalizaciones

La definición de conjuntos de Julia y Fatou se traslada fácilmente al caso de ciertos mapas cuya imagen contiene su dominio; en particular, las funciones meromorfas trascendentales y los mapas de tipo finito de Adam Epstein .

Los conjuntos de Julia también se definen comúnmente en el estudio de la dinámica en varias variables complejas.

Pseudocódigo

Las implementaciones en pseudocódigo que se muestran a continuación codifican las funciones para cada fractal. Considere implementar operaciones con números complejos para lograr un código más dinámico y reutilizable.

Pseudocódigo para conjuntos normales de Julia

F(z)=z2+do{\displaystyle f(z)=z^{2}+c}
R = radio de escape # elige R > 0 tal que R**2 - R >= sqrt(cx**2 + cy**2)para cada píxel ( x , y ) en la pantalla , hacer : { zx = coordenada x escalada del píxel ; # (escala para que esté entre -R y R) # zx representa la parte real de z. zy = coordenada y escalada del píxel ; # (escala para que esté entre -R y R) # zy representa la parte imaginaria de z.iteración = 0 ; iteración_máxima = 1000 ; mientras ( zx * zx + zy * zy < R ** 2 Y iteración < iteración_máxima ) { xtemp = zx * zx - zy * zy ; zy = 2 * zx * zy + cy ; zx = xtemp + cx ; iteración = iteración + 1 ; } si ( iteración == iteración_máxima ) devolver negro ; de lo contrario devolver iteración ; }

Pseudocódigo para conjuntos multi-Julia

F(z)=znorte+do{\displaystyle f(z)=z^{n}+c}
R = radio de escape # elige R > 0 tal que R**n - R >= sqrt(cx**2 + cy**2)para cada píxel ( x , y ) en la pantalla , hacer : { zx = coordenada x escalada del píxel ; # (escala para que esté entre -R y R) zy = coordenada y escalada del píxel ; # (escala para que esté entre -R y R) iteración = 0 ; iteración_máxima = 1001 ; mientras ( zx * zx + zy * zy < R ** 2 Y iteración < iteración_máxima ) { xtmp = ( zx * zx + zy * zy ) ^ ( n / 2 ) * cos ( n * atan2 ( zy , zx )) + cx ; zy = ( zx * zx + zy * zy ) ^ ( n / 2 ) * sin ( n * atan2 ( zy , zx )) + cy ; zx = xtmp ; iteración = iteración + 1 ; } si ( iteración == max_iteration ) devuelve negro ; de lo contrario, devuelve iteración ; }

Otra opción recomendada es reducir las bandas de color entre iteraciones utilizando una fórmula de renormalización para la iteración. [ 7 ]

Dicha fórmula se da como,

metro=k+1registroregistro|zk|registronorte{\displaystyle mu=k+1-{\frac {\log {\log {|z_{k}|}}}{\log {n}}}}
Fdo,norte(z)=znorte+términos de menor potencia+do{\displaystyle \forall f_{c,n}(z)=z^{n}+{\text{lower power terms}}+c}

dóndek{\displaystyle k}es la iteración de escape, limitada por algúnK{\displaystyle K}de tal manera que0k<K{\displaystyle 0\leq k<K}yKnorte{\displaystyle K\in \mathbb {N} }, y|zk|{\displaystyle |z_{k}|}es la magnitud de la última iteración antes de escapar.

Esto se puede implementar, de forma muy sencilla, de la siguiente manera:

# Simplemente reemplace las últimas 4 líneas de código del último ejemplo con estas líneas de código:if ( iteration == max_iteration ) return black ; else abs_z = zx * zx + zy * zy ; return iteration + 1 - log ( log ( abs_z )) / log ( n );

La diferencia se muestra a continuación con un conjunto de Julia definido comoFdo,2(z){\displaystyle f_{c,2}(z)}dóndedo=0,8350,321i{\displaystyle c=-0.835-0.321i}.

La función potencial y el número real de iteraciones

El set de Julia paraF(z)=z2{\displaystyle f(z)=z^{2}}es el círculo unitario, y en el dominio exterior de Fatou, la función potencial φ ( z ) se define por φ ( z ) = log| z |. Las líneas equipotenciales para esta función son círculos concéntricos. Como|F(z)|=|z|2{\displaystyle |f(z)|=|z|^{2}}tenemos

φ(z)=límitekregistro|zk|2k,{\displaystyle \varphi (z)=\lim _{k\to \infty }{\frac {\log |z_{k}|}{2^{k}}},}

dóndezk{\displaystyle z_{k}}es la secuencia de iteración generada por z . Para la iteración más generalF(z)=z2+do{\displaystyle f(z)=z^{2}+c}Se ha demostrado que si el conjunto de Julia es conexo (es decir, si c pertenece al conjunto de Mandelbrot (usual)), entonces existe una aplicación biholomorfa ψ entre el dominio exterior de Fatou y el exterior del círculo unitario tal que|ψ(F(z))|=|ψ(z)|2{\displaystyle |\psi (f(z))|=|\psi (z)|^{2}}. [ 8 ] Esto significa que la función potencial en el dominio de Fatou exterior definido por esta correspondencia viene dada por:

φ(z)=límitekregistro|zk|2k.{\displaystyle \varphi (z)=\lim _{k\to \infty }{\frac {\log |z_{k}|}{2^{k}}}.}

Esta fórmula también tiene sentido si el conjunto de Julia no es conexo, de modo que para todo c podemos definir la función potencial en el dominio de Fatou que contiene ∞ mediante esta fórmula. Para una función racional general f ( z ) tal que ∞ es un punto crítico y un punto fijo, es decir, tal que el grado m del numerador es al menos dos unidades mayor que el grado n del denominador, definimos la función potencial en el dominio de Fatou que contiene ∞ mediante:

φ(z)=límitekregistro|zk|dk,{\displaystyle \varphi (z)=\lim _{k\to \infty }{\frac {\log |z_{k}|}{d^{k}}},}

donde d = mn es el grado de la función racional. [ 9 ]

Si N es un número muy grande (por ejemplo, 10 100 ), y si k es el primer número de iteración tal que|zk|>norte{\displaystyle |z_{k}|>N}, tenemos eso

registro|zk|dk=registro(norte)dν(z),{\displaystyle {\frac {\log |z_{k}|}{d^{k}}}={\frac {\log(N)}{d^{\nu (z)}}},}

para algún número realν(z){\displaystyle \nu (z)}, que debe considerarse como el número de iteración real , y tenemos que:

ν(z)=kregistro(registro|zk|/registro(norte))registro(d),{\displaystyle \nu (z)=k-{\frac {\log(\log |z_{k}|/\log(N))}{\log(d)}},}

donde el último número está en el intervalo [0,  1).

Para la iteración hacia un ciclo atractor finito de orden r , tenemos que siz{\displaystyle z^{*}}es un punto del ciclo, entoncesF(F(...F(z)))=z{\displaystyle f(f(...f(z^{*})))=z^{*}}(la composición r -ésima), y el número

α=1|(d(F(F(F(z))))/dz)z=z|(>1){\displaystyle \alpha ={\frac {1}{\left|(d(f(f(\cdots f(z))))/dz)_{z=z^{*}}\right|}}\qquad (>1)}

es la atracción del ciclo. Si w es un punto muy cercanoz{\displaystyle z^{*}}y w es w iterado r veces, tenemos que

α=límitek|wz||wz|.{\displaystyle \alpha =\lim _{k\to \infty }{\frac {|w-z^{*}|}{|w'-z^{*}|}}.}

Por lo tanto, el número|zkrz|αk{\displaystyle |z_{kr}-z^{*}|\alpha ^{k}}es casi independiente de k . Definimos la función potencial en el dominio de Fatou mediante:

φ(z)=límitek1(|zkrz|αk).{\displaystyle \varphi (z)=\lim _{k\to \infty }{\frac {1}{(|z_{kr}-z^{*}|\alpha ^{k})}}.}

Si ε es un número muy pequeño y k es el primer número de iteración tal que|zkz|<ϵ{\displaystyle |z_{k}-z^{*}|<\epsilon }, tenemos eso

φ(z)=1(εαν(z)){\displaystyle \varphi (z)={\frac {1}{(\varepsilon \alpha ^{\nu (z)})}}}

para algún número realν(z){\displaystyle \nu (z)}, que debe considerarse como el número de iteración real, y tenemos que:

ν(z)=kregistro(ε/|zkz|)registro(α).{\displaystyle \nu (z)=k-{\frac {\log(\varepsilon /|z_{k}-z^{*}|)}{\log(\alpha )}}.}

Si la atracción es ∞, lo que significa que el ciclo es superatractivo , lo que nuevamente significa que uno de los puntos del ciclo es un punto crítico, debemos reemplazar α por

α=límitekregistro|wz|registro|wz|,{\displaystyle \alpha =\lim _{k\to \infty }{\frac {\log |w'-z^{*}|}{\log |w-z^{*}|}},}

donde w es w iterado r veces y la fórmula para φ ( z ) es:

φ(z)=límitekregistro(1/|zkrz|)αk.{\displaystyle \varphi (z)=\lim _{k\to \infty }{\frac {\log(1/|z_{kr}-z^{*}|)}{\alpha ^{k}}}.}

Y ahora el número de iteración real viene dado por:

ν(z)=kregistro(registro|zkz|/registro(ε))registro(α).{\displaystyle \nu (z)=k-{\frac {\log(\log |z_{k}-z^{*}|/\log(\varepsilon ))}{\log(\alpha )}}.}

Para la coloración debemos tener una escala cíclica de colores (construida matemáticamente, por ejemplo) que contenga H colores numerados del 0 al H −1 ( H = 500, por ejemplo). Multiplicamos el número realν(z){\displaystyle \nu (z)}mediante un número real fijo que determina la densidad de los colores en la imagen, y tomar la parte entera de este número módulo H.

La definición de la función potencial y nuestra forma de colorear presuponen que el ciclo es atractor, es decir, no neutro. Si el ciclo es neutro, no podemos colorear el dominio de Fatou de forma natural. Dado que el final de la iteración es un movimiento de rotación, podemos, por ejemplo, colorear según la distancia mínima al ciclo que queda fijo tras la iteración.

Líneas de campo

Las líneas equipotenciales para la iteración hacia el infinito
Líneas de campo para una iteración de la forma(1z3/6)(zz2/2)2+do{\displaystyle {\frac {(1-z^{3}/6)}{(z-z^{2}/2)^{2}}}+c}

En cada dominio de Fatou (que no es neutral) hay dos sistemas de líneas ortogonales entre sí: las líneas equipotenciales (para la función potencial o el número de iteración real) y las líneas de campo .

Si coloreamos el dominio de Fatou según el número de iteración (y no el número de iteración real)ν(z){\displaystyle \nu (z)}, como se define en la sección anterior), las bandas de iteración muestran el curso de las líneas equipotenciales. Si la iteración es hacia ∞ (como es el caso con el dominio de Fatou exterior para la iteración usualz2+do{\displaystyle z^{2}+c}), podemos mostrar fácilmente el curso de las líneas de campo, es decir, alterando el color según si el último punto en la secuencia de iteración está por encima o por debajo del eje x (primera imagen), pero en este caso (más precisamente: cuando el dominio de Fatou es superatractor) no podemos dibujar las líneas de campo de manera coherente, al menos no por el método que describimos aquí. En este caso, una línea de campo también se llama rayo externo .

Sea z un punto en el dominio de Fatou atractor. Si iteramos z un gran número de veces, el término de la secuencia de iteración es un ciclo finito C , y el dominio de Fatou es (por definición) el conjunto de puntos cuya secuencia de iteración converge hacia C. Las líneas de campo parten de los puntos de C y de los (infinitos) puntos que iteran hacia un punto de C. Y terminan en el conjunto de Julia en puntos que no son caóticos (es decir, que generan un ciclo finito). Sea r el orden del ciclo C (su número de puntos) y seaz{\displaystyle z^{*}}ser un punto en C. TenemosF(F(F(z)))=z{\displaystyle f(f(\dots f(z^{*})))=z^{*}}(la composición r-ésima), y definimos el número complejo α por

α=(d(F(F(F(z))))/dz)z=z.{\displaystyle \alpha =(d(f(f(\dots f(z))))/dz)_{z=z^{*}}.}

Si los puntos de C sonzi,i=1,,r(z1=z){\displaystyle z_{i},i=1,\dots ,r(z_{1}=z^{*})}, α es el producto de los números rF(zi){\displaystyle f'(z_{i})}El número real 1/|α| es la atracción del ciclo, y nuestra suposición de que el ciclo no es ni neutro ni superatractor, significa que 1 < 1 / | α | < . El puntoz{\displaystyle z^{*}}es un punto fijo paraF(F(F(z))){\displaystyle f(f(\dots f(z)))}y cerca de este punto el mapaF(F(F(z))){\displaystyle f(f(\dots f(z)))}tiene (en relación con las líneas de campo) carácter de una rotación con el argumento β de α (es decir,α=|α|miβi{\displaystyle \alpha =|\alpha |e^{\beta i}}).

Para colorear el dominio de Fatou, hemos elegido un número pequeño ε y establecido las secuencias de iteración.zk(k=0,1,2,,z0=z){\displaystyle z_{k}(k=0,1,2,\dots ,z_{0}=z)}parar cuando|zkz|<ϵ{\displaystyle |z_{k}-z^{*}|<\epsilon }y coloreamos el punto z según el número k (o el número de iteración real, si preferimos un coloreado suave). Si elegimos una dirección desdez{\displaystyle z^{*}}dada por un ángulo θ , la línea de campo que sale dez{\displaystyle z^{*}}en esta dirección consiste en los puntos z tales que el argumento ψ del númerozkz{\displaystyle z_{k}-z^{*}}satisface la condición de que

ψkβ=θmodπ.{\displaystyle \psi -k\beta =\theta \mod \pi .\,}

Porque si pasamos una banda de iteración en la dirección de las líneas de campo (y lejos del ciclo), el número de iteración k aumenta en 1 y el número ψ aumenta en β, por lo tanto el númeroψkβmodπ{\displaystyle \psi -k\beta \mod \pi }es constante a lo largo de la línea de campo.

Imágenes en las líneas del campo para una iteración de la formaz2+do{\displaystyle z^{2}+c}

Una coloración de las líneas de campo del dominio de Fatou significa que coloreamos los espacios entre pares de líneas de campo: elegimos una serie de direcciones situadas regularmente que parten dez{\displaystyle z^{*}}En cada una de estas direcciones, elegimos dos direcciones alrededor de esta dirección. Como puede ocurrir que las dos líneas de campo de un par no terminen en el mismo punto del conjunto de Julia, nuestras líneas de campo coloreadas pueden ramificarse (infinitamente) en su camino hacia el conjunto de Julia. Podemos colorear en función de la distancia a la línea central de la línea de campo, y podemos combinar este coloreado con el coloreado habitual. Tales imágenes pueden ser muy decorativas (segunda imagen).

Una línea de campo coloreada (el dominio entre dos líneas de campo) se divide mediante bandas de iteración, y dicha parte se puede correlacionar biunívocamente con el cuadrado unitario: una coordenada es (calculada a partir de) la distancia desde una de las líneas de campo delimitadoras, y la otra es (calculada a partir de) la distancia desde el interior de las bandas de iteración delimitadoras (este número es la parte no entera del número real de iteraciones). Por lo tanto, podemos representar gráficamente las líneas de campo (tercera imagen).

Trazando el conjunto de Julia

Descomposición binaria del interior en caso de ángulo interno 0

Métodos  :

  • Método de estimación de distancias para el conjunto de Julia (DEM/J)
  • Método de iteración inversa (IIM)

Utilizando la iteración inversa (IIM)

Un gráfico de conjuntos de Julia, generado mediante IIM aleatorio
Un gráfico de conjuntos de Julia, generado usando MIIM.

Como se mencionó anteriormente, el conjunto de Julia se puede encontrar como el conjunto de puntos límite del conjunto de preimágenes de (esencialmente) cualquier punto dado. Entonces podemos intentar graficar el conjunto de Julia de una función dada de la siguiente manera. Comenzamos con cualquier punto z que sabemos que está en el conjunto de Julia, como un punto periódico repulsivo, y calculamos todas las preimágenes de z bajo alguna iteración alta.Fnorte{\displaystyle f^{n}}apagado .

Desafortunadamente, dado que el número de preimágenes iteradas crece exponencialmente, esto no es factible computacionalmente. Sin embargo, podemos ajustar este método, de manera similar al método del "juego aleatorio" para sistemas de funciones iteradas . Es decir, en cada paso, elegimos al azar una de las imágenes inversas de f .

Por ejemplo, para el polinomio cuadrático f c , la iteración hacia atrás se describe mediante

znorte1=znortedo.{\displaystyle z_{n-1}={\sqrt {z_{n}-c}}.}

En cada paso, se selecciona al azar una de las dos raíces cuadradas.

Cabe señalar que ciertas partes del conjunto de Julia son bastante difíciles de acceder con el algoritmo inverso de Julia. Por este motivo, es necesario modificar IIM/J (un ejemplo de ello es MIIM/J [ 10 ] ) o utilizar otros métodos para obtener mejores imágenes.

Utilizando DEM/J

Dado que un conjunto de Julia es infinitamente delgado, no podemos dibujarlo eficazmente mediante iteración hacia atrás a partir de los píxeles. Aparecerá fragmentado debido a la impracticabilidad de examinar infinitos puntos de inicio. Como el número de iteraciones cambia drásticamente cerca del conjunto de Julia, una solución parcial consiste en inferir el contorno del conjunto a partir de los contornos de color más cercanos, pero el conjunto tenderá a verse borroso.

Una mejor manera de dibujar el conjunto de Julia en blanco y negro es estimar la distancia de los píxeles (DEM) al conjunto y colorear cada píxel cuyo centro esté cerca del conjunto. La fórmula para la estimación de la distancia se deriva de la fórmula para la función potencial φ ( z ). Cuando las líneas equipotenciales para φ ( z ) están cerca, el número|φ(z)|{\displaystyle |\varphi '(z)|}es grande y, a la inversa, por lo tanto, las líneas equipotenciales para la funciónδ(z)=φ(z)/|φ(z)|{\displaystyle \delta (z)=\varphi (z)/|\varphi '(z)|}debería ser aproximadamente regular. Se ha demostrado que el valor hallado por esta fórmula (salvo un factor constante) converge hacia la distancia real para z que converge hacia el conjunto de Julia. [ 9 ]

Suponemos que f ( z ) es racional, es decir,F(z)=pag(z)/q(z){\displaystyle f(z)=p(z)/q(z)}donde p ( z ) y q ( z ) son polinomios complejos de grados m y n , respectivamente, y tenemos que encontrar la derivada de las expresiones anteriores para φ ( z ). Y como es solozk{\displaystyle z_{k}}que varía, debemos calcular la derivadazk{\displaystyle z'_{k}}dezk{\displaystyle z_{k}}con respecto a z . Pero comozk=F(F(F(z))){\displaystyle z_{k}=f(f(\cdots f(z)))}(la composición k -ésima),zk{\displaystyle z'_{k}}es el producto de los númerosF(zk){\displaystyle f'(z_{k})}y esta secuencia se puede calcular recursivamente mediantezk+1=F(zk)zk{\displaystyle z'_{k+1}=f'(z_{k})z'_{k}}, comenzando conz0=1{\displaystyle z'_{0}=1}( antes del cálculo de la siguiente iteración)zk+1=F(zk){\displaystyle z_{k+1}=f(z_{k})}).

Para la iteración hacia ∞ (más precisamente cuando mn + 2 , de modo que ∞ es un punto fijo superatractor), tenemos

|φ(z)|=límitek|zk||zk|dk,{\displaystyle |\varphi '(z)|=\lim _{k\to \infty }{\frac {|z'_{k}|}{|z_{k}|d^{k}}},}

( d = mn ) y, en consecuencia:

δ(z)=φ(z)/|φ(z)|=límitekregistro|zk||zk|/|zk|.{\displaystyle \delta (z)=\varphi (z)/|\varphi '(z)|=\lim _{k\to \infty }\log |z_{k}||z_{k}|/|z'_{k}|.\,}

Para iterar hacia un ciclo atractor finito (que no sea superatractor) que contenga el punto z{\displaystyle z^{*}}y teniendo orden r , tenemos

|φ(z)|=límitek|zkr|/(|zkrz|2αk),{\displaystyle |\varphi '(z)|=\lim _{k\to \infty }|z'_{kr}|/(|z_{kr}-z^{*}|^{2}\alpha ^{k}),\,}

y, en consecuencia:

δ(z)=φ(z)/|φ(z)|=límitek|zkrz|/|zkr|.{\displaystyle \delta (z)=\varphi (z)/|\varphi '(z)|=\lim _{k\to \infty }|z_{kr}-z^{*}|/|z'_{kr}|.\,}

Para un ciclo súper atractivo, la fórmula es:

δ(z)=límitekregistro|zkrz|2/|zkr|.{\displaystyle \delta (z)=\lim _{k\to \infty }\log |z_{kr}-z^{*}|^{2}/|z'_{kr}|.\,}

Calculamos este número cuando finaliza la iteración. Cabe destacar que la estimación de la distancia es independiente de la atracción del ciclo. Esto significa que tiene sentido para funciones trascendentales de grado infinito (por ejemplo, sin( z ) y tan( z )).

Además de trazar el contorno, se puede introducir la función de distancia como una tercera dimensión para crear un paisaje fractal sólido.

Véase también

Notas

  1. Respecto a la notación: Para otras ramas de las matemáticas la notaciónJ(F),{\displaystyle \operatorname {J} (f),}También puede representar la matriz jacobiana de una aplicación de valor real f entre variedades diferenciables .

Referencias

  1. ^ Gaston Julia (1918) "Mémoire sur l'iteration des fonctions rationnelles", Journal de Mathématiques Pures et Appliquées , vol. 8, páginas 47–245.
  2. Pierre Fatou (1917) "Sur les substitutions rationnelles", Comptes Rendus de l'Académie des Sciences de Paris , vol. 164, páginas 806–808 y vol. 165, páginas 992–995.
  3. Beardon, Iteración de funciones racionales , Teorema 5.6.2.
  4. Beardon, Iteración de funciones racionales , Teorema 7.1.1.
  5. Beardon, Iteración de funciones racionales , Teorema 3.2.4.
  6. Tan Lei , "Similitud entre el conjunto de Mandelbrot y los conjuntos de Julia" , Communications in Mathematical Physics 134 (1990), pp. 587 617.
  7. Vepstas, Linas. "Renormalizando el escape de Mandelbrot" . linas.org . Creative Commons . Consultado el 5 de noviembre de 2023 .
  8. ^ Douady, Adrián; Hubbard, John H. (1984). "Etude dynamique des polynômes complexes". Prépublications mathémathiques d'Orsay . 2 ; "[ op.cit. ]". Prépublications mathémathiques d'Orsay . 4 . 1985.
  9. 1 2 Peitgen, Heinz-Otto; Richter Peter (1986). La belleza de los fractales . Heidelberg: Springer-Verlag. ISBN 0-387-15851-0.
  10. Saupe, Dietmar (octubre de 1987). "Cálculo eficiente de conjuntos de Julia y su dimensión fractal" . Physica D: Nonlinear Phenomena . 28 (3): 358–370 . doi : 10.1016/0167-2789(87)90024-8 . Consultado el 5 de febrero de 2026 .

Bibliografía

  • Carleson, Lennart ; Gamelin, Theodore W. (1993). Dinámica compleja . Springer.
  • Douady, Adrián ; Hubbard, John H. (1984). "Etude dynamique des polynômes complexes". Prépublications mathémathiques d'Orsay . 2 ; "[ op.cit. ]". Prépublications mathémathiques d'Orsay . 4 . 1985.
  • Milnor, JW (2006) [1990]. Dinámica en una variable compleja . Anales de estudios matemáticos. Vol.  160 (Tercera  ed.). Princeton University Press;Publicado originalmente como "Preimpresión del IMS de Stony Brook" . Archivado del original el 24 de abril de 2006.disponible como Milnor, John W. (1990). "Dinámica en una variable compleja: conferencias introductorias". arXiv : math.DS/9201272 .
  • Bogomolny, Alexander . "Conjunto de Mandelbrot e indexación de conjuntos de Julia" . cut-the-knot . Plan de estudios de álgebra.
  • Demidov, Evgeny (2003). "La anatomía de los conjuntos de Mandelbrot y Julia" .
  • Beardon, Alan F. (1991). Iteración de funciones racionales . Springer. ISBN 0-387-95151-2.
  • "Julia set" , Enciclopedia de Matemáticas , EMS Press, 2001 [1994]
  • Weisstein, Eric W. "Julia Set" . MundoMatemático .
  • Bourke, Paul. "Conjunto fractal Julia (2D)" (sitio web personal).
  • Sawyer, Jamie (6 de abril de 2007). "Julia establece" (blog).
  • McGoodwin, Michael. "Joyas Julia: Una exploración de los conjuntos Julia" (sitio web personal).
  • Pringle, Lucy. "Círculos en los cultivos, Julia Set" (sitio web personal).
  • Greig, Josh. "Applet interactivo para Julia" . Archivado del original el 26 de marzo de 2012.
  • Joyce, David E. "Julia y el explorador del conjunto de Mandelbrot" (sitio web personal académico). Universidad Clark.
  • "Un programa sencillo para generar conjuntos de Julia" . liazardie.com . Archivado del original el 17 de marzo de 2011.– Windows, 370  kB
  • "Una colección de applets" . SourceForge .– uno de los applets puede renderizar conjuntos de Julia, a través de sistemas de funciones iteradas.
  • "Julia se encuentra con HTML5" . Google Labs. Archivado del original el 18 de febrero de 2011. Generador de fractales HTML5 para tu navegador.
  • "Julia" . r-project.org . Paquete GNU R. 25 de noviembre de 2014. Genera un conjunto de Julia o Mandelbrot en una región y resolución determinadas.
  • "Julia establece" . Archivado del original el 16 de junio de 2017. Consultado el 24 de agosto de 2011 .– Una explicación visual de los conjuntos de Julia.
  • "FractalTS" . github.io .– Mandelbrot, el barco en llamas y el generador de conjuntos de Julia correspondiente.
  • "Imágenes del set de Julia, renderizado en línea" . finengin.net . Archivado del original el 22 de junio de 2019. Consultado el 24 de enero de 2017 .
  • "Comprendiendo a Julia y los conjuntos de Mandelbrot" .- Una explicación visual.