Korrekter Weg zur Lösung der Gleichung für einfache harmonische Bewegung

Ich betrachte die Gleichung für einfache harmonische Bewegung, die ist X ¨ + ω 2 X = 0 Um dies zu lösen, habe ich drei Ansätze gesehen. Dies ist verwirrend, da ich nicht weiß, welcher Ansatz physikalisch korrekt ist oder, falls es keinen korrekten Ansatz gibt, welche physikalische Bedeutung die drei verschiedenen Ansätze haben. Meine Vermutung ist, dass der erste Ansatz leicht redundant ist, da er nur das bekannte Ergebnis (von Kosinus- und Sin-Lösungen dieser Form einer Differentialgleichung zweiter Ordnung) verwendet, das aus der Verwendung von Methoden der vollständigen Lösung abgeleitet wurde e ich ϕ , wie im zweiten und dritten Ansatz. Ich weiß jedoch nicht, was davon richtig ist oder warum wir es auf zwei verschiedene Arten betrachten können.

1. Nehmen Sie die Sinus- und Cosinus-Ergebnisse an

Cosinus und Sinus sind beide Lösungen der obigen Gleichung, also ist die vollständige Lösung eine Linearkombination der beiden X = A cos ( ω T ) + B Sünde ( ω T ) was äquivalent ist X = A cos T ( ω T + ϕ ) , was die gemeinsame Form ergibt.

2. Löse allgemeiner nach einem komplexen x.

Ich habe gesehen, dass die Gleichung auch gelöst wird, indem die homogene Differentialgleichung für gelöst wird X = A e ich ω T + B e ich ω T . Dann werden die Konstanten verwendet, um die reelle Lösung zu geben, also die Randbedingung (das X ist real) ergibt das A = B . Hier haben wir also komplexe Konstanten und erhalten das Ergebnis X = ( A + B ) cos ( ω T ) + ich ( A B ) Sünde ( ω T ) = 2 ( A ) cos ( ω T ) 2 ( A ) Sünde ( ω T ) was sich, wie zuvor mit Konstanten, die sin- und cos-Terme multiplizieren, auf reduziert X = D cos ( ω T + ϕ )

3. Der letzte Ansatz, den ich gesehen habe, besteht darin, die Differentialgleichung für zu lösen z und nur Einstellung X = ( z ) .

Auflösen der Gleichung nach z ergibt X = A e ich ω T + B e ich ω T wo ich glaube (wenn ich die Methode richtig verstanden habe), sind die Konstanten hier REAL. Da reduziert sich das noch auf z = ( A + B ) cos ( ω T ) + ich ( A B ) Sünde ( ω T ) , aber diesmal für echtes A und B und daher z eine komplexe Zahl ist, können Sie sehen, dass dies äquivalent zu ist z = D e ich ϕ e ich ω T wo die komplexe Amplitudenkomponente e ich ϕ rotiert die komplexe Zahl e ich ω T in der komplexen Ebene so, dass das Verhältnis der Konstanten der Cosinus- und Sin-Terme von z ist korrekt und wie in der angegeben z = ( A + B ) cos ( ω T ) + ich ( A B ) Sünde ( ω T ) Eine Form von z . Dies reduziert sich auf z = e ich ( ω T + ϕ ) nachgeben X = ( z ) = D cos ( ω T + ϕ ) wie vorher.

Wie bereits erwähnt, bin ich mir nicht sicher, was der „richtige“ Ansatz ist, wenn überhaupt. Aber es muss zumindest eine physikalische Bedeutung haben, dies in den verschiedenen Methoden anzugehen, und ich wäre dankbar, wenn jemand es erklären könnte.

