Crisis, intermitencia, control y sincronización del caos

Cuando el atractor cambia de forma global

Una bifurcación local se detecta mediante multiplicadores o autovalores que cruzan una frontera de estabilidad. Una crisis es un cambio global y repentino de un atractor caótico causado por su colisión con un conjunto inestable, con frecuencia una órbita periódica y su variedad estable. La figura puede cambiar de tamaño, unir componentes o desaparecer, aunque los equilibrios visibles no sufran una bifurcación local en ese mismo parámetro.

Se distinguen tres mecanismos frecuentes:

• en una crisis interior, el atractor se expande al acceder a una región antes separada por una barrera inestable;
• en una crisis de fusión, dos o más componentes caóticas se comunican;
• en una crisis de frontera, el atractor choca con el borde de su cuenca y se destruye, dejando caos transitorio.

La clasificación exige seguir el conjunto inestable implicado. Un salto en el máximo de una serie o un hueco en un diagrama de bifurcación es una señal candidata, no la prueba del mecanismo. Se continúa una órbita periódica, se aproximan sus variedades, se barre el parámetro en ambas direcciones y se comprueba la persistencia al refinar tiempo y malla.

Ott desarrolla las crisis como colisiones entre atractores y conjuntos inestables y conecta su geometría con escalas estadísticas Ott, 2002 Edward Ott (2002) Chaos in Dynamical Systems 2.ª ed. · Cambridge University Press Ubicación consultada: cap. 8 . Esa conexión permite buscar algo más preciso que una imagen: tiempos medios entre estallidos o vidas medias que cambian con la distancia al parámetro crítico.

Intermitencia y episodios laminares

La intermitencia alterna intervalos de comportamiento casi regular, llamados fases laminares, con estallidos irregulares. No es simplemente “caos que a veces se calma”. Diferentes mecanismos producen leyes de escala distintas.

La intermitencia tipo I aparece cerca de una bifurcación silla-nodo de un mapa. Una forma normal para la región estrecha es

$$x_{n+1}=x_n+\varepsilon+x_n^2, \qquad \varepsilon>0.$$

Ya no existe el punto fijo que habría para el otro signo, pero cerca de $x=0$ el incremento es pequeño. La órbita atraviesa lentamente un canal y luego es reinyectada por la dinámica global. Aproximando la iteración por

$$\frac{dx}{dn}=\varepsilon+x^2,$$

el tiempo para cruzar de $-c$ a $c$ es

$$\ell(\varepsilon) \approx\int_{-c}^{c}\frac{dx}{\varepsilon+x^2} =\frac{2}{\sqrt\varepsilon} \arctan\frac{c}{\sqrt\varepsilon}.$$

Cuando $\varepsilon\ll c^2$,

$$\ell(\varepsilon)\sim\frac{\pi}{\sqrt\varepsilon}.$$

Esta es una predicción cuantitativa: el tiempo laminar medio crece con exponente $-1/2$ bajo hipótesis de reinyección apropiadas. En datos finitos se define antes el umbral que separa fase laminar y estallido, se mide la distribución completa de longitudes y se verifica la pendiente en más de una década de $\varepsilon$.

La intermitencia tipo II se asocia con la pérdida de estabilidad de un par complejo y la tipo III con un multiplicador que cruza $-1$. Sus leyes dependen de la forma normal y la reinyección. La intermitencia on–off aparece cuando una variedad invariante pierde estabilidad transversal: la distancia a la variedad pasa largos periodos pequeña y sufre estallidos. Etiquetar el tipo solo por la serie temporal puede ser ambiguo; se necesita el multiplicador o exponente transversal correspondiente.

Intermitencia inducida por una crisis

Después de una crisis interior, la órbita suele permanecer en la región del antiguo atractor y ocasionalmente visitar la nueva zona accesible. El tiempo medio entre visitas puede seguir

$$\langle\tau\rangle\sim |\mu-\mu_c|^{-\gamma},$$

con un exponente relacionado con los autovalores del conjunto inestable responsable y con la geometría de reinyección. Tras una crisis de frontera, la vida media del transitorio presenta una ley análoga al acercarse al umbral desde el lado sin atractor.

Alligood, Sauer y Yorke relacionan crisis, órbitas periódicas inestables y variedades mediante continuación geométrica Alligood, 1996 Kathleen T. Alligood, Tim D. Sauer y James A. Yorke (1996) Chaos: An Introduction to Dynamical Systems Springer Ubicación consultada: caps. 10–12 . En una auditoría numérica se estiman $\mu_c$ y $\gamma$ conjuntamente, se muestran pendientes locales y se repite cambiando moderadamente la definición de visita. Ajustar $\mu_c$ hasta obtener la pendiente deseada es circular.

Controlar una órbita periódica inestable

Un atractor caótico contiene muchas órbitas periódicas inestables. El método OGY aprovecha una de ellas: espera hasta que la trayectoria entra en su vecindad y aplica perturbaciones paramétricas pequeñas para estabilizarla. No elimina la no linealidad ni impone una referencia externa arbitraria; modifica localmente el multiplicador del mapa de retorno.

Sea

$$x_{n+1}=F(x_n,p)$$

y sea $(x^*,p_0)$ un punto fijo inestable. Con $\delta x_n=x_n-x^*$ y $\delta p_n=p_n-p_0$,

$$\delta x_{n+1}=A\delta x_n+B\delta p_n+O(\delta^2),$$

donde $A=F_x(x^*,p_0)$ y $B=F_p(x^*,p_0)$. Si se desea un multiplicador cerrado $\alpha$ con $|\alpha|<1$, se elige

$$\delta p_n=K_c\delta x_n, \qquad K_c=\frac{\alpha-A}{B}.$$

