La PizarraMatemática.

Máster y doctorado

Exponentes de Lyapunov

Si una perturbación se duplica en cada paso, su crecimiento es exponencial. El exponente de Lyapunov mide esa tasa usando logaritmos. Primero lo calcularemos con factores constantes; después distinguiremos dirección, promedio y lo que una simulación de tiempo finito puede asegurar.

1. El exponente sigue perturbaciones infinitesimales, no dos puntos para siempre

Para n fijo, la diferenciabilidad da Tⁿ(x+εv)−Tⁿ(x)=εDₓTⁿv+o(ε).

Primero se aproxima una perturbación pequeña con la derivada; después se pregunta cómo crece esa parte lineal cuando n aumenta. Intercambiar ambos límites sin justificación puede dar una interpretación incorrecta.

Si ||DₓTⁿv|| crece como e^{\lambda n}, su logaritmo crece aproximadamente como λn. Positivo significa amplificación exponencial infinitesimal, negativo contracción y cero ausencia de una tasa exponencial no nula. En un espacio compacto la distancia entre dos puntos no puede crecer indefinidamente sin límite.

Un mapa lineal donde la diferencia sí es exacta

T(x)=3x y una perturbación ε.

  1. Iteración

    Tⁿ(x+ε)−Tⁿ(x)=3ⁿε.
  2. Derivada

    (Tⁿ)′=3ⁿ: aquí no hay resto no lineal.
  3. Tasa

    El crecimiento relativo es 3ⁿ y log(3ⁿ)/n=log3.
Resultado

La tasa se mide por paso.

Lo comprobamos

En mapas no lineales la igualdad de diferencias finitas se sustituye inicialmente por una aproximación.

Una separación finita puede volver a cero

D(x)=2x mod1, con x=0 y y=1/16.

  1. Repetir

    La separación avanza por 1/16,1/8,1/4,1/2 y después ambos puntos llegan a cero.

  2. Derivada

    La derivada local de Dⁿ tiene factor 2ⁿ.

  3. Distinguir

    El exponente tangente es log2, aunque esta pareja finita acabe coincidiendo.

Resultado

La distancia global módulo uno y la perturbación infinitesimal son objetos distintos.

Lo comprobamos

Una pareja concreta no define por sí sola el exponente del mapa.

2. Tomar logaritmo y dividir por el tiempo extrae la tasa

Para v≠0 definimos \lambda_n(x,v)=\frac1n\log\frac{\|D_xT^n v\|}{\|v\|}. Si el límite existe, es λ(x,v). El logaritmo es natural. Si no existe, se indica si se usa limsup o liminf; una tasa finita no es automáticamente un límite demostrado.

En escala logarítmica las tasas exponenciales se vuelven pendientes
A=diag(2,1/2)Logaritmo de la norma de Aⁿvn06v=(1,0)v=(0,1)Pendientes: log 2 y −log 2La dirección elegida importa.
Se grafica el logaritmo natural de la norma euclídea de A elevado a n por un vector unitario, con A diagonal dos, un medio. Para el vector horizontal la recta tiene pendiente log dos; para el vertical menos log dos. El mismo eje temporal va de cero a seis y la escala vertical es común. No son logaritmos de coordenadas negativas: se toma la norma positiva. Los exponentes se obtienen al dividir el crecimiento logarítmico por n.

La tasa máxima se relaciona con log||DₓTⁿ||/n cuando ese límite existe. Una dirección particular puede crecer menos. Las normas de matrices cumplen ||AB||≤||A||||B||; no podemos reemplazar sin más el logaritmo de la norma del producto por la suma de logaritmos de normas de cada paso.

La longitud inicial no cambia el exponente

T(x)=x/2 y v=7.

  1. Derivada iterada

    D Tⁿv=7·2⁻ⁿ.
  2. Razón

    ||D Tⁿv||/||v||=2⁻ⁿ; el siete se cancela.
  3. Tasa

    λₙ=−log2 para cada n>0.
Resultado

Escalar el vector inicial no cambia la tasa.

Lo comprobamos

El vector cero se excluye porque la razón no estaría definida.

