Synchronisation verständlich erklärt

Viele Metronome, jedes in einem leicht anderen Tempo, werden so gekoppelt, dass jedes auf den Takt der anderen reagiert. Erst ab einer bestimmten Kopplungsstärke entsteht ein gemeinsamer Takt. Die Seite erklärt zuerst ohne Formeln, warum, und dann mit Modell und Messungen, was das Rechenbeispiel zeigt und was nicht.

Zuletzt aktualisiert: 15. November 2026

Ein Wort unbekannt? Das Begriffsverzeichnis erklärt es in wenigen Sätzen.

Viele Metronome, jedes in eigenem Tempo

Stellen Sie sich viele Metronome vor, die nebeneinander ticken. Jedes hat sein eigenes, leicht abweichendes Tempo (sein Eigentempo, so tickte es allein): eines etwas schneller, eines etwas langsamer. Wie weit ein Metronom gerade in seinem Takt ist, heißt Phase, wie die Stellung eines Zeigers auf dem Zifferblatt. Zu Beginn steht jedes an einer anderen Stelle. Das klingt wie Regen: kein gemeinsamer Rhythmus.

Nun verbinden Sie die Metronome so, dass jedes auf alle anderen hört und sich ein wenig zu deren gemeinsamem Takt hinziehen lässt. Diese Verbindung heißt Kopplung, und wie stark sie ist, lässt sich einstellen. Diese Stärke nennt die Seite K. Die Frage: Wie groß muss K sein, damit aus dem Regen ein gemeinsamer Takt wird?

Das Metronom-Bild dient nur zum Vorstellen. Die Seite beschreibt kein echtes Brett mit echten Metronomen, sondern ein Rechenmodell mit solchen Taktgebern. Man nennt sie allgemein Schwinger (Fachwort: Oszillatoren).

Unter der Schwelle: jedes bleibt bei sich

Ist K klein, gewinnen die Unterschiede im Tempo. Jedes Metronom wird zwar ein wenig gezogen, entfernt sich aber von den anderen schneller, als es herangezogen wird. Es entsteht kein gemeinsamer Takt.

Zu dieser Seite gehört ein Programm zum Ausprobieren im Browser, die Demo. Sie misst den Gleichlauf mit einer Zahl zwischen 0 und 1: 0 heißt, die Stellungen sind gleichmäßig verteilt, 1 heißt, alle stehen im Gleichschritt (Fachwort: Ordnungsparameter). Unter der Schwelle liegt sie bei etwa 0.

Über der Schwelle: der Takt verstärkt sich selbst

Dreht man K hoch, passiert ab einem bestimmten Wert etwas Neues. Zuerst machen die Metronome mit, deren Tempo dem Durchschnitt aller Tempi am nächsten liegt. Sie bilden einen gemeinsamen Takt, und der zieht stärker an den übrigen. Je mehr mitmachen, desto stärker der Zug auf alle anderen.

Dieses Aufschaukeln (Fachwort: Rückkopplung) erklärt, warum es eine Kopplungsschwelle gibt. Darunter reicht der Zug nicht, um damit anzufangen. Darüber verstärkt er sich selbst und der Gleichlauf wächst, zunächst steil und dann langsamer.

Beispiel für eine Zahl: In einer der Einstellungen der Demo streuen die Tempi so um ihren Durchschnitt, dass die meisten nahe daran liegen und nur wenige weit weg. Dort liegt die Schwelle bei K gleich 2. Darunter bleibt der Gleichlauf bei etwa 0, darüber steigt er.

Wovon hängt die Schwelle ab?

Davon, wie weit die Tempi um ihren Durchschnitt streuen. Liegen viele Metronome nahe am Durchschnitt, genügt ein kleineres K. Streuen die Tempi weit, braucht es ein größeres. Die Demo hat zwei Arten, die Tempi zu streuen (der zweite Teil nennt sie), und jede hat ihre eigene Schwelle.

Auch unter der Schwelle bleibt ein kleiner Rest

