Lissajous Figuren in Gnuplot

Da nicht jeder das nötige Kleingeld für ein Oszilloskop und Funktionsgenerator hat, aber jeder gerne eine Lissajous-Figur laufen haben möchte, liefere ich hier den entsprechenden Gnuplot Code.

reset
set term gif animate optimize
set output "lissajous.gif"
n=6250

set xr [-1:1]
set yr [-1:1]

set parametric
unset border
unset xtics
unset ytics

fx(t) = sin(t)
fy(t) = sin(2.999*t)

i=0
load "animateLissajou.gp"
set output

Die Datei „animateLissajou.gp“ sieht dann so aus:

set trange [i:i+2*pi]
plot fx(t),fy(t) lc rgb 'black' notitle

i=i+2*pi*10
if (i < n) reread

Stark angelehnt an diesen Blogeintrag. Das Ergenis sieht dann so aus.

Lissajous Figur

Bootstrapping

Wer kennt das nicht: Man hat sich ein Python Skript geschrieben, um seine Daten per Bootstrap Resampling auszuwerten und stellt fest, dass das Konstrukt zur Bildung des „Samples mit Ersetzungen“

    import random
    x = [1,2,3]

    bootstrapSample = [random.choice(x) for _ in x]

einfach nicht schnell genug ist.
Aber glücklicherweise gibt es numpy!

    import numpy
    x = [1,2,3]

    bootstrapSample = list(numpy.random.choice(x, len(x)))

Das ist — zumindest in meinem Anwendungsfall — spürbar schneller. Ich werde in Zukunft also immer optimale Fehlerbalken erzeugen.

Seltsamer Attraktor

Zuvor habe ich bereits den Schmetterlingseffekt erwähnt. Um den Zusammenhang mit Chaos zu zeigen, betrachten wir folgendes Video von der Projektion in die y-z-Ebene von 13 Teilchen, die den Attraktor durchlaufen.

Alle Teilchen starten auf fast dem selben Punkt, aber nehmen sehr verschiedene Wege. Nach kurzer Zeit kann man den einzelnen Teilchen nicht mehr ansehen, dass sie fast die gleichen Anfangsbedingungen hatten.

Lorenz war Meteorologe und sein Differentialgleichungssystem \begin{align} \dot{X} &= a(Y - X) \ \dot{Y} &= X(b - Z) - Y \ \dot{Z} &= XY - cZ, \ \end{align} das dieses chaotische Verhalten zeigt, sollte die Atmosphäre modellieren.

Jetzt kann man verstehen, was es mit dem Schmetterling aus Jurassic Park auf sich hat.

Er bewegt in Peking die Flügel, und im Central Park gibt’s Regen statt Sonne.

Dr. Ian Malcolm (1993)

Sein Flügelschlag ändert den Zustand eines chaotischen Systems, dem Wetter, ein wenig und nach einiger Zeit hat das System einen grundlegend anderen Weg eingeschlagen, als ohne diesen Flügelschlag.

Dennoch sieht das Video irgendwie geordnet aus. Fast schon vorhersagbar. Seltsam.