Una dirección especial puede ocultar la máxima

A=diag(2,1/2).

  1. Horizontal

    v=(1,0) crece como 2ⁿ: tasa log2.
  2. Vertical

    v=(0,1) decrece como 2⁻ⁿ: tasa −log2.
  3. Máxima

    La norma de operador de Aⁿ es 2ⁿ; su tasa máxima es log2.

Resultado

La dirección de la medición forma parte de la pregunta.

Lo comprobamos

No se debe reportar la tasa vertical como si fuera la máxima del sistema.

3. La regla de la cadena produce un producto a lo largo de la órbita

La derivada iterada es DT(Tⁿ⁻¹x)···DT(Tx)DT(x). El orden importa: primero actúa la matriz de la derecha. En dimensión uno, si las derivadas no se anulan, el logaritmo del módulo convierte ese producto escalar en la suma de log|T′(Tʲx)|.

Por tanto λₙ es el promedio de esos logaritmos. Con medida invariante ergódica y log|T′| integrable, Birkhoff da para casi todo x el límite ∫log|T′|dμ. Una derivada cero puede producir −∞ con la convención log0=−∞; no se incluye sin más en un teorema invertible.

Un producto periódico de dos factores

Un cociclo escalar alterna factores cuatro y un medio.

  1. Dos pasos

    El producto es 4·(1/2)=2.
  2. Muchos pasos

    En 2m pasos el producto es 2ᵐ.

  3. Por paso

    log(2ᵐ)/(2m)=(log2)/2, igual al promedio de log4 y log(1/2).

Resultado

Promediamos logaritmos, no los propios factores.

Lo comprobamos

Es un cociclo escalar sobre un ciclo de dos estados; no se ha supuesto una derivada constante.

Las matrices no se reordenan como números

A=diag(2,1/2) y B intercambia las coordenadas.

  1. Primero A

    BA(1,0)=B(2,0)=(0,2).
  2. Primero B

    AB(1,0)=A(0,1)=(0,1/2).
  3. Comparar

    BA y AB no coinciden; cambiar el orden altera el crecimiento del vector.

Resultado

La regla de la cadena conserva el orden de la órbita.

Lo comprobamos

La suma escalar de logaritmos no se transfiere a matrices reemplazando cada una por su norma.

4. Los factores polinómicos no cambian la tasa exponencial límite

Para una matriz diagonal invertible, cada eje tiene exponente log|a|. Un vector con componentes en varios ejes tiene la tasa mayor entre los ejes que realmente aparecen. Si la matriz es constante pero no diagonalizable, los bloques de Jordan añaden factores nʲ.

Como log(nʲ)/n=j log n/n→0, los exponentes siguen siendo logaritmos de módulos de autovalores, con multiplicidad algebraica. Los valores singulares de una sola iteración describen deformación instantánea y pueden ser muy diferentes, sobre todo en matrices no normales.

Un vector mixto termina dominado por la expansión

A=diag(2,1/2), v=(1,1).

  1. Iterar

    Aⁿv=(2ⁿ,2⁻ⁿ).
  2. Norma

    Es 2ⁿ√(1+2⁻⁴ⁿ); la norma inicial es √2.

  3. Logaritmo

    Al dividir entre n, el término principal es log2 y los factores acotados aportan una cantidad que tiende a cero.

Resultado

El exponente del vector mixto es log2.

Lo comprobamos

La tasa vertical especial sigue siendo distinta.

Una gran cizalla instantánea no da exponente positivo

J tiene filas (1,100) y (0,1).

  1. Iterar

    Jⁿ tiene filas (1,100n) y (0,1).

  2. Vector

    Jⁿ(0,1)=(100n,1), con norma grande para n grande.
  3. Tasa

    log√(10000n²+1)/n→0.
Resultado

La tasa es cero aunque el crecimiento y la deformación de un paso sean grandes.

Lo comprobamos

El logaritmo de un valor singular de J en un solo paso no es por sí solo el exponente asintótico.

5. Oseledets da tasas y subespacios típicos bajo hipótesis integrables