Die Zahl der Metronome ist begrenzt. Stehen die Zeiger zufällig verteilt, liegen von endlich vielen immer ein paar zufällig nahe beieinander, auch wenn keines das andere zieht. Dadurch ist der Gleichlauf unter der Schwelle nicht genau 0. Vergleichen Sie eine Münze: Nach vielen Würfen gleichen sich Kopf und Zahl ungefähr aus, aber nicht genau. Das Bild erklärt die Richtung; die Größe des Rests misst die Demo.

In den Messungen der Demo wurde der Rest kleiner, wenn mehr Schwinger beteiligt waren. Das passt dazu, dass er ungefähr mit der Wurzel aus ihrer Zahl schrumpft. Beispiel: Bei 100 statt 25 Schwingern wäre er etwa halb so groß. Bewiesen ist das nicht, und für sehr viele Schwinger ist es nicht gemessen.

Ein Test am Computer

Der Knopf „Demo ausprobieren“ am Ende dieses Teils öffnet das Programm, der Knopf daneben führt zum zweiten Teil für Fachleute. Zum Verstehen brauchen Sie die Demo nicht.

In der Demo rechnet der Computer die Schwinger nach. Die Tempi werden nach einer von zwei festen Arten gewählt, nur der Start der Stellungen ist Zufall. Dieser Zufall hängt an einer Zahl, dem Zufalls-Seed: Mit derselben Zahl entsteht derselbe Start, der Lauf lässt sich wiederholen. Drei Ergebnisse, jeweils für dieses Modell, in dem jeder mit jedem gleich stark gekoppelt ist:

• Im Modell stimmt die Formel. Über der Schwelle sagt eine Formel (sie steht im zweiten Teil) den Gleichlauf voraus. Der Computerwert liegt nahe daran. Beispiel: Bei 4000 Schwingern, gemittelt über 3 Läufe mit verschiedenem Zufalls-Seed, weicht er um weniger als 1 % vom Wert der Formel ab.
• Auch das Rechenverfahren kann täuschen. Der Computer rechnet in kleinen Zeitschritten. Bei gleich großen Zeitschritten lag ein einfacheres Verfahren (Euler) knapp über der Schwelle zu niedrig, das genauere der Demo (Heun) auf der Formel. Der Unterschied kam also vom Rechnen, nicht von den Schwingern.
• Der Rest ist messbar. Unter der Schwelle bleibt bei 4000 Schwingern ein kleiner Gleichlauf, nicht genau 0.

Der Haken: ein Modell mit gewählten Tempi

Die Tempi sind vom Computer gewählt, nicht an echten Systemen gemessen. Die Aussagen gelten nur, solange jeder mit jedem gleich stark gekoppelt ist. Hört jeder nur auf wenige Nachbarn, gilt die Formel nicht. Die Seite sagt nichts darüber, wann echte Uhren, Zellen oder Menschen im Takt laufen.

Die Demo enthält außerdem ein eigenes, selbst gewähltes Modell, in dem auch Dreiergruppen von Schwingern aufeinander wirken. Dort springt der Gleichlauf beim Hochdrehen von K plötzlich an und bleibt beim Zurückdrehen länger bestehen. Das gilt nur in diesem Modell.

Häufige Missverständnisse

• „Synchronisation heißt, alle ticken gleich.“ Nein. Die Eigentempi bleiben verschieden. In den Messungen der Demo ist der Gleichlauf bei K gleich 6 noch kleiner als 1.
• „Unter der Schwelle ist der Gleichlauf null.“ Nein. Bei begrenzter Zahl von Schwingern bleibt ein kleiner Rest.
• „An der Schwelle springt der Gleichlauf an.“ Nicht im Hauptmodell: Dort wächst er ab der Schwelle ohne Sprung an, nur steil. Ein Sprung zeigte sich nur im Modell mit Dreiergruppen.

Fragen an eine Behauptung

Hören Sie, etwas „synchronisiere von selbst“, helfen diese Fragen:

