Die Jacobi-Matrix verstehen

Am Beispiel eines zweidimensionalen Systems, beschrieben durch die folgenden ODEs:

D X 1 D T = F 1 ( X 1 , X 2 ) D X 2 D T = F 2 ( X 1 , X 2 )

Die Jacobi-Matrix JM ist dann gegeben durch:

J M = ( F 1 X 1 F 1 X 2 F 2 X 1 F 2 X 2 )

Jetzt Zitat aus Scholarpedia :

Die Stabilität typischer Gleichgewichte glatter ODEs wird durch das Vorzeichen des Realteils der Eigenwerte der Jacobi-Matrix bestimmt. Diese Eigenwerte werden oft als „Eigenwerte des Gleichgewichts“ bezeichnet. Die Jacobi-Matrix eines Systems glatter ODEs ist die Matrix der partiellen Ableitungen der rechten Seite nach Zustandsvariablen, wobei alle Ableitungen am Gleichgewichtspunkt x=xe ausgewertet werden. Seine Eigenwerte bestimmen die linearen Stabilitätseigenschaften des Gleichgewichts.

Ein Gleichgewicht ist asymptotisch stabil, wenn alle Eigenwerte negative Realteile haben; es ist instabil, wenn mindestens ein Eigenwert einen positiven Realteil hat.

  • Gibt es eine intuitive Art zu verstehen, warum das Vorzeichen der Eigenwerte von JM den Stabilitätszustand des Systems impliziert?
Hinweis: Denken Sie an den eindimensionalen Fall und wie der Graph von F betrachtet den Gleichgewichtspunkt, wenn die Ableitung (oder 1D Jacobi) F ' positiv ist und wie es aussieht, wenn es negativ ist. Die mehrdimensionale Aussage ist ähnlich.
Die Intuition ist, dass der Realteil der Eigenwerte einem exponentiellen Wachstum (Abfall) für positives (negatives) Vorzeichen entspricht, was zu einem instabilen (stabilen) Verhalten führt.

Antworten (2)

Betrachten wir ein eindimensionales Beispiel:

Gleichgewichtsdiagramm

Erinnere dich daran F ( X ) = X ˙ , So F kodiert die zeitliche Entwicklung von X . Wenn F < 0 , Dann X wird sich nach links bewegen. Wenn F > 0 , Dann X wird sich nach rechts bewegen. Wenn F = 0 , X wird sich überhaupt nicht bewegen, deshalb F ( X 0 ) = 0 ist die Gleichgewichtsbedingung.

Sehen Sie sich nun an, was passiert, wenn Sie die Gleichgewichte stören X ich leicht nach rechts: Wenn F ' ( X ich ) < 0 , dann, in einer kleinen Umgebung von X ich , F ( X ) < 0 für X > X 0 in dieser Umgebung, dh rechts von einem Gleichgewicht mit F ' ( X ich ) < 0 , X bewegt sich nach links - kehrt zum Gleichgewicht zurück! Umgekehrt, wenn F ' ( X ich ) > 0 , Dann F ( X ) > 0 für alle X in der kleinen Umgebung nach rechts, d. h. nach einem kleinen Schubs nach rechts X wird sich noch weiter nach rechts bewegen und das Gleichgewicht verlassen!

Die gleiche Argumentationslinie kann auf Störungen nach links angewendet werden, was insgesamt zeigt, dass, wenn F ' ( X ich ) < 0 , dann eine kleine Störung herum X ich wird immer wieder einziehen X ich , und wenn F ' ( X ich ) > 0 , dann wird eine kleine Störung immer größer und entfernt sich von X ich .

Der mehrdimensionale Fall ist weniger anschaulich, aber die Intuition ist die gleiche – negative Eigenwerte der Jacobi-Kurve bedeuten, dass die Zeitentwicklung zurück ins Gleichgewicht zeigt, positive Eigenwerte, dass sie vom Gleichgewicht wegzeigt .