Sea T invertible y bimedible en un espacio de probabilidad estándar, preservando μ, y A(x)∈GL(d,R) medible. Pedimos integrabilidad de log⁺||A|| y log⁺||A⁻¹||. Oseledets proporciona, en un conjunto invariante de medida total, exponentes reales y un desdoblamiento medible E₁(x)⊕···⊕Eₖ(x).

A(x) lleva Eᵢ(x) a Eᵢ(Tx) y los vectores no nulos de cada subespacio tienen su tasa λᵢ. Si la base es ergódica, las tasas y multiplicidades son constantes casi en todas partes. No se afirma continuidad de los subespacios, cotas uniformes ni existencia en absolutamente cada punto.

Para un difeomorfismo C¹ de una variedad compacta, las derivadas directa e inversa están acotadas, de modo que las condiciones logarítmicas valen para cualquier probabilidad invariante. En espacios no compactos se verifican, aunque cada matriz individual sea finita.

Sin ergodicidad las tasas pueden depender de la componente

Base de dos puntos fijos, cada uno con medida un medio; A vale dos en uno y un medio en el otro.

  1. Órbitas

    Cada punto permanece en su componente.

  2. Tasas

    En la primera el exponente es log2; en la segunda −log2.

  3. Medida

    Las condiciones integrables se cumplen, pero la base no es ergódica.

Resultado

No existe una única tasa común a casi todo el espacio de esta base.

Lo comprobamos

El teorema sin ergodicidad permite funciones de exponentes invariantes.

Matrices finitas no garantizan integrabilidad

En puntos j=1,2,…, asignar probabilidad 1/[j(j+1)] y A(j)=eʲ.

  1. Probabilidad

    Las probabilidades suman uno por telescopado.

  2. Logaritmo

    log⁺|A(j)|=j.
  3. Integral

    La media es Σ1/(j+1), que diverge.

Resultado

No podemos invocar esta versión de Oseledets con esa hipótesis faltante.

Lo comprobamos

Esto no prueba que todos los límites individuales fallen; prueba que finitud punto a punto no es integrabilidad.

6. El determinante suma las tasas contando sus multiplicidades

El módulo del determinante mide cambio de volumen. Para productos, log|det A⁽ⁿ⁾| es la suma de log|det A| a lo largo de la órbita. Bajo las hipótesis integrables y ergódicas, la suma de exponentes con multiplicidad es ∫log|det A|dμ.

Preservar volumen de una métrica fija da suma cero, no cada exponente cero. En un flujo suave de divergencia cero, el jacobiano del mapa temporal es uno mientras el flujo exista en el marco considerado. La orientación puede invertirse en un mapa: por eso usamos módulo del determinante.

Una suma cero con tres tasas no nulas

A=diag(2,2,1/4).

  1. Volumen

    det A=2·2·(1/4)=1.
  2. Exponentes

    log2 aparece dos veces; log(1/4)=−2log2 aparece una.

  3. Sumar

    log2+log2−2log2=0.
Resultado

Contar multiplicidades es necesario para la suma.

Lo comprobamos

Dos direcciones se expanden mientras una compensa contrayendo volumen.

Invertir orientación sin alterar el balance absoluto

A=diag(−2,1/2).

  1. Determinante

    det A=−1; su módulo es uno.
  2. Tasas

    Los módulos de autovalores son dos y un medio.

  3. Suma

    log2−log2=0.
Resultado

El volumen absoluto se conserva aunque cambie la orientación.

Lo comprobamos

El logaritmo usado aquí es log|det A|, no el logaritmo real de un número negativo.

7. Las coordenadas deben ser comparables y la unidad de tiempo debe declararse

Si c||v||≤||v||′≤C||v|| con constantes positivas uniformes a lo largo de la órbita, la diferencia entre tasas queda acotada por log(C/c)/n y tiende a cero. Un cambio suave de coordenadas con derivada e inversa uniformemente acotadas también preserva el exponente. En una variedad compacta, métricas riemannianas suaves son comparables de ese modo.

Sin uniformidad en un espacio no compacto puede cambiar.

Para T(x)=x+1 en R, la derivada es uno y la tasa cero.

