Differentialgleichungen · Phasenporträt

Phasenporträt

Ein Phasenporträt stellt jede Lösung eines Differentialgleichungssystems als Kurve in der Ebene seiner Variablen dar. Für ein lineares System x′ = Ax liefern Spur und Determinante von A die Eigenwerte. Solange diese nicht zusammenfallen, bestimmen sie den Typ des Bildes.

Phasenporträts in dieser Visualisierung

Ein Phasenporträt zeichnet die Lösungen eines Differentialgleichungssystems als Kurven in der Ebene seiner Variablen, eine durch jeden Startpunkt. Diese Seite zeigt das lineare System x′=Axx' = Ax für eine 2×22 \times 2-Matrix AA. Bewegst du den Punkt (τ,Δ)(\tau, \Delta) im Spur-Determinanten-Diagramm, wechselt das Bild zwischen Sätteln, Knoten, Foki und Zentren. Die Eigenwertebene daneben zeigt die beiden Eigenwerte hinter jedem Bild.

Das System, seine Spur und seine Determinante
(xy)′=(abcd)(xy),τ=a+d,Δ=ad−bc\begin{pmatrix} x \\ y \end{pmatrix}' = \begin{pmatrix} a & b \\ c & d \end{pmatrix} \begin{pmatrix} x \\ y \end{pmatrix}, \qquad \tau = a + d, \quad \Delta = ad - bc

Jede Kurve folgt dem Richtungsfeld. Dessen Pfeile sind gleich lang und zeigen nur die Richtung. Die Geschwindigkeit erkennt man an den Markierungen in gleichen Zeitabständen: Wo die Bewegung langsam ist, liegen sie dicht zusammen. Manche Lösungen verlaufen auf einer Geraden durch den Ursprung. Gilt Av=λvAv = \lambda v, dann ist x(t)=eλtvx(t) = e^{\lambda t} v eine Lösung, denn ihre Ableitung λeλtv\lambda e^{\lambda t} v ist genau das Ergebnis der Anwendung von AA auf sie. Die Geraden sind also die Eigenrichtungen von AA, und λ\lambda gibt die Wachstumsrate entlang jeder an: Bei λ>0\lambda > 0 geht die Bewegung nach außen, bei λ<0\lambda < 0 nach innen und bei λ=0\lambda = 0 steht sie still.

Die Eigenwerte lösen λ2−τλ+Δ=0\lambda^2 - \tau\lambda + \Delta = 0. Daher ist τ\tau ihre Summe, Δ\Delta ihr Produkt und τ2−4Δ\tau^2 - 4\Delta die Diskriminante. Jede Grenze im Diagramm entspricht einem Ereignis in der Eigenwertebene. Unterhalb von Δ=0\Delta = 0 haben die Eigenwerte entgegengesetzte Vorzeichen: Das Bild ist ein Sattel. Auf dieser Linie ist einer der Eigenwerte 0, und eine ganze Gerade aus Punkten bleibt in Ruhe. Auf der Parabel τ2=4Δ\tau^2 = 4\Delta treffen sich die Eigenwerte; darüber bilden sie ein komplex konjugiertes Paar. Bei Δ>0\Delta > 0 trennt das Vorzeichen von τ\tau Abklingen und Anwachsen. Auf τ=0\tau = 0 liegt das Paar auf der imaginären Achse. Somit legen τ\tau und Δ\Delta die Eigenwerte und den Bildtyp fest — außer bei doppelten Eigenwerten, wo der Typ auch davon abhängt, ob AA ein Vielfaches der Einheitsmatrix ist.

Für komplexe Eigenwerte α±βi\alpha \pm \beta i mit β>0\beta > 0 umläuft jede von null verschiedene Lösung den Ursprung einmal in der Zeit 2π/β2\pi/\beta und wird dabei mit e2πα/βe^{2\pi\alpha/\beta} skaliert. Deshalb liegen Markierungen im Abstand eines Umlaufs auf demselben Strahl vom Ursprung, jeweils mit demselben Faktor skaliert. Bei α=0\alpha = 0 ist der Faktor 1, und jede von null verschiedene Bahn schließt sich zu einer Ellipse: ein Zentrum. Mit Δ=1\Delta = 1 und A=(01−1τ)A = \begin{pmatrix} 0 & 1 \\ -1 & \tau \end{pmatrix} ist das System der gedämpfte Oszillator x′′+2ζx′+x=0x'' + 2\zeta x' + x = 0 mit ζ=−τ/2\zeta = -\tau/2.

Beim Überqueren der Parabel schließen sich die beiden Geraden wie eine Schere zu einer und verschwinden dann, während das Bild zur Spirale wird. Abseits des Scheitels hat eine typische Matrix auf der Parabel nur diese eine Gerade. Jede andere Lösung nähert sich dem Ursprung tangential zu ihr oder verlässt ihn so. Bei (τ,Δ)=(−2,1)(\tau, \Delta) = (-2, 1) hat die Matrix −I-I dieselben Eigenwerte, doch jede Richtung ist geradlinig. Der Basiswechsel A↦PAP−1A \mapsto PAP^{-1} verformt das Bild, erhält aber τ\tau, Δ\Delta, die Eigenwerte und den Typ. Ein Vielfaches der Einheitsmatrix bleibt unverändert.

Worauf du achten kannst

  • Die beiden Geraden eines Knotens stehen meist nicht senkrecht aufeinander. Eine Lösung abseits der Geraden nähert sich dem Ursprung entlang der langsameren oder verlässt ihn entlang dieser; weit draußen biegt sie zur schnelleren ab.
  • Eine blaue Kurve ist die gesamte Lösung durch einen Startpunkt. Der gelbe Pfad ist der Abschnitt seit dem Setzen dieses Punktes. Nach einer Änderung von AA zeigt er daher die neue Lösung, keine Aufzeichnung der alten.

Bedienung

  • Drücke auf eine freie Stelle der Phasenebene, um einen Startpunkt hinzuzufügen. Nahe einer geradlinigen Lösung rastet er auf der Geraden ein.

Verwandte Themen

Weiterführende Literatur: Wikipedia: Phase portrait (Englisch).