Cómo encontrar el exponente de Lyapunov para el sistema acoplado

Answer proporciona un software para calcular el exponente condicional de Lyapunov (CLE) para osciladores acoplados en sincronización de caos. Sin embargo, es difícil de seguir y no hay salida gráfica de la trama (y está en C, más complejo). ¿Alguien sabe cómo modificar LET Toolbox, que es muy bueno para el sistema desacoplado, pero no entiendo cómo trabajar con un sistema sincronizado para acomodar el CLE.

1. Tengo una confusión sobre cómo incluir la señal del controlador al calcular la matriz jacobiana para CLE, ya que la teoría dice que se encuentra CLE para el sistema de respuesta, por lo que no tenemos que encontrar el jacobiano para el controlador y el sistema de respuesta para osciladores similares (accionamiento y respuesta). ¿O deberíamos considerar tanto la unidad como el sistema de respuesta en el software y proceder como si fuera un sistema único?
2. Cómo acomodar un forzamiento externo como un proceso aleatorio en una ecuación de estado, si hay alguno en CLE.
3. ¿Hay alguna otra implementación para CLE?

gracias

Aquí hay un enlace a un documento antiguo en J Comp Phys que calcula la LE para sistemas dinámicos aleatorios acoplados http://mesoscopic.mines.edu/acoustics-old/preprints/lyapunov.pdf Quizás lo encuentre más legible. Es un algoritmo estable y bien probado. El código fue desarrollado por mi coautor que ahora está en KU (en caso de que desee el código).