Fyskode Learning

Caos y diagnósticos · Universitario · 7 horas

Dinámica simbólica, particiones de Markov y entropía

Codificación de órbitas mediante símbolos, matrices de transición, subdesplazamientos de tipo finito y cálculo de entropía topológica.

Sustituir coordenadas por itinerarios

Una órbita numérica contiene coordenadas con muchos decimales. La dinámica simbólica conserva otra información: el orden en que la órbita visita regiones previamente definidas. Sea F:XXF:X\to X un mapa y sea

P={P0,P1,,Pk1}\mathcal P=\{P_0,P_1,\ldots,P_{k-1}\}

una partición del dominio relevante. A un punto xx se le asigna el itinerario

s(x)=s0s1s2,sn=jsiFn(x)Pj.s(x)=s_0s_1s_2\ldots, \qquad s_n=j \quad\text{si}\quad F^n(x)\in P_j.

El desplazamiento a la izquierda,

σ(s0s1s2)=s1s2s3,\sigma(s_0s_1s_2\ldots)=s_1s_2s_3\ldots,

representa una iteración de FF. Si la codificación está bien construida, se cumple s(F(x))=σ(s(x))s(F(x))=\sigma(s(x)). Esta identidad es una semiconjugación: la descripción simbólica reproduce la evolución, aunque distintos puntos pueden compartir un itinerario, sobre todo en fronteras de la partición. Una conjugación exige además correspondencia uno a uno y continuidad en ambas direcciones.

No cualquier división coloreada produce un modelo simbólico útil. Si una frontera atraviesa repetidamente las órbitas o si el siguiente símbolo depende de una historia arbitrariamente larga, la secuencia pierde la propiedad de memoria finita. Devaney presenta el desplazamiento como un sistema dinámico por derecho propio y usa conjugaciones para transferir propiedades como periodicidad y sensibilidad Devaney, 2022 Robert L. Devaney (2022) An Introduction to Chaotic Dynamical Systems 3.ª ed. · CRC Press Ubicación consultada: caps. 1–3 .

Ejemplo trabajado: el mapa tienda

Considérese el mapa tienda completo

