Entender la ecuación del exponente de Lyapunov

Estoy tratando de calcular el exponente de Lyapunov para un sistema dinámico simple, pero creo que no entendí bien la ecuación. Mis cálculos conducen constantemente a cero, aunque estoy variando las condiciones iniciales y otros parámetros y el sistema definitivamente no es estable sino caótico cuando lo estoy trazando. Parece que no puedo entender qué está mal, especialmente porque cuando miro la ecuación y solo me doy cuenta de que el exponente de Lyapunov siempre debe ser cero. ¡Obviamente mi entendimiento está faltando!

La ecuación que he usado es

λ = límite norte 1 norte i = 0 norte 1 registro | F ( X i ) |

Cuando norte Creo que el exponente de Lyapunov, λ , debe ser cero debido a la 1 / norte factor, entonces, ¿qué es lo que no estoy entendiendo? También me siento un poco confundido acerca de cuál es la variable del iterador, i , es. ¿Es el paso del tiempo? Pero, ¿no debería el exponente de Lyapunov comparar trayectorias entre sí? ¿Cómo se incluyen en la ecuación? Siento que no he entendido la definición del exponente de Lyapunov, ¡así que espero que alguien pueda explicarme un poco!

tenga en cuenta que su suma también depende de n. Si la suma aumenta más rápido que ~n, el exponente no será cero.
¡Gracias! Me di cuenta de que mis cálculos en realidad no estaban equivocados (espero). Obtuve la suma 89 y la dividí con el número de iteraciones 400, pero mi programa estimó que 89/400 era 0,0 en lugar de 0,22. Solo tuve que convertir el número de iteración en un flotante... suspiro
Primera ley de la física numérica: El código siempre se rompe :)

Respuestas (1)

λ = límite norte 1 norte i = 0 norte 1 registro | F ( X i ) |

Cuando norte Creo que el exponente de Lyapunov, λ , debe ser cero debido a la 1 / norte factor, entonces, ¿qué es lo que no estoy entendiendo?

No, el número de sumandos en tu suma va hacia tan rápido como tu 1 norte el factor va hacia 0 . por finito norte (es decir, en la aplicación), la ecuación significa que obtienes tantos valores de | F ( X i ) | como puedas y promediarlos.

También me siento un poco confundido acerca de cuál es la variable del iterador, i , es. ¿Es el paso del tiempo?

Si por paso de tiempo se refiere al paso de tiempo de integrar un sistema de tiempo continuo (es decir, una ecuación diferencial), entonces está aplicando la fórmula incorrecta. (Aún así, la aplicación de esta fórmula debería darle algo diferente de 0 , pero no el exponente de Lyapunov.)

La fórmula anterior solo funciona para sistemas de tiempo discreto (es decir, mapas), cuyas iteraciones puede llamar pasos de tiempo.

Pero, ¿no debería el exponente de Lyapunov comparar trayectorias entre sí? ¿Cómo se incluyen en la ecuación?

Como se dijo anteriormente, la ecuación que está utilizando solo funciona para mapas que no tienen trayectorias. Supongamos que el estado actual de nuestro sistema es X i y el estado perturbado es X i + ε (con ε 0 por un cambio). La separación de su sistema perturbado de su sistema no perturbado en el próximo paso de tiempo es F ( X i + ε ) F ( X i ) . Y el exponente "local" de Lyapunov λ i – es decir, el punto de Lyapunov para el mismo punto X i  - es:

λ i = límite ε 0 registro | F ( X i + ε ) F ( X i ) ε | = registro | límite ε 0 F ( X i + ε ) F ( X i ) ε | = registro | F ( X i ) |

En el último paso, usamos la definición misma de F . (Para obtener el exponente de Lyapunov "global", solo promedie los locales).