• Wie stark ist die Kopplung, und liegt sie über einer Schwelle?
• Wie verschieden sind die Eigentempi, wie weit streuen sie um den Durchschnitt? Je weiter, desto mehr Kopplung braucht es.
• Hört wirklich jeder auf jeden, oder nur auf wenige Nachbarn? Dann gilt die Formel der Demo nicht.
• Wie viele Beteiligte gibt es? Bei begrenzter Zahl bleibt immer ein Rest.

Wie ernst ist das Modell zu nehmen?

Es ist ein Rechenbeispiel mit gewählten Tempi, kein Befund über echte Systeme. Es zeigt, was in diesem Modell gilt, und liefert Fragen, keine Antworten für die Praxis.

Was Sie mitnehmen:

• Stellt man die Kopplungsstärke K hoch, wächst ab einer Schwelle ein gemeinsamer Takt aus verschiedenen Tempi, weil sich der Zug selbst verstärkt.
• Darunter bleibt nur ein kleiner Rest, der mit mehr Schwingern schrumpft.
• Die Aussagen gelten für gewählte Tempi und Kopplung jeder mit jedem.

Hinweis: Alle Demos rechnen mit synthetischen oder frei verfügbaren Beispieldaten und laufen auf einem externen Anbieter (Datenschutz). Bitte geben Sie darin keine vertraulichen oder personenbezogenen Daten ein. Die Ergebnisse gelten für die jeweilige Demo-Instanz und lassen sich nicht ungeprüft auf reale Fälle übertragen; für darauf gestützte Entscheidungen übernehme ich keine Gewähr.

Modell und Annahmen

Dieser Teil ist eigenständig lesbar. Er fasst die formale Beschreibung der Demo zusammen; die vollständige Fassung steht im Expander „📐 Mathematische Formulierung“ der Demo.

• Zustand: \(N\) Oszillatoren mit Phase \(\theta_i\) und Eigenfrequenz \(\omega_i\), jeder mit jedem gekoppelt, Kopplung \(K/N\) je Paar (die Konvention ist \(K/N\), nicht \(K\)). Der Ordnungsparameter \(r \in [0, 1]\) ist der Betrag von \(z = \frac1N\sum_j e^{i\theta_j}\): 0 bei gleichmäßig verteilten Phasen, 1 bei gleichen Phasen.
• Eigenfrequenzen: deterministische Quantile einer Verteilung mit Dichte \(g\), Lorentz mit Breite 1 oder Gauß mit Standardabweichung 1. Die äußersten Lorentz-Quantile (bei \(N = 4000\) bis ±2546) sind auf ±200 begrenzt; das ändert \(r\) um höchstens 0,001. Zufall steckt nur in den Startphasen (Zufalls-Seed, Zufallszahlengenerator SplitMix64).
• Rechenverfahren: Heun-Verfahren mit \(dt = 0{,}05\) (Zeitverlauf: 0,02). \(r\) ist das Mittel über die zweite Hälfte eines Laufs der Dauer 200 und, in der Demo, über 3 Zufalls-Seeds (der gewählte und die zwei folgenden).
• Dreier-Kopplung: ein eigenes Modell mit zusätzlichem Term \(K_2\), unten als Mittelfeld-Form angegeben. Es ist eine eigene Modellwahl nach Art der Dreier-Kopplung, keine gelesene Quelle.

Die Formeln

Das Modell in Mittelfeld-Form lautet (mit \(K_2 = 0\) ist es das gewöhnliche Kuramoto-Modell):

\[\dot\theta_i = \omega_i + K_1\,r\,\sin(\psi - \theta_i) + K_2\,r^2\,\sin\bigl(2(\psi - \theta_i)\bigr), \qquad z = \frac1N\sum_{j=1}^{N} e^{i\theta_j}.\]

Der \(K_2\)-Term ist die Mittelfeld-Form von \(\frac{K_2}{N^2}\sum_{j,k}\sin(\theta_j + \theta_k - 2\theta_i)\). Jeder Oszillator wird zur mittleren Phase \(\psi\) gezogen, umso stärker, je größer \(r\) ist: eine Rückkopplung.

