LECCION 23 | ~35 min

Diagramas de bifurcacion

El mapa completo de la dinamica

¿Que muestra un diagrama de bifurcacion?

Un diagrama de bifurcacion es quizas la herramienta visual mas poderosa de toda la teoria del caos. En una sola imagen comprime toda la informacion sobre como cambia el comportamiento de un sistema dinamico cuando variamos un parametro de control. Es, literalmente, el mapa completo de todas las dinamicas posibles.

La idea fundamental es elegante. Tomamos un sistema que depende de un parametro (llamemoslo $\rho$ para el Lorenz, o $r$ para el mapa logistico). Para cada valor del parametro, simulamos el sistema hasta que alcanza su comportamiento estacionario (eliminando el transitorio), y registramos los estados que visita a largo plazo. Graficamos estos estados en el eje vertical, y el parametro en el horizontal.

El resultado revela estructura donde parecia no haberla. Para valores del parametro donde el sistema es periodico, vemos puntos discretos: un punto para un punto fijo, dos para periodo 2, cuatro para periodo 4. Para valores donde el sistema es caotico, vemos nubes densas de puntos que, sin embargo, no llenan todo el espacio sino que forman estructuras con geometria fractal. Las transiciones entre estas regiones son las bifurcaciones.

Construccion algoritmica

La construccion de un diagrama de bifurcacion para un sistema de flujo continuo como el Lorenz sigue un procedimiento sistematico. No podemos iterar un mapa directamente, asi que usamos los maximos locales de una de las variables (tipicamente $z(t)$) como proxy del comportamiento estacionario:

Para cada $\rho$ en el rango deseado (por ejemplo, $\rho \in [0, 200]$), elegir una condicion inicial (por ejemplo, $[1, 1, 1]$).
Integrar el sistema de Lorenz durante un periodo transitorio (al menos 500 unidades de tiempo) para que la trayectoria se asiente en el atractor.
Continuar la integracion y detectar los maximos locales de $z(t)$: puntos donde $z$ crece, alcanza un pico, y luego decrece.
Registrar los valores de estos maximos (tipicamente 100-300 maximos).
Graficar un punto por cada $(\rho, z_{\max})$ encontrado.
Avanzar al siguiente valor de $\rho$ y repetir.

La eleccion de maximos locales no es arbitraria. Los maximos de $z$ corresponden a los puntos donde la trayectoria cruza una seccion de Poincare particular. Si la orbita es periodica de periodo $p$, habra exactamente $p$ maximos distintos (o un submultiplo). Si la orbita es caotica, habra infinitos maximos distintos que llenan una region con estructura fractal.

Los puntos densos indican caos, los puntos aislados indican periodicidad. Las transiciones entre estas regiones son las bifurcaciones, y su estudio revela la estructura universal del camino hacia el caos.

Tipos de bifurcacion

Cada transicion entre comportamientos dinamicos tiene un nombre y una estructura matematica precisa. Los tipos principales de bifurcacion son:

Bifurcacion silla-nodo (fold)

Dos puntos fijos (uno estable, uno inestable) colisionan y se aniquilan. La condicion es que un multiplicador de Floquet cruce el valor 1. Aparece como la aparicion o desaparicion subita de una rama en el diagrama.

Bifurcacion pitchfork

Un punto fijo se divide en tres: uno inestable y dos estables (supercritica) o uno estable y dos inestables (subcritica). Ocurre en sistemas con simetria, como el Lorenz en $\rho = 1$.

Bifurcacion de Hopf

Un punto fijo pierde estabilidad y nace un ciclo limite. La condicion es que un par de autovalores complejos cruce el eje imaginario. Aparece como la transicion de un punto a un anillo en el diagrama.

Duplicacion de periodo

Una orbita periodica se vuelve inestable y nace una orbita de periodo doble. Un multiplicador de Floquet cruza el valor -1. Es la ruta clasica al caos: periodo 1 $\to$ 2 $\to$ 4 $\to$ ... $\to$ caos.

