Verallgemeinerte Koordinaten finden, wenn der Satz über implizite Funktionen fehlschlägt

Gegeben einige Koordinaten ( x 1 , , x N ) und h holonomen Nebenbedingungen sollte es immer möglich sein, die Koordinaten auf diese zu reduzieren n = N h verallgemeinerte Koordinaten ( q 1 , , q n ) . Dies wird durch den Satz über implizite Funktionen gewährleistet (seine Anwendung zur Reduktion von Koordinaten wird hier kurz erwähnt ) .

Für eine Einschränkung F ( x , j ) = 0 In Bezug auf nur zwei Koordinaten besagt der implizite Funktionssatz, dass eine Funktion g warum j = g ( x ) existiert in einer Umgebung eines Punktes ( x 0 , j 0 ) die Bedingung nur erfüllen , wenn

  1. F ist stetig differenzierbar und
  2. F j ( x 0 , j 0 ) 0 .

Dies ist jedoch nicht immer der Fall.

Beispiel 1: Einheitskreis

Wie Emilio Pisanty darauf hingewiesen hat, ist ein einfaches Beispiel ein Teilchen, das sich auf dem Einheitskreis in bewegen muss ( x , j ) Flugzeug: F ( x , j ) = x 2 + j 2 1 = 0 mit Ableitung F j ( x , j ) = 2 j . Beim Versuch zu eliminieren j als Koordinate gibt es zwei lokale Lösungen: entweder j = 1 x 2 oder j = 1 x 2 , man kann also die Bewegung des Teilchens lokal leicht mit der Lagrange-Funktion beschreiben. Leider gibt es ein Problem mit zwei Punkten, ( 1 , 0 ) und ( 1 , 0 ) , weil die Ableitung F j Null wird, also die Punkte, die die Lösungen von verbinden j die zweite Bedingung verletzen und wir sind nicht in der Lage, die globale Bewegung des Teilchens zu beschreiben (das Teilchen könnte von einer Lösung zur anderen „hüpfen“)!

Natürlich lässt sich das Problem im obigen Beispiel leicht lösen, indem man andere Anfangskoordinaten wählt, nämlich Polarkoordinaten ( r , θ ) . Die Einschränkung wird F ( r , θ ) = r 1 = 0 , die beide Bedingungen erfüllt.

Beispiel 2: Viereck

Stellen Sie sich zwei Pendel vor, die an einer festen Decke hängen und durch eine starre Schnur zusammengebunden sind:

Beispiel zweier miteinander verbundener Pendel

Wählen a 1 und a 2 Als Anfangskoordinaten gibt es nur eine Einschränkung: die d -Saite, die l -Saiten und die Decke sollten immer ein Viereck mit vorgegebenen Längen für alle Seiten bilden.

Mit etwas uninteressanter Geometrie kann man eine Beziehung zwischen finden a 1 und a 2 , wie hier in einer Überarbeitung durchgeführt. Der Punkt ist, dass, wenn etwas gegeben wird a 1 , gibt es wieder zwei mögliche Lösungen für a 2 : a 2 = ϕ 1 ( a 1 ) + ϕ 2 ( a 1 ) oder a 2 = ϕ 1 ( a 1 ) ϕ 2 ( a 1 ) , wo ϕ 1 und ϕ 2 sind wie in der Revision definiert, aber ohne Bedeutung. Hier ist eine Veranschaulichung der beiden Möglichkeiten:

zwei Möglichkeiten dargestellt

Die Einschränkung kann geschrieben werden als F ( a 1 , a 2 ) = | a 2 ϕ 1 ( a 1 ) | ϕ 2 ( a 1 ) = 0 , was die erste Bedingung verletzt, oder es kann geschrieben werden als F ( a 1 , a 2 ) = ( a 2 ϕ 1 ( a 1 ) ) 2 ϕ 2 ( a 1 ) 2 = 0 , was die zweite Bedingung verletzt. Auch hier bin ich nicht in der Lage, die globale Bewegung des Systems zu beschreiben.

Im vorherigen Beispiel könnte dies durch die Wahl anderer Anfangskoordinaten gelöst werden, aber ich weiß nicht, ob dies hier möglich ist. Damit komme ich zu meinen Fragen.

Fragen

  • Ist es immer möglich, Anfangskoordinaten so zu wählen, dass alle Nebenbedingungen die Bedingungen des Satzes über implizite Funktionen erfüllen?
  • Wenn ja, gibt es einen systematischen Weg, sie zu finden?
  • Wenn nein, ist es dann noch möglich, die globale Bewegung des Systems mit der Lagrange-Funktion zu beschreiben?