Für zwei Oszillatoren ist mit \(\varphi = \theta_1 - \theta_2\) die Gleichung \(\dot\varphi = \Delta\omega - K\sin\varphi\). Sie hat genau dann einen Fixpunkt (die beiden verriegeln), wenn \(\lvert\Delta\omega\rvert \le K\) ist; sonst wandert \(\varphi\) mit der Periode \(2\pi/\sqrt{\Delta\omega^2 - K^2}\).

Für große \(N\) und \(K_2 = 0\) gilt: Eine Lösung \(r > 0\) erfüllt die Selbstkonsistenz

\[r = K\,r\int_{-\pi/2}^{\pi/2}\cos^2\theta\; g(K\,r\sin\theta)\,d\theta .\]

Die Schwelle ist \(K_c = 2/(\pi g(0))\). Bei der Lorentz-Verteilung mit Breite 1 ist \(K_c = 2\) und \(r = \sqrt{1 - 2/K}\) für \(K > K_c\); bei der Gauß-Verteilung mit Standardabweichung 1 ist \(K_c = 1{,}596\), und \(r\) folgt aus der Selbstkonsistenz.

Nur für Lorentz gibt es auf der Mannigfaltigkeit der Poisson-Kerne (Ott–Antonsen) die geschlossene Zeitlösung. Mit \(s = r^2\) und \(s^* = 1 - 2\gamma/K\) (Lorentz-Breite \(\gamma\), in der Demo \(\gamma = 1\)) lautet sie:

\[s(t) = \frac{s^*}{1 + (s^*/s_0 - 1)\,e^{-2at}}, \qquad a = \frac K2 - \gamma .\]

Zum Verfahren (Heun): \(\theta^{n+1} = \theta^n + \frac{dt}{2}(v_1 + v_2)\) mit \(v_1 = f(\theta^n)\) und \(v_2 = f(\theta^n + dt\,v_1)\). Die Eigenfrequenzen sind \(\omega_i = F^{-1}\bigl((i - \tfrac12)/N\bigr)\) mit der Verteilungsfunktion \(F\).

Woher die Schwelle kommt

Skizze (Standardwissen, Rechnung dieser Seite, keine Messung der Demo und kein Beweis): Für kleines \(r\) ist \(g(Kr\sin\theta) \approx g(0)\), und \(\int_{-\pi/2}^{\pi/2}\cos^2\theta\,d\theta = \pi/2\). Die rechte Seite der Selbstkonsistenz wird dann zu \((\pi/2)\,K\,g(0)\,r\). Eine Lösung \(r > 0\) nahe null kann nur entstehen, wenn dieser Faktor 1 erreicht, also bei \(K = 2/(\pi g(0)) = K_c\). Ob sich dort tatsächlich eine positive Lösung abzweigt und wie sie wächst, folgt aus dieser Skizze nicht; für Lorentz und Gauß belegen es die Messungen unten.

Für Lorentz mit Breite \(\gamma\) ist \(g(0) = 1/(\pi\gamma)\), also \(K_c = 2\gamma\). Die Demo rechnet nur \(\gamma = 1\). Dazu passt die Wachstumsrate \(a = K/2 - \gamma\) der Zeitlösung: Sie wechselt bei \(K = 2\gamma\) das Vorzeichen.

Weil \(g(0)\) im Nenner steht, senkt eine dichter um die Mitte liegende Frequenzverteilung die Schwelle. Das ist die Lesart von Teil 1: Tempi nahe der Mitte machen zuerst mit.

Befunde

Alle Zahlen stammen aus der Messreihe der Demo (jede steht in deren Tests als Band oder Zählung) und gelten für dieses Modell: alle mit allen, Kopplung \(K/N\), \(N\) Quantile der Verteilung, Heun mit \(dt = 0{,}05\). „Zufalls-Seeds“ sind die Seeds der Startphasen.