La condicion matematica general para una bifurcacion es que al menos un multiplicador de Floquet (autovalor de la aplicacion de Poincare) cruce el circulo unidad en el plano complejo. El punto por donde cruza determina el tipo de bifurcacion:

$$|\mu_i(\rho_c)| = 1 \quad \Rightarrow \quad \text{bifurcacion en } \rho = \rho_c$$

Lorenz vs logistico

El diagrama de bifurcacion del sistema de Lorenz comparte muchas caracteristicas con el del mapa logistico, pero es considerablemente mas complejo debido a la naturaleza tridimensional del flujo. Ambos muestran cascadas de duplicacion de periodo, ventanas de periodicidad dentro del caos, y estructura fractal. Pero el Lorenz tiene peculiaridades propias.

Rango de ρ	Comportamiento	En el diagrama
0 < ρ < 1	Punto fijo en el origen	Sin maximos (z decae a 0)
1 < ρ < 24.74	Puntos fijos C+, C- estables	Un unico maximo por rama
ρ ≈ 24.74	Inicio del caos sostenido	Nube densa de puntos
ρ ≈ 99.5	Ventana periodica (knot)	Puntos discretos en la nube
ρ ≈ 148	Otra ventana periodica	Periodo 3 visible
ρ > 166	Caos robusto	Nube sin ventanas apreciables

Las ventanas de periodicidad en el Lorenz son analogas a las del mapa logistico. Dentro de cada ventana se repite la misma cascada de duplicaciones que en la region principal. Es una manifestacion de la autosimilitud del diagrama: la estructura fractal se repite a todas las escalas. El celebre teorema de Sharkovskii garantiza que si existe una orbita de periodo 3, existen orbitas de todos los periodos.

"El diagrama de bifurcacion es la radiografia definitiva de un sistema dinamico. Una sola imagen condensa un infinito de comportamientos posibles, desde la calma del punto fijo hasta la tormenta del caos, y revela la estructura fractal que los conecta."

Los diagramas de bifurcacion no son herramientas exclusivamente teoricas. Se usan en ingenieria para identificar regiones de operacion estable de motores y reactores quimicos, en cardiologia para detectar transiciones a arritmias, en ecologia para predecir colapsos de poblaciones, y en electronica para disenar circuitos osciladores. Donde hay un parametro que puede variar, hay un diagrama de bifurcacion esperando revelar la estructura oculta.

Ejercicios

Lectura del diagrama
En el LAB de abajo, identifica visualmente tres regiones distintas: (a) un rango de $\rho$ donde el sistema tiene un punto fijo estable, (b) un rango donde hay caos, y (c) una ventana de periodicidad dentro del caos. Registra los valores aproximados de $\rho$ para cada region.
Identificacion de bifurcaciones
Haz clic en $\rho \approx 24.74$ y luego en $\rho \approx 25$. ¿El atractor cambia gradualmente o hay una transicion abrupta? ¿Que tipo de bifurcacion crees que es? Compara con la tabla de tipos de bifurcacion.
Ventanas de periodicidad
Busca la ventana de periodicidad cerca de $\rho \approx 100$. Haz clic alli y observa el atractor en la esquina. ¿Cuantas orbitas distintas ves? ¿Que periodo tiene la orbita? Ahora busca la ventana de periodo 3 y repite el analisis.

LAB: Bifurcacion interactiva

Diagrama de bifurcacion del sistema de Lorenz, barriendo $\rho$ de 0 a 200. Se muestran los maximos locales de $z(t)$ para cada valor de $\rho$. Los puntos se colorean por densidad (mas brillante = mas visitado). Haz clic en cualquier punto del diagrama para ver el atractor correspondiente en la esquina inferior derecha (proyeccion x-z). El calculo se realiza al cargar la pagina.

bifurcacion-lorenz.js

Calculando...