Einige der Gedanken, die ich bisher hatte, sind:

  • Die erste Sin-and-Cosinus-Lösung ist nur eine Abkürzung unter Verwendung des bekannten Ergebnisses, das aus der formelleren Lösung abgeleitet wird e ich ϕ . Die Sin/Cos-Form gilt jedoch NUR, wenn die Variable in der Differentialgleichung REAL ist. Zum Beispiel gibt es die richtige Lösung für real X , aber wenn hatte einen Komplex z dort, von dem ich wusste, dass es komplex sein sollte, dann wäre die Annahme, dass die sin / cos-Form falsch wäre und nicht die vollständige Lösung / das vollständige Bild liefern würde - es würde nur die reale Komponente von geben z .
  • Der Unterschied zwischen dem zweiten und dem dritten Ansatz scheint darin zu bestehen, dass wir in einem Fall wissen, dass wir eine komplexe Variable haben z . In diesem Fall stellen wir tatsächlich fest, dass die Konstanten real sein können (obwohl ich denke, dass sie auch nicht real sein können. Vorausgesetzt, das A B , dies ergibt immer noch einen Komplex z da sich die komplexen Teile nicht aufheben). Und im zweiten Fall finden wir, dass die Konstanten reell sein müssen und so A = B so dass wir ein echtes bekommen X , wie es die Randbedingungen vorgeben. In Anbetracht dessen kann ich also nicht ganz erkennen, wie die obigen Ansätze 2 und 3 auf zwei verschiedene Arten dieselbe x-Form ergeben würden: Verwenden von Randbedingungen für einen Komplex A und B, sodass sich komplexe Teile aufheben, oder Zulassen, dass z komplex ist und somit A B und x zum reellen Teil dieser komplexen Lösung zu machen.

Entschuldigung für den langen Beitrag. Ich hoffe, ich habe meine Verwirrung und meine Gedanken hinreichend deutlich gemacht!

Selbst wenn x komplex wäre, würden die sin/cos-Terme funktionieren, aber mit komplexen Amplituden.
Ich habe Ihren Beitrag mit der Notation für Real- und Imaginärteile leicht umformatiert . Hoffe, das ist in Ordnung.
@RazimanT.V. Ah ja ich verstehe! Die Verwendung von sin/cos und gerade ist also im Wesentlichen äquivalent, solange Sie komplexe Konstanten zulassen, um einen Umrechnungsfaktor zwischen den beiden bereitzustellen.
Es gibt eine eindeutige Lösung für das Anfangswertproblem. Gibt Ihnen Ihre Methode eine Lösung? Dann ist es richtig. Zeitraum. Keine Fragen mehr gestellt. (Und außerdem gibt es keine Methoden, die "richtiger" sind als andere - sie geben entweder eine Lösung oder nicht.) Deshalb sind die Existenz- und Eindeutigkeitstheoreme so stark.

Antworten (4)

Ich glaube, du machst dir zu viele Sorgen. Dies ist der richtige Ansatz (ich werde etwas leichtfertig sein, also nimm diesen ersten Absatz beim ersten Lesen nicht zu ernst :)):

  1. Schritt 1: Die Bedeutung des Satzes von Picard-Lindelöf verstehen ;
  2. Schritt 2: Verstehen Sie, dass Sie nacharbeiten können, indem Sie Zustandsvariablen allen außer der Ableitung höchster Ordnung zuweisen X ¨ + ω 2 X = 0 in eine Vektorversion der Standardform u ˙ = F ( u ) vom PL-Theorem angesprochen und dass in diesem Fall die F ( u ) erfüllt die Bedingungen des PL-Theorems (es ist Lipschitz-stetig)
  3. Schritt 3: Wählen Sie Ihre bevorzugte Methode, um eine Lösung für die DE und die Randbedingungen zu finden - Tricks, die Sie in Differentialgleichungen 101 lernen, Versuch und Irrtum, Vermutungen einfüllen und sehen, was passiert ... alles! .... und dann GO FOR IT!

Okay, das ist ein bisschen leichtsinnig, aber der Punkt ist, dass Sie aus grundlegenden theoretischen Überlegungen wissen, dass es eine Lösung geben muss, und wie auch immer Sie die Gleichung lösen, wenn Sie eine Lösung finden können, die zu der Gleichung und den Randbedingungen passt, müssen Sie einfach die haben richtige und einzige Lösung , egal wie Sie sie ableiten .