Frage	Befund
Zwei Oszillatoren	Δω = 1, K = 2: φ* = π/6 = 0,5236, r = cos(π/12) = 0,96593. Wandernd: Periode 2π/√(Δω² − K²) = 3,6276 (Δω = 2, K = 1); die Integration mit scipy gibt 3,6276.
Lorentz über der Schwelle	N = 4000, 3 Zufalls-Seeds, Dauer 200: r = 0,218; 0,446; 0,577; 0,707; 0,816 bei K = 2,1; 2,5; 3,0; 4,0; 6,0, Theorie √(1 − 2/K) = 0,218; 0,447; 0,577; 0,707; 0,816. Abweichung des Mittels über 3 Zufalls-Seeds höchstens 0,3 % (gemessen höchstens 0,2 %); einzelne Zufalls-Seeds weichen bis 0,5 % ab.
Wie viele Oszillatoren braucht der Befund?	Bei N = 100 und K = 2,1 liegt das Mittel der Zufalls-Seeds 1 bis 3 um 1,2 % neben der Theorie, ab N = 250 (K = 2,1; 2,5; 3,0) um höchstens 0,6 %; das gilt für diese Zufalls-Seeds. Bei 20 weiteren Dreiergruppen von Zufalls-Seeds (4 bis 63, K = 2,1) weicht das Mittel bei N = 100 im Mittel des Betrags um 2,8 % ab (bis 7,3 %), bei N = 250 um 0,69 % (bis 2,1 %).
Gauß	N = 4000, 3 Zufalls-Seeds: r = 0,426; 0,715; 0,925; 0,964 bei K = 1,7; 2,0; 3,0; 4,0; Unterschied zur Selbstkonsistenz höchstens 0,0005 (Schwelle K_c = 1,5958).
Unter der Schwelle	N = 4000, 3 Zufalls-Seeds: Lorentz r = 0,020; 0,024; 0,037 bei K = 1,0; 1,5; 1,9 (r·√N = 1,25; 1,54; 2,34); Gauß r = 0,029; 0,042 bei K = 1,2; 1,5. Unter K_c ist r kein Nullwert; der Faktor vor 1/√N wächst zur Schwelle hin.
Endlichkeit (vorgerechnet)	Lorentz, K = 1,0, Dauer 200, 5 Zufalls-Seeds: r·√N = 0,89 (N = 250), 0,97 (500), 1,10 (1000), 1,22 (2000), 1,22 (4000), 1,19 (8000), 1,25 (16 000). Ab N = 2000 liegen die 4 Mittel über 5 Zufalls-Seeds zwischen 1,19 und 1,25 (je Zufalls-Seed streut r·√N bei N = 2000 bis 16 000 zwischen 1,09 und 1,38). Das ist mit einem konstanten r·√N verträglich, aber nicht bewiesen; darunter steigt r·√N mit N.
Euler gegen Heun nahe der Schwelle	K = 2,1, N = 4000, Zufalls-Seed 1: Euler mit dt = 0,02 und Dauer 300 gibt r = 0,200, Euler mit dt = 0,05 gibt 0,167, Heun mit dt = 0,05 gibt 0,218 gegen die Theorie 0,218 (3 von 3 Zufalls-Seeds innerhalb 0,002). Die Abweichung von 8 % einer Vormessung bei K = 2,1 war damit ein Fehler des Euler-Verfahrens, kein Endlichkeitseffekt.
Noch näher an der Schwelle	Bei K = 2,05 und N = 4000 liegt r (Mittel der Zufalls-Seeds 35 bis 37) 3,4 % unter der Theorie, bei N = 16 000 0,4 % darüber; einzelne Zufalls-Seeds weichen bei N = 4000 bis 8,3 % ab. Endlichkeit zeigt sich also unter der Schwelle, bei sehr kleinem N und nahe an der Schwelle.
Schwelle auf der Rampe	N = 500, Gitter 0,25, Schrittzeit 30, 3 Zufalls-Seeds: Der erste K-Wert mit r > 0,1 ist bei Lorentz 2,00 (3 von 3; K_c = 2) und bei Gauß 1,75 (3 von 3; K_c = 1,596). Bei den Zufalls-Seeds 1 bis 10 ist er bei Lorentz in 7 von 10 Läufen 2,00 (sonst 2,25) und bei Gauß in 4 von 10 Läufen 1,75 (sonst 0,75 bis 1,50). Das ist eine Gitterangabe, kein Schwellenwert.
Dreier-Kopplung (vorgerechnet)	N = 2000, Lorentz, Schrittzeit 100, K₁-Gitter 0,05 (bis 3), 2 Zufalls-Seeds. Sprung hinauf (erster K₁ mit r > 0,2): in 10 von 10 Läufen (K₂ = 0,5; 1; 2; 3; 4) zwischen 1,90 und 2,05; bei K₂ = 0 bei 2,10. Der Rückweg bleibt über 0,2 bis K₁ = 2,10 (K₂ = 0), 1,95 (0,5), 1,80 (1), 1,30 (2), 0,65 (3), 0,15 (4), je Zufalls-Seed gleich. Die Breite der Hysterese (Sprung minus Rückweg) wächst mit K₂ in 5 von 5 Schritten. Bei K₂ = 0 folgt r der Formel √(1 − 2/K₁) ab 2,25 in allen 28 Gitterpunkten je Seed auf 0,003, hinauf und hinunter.
Dreier-Kopplung live	N = 500, Gitter 0,25, Schrittzeit 40, 3 Zufalls-Seeds: Bei K₂ = 3 springt r über 0,2 bei K₁ = 2,00, der Rückweg bleibt bis 0,75 (3 von 3 bei den Zufalls-Seeds 1 bis 3). Bei den Zufalls-Seeds 1 bis 10 springt r in 6 von 10 Läufen bei K₁ = 2,00 und in 4 von 10 schon bei 1,75 (unter K_c); der Rückweg bleibt in 10 von 10 Läufen bis 0,75.
Gegenproben der Dreier-Kopplung	Mittelfeld-Form gegen Doppelsumme (N = 6, Zufallsphasen): 1e-12; Eigenwerte der Jacobi-Matrix am Gleichlauf (N = 6, ω = 0): 0 und −(K₁ + 2K₂) auf 1e-4; Wachstumsrate der Gleichverteilung bei K₁ = 1,3 gleich 0,650 (Soll K₁/2); bei K₁ = −1, K₂ = 2 sind Gleichlauf (r > 0,999) und Gleichverteilung (r < 0,001) beide stabil, bei K₂ = 0,3 zerfällt der Gleichlauf.