T(x)={2x,0x12,2(1x),12<x1.T(x)= \begin{cases} 2x, & 0\le x\le \tfrac12,\\ 2(1-x), & \tfrac12<x\le1. \end{cases}

Dividimos [0,1][0,1] en P0=[0,1/2]P_0=[0,1/2] y P1=(1/2,1]P_1=(1/2,1]. Cada rama cubre todo el intervalo: T(P0)=T(P1)=[0,1]T(P_0)=T(P_1)=[0,1]. Por ello, después de cualquiera de los símbolos pueden aparecer 00 o 11. Hay 2n2^n palabras admisibles de longitud nn, salvo las convenciones duplicadas en preimágenes de la frontera.

Calculemos un itinerario. Con x0=0.2x_0=0.2,

x1=0.4,x2=0.8,x3=0.4,x4=0.8,x_1=0.4, \quad x_2=0.8, \quad x_3=0.4, \quad x_4=0.8, \ldots

y la secuencia es 0010101000101010\ldots. Después de dos iteraciones entra en un ciclo de periodo dos. En cambio, especificar una palabra finita permite localizar un intervalo de condiciones iniciales. Para la palabra 010010, exigimos

xP0,T(x)P1,T2(x)P0.x\in P_0,\qquad T(x)\in P_1,\qquad T^2(x)\in P_0.

La primera condición da x[0,1/2]x\in[0,1/2]. En esa rama T(x)=2xT(x)=2x, por lo que la segunda exige x(1/4,1/2]x\in(1/4,1/2]. Allí T(x)>1/2T(x)>1/2 y

T2(x)=2(12x)=24x.T^2(x)=2(1-2x)=2-4x.

Pedir 024x1/20\le2-4x\le1/2 produce 3/8x1/23/8\le x\le1/2. Al resolver cada condición en orden, la palabra queda asociada con un intervalo de longitud 232^{-3}, salvo extremos compartidos. Una gráfica de telaraña enseña las coordenadas; el itinerario enseña la combinatoria de las ramas.

Telarañas de un mapa unimodal sobre su gráfica y la diagonal
La partición natural separa las ramas a ambos lados del punto crítico. Cada tramo vertical determina el siguiente símbolo; los puntos que caen exactamente en la frontera requieren una convención explícita.

Atlas reproducible de sistemas dinámicos Función: plot_logistic_cobweb

Descargar .py

Matrices de transición y palabras admisibles

Cuando no todas las transiciones son posibles, se usa una matriz AA de ceros y unos:

Aij=1se permite pasar de Pi a Pj.A_{ij}=1 \quad\Longleftrightarrow\quad \text{se permite pasar de }P_i\text{ a }P_j.

El espacio simbólico asociado contiene las secuencias para las que Asn,sn+1=1A_{s_n,s_{n+1}}=1 en todo nn. Se llama subdesplazamiento de tipo finito. El número de caminos de longitud nn se obtiene sumando las entradas de AnA^n; la entrada (An)ij(A^n)_{ij} cuenta caminos de ii a jj con nn transiciones. La traza tr(An)\operatorname{tr}(A^n) cuenta caminos cerrados con un origen marcado y se relaciona con puntos periódicos simbólicos.

Tomemos la regla en la que nunca aparecen dos unos consecutivos:

A=(1110).A= \begin{pmatrix} 1&1\\ 1&0 \end{pmatrix}.

Desde 00 puede seguir 00 o 11; desde 11 solo puede seguir 00. Para longitud uno hay dos palabras. Para longitud dos son 00,01,1000,01,10; para longitud tres, 000,001,010,100,101000,001,010,100,101. Los conteos siguen la recurrencia de Fibonacci. El polinomio característico es

λ2λ1=0,\lambda^2-\lambda-1=0,

y el radio espectral es la razón áurea φ=(1+5)/2\varphi=(1+\sqrt5)/2. Así, el número de palabras crece aproximadamente como CφnC\varphi^n, menor que las 2n2^n palabras del desplazamiento completo.

Alligood, Sauer y Yorke conectan las bandas de una herradura con secuencias permitidas y muestran cómo el conteo simbólico traduce geometría en combinatoria Alligood, 1996 Kathleen T. Alligood, Tim D. Sauer y James A. Yorke (1996) Chaos: An Introduction to Dynamical Systems Springer Ubicación consultada: caps. 5 y 7 . La matriz no debe inferirse solo porque dos colores aparecen consecutivamente en una órbita corta. Hay que verificar imágenes de regiones: F(Pi)PjF(P_i)\cap P_j\ne\varnothing, controlar fronteras y comprobar que las transiciones puedan concatenarse dentro del conjunto invariante.

Particiones de Markov

En un sistema hiperbólico, una partición de Markov alinea las fronteras de sus rectángulos con direcciones estables e inestables. La imagen de una franja estable atraviesa los rectángulos permitidos de forma compatible, y la preimagen hace lo correspondiente con las franjas inestables. Esta propiedad permite determinar itinerarios mediante transiciones locales y reconstruir puntos a partir de secuencias, excepto ambigüedades de frontera.

La palabra “rectángulo” es dinámica, no necesariamente euclidiana: significa que, dados dos puntos cercanos del conjunto, la intersección entre la hoja estable de uno y la inestable del otro permanece en la pieza. En una herradura ideal las dos tiras forman una partición evidente. En datos o en un atractor curvo, dibujar cajas uniformes rara vez produce una partición de Markov exacta.

Un procedimiento numérico razonable comienza con una sección de Poincaré o un mapa conocido, aproxima variedades, propone fronteras siguiendo esas curvas y construye la tabla de intersecciones. Después prueba palabras: para cada secuencia permitida de longitud nn, busca una condición inicial que la realice y refina la malla. Las palabras espurias indican una transición demasiado permisiva; las palabras observadas pero prohibidas indican una frontera o una tabla incorrecta. El resultado sigue siendo aproximado hasta disponer de cotas geométricas.

Región estirada y plegada que regresa en dos tiras
Etiquetar las dos tiras como $0$ y $1$ produce secuencias, pero una partición de Markov requiere compatibilidad con las direcciones estables e inestables y no solo una división visual.

Atlas reproducible de sistemas dinámicos Funciones: plot_henon_folding, animate_horseshoe_mechanism

Descargar .py

Entropía topológica

La entropía topológica mide la tasa exponencial con la que crece el número de historias orbitales distinguibles al aumentar el tiempo y refinar la resolución. Para un subdesplazamiento de tipo finito con matriz irreducible AA,

htop=logρ(A),h_{\mathrm{top}}=\log\rho(A),

donde ρ(A)\rho(A) es el radio espectral. En el desplazamiento completo con kk símbolos, ρ(A)=k\rho(A)=k y htop=logkh_{\mathrm{top}}=\log k. En el ejemplo sin unos consecutivos,

htop=logφ0.4812,h_{\mathrm{top}}=\log\varphi\approx0.4812,

mientras que el desplazamiento binario completo tiene log20.6931\log2\approx0.6931. La diferencia cuantifica una restricción combinatoria: ambos sistemas usan dos símbolos, pero el segundo admite más historias.

La base del logaritmo fija la unidad. Con logaritmo natural se obtienen nats por iteración; con base dos, bits por iteración. También importa la escala temporal. Si FF es un mapa de retorno cuyo retorno medio dura Tˉ\bar T, dividir informalmente por Tˉ\bar T solo produce una tasa física aproximada cuando la variabilidad de retornos está controlada.

Ott distingue la entropía topológica, que considera todas las órbitas distinguibles, de las entropías dependientes de una medida, que ponderan lo que una órbita típica visita Ott, 2002 Edward Ott (2002) Chaos in Dynamical Systems 2.ª ed. · Cambridge University Press Ubicación consultada: cap. 4 . Contar cajas ocupadas en una serie finita no entrega automáticamente ninguna de las dos. Una estimación debe declarar partición, longitud de palabras, corrección por muestreo y región donde el crecimiento parece exponencial.

Procedimiento de verificación

Para pasar de una trayectoria a un modelo simbólico se propone la siguiente secuencia:

  1. Define las piezas PiP_i mediante fronteras reproducibles, no a ojo después de observar la respuesta.
  2. Decide cómo clasificar puntos en las fronteras y registra los casos ambiguos.
  3. Calcula imágenes y preimágenes de cada pieza para proponer la matriz AA.
  4. Enumera todas las palabras admitidas de longitudes pequeñas y busca realizaciones dinámicas.
  5. Compara conteos observados con ij(An1)ij\sum_{ij}(A^{n-1})_{ij} sin confundir ausencia en una muestra con prohibición matemática.
  6. Estima ρ(A)\rho(A) y documenta la base del logaritmo y la unidad temporal.

Una matriz estable al refinar la partición es evidencia más sólida que una secuencia larga. Aun así, afirmar conjugación requiere justificar unicidad de la codificación y continuidad. Si solo se verifica que todo itinerario dinámico respeta la matriz, la conclusión adecuada es una codificación o semiconjugación candidata.

Ejercicios de codificación y entropía

  1. Para el mapa tienda, determina los intervalos de condiciones iniciales asociados con 000000, 011011 y 10111011. Verifica una condición de cada intervalo por iteración directa y trata por separado los extremos.
  2. Construye todas las palabras de longitud hasta seis que no contengan 1111. Comprueba la recurrencia de Fibonacci y compárala con las potencias de la matriz AA.
  3. Considera A=(110011100)A=\begin{pmatrix}1&1&0\\0&1&1\\1&0&0\end{pmatrix}. Dibuja su grafo, enumera los ciclos de longitud tres y calcula numéricamente logρ(A)\log\rho(A).
  4. Codifica una órbita logística usando la frontera x=1/2x=1/2. Repite desplazando la frontera a 0.490.49 y 0.510.51. Mide cuántos símbolos cambian y explica por qué esa sensibilidad no prueba por sí sola una partición de Markov.
  5. A partir de la animación de herradura, propón dos piezas y una matriz de transición. Enumera qué propiedades geométricas no pueden verificarse solo con la animación.
  6. Compara htoph_{\mathrm{top}} en nats y bits para el desplazamiento completo de dos y tres símbolos. Interpreta la razón entre los números de palabras después de diez iteraciones.
  7. Diseña una prueba que intente refutar una matriz simbólica propuesta: incluye palabras prohibidas, refinamiento de fronteras, múltiples condiciones iniciales y un criterio para registrar ambigüedades.
  8. Explica con un contraejemplo conceptual por qué observar todos los pares ijij permitidos no garantiza que toda palabra larga formada con esos pares sea realizable en el sistema original.

Fuentes consultadas

Obras citadas en el desarrollo; los localizadores indican los capítulos o secciones consultados.
  1. Robert L. Devaney (2022). An Introduction to Chaotic Dynamical Systems. 3.ª ed. CRC Press. ISBN 978-1-032-15046-8.

    caps. 1–3 · Conjugación, desplazamiento simbólico y definiciones de caos.
  2. Kathleen T. Alligood, Tim D. Sauer y James A. Yorke (1996). Chaos: An Introduction to Dynamical Systems. Springer. ISBN 0-387-94677-2.

    caps. 1, 5 y 7 · Itinerarios, herradura, matrices de transición y conteo de órbitas.
  3. Edward Ott (2002). Chaos in Dynamical Systems. 2.ª ed. Cambridge University Press. ISBN 978-0-521-01084-9.

    caps. 2 y 4 · Entropía topológica, particiones y descripciones simbólicas.