Comment trouver l'exposant de Lyapunov pour le système couplé

Answer donne un logiciel pour calculer l' exposant de Lyapunov conditionnel (CLE) pour les oscillateurs couplés dans la synchronisation du chaos. Cependant, c'est difficile à suivre et il n'y a pas de sortie graphique du tracé (et c'est en C, plus complexe). Est-ce que quelqu'un sait comment modifier la boîte à outils LET qui est très bonne pour un système non couplé mais je ne comprends pas comment travailler avec un système synchronisé de manière à accueillir le CLE.

1. J'ai une confusion quant à la façon d'inclure le signal du pilote lors du calcul de la matrice jacobienne pour CLE car la théorie dit que CLE est trouvé pour le système de réponse, donc n'avons-nous pas à trouver le jacobien pour le pilote ainsi que le système de réponse pour oscillateurs similaires (entraînement et réponse). Ou devrions-nous considérer à la fois le système d'entraînement et de réponse dans le logiciel et procéder comme s'il s'agissait d'un seul système?
2. Comment gérer un forçage externe comme un processus aléatoire dans une équation d'état le cas échéant dans CLE.
3. Existe-t-il une autre implémentation pour CLE?

Merci

Voici un lien vers un ancien article de J Comp Phys qui calcule le LE pour les systèmes dynamiques aléatoires couplés http://mesoscopic.mines.edu/acoustics-old/preprints/lyapunov.pdf Peut-être le trouverez-vous plus lisible. C'est un algorithme stable et bien testé. Le code a été développé par mon co-auteur qui est maintenant à KU (au cas où vous voudriez le code).