La acción se satura a $|\delta p|\le p_{\max}$ y solo se activa dentro de una vecindad donde la linealización sea defendible.

Ejemplo calculado en el mapa logístico

Para $F(x,r)=rx(1-x)$, el punto fijo no nulo es

$$x^*=1-\frac1r.$$

Con $r_0=3.9$,

$$x^*\approx0.743590, \qquad A=F_x(x^*,r_0)=2-r_0=-1.9.$$

Además,

$$B=F_r(x^*,r_0)=x^*(1-x^*) =\frac{r_0-1}{r_0^2} \approx0.190664.$$

Para colocar el multiplicador lineal en $\alpha=0$, se obtiene

$$K_c=\frac{0-(-1.9)}{0.190664} \approx9.9655.$$

Si $x_n-x^*=10^{-4}$, la corrección es $\delta r\approx9.97\times10^{-4}$, pequeña aunque la ganancia parezca grande. Si la desviación es $0.1$, la corrección ya no es pequeña y la linealización no debe usarse. El experimento compara control apagado, control activo sin saturación y control activo con límites físicos; informa tiempo de captura, esfuerzo y robustez ante ruido.

Sincronización de sistemas caóticos

Dos sistemas sensibles pueden sincronizarse si el acoplamiento estabiliza su diferencia. Para sistemas idénticos con acoplamiento difusivo,

$$\dot x=f(x), \qquad \dot y=f(y)+kC(x-y),$$

la variedad $x=y$ es invariante. Definiendo $e=y-x$ y linealizando sobre la trayectoria sincronizada $s(t)$,

$$\dot e=\bigl[Df(s(t))-kC\bigr]e.$$

La sincronización completa es transversalmente estable si el mayor exponente condicional de esta ecuación es negativo. Reducir temporalmente $\|e(t)\|$ no basta: puede existir un estallido posterior por inestabilidad transversal, desajuste de parámetros o ruido.

La sincronización completa exige estados de igual dimensión y parámetros compatibles. La sincronización de fase solo fija una relación entre fases mientras las amplitudes difieren. La sincronización generalizada propone $y=\Phi(x)$ para alguna relación funcional. La sincronización inducida por ruido o por forzamiento común requiere controles adicionales: dos respuestas parecidas pueden seguir la misma entrada sin acoplarse entre sí.

Slotine y Li proporcionan herramientas de estabilidad y realimentación para derivar condiciones sobre errores en sistemas no lineales Slotine, 1991 Jean-Jacques E. Slotine y Weiping Li (1991) Applied Nonlinear Control Prentice Hall Ubicación consultada: caps. 3 y 8 . Para un acoplamiento concreto se construye una función de Lyapunov o se calcula el sistema variacional transversal; no se infiere estabilidad de una correlación alta.

Validación experimental y límites

Un protocolo de crisis barre el parámetro con reinicios y continuación, conserva transitorios, localiza órbitas inestables y mide tiempos. Un protocolo de control registra perturbaciones y restricciones. Un protocolo de sincronización informa error completo, exponentes condicionales, retraso, ruido y desajuste. Son preguntas relacionadas, pero sus evidencias no se sustituyen.

Para sincronización se comparan al menos cuatro casos: sistemas idénticos sin acoplar, acoplados, acoplados con desajuste y respuestas independientes bajo forzamiento común. Se mide

$$E_{\mathrm{RMS}}= \sqrt{\frac1T\int_{t_0}^{t_0+T}\|x(t)-y(t)\|^2dt},$$

pero también máximos por bloque y distribución de estallidos. Tél y Gruiz sitúan control y manifestaciones físicas del caos dentro de experimentos donde transitorios y ruido cambian la interpretación Tél, 2006 Tamás Tél y Márton Gruiz (2006) Chaotic Dynamics: An Introduction Based on Classical Mechanics Cambridge University Press Ubicación consultada: caps. 5, 6 y 9 .

Ejercicios de crisis, control y acoplamiento

1. Simula la forma normal de intermitencia tipo I para diez valores de $\varepsilon$. Define una región laminar antes de calcular y contrasta $\langle\ell\rangle\propto\varepsilon^{-1/2}$ con pendientes locales.
2. Construye datos sintéticos con reinyección no uniforme. Determina cómo cambia la distribución de fases laminares y explica qué parte de la ley depende del mecanismo local.
3. En un diagrama logístico, selecciona un cambio abrupto candidato a crisis. Diseña pruebas con barrido bidireccional, órbitas periódicas y estadísticas que puedan refutar esa clasificación.
4. Reproduce el cálculo OGY para $r_0=3.8$ y objetivos $\alpha=0$ y $\alpha=0.5$. Determina la ganancia y el rango de estados compatible con una saturación $|\delta r|\le0.01$.
5. Implementa control OGY del mapa logístico con ruido de medición. Informa tiempo de captura, fracción de pasos controlados, esfuerzo cuadrático y pérdidas de control.
6. Deriva la ecuación transversal para dos sistemas lineales acoplados y encuentra analíticamente el umbral de $k$. Usa ese caso para verificar el código antes de aplicar Lorenz.
7. Acopla dos sistemas de Lorenz por una sola variable. Barre $k$, calcula $E_{\mathrm{RMS}}$ y el mayor exponente condicional, y compara los umbrales obtenidos.
8. Diseña un control de forzamiento común: integra dos respuestas no acopladas sometidas a la misma entrada. Explica cómo distinguir coincidencia inducida por la entrada de sincronización transversal estable.
9. Redacta una tabla que separe evidencia necesaria para crisis interior, crisis de frontera, intermitencia tipo I, control exitoso y sincronización completa.