Durchgerechnet: Zeitverlauf bei K = 3

Die Zeilen stammen aus der Messreihe der Demo: Lorentz, K = 3, Start auf der Ott–Antonsen-Mannigfaltigkeit (Poisson-Kern) mit r₀ = 0,05, zufällig den Frequenzen zugeordnet. Die geschlossene Lösung \(r(t)\) steht in der zweiten Spalte. Die beiden rechten Spalten geben den Faktor Simulation durch Theorie an.

t	Geschlossene Lösung r(t)	N = 4000, 20 Zufalls-Seeds	N = 40 000, 5 Zufalls-Seeds
0	0,050	–	–
2	0,133	im Mittel 1,04, Streuung 0,31, Spanne 0,60 bis 1,65	0,89 bis 1,16
4	0,312	0,57 bis 1,46	0,88 bis 1,11
6	0,501	0,71 bis 1,10	0,95 bis 1,03
8	0,565	0,92 bis 1,02	1,00 bis 1,01
10	0,576	0,98 bis 1,01	1,00

Lesart: Der Start ist weit von der Endlage entfernt, und bei \(N = 4000\) ist die Anfangsunschärfe von der Größe \(1/\sqrt{N}\) nicht klein gegen \(r_0\) (Deutung dieser Seite). Der Zeitverlauf hängt darum vom Zufalls-Seed ab. Eine frühere Vormessung mit einer einzelnen Zuordnung von Phasen zu Frequenzen hatte bei \(N = 4000\) den Faktor 1,8 bei \(t = 2\) gefunden. Über 20 Zufalls-Seeds ist der Faktor im Mittel 1,04 bei der Spanne 0,60 bis 1,65.

