Como encontrar o expoente Lyapunov para o sistema acoplado

O Answer fornece um software para calcular o expoente Lyapunov condicional (CLE) para osciladores acoplados na sincronização do caos. No entanto, é difícil de seguir e não há saída gráfica do gráfico (e é em C, mais complexo). Alguém sabe como modificar a LET Toolbox, que é muito boa para o sistema desacoplado, mas eu não entendo como trabalhar com o sistema sincronizado para acomodar o CLE.

1. Eu tenho uma confusão sobre como incluir o sinal do driver ao calcular a matriz jacobiana para o CLE, uma vez que a teoria diz que o CLE é encontrado para o sistema de resposta, portanto, não precisamos encontrar o jacobiano para o motorista, bem como o sistema de resposta para o CLE. osciladores similares (drive e resposta). Ou devemos considerar o sistema de acionamento e resposta no software e proceder como se fosse um sistema único?
2. Como acomodar uma força externa como um processo aleatório em uma equação de estado, se houver no CLE.
3. Existe alguma outra implementação para o CLE?

obrigado

Aqui está um link para um artigo antigo no J Comp Phys que calcula o LE para sistemas dinâmicos aleatórios acoplados http://mesoscopic.mines.edu/acoustics-old/preprints/lyapunov.pdf Talvez você o ache mais legível. É um algoritmo estável e bem testado. O código foi desenvolvido pelo meu co-autor, que agora está na KU (caso você queira o código).