Synchronisation schwach gekoppelter Oszillatoren
Teil 1: Theoretische Grundlagen
Seminar: Physik in der Biologie
Raphael Engesser
In der Biologie: Oszillatoren von grundlegender Bedeutung:
• Herzschlag
• Neuronen
• Parkinson
• Lotka – Volterra
• Glühwürmchen
• …
Ein Oszillator ist ein dynamisches System mit
einem beschränkten periodischen Attraktor
Biologie: immer Dissipation und Fluktuation vorhanden
=> Es müssen aktive System sein (zB van der Pol) Hamiltonsche Systeme:
• klingen ab oder
• laufen aus dem Ruder
Grenzzyklen
Energie
x
Dissipation Von außen
zugeführte Energie
• Amplitude unempfindlich gg Störungen
Von Interessere:
• nicht die Ursache einer Oszillation
• sondern Wechselwirkungen (Kopplungen) zwischen einzelnen Oszillatoren
Mögliche Effekte:
• Schwebungen
• Chaos
Synchronisation
Anpassung der Frequenzen von periodisch
schwingenden, selbständigen Systemen (Oszillatoren) aufgrund einer schwachen Wechselwirkung
• frequency entrainment
• phase locking
gleichphasig gegenphasig
Konstante Phasendifferenz keine Synchronistation
Beispiel: Millennium Bridge in London
Synchronisation in der Biologie
• Herz
• Neuronen
• Glühwürmchen
• Tausendfüssler
• Grillen
• …
Entdeckung durch Christian Huygens (1629 – 1695)
Arten von Kopplungen:
a) Unidirektionale Kopplung
Bsp: getriebener linearer Oszillator Jahreszyklus der Bäume
b) Bidirektionale Kopplung
Bsp: Gekoppeltes Pendel (siehe AP I)
Kopplung von linearen Oszillatoren:
Beispiel:
Gekoppelte Federpendel (lineare Näherung)
Allg. Lösung: Überlagerung der Normalschwingungen Фgleich und Фgegen
X1(t) = Фgleich + Фgegen X2(t) = Фgleich - Фgegen
• Schwebungen
• Maxima versetzt
Kopplung von nichtlinearen Oszillatoren
Beispiel:
Van-der-Pol Oszillator
• periodisches Störsignal
• unidirektionale Kopplung
) sin(
) 1
(
12 2 2 12
2 1
t x
x x
x
x x
Störsignal
Van-der-Pol ohne Störsignal
mit μ = 3
(a) ε = 0, d.h. ohne Kopplung (b) ε = 0.24
Synchronisation eines periodisch getriebenen van-der-Pols
Das ganze bisschen mathematischer:
• Ein Oszillator ist ein dynamisches System
• Mit einem beschränktem periodischem Attraktor
• Periode T>0: kleinstes T für das gilt
m
f x x x ( ),
m
Phasenbeschreibung
• Beschreibung eines Oszillators durch nur eine Variable
• definiere Transformation
• Θ bildet Lösungen x(t) € R auf Ф(t) € S1 ab
• Entspricht Parametrisierung des Grenzzyklus
S
1:
(x(t))
Eigenschaften von Φ(t):
• Koordinate entlang des Grenzzyklus
• steigt monoton an
• bei einem Umlauf um den Grenzzyklus um 2π
• gleichförmige Bewegung gemäß:
T dt
d 2
0
Phasenbeschreibung sinnvoll da:
• Störungen wirken sich nur auf Phase aus
• Grenzzyklus: Amplitude ist stabil
• System nur eindimensional
) ,
(
) ,
(
1 2
2 2
2 2
2 1
1 1
1 1
x x
p x
x f
x x
p x
x f
ε )
dt ( d
ε )
dt ( d
Betrachte zwei miteinander gekoppelte Oszillatoren:
Frage: Wie sieht Phasenbeschreibung aus?
Wegen Störungen muss man die Phase auch auf einer Umgebung des Attraktors definieren
1 1 1
( )
dt
d x
1 2
1 1 1
1
( ) ( , )
x
x f x x
dt d
Ungestörter Oszillator auf Umgebung des Attraktors
Kettenregel
) ,
(
) ,
) ( (
2 1
1 1
1
2 1
1 1
1 1
1 1
1 1
x x
p
x x
p x
f x x
ε
ε )
dt (
d
( ), ( )
) ,
(
1 2 1 1 1 1 2 21
p x x
h
mit Kopplung
definiere 2π-periodische Funktionen h1,2
) ,
(
) ,
(
2 1
2 2
2 2
2 1
1 1
1 1
x x p
x x f
x x p x
x f
ε )
dt ( d
ε )
dt ( d
) ,
(
) ,
(
1 2 2 2 2
2 1 1 1 1
dt εh d
dt εh d
Dynamische System:
lässt sich überführen in:
betrachte Störung auf dem Grenzzyklus:
) ,
(
1 21
1
1
εh
dt
d
1
1
12(
1
2) dt H
d
) (
) (
1 2
21 2
(2) 2
2 1
12 1
(1) 1
H H
) (
) (
) (
und
) (
12 21
1
2
H H
H mit
H
(2) – (1) ergibt Phasendifferenz ΔФ = Ф2 - Ф1 man erhält neue Koordinate ΔФ:
Fixpunkte ΔФ´ = 0:
) (
) (
0
H
H
Annahme: identischen Oszillatoren und WW )
( )
H(
und
1
0
2
H
Stabilitätsanalyse:
• System: ΔФ´=εH(ΔФ)
• Fixpunkt ΔФ*
• Stabil wenn H´(ΔФ
*) < 0
Beispiel für H(Δφ) und H
12(Δφ ) bzw. H
21(Δφ)
Fixpunkt bei ΔФ = 0 stabil - gleichphasig
Adler Gleichung
Zur Veranschaulichung:
wähle für H(ΔФ) = sin(ΔФ)
Adlergleichung:
)
(
H
) sin(
)
(t
Adlergleichung – Lösungen für verschiedene ε
) sin(
) (
H
„Washboard“ - Potential
K x
V d H
dV
( ) , v gl. klass. Mech : ( )
) cos(
) )(
( )
( )
)(
( )
(
0
H x dx
V
Gleichung für Phasendifferenz
Rechte Seite als Potential:
V(ΔФ) ergibt sich mit H(ΔФ) = sin(ΔФ) als:
ΔФ
Untersuchung der Potentialgleichung:
) cos(
) )(
( )
( V
ΔФ
ΔФ