Grenzen

• Ein Modell: alle mit allen, gleiche Kopplung, \(N\) Quantile der Verteilung. Auf einem Netz mit wenigen Nachbarn gilt die Formel nicht. Gezeigt sind nur Lorentz und Gauß (symmetrisch, einhöckrig); zweihöckrige Verteilungen sind nicht gemessen.
• Eigenfrequenzen gewählt: Die Frequenzen sind Quantile einer Verteilung, nicht an einem echten System gemessen. Die Seite sagt nichts darüber, wann echte Systeme synchronisieren.
• Dreier-Kopplung ist ein eigenes Modell (Mittelfeld-Form, Lorentz mit Breite 1, ein \(N\) und eine Schrittzeit je Reihe). Die Hysterese ist eine Messung in diesem Modell, keine Aussage über die Literatur zur Kopplung höherer Ordnung. Dass die Ott–Antonsen-Reduktion für die zweite Harmonische nicht ohne Weiteres geschlossen bleibt, ist nicht verfolgt. Die Schwelle \(K_c = 2\) beim Hinaufdrehen folgt aus der Linearisierung der Gleichverteilung (der \(K_2\)-Term kommt dort nur in zweiter Ordnung vor), die Lage des Rückwegs nicht.
• Der Sprung hinauf bei \(K_2 \ne 0\) liegt in allen 10 gemessenen Läufen bei 1,90 bis 2,05, also bei oder knapp unter \(K_c = 2\). Ob eine Unterkritikalität dort durch Rauschen früher springt, ist bei dieser Auflösung nicht getrennt.
• Die Schwelle auf der Rampe ist eine Gitterangabe (Schrittweite 0,25 in \(K\), begrenzte Schrittzeit); nahe der Schwelle stellt sich das Gleichgewicht langsam ein.
• Endlichkeit: Das vorgerechnete Bild hat 5 Zufalls-Seeds je \(N\) und eine Dauer von 200. Ob \(r\sqrt{N}\) für \(N\) über 16 000 konstant bleibt, ist nicht gemessen.
• Rechenverfahren: Nahe der Schwelle hängt das Ergebnis am Verfahren (siehe Euler gegen Heun); ein gröberes Verfahren verschiebt die Aussage.
• Beweise: Die Schwellenformel, die Ott–Antonsen-Reduktion und die Hysterese sind hier nicht bewiesen; die Messungen prüfen sie nur für die genannten Einstellungen.
• Nicht mit echtem Mausklick im Browser getestet: Die Diagramme der Demo haben keine Auswahl; die App wurde mit Streamlits AppTest und einem Start auf Skriptfehler geprüft.

Quellen

• Quelltext, Messreihe und Tests der Demo: github.com/sebastian-hanisch/kuramoto-demo. Die Orakel-Tests rechnen die Selbstkonsistenz mit scipy.integrate.quad und brentq nach und prüfen die geschlossene Ott–Antonsen-Lösung gegen die Simulation.
• Kuramoto: Modell gekoppelter Oszillatoren mit mittlerem Feld (Stufe: aus dem Gedächtnis, Originalarbeit nicht gelesen).
• Ott und Antonsen: Reduktion auf die Mannigfaltigkeit der Poisson-Kerne (Stufe: aus dem Gedächtnis, Originalarbeit nicht gelesen).
• Literatur zur Kopplung höherer Ordnung und zu „explosiven Übergängen“: nicht gelesen und deshalb nicht zitiert.

Die Literaturhinweise stehen nicht in der README der Demo und wurden für diese Seite nicht nachgeschlagen; sie belegen hier nichts. Die Skizze zur Schwelle ist Standardwissen und eine Rechnung dieser Seite, keine Messung der Demo. Alle Zahlen stammen aus der README der Demo.