Insbesondere gelten die obigen theoretischen Überlegungen, ob die Variablen reell oder komplex sind. Wenn Sie also eine Lösung mit komplexen Variablen finden und sie zu den realen Randbedingungen passen, dann muss die Lösung die gleiche sein wie die, die durch Kleben gefunden werden soll mit reeller Variablenschreibweise. In der Tat kann man die Begriffe von definieren Sünde Und cos durch die Lösungen von X ¨ + ω 2 X = 0 und sie müssen aufgrund der obigen Überlegungen zum PL-Theorem äquivalent zu komplexen Exponentiallösungen sein. Diese erzwungene Gleichwertigkeit können Sie sich dann als Grund für Ihre eigene schön formulierte Erkenntnis vorstellen, die Sie sich erarbeitet haben:

"Die Verwendung von sin / cos und gerade ist also im Wesentlichen äquivalent, solange Sie komplexe Konstanten zulassen, um einen Umrechnungsfaktor zwischen den beiden bereitzustellen."

Lassen Sie das Wort "im Wesentlichen" fallen und Sie haben alles sortiert!

Lassen Sie uns eigentlich zu Schritt 2 in meiner "augenzwinkernden" (aber insgesamt theoretisch fundierten) Antwort zurückkehren, da sie uns zeigt, wie wir all diese Ansätze vereinen und die Physik gut einbringen können. Brechen Sie die Gleichung in ein gekoppeltes Paar von Gleichungen erster Ordnung auf, indem Sie schreiben:

X ˙ = ω v ; v ˙ = ω X

und jetzt können wir die Dinge kurz und bündig als Matrixgleichung schreiben:

(1) X ˙ = ich ω X ; ich = D e F ( 0 1 1 0 )  Und  X = ( X v )

deren eindeutige Lösung die Matrixgleichung ist X = exp ( ich ω T ) X ( 0 ) . Hier exp ist die Exponentialmatrix. Beachten Sie auch, dass als reelle Koeffizientenmatrix ich 2 = ich D . Nun, Sie wissen vielleicht, dass eine absolut gute Möglichkeit, komplexe Zahlen darzustellen, die folgende ist: das Feld ( C , + , ) ist isomorph zum Kommutativkörper von Matrizen der Form:

(2) ( X j j X ) ; X , j R

zusammen mit Matrixmultiplikation und -addition. Für Matrizen dieser speziellen Form ist die Matrizenmultiplikation kommutativ (obwohl es natürlich nicht allgemein so ist) und der Isomorphismus wird durch die Bijektion gezeigt

(3) z C ( R e ( z ) ICH M ( z ) ICH M ( z ) R e ( z ) )

Also, wenn wir jetzt lassen Z sei ein 2 × 2 Matrix dieser Form, dann können wir (1) lösen, indem wir den Zustandsvektor abbilden X = ( X v ) bijektiv zum 2 × 2 Matrix Z = ( X v v X ) , Lösung der Gleichung Z ˙ = ich ω Z , dh Z ( T ) = exp ( ich ω T ) Z ( 0 ) , Wo Z ( 0 ) ist der 2 × 2 Matrix der Form (2) mit den korrekten Werten von X ( 0 ) Und v ( 0 ) die die Randbedingungen erfüllen, und dann nur die erste Spalte des Ergebnisses nehmen 2 × 2 Matrixlösung Z ( T ) zu bekommen X ( T ) .

Dies entspricht genau der komplexen Notationsmethode, die Sie verwendet haben, wie Sie hoffentlich sehen werden, wenn Sie das Obige ein wenig untersuchen. Die Phasenwinkel werden durch die Phase der kodiert 2 × 2 Matrix Z , gedacht als komplexe Zahl durch den oben beschriebenen Isomorphismus.

Außerdem gibt es hier einige schöne Physik. Betrachten Sie die quadratische Norm des Zustandsvektors X ; es ist E = 1 2 X , X = 1 2 ( X 2 + v 2 ) und das können Sie sofort aus (1) ableiten

