Come trovare l'esponente di Lyapunov per il sistema accoppiato

Answer fornisce un software per il calcolo dell'esponente condizionale di Lyapunov (CLE) per oscillatori accoppiati nella sincronizzazione del caos. Tuttavia, è difficile da seguire e non esiste un output grafico della trama (ed è in C, più complesso). Qualcuno sa come modificare il LET Toolbox che è molto buono per il sistema non accoppiato ma non capisco come lavorare con il sistema sincronizzato in modo da adattarsi al CLE.

1. Ho una confusione su come includere il segnale del driver durante il calcolo della matrice Jacobiana per CLE poiché la teoria dice che CLE viene trovato per il sistema di risposta, quindi non dobbiamo trovare il Jacobiano per il driver così come il sistema di risposta per oscillatori simili (azionamento e risposta). O dovremmo considerare sia il drive che il sistema di risposta nel software e procedere come se fosse un singolo sistema?
2. Come accogliere una forzatura esterna come un processo casuale in un'equazione di stato se presente in CLE.
3. C'è qualche altra implementazione per CLE?

grazie

Ecco un link ad un vecchio documento in J Comp Phys che calcola il LE per sistemi dinamici casuali accoppiati http://mesoscopic.mines.edu/acoustics-old/preprints/lyapunov.pdf Forse lo troverai più leggibile. È un algoritmo stabile e ben testato. Il codice è stato sviluppato dal mio co-autore, che è ora alla KU (nel caso tu voglia il codice).