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.
Si necesitas repasar:
1. El exponente sigue perturbaciones infinitesimales, no dos puntos para siempre
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 , 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 ε.
Iteración
Tⁿ(x+ε)−Tⁿ(x)=3ⁿε.Derivada
(Tⁿ)′=3ⁿ: aquí no hay resto no lineal.Tasa
El crecimiento relativo es 3ⁿ y log(3ⁿ)/n=log3.
La tasa se mide por paso.
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.
Repetir
La separación avanza por 1/16,1/8,1/4,1/2 y después ambos puntos llegan a cero.
Derivada
La derivada local de Dⁿ tiene factor 2ⁿ.
Distinguir
El exponente tangente es log2, aunque esta pareja finita acabe coincidiendo.
La distancia global módulo uno y la perturbación infinitesimal son objetos distintos.
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 . 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.
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.
Derivada iterada
D Tⁿv=7·2⁻ⁿ.Razón
||D Tⁿv||/||v||=2⁻ⁿ; el siete se cancela.Tasa
λₙ=−log2 para cada n>0.
Escalar el vector inicial no cambia la tasa.
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).
Horizontal
v=(1,0) crece como 2ⁿ: tasa log2.Vertical
v=(0,1) decrece como 2⁻ⁿ: tasa −log2.Máxima
La norma de operador de Aⁿ es 2ⁿ; su tasa máxima es log2.
La dirección de la medición forma parte de la pregunta.
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.
Dos pasos
El producto es 4·(1/2)=2.Muchos pasos
En 2m pasos el producto es 2ᵐ.
Por paso
log(2ᵐ)/(2m)=(log2)/2, igual al promedio de log4 y log(1/2).
Promediamos logaritmos, no los propios factores.
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.
Primero A
BA(1,0)=B(2,0)=(0,2).Primero B
AB(1,0)=A(0,1)=(0,1/2).Comparar
BA y AB no coinciden; cambiar el orden altera el crecimiento del vector.
La regla de la cadena conserva el orden de la órbita.
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).
Iterar
Aⁿv=(2ⁿ,2⁻ⁿ).Norma
Es 2ⁿ√(1+2⁻⁴ⁿ); la norma inicial es √2.
Logaritmo
Al dividir entre n, el término principal es log2 y los factores acotados aportan una cantidad que tiende a cero.
El exponente del vector mixto es log2.
La tasa vertical especial sigue siendo distinta.
Una gran cizalla instantánea no da exponente positivo
J tiene filas (1,100) y (0,1).
Iterar
Jⁿ tiene filas (1,100n) y (0,1).
Vector
Jⁿ(0,1)=(100n,1), con norma grande para n grande.Tasa
log√(10000n²+1)/n→0.
La tasa es cero aunque el crecimiento y la deformación de un paso sean grandes.
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.
Órbitas
Cada punto permanece en su componente.
Tasas
En la primera el exponente es log2; en la segunda −log2.
Medida
Las condiciones integrables se cumplen, pero la base no es ergódica.
No existe una única tasa común a casi todo el espacio de esta base.
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ʲ.
Probabilidad
Las probabilidades suman uno por telescopado.
Logaritmo
log⁺|A(j)|=j.Integral
La media es Σ1/(j+1), que diverge.
No podemos invocar esta versión de Oseledets con esa hipótesis faltante.
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).
Volumen
det A=2·2·(1/4)=1.Exponentes
log2 aparece dos veces; log(1/4)=−2log2 aparece una.
Sumar
log2+log2−2log2=0.
Contar multiplicidades es necesario para la suma.
Dos direcciones se expanden mientras una compensa contrayendo volumen.
Invertir orientación sin alterar el balance absoluto
A=diag(−2,1/2).
Determinante
det A=−1; su módulo es uno.Tasas
Los módulos de autovalores son dos y un medio.
Suma
log2−log2=0.
El volumen absoluto se conserva aunque cambie la orientación.
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.
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.
Comparar
||v||∞≤||v||₂≤√2||v||∞.Cocientes
El cambio en el logaritmo de la razón está acotado por log√2.
Dividir
Ese error dividido por n tiende a cero.
Las dos normas dan el mismo exponente.
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.
Una muestra
La contracción lleva factor e⁻³ᐟ².
Por muestra
El exponente es −3/2.
Por tiempo
Dividir por h recupera (−3/2)/(1/2)=−3.
Un cambio de unidad no cambia el fenómeno, pero sí el número reportado por muestra.
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
La norma crece sin cota pero su tasa logarítmica tiende a cero. Un exponente cero no descarta crecimiento polinómico ni decide estabilidad.
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).
Norma
||Jⁿv||=√(n²+1).Crecimiento
Tiende a infinito aproximadamente como n.
Exponente
log√(n²+1)/n→0.
Ausencia de crecimiento exponencial no significa perturbación acotada.
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 μ=δ₀.
Invariante
T(0)=0, así que la masa puntual se conserva.Órbita típica para μ
Solo se observa el punto cero.
Derivada
La tasa tangente allí es log3, positiva.
Positividad tangente y complejidad estadística no son sinónimos.
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.
Cada paso
La norma del vector propagado es dos; al normalizar vuelve al vector horizontal unitario.
Acumular
La suma logarítmica es 2000log2, sin construir 2²⁰⁰⁰.
Promediar
Dividir entre dos mil da log2.
La renormalización conserva la tasa y evita el desbordamiento.
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₂.
Propagar
A convierte las columnas en 2e₁ y e₂/2.
QR
Q vuelve a ser la identidad y la diagonal de R es dos, un medio.
Acumular
Los promedios de logaritmos dan log2 y −log2, con suma cero.
Ortogonalizar permite registrar también la dirección que se contrae.
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.
Mapa
T(x)=3x.Derivada
T′(x)=3 en todo punto.Iterada
(Tⁿ)′(x)=3ⁿ.Razón
Para v≠0, ||D Tⁿv||/||v||=3ⁿ.Logaritmo
log(3ⁿ)/n=log3.Límite
La expresión ya es constante en n: el exponente es log3, positivo y medido por iteración.
El exponente es log3.
Es positivo porque 3>1.
EJEMPLO RESUELTO 2
Contracción escalar
Calcula el exponente de T(x)=x/2.
Mapa
T(x)=x/2.Derivada
T′=1/2.Iteración
(Tⁿ)′=2⁻ⁿ.Norma
La razón de longitudes es 2⁻ⁿ para todo vector inicial no nulo.
Tasa
log(2⁻ⁿ)/n=−log2.Límite
El exponente es −log2. El signo negativo indica contracción exponencial, no un logaritmo de una longitud negativa.
El exponente es −log2.
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 →