Con y=eˣ, el mapa se transforma en S(y)=ey y su tasa es uno; el cambio de coordenadas crece a lo largo de la órbita.

Un paso de Tᵐ equivale a m pasos de T, por lo que la tasa por nueva iteración es mλ. En un flujo se divide por tiempo físico; si se muestrea cada h unidades, la tasa por muestra es hλ.

Cambiar la norma sin cambiar la tasa

En R², usar norma euclídea o norma máxima.

  1. Comparar

    ||v||∞≤||v||₂≤√2||v||∞.
  2. Cocientes

    El cambio en el logaritmo de la razón está acotado por log√2.

  3. Dividir

    Ese error dividido por n tiende a cero.

Resultado

Las dos normas dan el mismo exponente.

Lo comprobamos

El argumento depende de constantes independientes del tiempo.

Comparar tasas con unidades distintas

Un flujo tiene tasa física −3 por unidad de tiempo y se observa cada h=1/2.

  1. Una muestra

    La contracción lleva factor e⁻³ᐟ².

  2. Por muestra

    El exponente es −3/2.

  3. Por tiempo

    Dividir por h recupera (−3/2)/(1/2)=−3.
Resultado

Un cambio de unidad no cambia el fenómeno, pero sí el número reportado por muestra.

Lo comprobamos

Para T⁴ la tasa por iteración del nuevo mapa sería cuatro veces la original.

8. Cero no equivale a estabilidad y positivo no prueba caos por sí solo

Para J con filas (1,1),(0,1), Jⁿ(0,1)=(n,1).

La norma crece sin cota pero su tasa logarítmica tiende a cero. Un exponente cero no descarta crecimiento polinómico ni decide estabilidad.

Un exponente cero permite crecimiento que no es exponencial
Jⁿ(0,1)=(n,1)Tasa finita: log(√(n²+1))/nlog 2tasa de J064La norma crece como n.El exponente límite es cero.
Para la cizalla J con filas uno, uno y cero, uno, el vector vertical se transforma tras n pasos en n, uno. Su norma es raíz de n al cuadrado más uno y crece sin cota. La curva verde grafica exactamente su tasa logarítmica dividida por n para n de uno a sesenta y cuatro; se aproxima a cero. La recta azul log dos sirve como referencia de una expansión exponencial distinta. No se calcula una tasa en n igual a cero, donde la división no está definida.

T(x)=3x en R tiene exponente positivo y órbitas muy simples que se alejan; incluso una medida puntual en un punto fijo repulsivo puede registrar una tasa positiva. Si todos los exponentes de una medida diferenciable son no nulos, se habla de medida hiperbólica. La teoría de Pesin necesita regularidad adicional como C¹⁺ᵅ y da variedades en puntos típicos, con tamaños y tasas que pueden depender del punto.

Una perturbación no acotada con tasa cero

J con filas (1,1),(0,1), v=(0,1).

  1. Norma

    ||Jⁿv||=√(n²+1).
  2. Crecimiento

    Tiende a infinito aproximadamente como n.

  3. Exponente

    log√(n²+1)/n→0.
Resultado

Ausencia de crecimiento exponencial no significa perturbación acotada.

Lo comprobamos

La gráfica muestra tasas de tiempo finito que se acercan a cero.

Una medida que vive solo en el punto repulsivo

T(x)=3x y μ=δ₀.

  1. Invariante

    T(0)=0, así que la masa puntual se conserva.
  2. Órbita típica para μ

    Solo se observa el punto cero.

  3. Derivada

    La tasa tangente allí es log3, positiva.

Resultado

Positividad tangente y complejidad estadística no son sinónimos.

Lo comprobamos

La medida elegida forma parte de cualquier afirmación sobre puntos típicos.

9. Renormalizar evita perder la cuenta y ortogonalizar evita perder direcciones

Para estimar la tasa máxima, propagamos un vector tangente, anotamos el logaritmo de su norma y lo normalizamos antes del siguiente paso. Sumar esos logaritmos conserva el crecimiento total sin formar números enormes. Un vector en una dirección excepcional puede medir una tasa menor.