Danke jetzt ganz klar. Um es zusammenzufassen: Wenn wir nach stabilen Gleichgewichtspunkten suchen, müssen wir sie finden X 0 'S ( F ( X 0 ) = 0 ), für die die 2. Ableitung von X rechts und links ist beides negativ, oder? (bedeutet also X 0 sind das Maxima unserer Flugbahn? wie kann ein Maximum der Stabilitätspunkt sein)
@user929304: X hat keine Ableitungen "rechts und links" ("links" und "rechts", wie ich sie verwendet habe, beziehen sich auf "kleiner". X “ und „größer X "). X hat nur Zeitableitungen. Ein stabiles Gleichgewicht ist eines, wo F ' ( X 0 ) < 0 , so dass F ist positiv nach links vom Gleichgewicht und negativ nach rechts - weil es sich dann nach links bewegt, wenn es nach rechts gestoßen wird, und nach rechts, wenn es nach links gestoßen wird.
  1. F = k X ist dagegen stabil F = k X ist instabil.
  2. Normalerweise können Sie eine Matrix umschreiben A als A = P D P 1 Wo P ist eine Matrix von Eigenvektoren und D ist eine Diagonalmatrix von Eigenwerten.
  3. Wenn F = A X , dann durch das obige, ( P 1 F ) = D ( P 1 X ) . Jetzt hast du N unabhängige Gleichungen genau von der Form F = k X oder F = k X . Wenn einer von ihnen ist wie F = k X , wird die Lösung bei einer winzigen Störung ins Unendliche explodieren.

Zum Beispiel: ( 2 2 1 1 ) ( 1 0 0 1 ) ( 2 2 1 1 ) 1 = ( 0 2 1 2 0 ) = A

Nur vom Anschauen der P D P 1 Diagonalisierung, kann ich sagen - wenn F = A X , dann irgendwelche X in die Richtung gehen ( 2 , 1 ) wird dazu neigen, vom Ursprung weggeschoben zu werden, und alle X in die Richtung gehen ( 2 , 1 ) wird tendenziell zum Ursprung zurückgedrängt. Du solltest verstehen P 1 X als Basiswechsel in ein Koordinatensystem, in dem die Achsen einige sind N willkürliche Richtungen, die als Kraftlinien wirken.

Hier ist eine Darstellung dieses Systems F = ( 0 2 1 2 0 ) X .

vom Ursprung wegdrängen

Mir scheint, dass Sie in Ihrer Antwort davon ausgehen J M kann immer diagonalisiert werden. Es ist auch in falsch C ...Eigentlich hast du "normalerweise" geschrieben. Was passiert, wenn dies nicht der Fall ist?
@ValterMoretti Fair Point, könnten Sie dies allgemein erläutern? (dh ob es diagonalisiert werden kann oder nicht). Es wäre auf jeden Fall schön, Ihre Meinung zu diesem Beitrag zu erfahren.
@ user929304 Sie sollten einfach (a) die Nicht- C -diagonalisierbares Real 2 × 2 Matrizen A , (b) Untersuchen Sie das Verhalten der Lösungen von D X / D T = A X um den kritischen Punkt X = 0 und schließlich (c) Ausnutzen eines (nicht trivialen) Arguments, das das Verhalten der Lösung Ihres dynamischen Systems mit dem des zugehörigen linearisierten Systems in Beziehung setzt. Letzteres wird durch Auswahl definiert A = J M X 0 am kritischen Punkt berechnet X 0 so dass F ich ( X 0 ) = 0 , ich = 1 , 2 . Natürlich müssen die Koordinaten zentriert werden X 0 .
@ValterMoretti Ich konnte keine Intuition dahinter erkennen! Ich denke, Sie haben immer noch die Eigenwerte von A negativ sein (außer bei den schlimmsten gestellten Problemen, mit denen Sie immer noch eine charakteristische Gleichung haben N Lösungen zählen Multiplizität), aber es sieht aus wie in der X = A X Fall können Sie Ihren Eigenwert negativ haben, aber Ihre Lösung explodiert wie unendlich X Sünde ( X ) ! Aber all das ist wirklich weit von einem intuitiven Verständnis entfernt, also denke ich, dass es nicht in die Antwort gehört. Ich würde es lieber einfach unter den Teppich kehren.
Tatsächlich wollte ich keine Intuition geben. Was ich sagen wollte ist, dass der linearisierte Fall das gleiche qualitative Verhalten hat wie der nicht-linearisierte Fall. Dies ist ein (schwieriger) Satz. Sobald Sie es jedoch wissen, können Sie sich nur auf den linearen Fall konzentrieren und die verschiedenen Fälle klassifizieren, da in diesem Fall explizite Lösungen für 2x2 reelle Matrizen aufgeschrieben werden können. Es gibt nur wenige mögliche Typen. Siehe zum Beispiel en.wikipedia.org/wiki/Phase_plane