Kommentare sind nicht für längere Diskussionen gedacht; Diese Konversation wurde in den Chat verschoben .
Die Seite „Diese Konversation wurde in den Chat verschoben“ wurde nicht gefunden

Antworten (1)

Im Allgemeinen, wenn ein Satz von Koordinaten gegeben ist x 1 , , x N unter einer Reihe von h = N n holonome Beschränkungen der Form F j ( x 1 , , x N ) = 0 , können Sie keine Teilmenge finden x 1 , , x n Ihrer ursprünglichen Koordinaten, die global als verallgemeinerte Koordinaten funktionieren: Der Satz über implizite Funktionen sagt Ihnen, dass Ihnen lokal eine solche Teilmenge garantiert ist, aber im Allgemeinen wird diese Teilmenge nicht überall funktionieren.

Das ist sozusagen der Hauptgrund, warum wir mit verallgemeinerten Koordinaten arbeiten q 1 , , q n , weil sie es uns ermöglichen, den verfügbaren Konfigurationsraum so neu zu parametrisieren, dass weniger Koordinatendiagramme erforderlich sind. Aber selbst diese sind keine vollständige Lösung, da die Beschränkungen im Allgemeinen den verfügbaren Koordinatenraum als eine Mannigfaltigkeit definieren , die möglicherweise zwei oder mehr Diagramme in ihrem minimalen Atlas erfordert , dh Ihnen wird niemals die Existenz eines Satzes von Koordinaten garantiert, der dies tun wird global arbeiten. Einige Beispiele:

  • Ein Teilchen in einer Ebene, die auf den Einheitskreis beschränkt ist, wird am besten über den Polarwinkel parametrisiert, aber selbst das ist nicht perfekt (es kann die Periodizität wörtlich genommen nicht reproduzieren).
  • Für etwas, bei dem klarer ist, dass Koordinatenpatches einfach nicht funktionieren, betrachten Sie zB ein Partikel in 3D, das auf die Oberfläche eines Doppeltorus beschränkt ist : Ein einzelnes Koordinatenpatch wird einfach nicht in der Lage sein, mit der Nicht-Null - Gattung der Oberfläche fertig zu werden.