(4) E ˙ = X , X ˙ = X T X ˙ = ω X T ich X = 0

Dies hat zwei Interpretationen. Zuerst, E ist die Gesamtenergie des Systems, aufgeteilt in potentielle Energie 1 2 X 2 und kinetisch 1 2 v 2 . Zweitens zeigt (4), dass der Zustandsvektor, geschrieben als kartesische Komponenten, dem Kreis folgt X 2 + v 2 = 2 E und tatsächlich ist diese Bewegung eine gleichmäßige kreisförmige Bewegung ω Radianten pro Zeiteinheit. Diese einfache harmonische Bewegung ist also die Bewegung einer beliebigen kartesischen Komponente einer gleichförmigen Kreisbewegung.

Sie könnten das Problem auch lösen, indem Sie mit (1) beginnen, (4) ableiten und dann die Substitution vornehmen

(5) X = 2 E cos ( θ ( T ) ) ; v = 2 E Sünde ( θ ( T ) )

was durch das Denkmalschutzgesetz bestätigt wird X 2 + v 2 = 2 E mit E ˙ = 0 . Dann ersetzen X zurück in die ursprüngliche SHM-Gleichung, um dies abzuleiten

(6) θ ( T ) = ± ω T + θ ( 0 )

Brillant! Wenn ich fragen darf, haben Sie einen mathematischen Hintergrund? Ich wünschte, ich hätte das gleiche Verständnis für diese Themen, wie Sie in Ihrem Beitrag demonstriert haben! Haben Sie Vorschläge für Kurse oder Bücher, die Sie darüber lesen können?
@21joanna12 Schön, dass es dir gefallen hat. Wo soll ich anfangen: Das ist der Hintergrund, den ich mir über lange Zeit angeeignet habe. Wenn Sie die ersten Kapitel von „Naive Lügentheorie“ von einem meiner sehr geliebten Lehrer, Prof. John Stillwell , lesen, werden Sie eine Menge Zeug dieser Richtung bekommen. Außerdem können Sie Kapitel 22 von Band 1 der Feynman-Vorlesungen genießen. Zufälligerweise steht dies direkt nach Kapitel 21 über SHM, das Sie ebenfalls lesen sollten.

Der einfachste Einblick in eine Differentialgleichung wie

X ¨ + A X ˙ = B X
ist zu beachten, dass die linke Seite Ihnen die Ableitungen der Funktion anzeigt X ( T ) muss ein Vielfaches dieser Funktion sein, da die rechte Seite proportional zu ist X ( T ) .

Was sind denn die (einfachen) Funktionen, die bei Differentiation ein Vielfaches von sich zurückgeben? Da ist nur A e λ T (mit λ Und A 0 Konstanten) oder die triviale Funktion X = Konstante. Keine andere Funktion gibt unter Differenzierung ein Vielfaches von sich selbst zurück.

Es folgt dem X muss die Form haben X = A e λ T . Wenn Sie dies in die Differentialgleichung einsetzen, wandeln Sie das Problem vom Finden der Funktion zum Finden der Konstanten um λ Und A .

Tatsächlich erhalten wir die Hilfsgleichung λ 2 + A λ = B mit dem Faktor A abbrechen. λ kann dann durch Lösen der Quadrate gefunden werden. Da es (normalerweise) zwei Wurzeln gibt λ ± die allgemeine Lösung wird sein

X ( T ) = A + e λ + T + A e λ T .
Denn die Stellung X ( T ) eine reelle Zahl ist, muss es sein, dass die Konstante A ± , die durch die Anfangsbedingungen des Problems bestimmt werden, verbinden sich zu reellen Funktionen; wenn die Wurzeln λ ± Imaginäre Teile haben, muss es sein, dass bei Verwendung des Satzes von Euler X ( T ) nimmt die Gestalt an
X ( T ) = e γ T ( B + cos ( ω T ) + B Sünde ( ω T ) )
Hier γ ist der Realteil und ω ist der Imaginärteil von λ ± = γ ± ich ω .