Para varias tasas se propaga un marco y se repite QR, acumulando logaritmos de las diagonales de R con una convención consistente. Sin ortogonalizar, los vectores tienden a alinearse con la dirección dominante. Transitorios, precisión y tiempo total afectan al cálculo: se necesitan comprobaciones de convergencia y no se demuestra un límite solo con una tabla larga.

Contar un producto que desbordaría

A=diag(2,1/2), vector inicial horizontal unitario y dos mil pasos.

  1. Cada paso

    La norma del vector propagado es dos; al normalizar vuelve al vector horizontal unitario.

  2. Acumular

    La suma logarítmica es 2000log2, sin construir 2²⁰⁰⁰.

  3. Promediar

    Dividir entre dos mil da log2.

Resultado

La renormalización conserva la tasa y evita el desbordamiento.

Lo comprobamos

Con vector vertical el mismo procedimiento daría −log2.

QR distingue las dos direcciones en el modelo diagonal

Marco inicial con columnas e₁ y e₂.

  1. Propagar

    A convierte las columnas en 2e₁ y e₂/2.

  2. QR

    Q vuelve a ser la identidad y la diagonal de R es dos, un medio.

  3. Acumular

    Los promedios de logaritmos dan log2 y −log2, con suma cero.

Resultado

Ortogonalizar permite registrar también la dirección que se contrae.

Lo comprobamos

En problemas generales se comprueba estabilidad numérica y convergencia de las estimaciones.

Ahora lo resolvemos de principio a fin.

EJEMPLO RESUELTO 1

Expansión escalar

Calcula el exponente de T(x)=3x.

  1. Mapa

    T(x)=3x.
  2. Derivada

    T′(x)=3 en todo punto.
  3. Iterada

    (Tⁿ)′(x)=3ⁿ.
  4. Razón

    Para v≠0, ||D Tⁿv||/||v||=3ⁿ.
  5. Logaritmo

    log(3ⁿ)/n=log3.
  6. Límite

    La expresión ya es constante en n: el exponente es log3, positivo y medido por iteración.

Resultado

El exponente es log3.

Lo comprobamos

Es positivo porque 3>1.

EJEMPLO RESUELTO 2

Contracción escalar

Calcula el exponente de T(x)=x/2.

  1. Mapa

    T(x)=x/2.
  2. Derivada

    T′=1/2.
  3. Iteración

    (Tⁿ)′=2⁻ⁿ.
  4. Norma

    La razón de longitudes es 2⁻ⁿ para todo vector inicial no nulo.

  5. Tasa

    log(2⁻ⁿ)/n=−log2.
  6. Límite

    El exponente es −log2. El signo negativo indica contracción exponencial, no un logaritmo de una longitud negativa.

Resultado

El exponente es −log2.

Lo comprobamos

La distancia decrece exponencialmente.

Los símbolos que hemos usado

T y Tⁿx
Aplicación diferenciable del espacio de estados y órbita de x tras n iteraciones.
v≠0
Vector tangente inicial que representa una dirección infinitesimal de perturbación.
DₓTⁿ
Derivada de la aplicación iterada en x, producto ordenado de derivadas a lo largo de la órbita.
λ(x,v)
Límite, si existe, de n⁻¹log(||DₓTⁿv||/||v||), usando logaritmo natural.
λₙ(x,v)
Tasa de tiempo finito calculada con n iteraciones; no es por sí misma el límite asintótico.
μ
Probabilidad invariante para la transformación base: μ(T⁻¹B)=μ(B) para todo conjunto medible B.
A⁽ⁿ⁾(x)
Producto A(Tⁿ⁻¹x)···A(Tx)A(x) de un cociclo de matrices.
log⁺r=max(log r,0)
Parte positiva del logaritmo, usada para controlar la integrabilidad de matrices y sus inversas.
Eᵢ(x)
Subespacios medibles invariantes asociados a exponentes distintos en el teorema de Oseledets.

Ahora te toca a ti.

Tienes cinco ejercicios. Inténtalos; las pistas te ayudan si te atascas.

Practicar este tema →

Para ampliar