Aus diesem Grund arbeitet die analytische Mechanik, wenn sie richtig gemacht wird (wie zB im Stil von VI Arnol'd ), abstrakt mit dem Konfigurationsraum als Mannigfaltigkeit und erkennt an, dass Koordinatenflecken immer nur lokal sein können – und das ist in Ordnung. Das bedeutet, dass Sie überall dort, wo Ihr Konfigurationsraum eine regelmäßige Mannigfaltigkeit ist, immer ein Diagramm finden können, das funktioniert, und Sie können die Dynamik dort explizit koordinatenabhängig lösen. Wenn die Lösung über den Rand des Diagramms hinauswandert, können Sie immer ein separates Diagramm verwenden, in dem die Dinge in Ordnung sind.

Bei Ihrer Konfiguration müssen Sie jedoch nicht besonders ausgefallen sein, da die Kerntopologie des Problems im Allgemeinen dieselbe ist wie die des oben erwähnten Einheitskreises. Dies kann durch Umformulieren der Einschränkung im Formular angezeigt werden F ( a 1 , a 2 ) = d wo

F ( a 1 , a 2 ) = ( D l cos ( a 2 ) , Sünde ( a 2 ) ) ( l cos ( a 1 ) , l Sünde ( a 1 ) ) .
Wie analysieren Sie das? Nun, werfen Sie einen Blick auf ein Konturdiagramm von F ( a 1 , a 2 ) auf der a 1 , a 2 Flugzeug (nehmen l = D / 2 zur Bestimmtheit):

Mathematica-Grafiken

Je nachdem was d Das heißt, Ihr System ist auf eine der Konturlinien im Diagramm beschränkt, und zum größten Teil sind dies topologische Kreise der langweiligsten, banalsten Art. Das bedeutet, dass Sie eine einzelne verallgemeinerte Winkelkoordinate finden können, die die Gesamtheit der Konfigurationsraum-Mannigfaltigkeit modulo der gleichen Periodizitätsprobleme parametrisiert, die Sie für ein auf den Einheitskreis beschränktes Teilchen erhalten.

Das soll nicht heißen, dass es einfach sein wird, diese Koordinate zu finden, und tatsächlich ist (i) ihre analytische Form wahrscheinlich ziemlich unordentlich und (ii) diese Unordnung wird in den Euler-Lagrange-Bewegungsgleichungen widergespiegelt und verstärkt. In diesem speziellen Fall glaube ich nicht, dass es etwas zu gewinnen gibt, wenn man versucht, eine solche Koordinate zu finden (im Gegensatz zur einfachen Verwendung von a 1 und a 2 als unabhängige Variablen in den Diagrammen, in denen sie arbeiten, das Erhalten und Lösen der Euler-Lagrange-EOMs für diese Parametrisierung und nur das Wissen, wann man möglicherweise die Diagramme wechseln muss), aber letztendlich ist es eine persönliche Entscheidung.


Nun, das heißt, es gibt einige (sehr spezifische) Situationen, in denen der Koordinatenraum möglicherweise keine Mannigfaltigkeit ist, im Wesentlichen weil F ( a 1 , a 2 ) hat keinen Gradienten ungleich Null. Bei deinem Problem tritt dies zB bei der auf d = 0,9 D Konturlinie wann l = 0,95 D , das sieht so aus,

Mathematica-Grafiken

dh es hat eine Kreuzung bei a 1 = a 2 = 0 , wo die beiden Arme „gekreuzt“ sind, und Sie können aus der Konfiguration herauskommen, indem Sie einen Arm wegziehen und den anderen folgen lassen, dh über zwei unterschiedliche Wege. Diese Art von Konfiguration ist die schlimmste, die Sie auf den von Ihnen erwähnten Systemen finden werden – aber Sie werden ihr mit ziemlicher Sicherheit nie begegnen, weshalb es im Allgemeinen in Ordnung ist, die Möglichkeit als pathologischen Sonderfall abzutun, über den sich die Mathematiker Sorgen machen müssen um.

Wenn Sie darauf bestehen, damit umzugehen, dann wird aus der Sicht eines Physikers im Allgemeinen passieren, dass das System dazu neigt, eine gewisse Trägheit zu haben, und das wird dazu neigen, es auf demselben „Zweig“ zu tragen, auf dem es läuft. (Oder mit anderen Worten, Sie parametrisieren die Achterkurve auf reguläre Weise, mit einem Selbstdurchgang, den Sie dann ignorieren.) Dies funktioniert in allen Fällen, außer in den Fällen, in denen sich das System am Schnittpunkt mit Nullgeschwindigkeit befindet , in diesem Fall bleibt es dort oder kommt nur asymptotisch an (aber, wenn Sie die Dinge wirklich vorantreiben, ja: Die Lagrange-Mechanik ist undefiniert, nachdem Sie einen solchen pathologischen Punkt erreicht haben).

Oder anders gesagt, diese Art von Problem ist nicht wirklich ein Problem.

Gute Antwort. Ich kann es nur bewundern.
Obwohl ich einige Bemerkungen zu Kleinigkeiten habe, nur um zu überprüfen, ob ich etwas falsch verstanden habe. (1) "im Gegensatz zur bloßen Verwendung a 1 und a 2 als unabhängige Variablen" Ich vermute, Sie meinen abhängige Variablen, da es um ihre Abhängigkeit geht. (2) "Dies geschieht z. B. bei der d = 0,9 D Konturlinie wann l = 0,95 d Ich glaube du meintest l = 0,95 D , Andernfalls a 1 = a 2 = 0 ist keine gültige Konfiguration.
(3) "dass Koordinatenpatches immer nur lokal sein können - und das ist in Ordnung" Es fällt mir schwer zu verstehen, wie das praktisch in Ordnung ist. Für den Fall des Einheitskreises wo j = ± ( 1 x 2 ) , hat jede Bewegungsgleichung einen singulären Punkt für x = ± 1 , wie werden Sie also die Diagramme wechseln, wenn die Bewegung dazwischen undefiniert ist?
(Für (3) muss ich gestehen, dass ich ein Anfänger bin, wenn es um Topologie geht. Es könnte sein, dass ich die genaue Bedeutung von Mannigfaltigkeiten, Koordinatenfeldern oder Diagrammen falsch verstanden habe, also hoffe ich, dass die Frage sinnvoll ist.)
Zu (2): Du hast vollkommen Recht, das war ein Tippfehler.
Bei (1) meine ich die Verwendung mehrerer verschiedener Diagramme: zwei davon a 1 ist unabhängig u a 2 abhängig ist, und zwei in denen a 2 ist die unabhängige Variable. Dies ist ähnlich wie die Antwort auf (3): wenn die Bewegung nahe kommt ( 1 , 0 ) , wechseln Sie einfach in das Diagramm x = 1 j 2 und verwenden j als Ihre unabhängige Variable dort. Sobald Sie diesen Punkt verlassen haben, können Sie zum anderen Zeichen von wechseln j = ± 1 x 2 wenn Sie es wünschen.
Mit anderen Worten, die Idee von Diagrammen besteht nicht darin, dass ihre Grenzen genau übereinstimmen, sondern dass sich ihre Domänen überlappen, sodass es für jeden bestimmten Punkt ein Diagramm gibt, bei dem sich der Punkt im Inneren des Diagramms befindet.