Das gleiche Argument gilt, wenn A = 0 : Welche Funktionen sind so, dass ihre zweiten Ableitungen ein Vielfaches von sich selbst sind? Nochmal e ± λ T aber jetzt auch Sünde ( ω T ) Und cos ( ω ) . Ersteres verwandelt sich in ein positives Vielfaches von sich selbst und die beiden letzteren in negative Vielfache von sich selbst. So können Sie die Form der Lösung aus dem Vorzeichen von spezialisieren B in der Differentialgleichung. Für reine harmonische Bewegung, B = ω > 0 .

Wie Sie erraten haben, gibt es also eine Redundanz in den Methoden. Ein wenig Einblick in das Verhalten von Funktionen unter Differenzierung, kombiniert mit einigen physikalischen Anforderungen, reicht aus, um die Form der Lösung schnell zu isolieren, wobei unbekannte Koeffizienten mit Ihrem spezifischen Problem abgeglichen werden müssen.

Da alle drei mehr oder weniger gleichwertig sind, verwenden Sie, was immer Sie für richtig halten; Einige sind jedoch je nach Situation nützlicher als andere.

Komplexe Exponentiale zu verwenden und am Ende den Realteil zu nehmen, ist nützlich, wenn Sie kompliziertere Probleme lösen, z. B. in erzwungenen einfachen harmonischen Schwingungen mit Dämpfung:

Erzwungener Oszillator

X ¨ + γ X ˙ + ω 0 2 X = F 0 M cos ( ω T )

Wir streben eine Steady-State-Lösung an. In der komplexen Ebene lautet die Bewegungsgleichung

X ~ ¨ + γ X ~ ˙ + ω 0 2 X ~ = F 0 M e ich ω T

Wir lösen, indem wir die Versuchslösung einsetzen

X ~ = A e ich ( ω T + ϕ )

Das ist eine komplexe Lösung, und wir würden bis zum Ende warten (nachdem alle Differenzierungen und Ersetzungen vorgenommen wurden), um den Realteil zu nehmen, um die physikalische Lösung zu erhalten.


Aber für einfachere Situationen wie

X ¨ + ω 2 X = 0

Es ist absolut in Ordnung, Methode 1 zu verwenden:

X = A cos ( ω T ) + B Sünde ( ω T )

Aber es muss zumindest eine physikalische Bedeutung haben, dies in den verschiedenen Methoden anzugehen, und ich wäre dankbar, wenn jemand es erklären könnte.

Es gibt viele Probleme, die mit (scheinbar) unterschiedlichen mathematischen Werkzeugen gelöst werden können. Es ist ein Fehler zu glauben, dass diese unterschiedlichen Ansätze eine unterschiedliche Bedeutung haben.

Sie versuchen, mehr als nötig in den Prozess einzulesen.

Mathematisch erhalten Sie das gleiche Ergebnis, unabhängig davon, welchen Ansatz Sie wählen. Die drei Ergebnisse sind verschiedene Formen desselben Ausdrucks.

Jede Bedeutung ergibt sich aus dem Kontext, in dem Ihr Ergebnis verwendet wird.

Bei manchen Problemen sind die Begriffe Phase, Frequenz und Amplitude sinnvoller. In manchen Kontexten ist die Phase nicht unbedingt hilfreich, aber die Form, die eine Überlagerung von zwei Funktionen ist, ist nützlicher. In manchen Situationen würde die Verwendung komplexer Formulare praktischen Notwendigkeiten einfach im Wege stehen.

Es gibt keine feste Regel, was funktioniert, und die Verwendung dieses speziellen Problems ist insofern irreführend, als echte Probleme selten so sauber funktionieren. Bei der Suche nach Näherungslösungen für einige Probleme kann es wichtig sein, die erwartete Form der Lösung in eine im Kontext des Problems sinnvolle Form zu gießen. Dies ist eine Art Kunst und wird mit Übung und Vertrautheit mit dem spezifischen Problem natürlicher.