Synchronisierung von Lorenz-Systemen

vorhergehende Artikel in: Chaos Links Numerik
27.09.2020

Nachdem ich das Buch und die Videos von Steven Strogatz durch hatte, habe ich ein Experiment nachvollziehen wollen, zu dem die Idee von einem seiner Studenten gekommen war:

Ich habe zunächst das traditionelle Lorenz-System mit

\begin{matrix}\dot x&=&\sigma (y-x)\hfill\cr \dot y&=&x(\rho -z)-y\hfill\cr \dot z&=&xy-\beta z\hfill\end{matrix}
benutzt. Dieses System löste ich numerisch mit je zwei unterschiedlichen Startbedingungen für eine gewisse Zeit bevor ich ein der beiden Systeme mit einer der Zustandsvariablen des anderen Systems beaufschlagte. Ich benutzte dazu den Parametersatz

\sigma=10\;\; \rho=28\;\; \beta=8\over 3
Am Beispiel der Übertragung der Zustandsvariable x sieht das wie folgt aus:

Driver:

\begin{matrix}\dot x_1&=&\sigma (y_1-x_1)\hfill\cr \dot y_1&=&x_1(\rho -z_1)-y\hfill\cr \dot z_1&=&x_1y_1-\beta z_1\hfill\end{matrix}

Driven:

\begin{matrix}\dot x_2&=&\sigma (y_2-x_1)\hfill\cr \dot y_2&=&x_1(\rho -z_2)-y_2\hfill\cr \dot z_2&=&x_1y_2-\beta z_2\hfill\end{matrix}

Das Ergebnis sieht man hier: in der ersten Zeile sind beide Trajektorien zum Zeitpunkt unmittelbar bevor die Kopplung erfolgt dargestellt. Man erkennt deutlich, dass die Zustände und damit die aktuelle Position im Zustandsraum beider Systeme stark voneinander abweicht. In der zweiten Zeile ist das Verhalten beider Systeme um diesen Zeitraum herum dargestellt. Man erkennt noch Unterschiede in Lage und Form der Trajektorie - schaut man sich jedoch die Position am Ende des dargestellten Ausschnitts an, kann man bereits erahnen, dass die Unterschiede immer kleiner werden. In der letzten Zeile schließlich erkennt man keine Unterschiede mehr zwischen beiden Trajektorien - Das zweite System wurde mit dem ersten nur durch die Übertragung einer Statusvariablen synchronisiert!

Screenshot Untersuchung des Verhaltens zweier Lorenz-Systeme bei Übertragung einer Statusvariablen (x)

Bei Übertragung der Statusvariablen y erkennt man exakt dasselbe Verhalten:

Driver:

\begin{matrix}\dot x_1&=&\sigma (y_1-x_1)\hfill\cr \dot y_1&=&x_1(\rho -z_1)-y\hfill\cr \dot z_1&=&x_1y_1-\beta z_1\hfill\end{matrix}

Driven:

\begin{matrix}\dot x_2&=&\sigma (y_1-x_2)\hfill\cr \dot y_2&=&x_2(\rho -z_2)-y_1\hfill\cr \dot z_2&=&x_2y_1-\beta z_2\hfill\end{matrix}

Screenshot Untersuchung des Verhaltens zweier Lorenz-Systeme bei Übertragung einer Statusvariablen (y)

Wie Strogatz et. al. bereits in ihren ursprünglichen Veröffentlichungen feststellten, gelingt die Synchronisation nicht mittels der Übertragung der Statusvariablen z: Man erreicht zwar eine qualitative Annäherung der Trajektorien, allerdings weichen sie qualitativ (man beachte die Skalenbeschriftungen!) stark voneinander ab.

Driver:

\begin{matrix}\dot x_1&=&\sigma (y_1-x_1)\hfill\cr \dot y_1&=&x_1(\rho -z_1)-y\hfill\cr \dot z_1&=&x_1y_1-\beta z_1\hfill\end{matrix}

Driven:

\begin{matrix}\dot x_2&=&\sigma (y_2-x_2)\hfill\cr \dot y_2&=&x_2(\rho -z_1)-y_2\hfill\cr \dot z_2&=&x_2y_2-\beta z_1\hfill\end{matrix}

Screenshot Untersuchung des Verhaltens zweier Lorenz-Systeme bei Übertragung einer Statusvariablen (z)

Ich habe die Untersuchung für verschiedene Parameter wiederholt:

\beta \in \{0.9,8.0/3.0,3.4\}

Die Ergebnisse sind im folgenden über ihre Thumbnails zu erreichen. Interessant war das Ergebnis beim Wert 0.9 und Synchronisierung über z: hier zeigte sich, dass die Synchronisation über z bei entsprechender Wahl der Parameter des Systems ebenfalls möglich ist:

Screenshot Untersuchung des Verhaltens zweier Lorenz-Systeme bei Übertragung einer Statusvariablen (z) und Parametrierung mit β=0.9

Artikel, die hierher verlinken

Synchronisierung von Lorenz-Systemen II

09.10.2020

Ich habe in einem vorhergehenden Artikel ein Paper zur Synchronisierung chaotischer Systeme nachvollzogen. Dort hatte ich gezeigt, dass - anders als im ursprünglichen Paper - eine Synchronisierung zweier gleich parametrierter Lorenz-Systeme bei geeigneter Parameterwahl auch über die Zustandsvariable z möglich ist.

Alle Artikel rss Wochenübersicht Monatsübersicht Github Repositories Gitlab Repositories Mastodon Über mich home xmpp


Vor 5 Jahren hier im Blog

  • Virtuelles Netzwerklabor

    11.01.2017

    Nachdem ich hin und wieder vor der Herausforderung stehe, Anwendungen unter realen Netzwerkbedingungen zu testen, habe ich bereits vor längerer Zeit begonnen, ein virtuelles Netzwerklabor aufzubauen...

    Weiterlesen...

Neueste Artikel

  • Stratum-1-NTP-Server Links

    Ich habe meinen eigenen Stratum-1-NTP-Server mittels eines GPS-Empfängers, einer Adapterschaltung und einem Raspi gebaut. Hier fasse ich einige nützliche Links zu diesem Themengebiet zusammen

    Weiterlesen...
  • EBCMS threaded und mit mehr Markdown-Unterstützung

    Mein eigener Static Site Generator hat in den letzten Wochen einige größere Umbauten erfahren

    Weiterlesen...
  • Fork der BeanShell wegen Trojan Source

    Es gibt inzwischen einen von mir erstellten Fork des originalen Repository, in dem ich die Komponente zur Darstellung der Konsole gegen die ausgetauscht habe, die in der sQLshell in den Plugins MDIJavaEditor und MDISqlEditor zum Einsatz kommt - dadurch wird wenigstens durch das Syntax-Highlighting auf problematische Stellen im Code hingewiesen.

    Weiterlesen...

Manche nennen es Blog, manche Web-Seite - ich schreibe hier hin und wieder über meine Erlebnisse, Rückschläge und Erleuchtungen bei meinen Hobbies.

Wer daran teilhaben und eventuell sogar davon profitieren möchte, muß damit leben, daß ich hin und wieder kleine Ausflüge in Bereiche mache, die nichts mit IT, Administration oder Softwareentwicklung zu tun haben.

Ich wünsche allen Lesern viel Spaß und hin und wieder einen kleinen AHA!-Effekt...

PS: Meine öffentlichen GitHub-Repositories findet man hier.