Periode auf der Phasenebene (kleine Schwingungen)

Ich habe diese Formel, um die Periode einer Bewegung im Phasenraum (in diesem Fall Plan) entlang einer Phasenkurve zu berechnen.

T ( E ) = X 1 X 2 D X 2 ( E U ( X ) )
Wo E ist die Gesamtenergie des Systems, U ( X ) ist die potentielle Energie und X 1 Und X 2 sind die Punkte meiner Phasenkurve, die die abfangen X -Achse des Phasenplans (Inversionspunkte). Nun soll ich die Periode der kleinen Schwingungen in der Umgebung des Punktes berechnen X 0 , das ist der Punkt des potentiellen Energieminimums. Also muss ich im Grunde finden:
T 0 = lim E E 0 T ( E )
Wo E 0 = U ( X 0 ) .

Mein Versuch: Ich habe Taylor verwendet, um U(x) in einer Umgebung von auszuwerten X 0 . Also, wenn man bedenkt, dass E dazu neigt U ( X 0 ) Ich habe:

T ( E ) = X 1 X 2 D X 2 ( U ( X 0 ) U ( X 0 ) U ' ( X 0 ) ( X X 0 ) U ( X 0 ) ( X X 0 ) 2 2
und da U ' ( X 0 ) = 0 (Weil X 0 ist ein Minimum für U ( X ) )
T ( E ) = X 1 X 2 D X U ( X 0 ) ( X X 0 ) 2
Also an dieser Stelle bin ich verwirrt, weil U ( X 0 ) ist sicher positiv, und das lässt mich mit einer negativen Zahl unter einer Quadratwurzel zurück. Ich bin mir sicher, dass ich irgendwo einen Satz angewendet habe, wo ich es nicht hätte tun sollen. Ich habe das Gefühl, dass es ein mathematisches Problem ist.

Mein Buch (Arnold) sagt, dass die Lösung ist 2 π / U ( X 0 ) aber ich bin hier total hängengeblieben. Kannst du mir helfen? Wo ist mein Fehler?

PS Entschuldigung für meinen Mangel an Formalismus, wie das Weglassen der Little-O-Notation usw., aber ich denke, Sie haben trotzdem meinen Punkt verstanden :)

Antworten (1)

Ihr Ausdruck ist für die halbe Periode, also verdoppeln Sie ihn, um zurück zu oszillieren. Sie haben Recht, das Potential soll konkav sein, also positive zweite Ableitung, aber Sie haben nicht berücksichtigt, dass, wenn die Energie des Systems unten ist, keine kinetische Energie, es sich nicht bewegt: Sie brauchen also wirklich eine kleine Menge von Energie oben E 0 , und wie jedes Pendel zeigt, spielt es keine Rolle, wie viel , solange es nicht Null ist! So nimm E = U ( X 0 ) + ϵ / 2 , den Atem anhalten, dass die Größe von ϵ ist irrelevant - Spoiler-Alarm.

Der Rest deiner Bewertung ist in Ordnung,

T ( E ) = X 1 X 2 D X 2 ( U ( X 0 ) + ϵ / 2 U ( X 0 ) U ' ( X 0 ) ( X X 0 ) U ( X 0 ) ( X X 0 ) 2 2
somit
T ( E ) = X 1 X 2 D X ϵ U ( X 0 ) ( X X 0 ) 2 2 = X 1 X 2 D X / ϵ 1 U ( X 0 ) ( X X 0 ) 2 2 ϵ

Nun sehen Sie, dass die Integration variabel ist und Umkehrpunkte verschoben werden können X 0 , Und ϵ in die Normalisierung der neuen Variablen absorbiert werden und verschwinden, solange sie nicht 0 ist. Außerdem kann die neu definierte Dummy-Variable durch ähnliches Absorbieren noch neu definiert werden U ( 0 ) darin, die nicht verschwindet, und dann weiter neu definiert als = Sünde θ , mit der offensichtlichen Übersetzung der Umkehrpunkte der Halbschwingung:

T ( E ) = 1 U ( 0 ) π / 2 π / 2 cos θ D θ 1 Sünde 2 θ = π U ( 0 )
für die Halbzeit und das Doppelte für die Vollzeit.

Danke Kosmas! Jetzt kann ich endlich verstehen, wo mein Fehler lag! Du hast Recht, ich habe einfach E=U(x_0) gesetzt, nicht berücksichtigt, dass in diesem Fall TENDING TO und EQUAL TO zwei völlig verschiedene Dinge sind! Vielen Dank, du hast alles klarer gemacht!
Wie sind die Extrema des Integrals von T ( E ) Konstante? sie hängen eindeutig davon ab E
@Pietro Nein, tun sie nicht! Nehmen Sie die Skalierung explizit vor. Wenn Sie das Verschwinden von verfolgt ϵ , ein Proxy für E , in der Skalierung sehen Sie, dass die skalierten dies nicht tun.
@Pierre Man sollte an die Inversionspunkte denken X 1 Und X 2 als Phasenpunkte ( X 0 ± ϵ 2 , 0 ) , was dem Zeitpunkt entspricht, an dem der gestörte Oszillator auf seiner Spitze und seinem Tiefpunkt ist. Dann wenden wir die Änderung der Variablen an Sünde θ = ( X X 0 ) 2 ϵ . Wir bekommen dann die ± π 2 als Integrationsgrenzen.