Eigenschwingung einer eingespannten Saite (DGL + Eigenwertproblem) |
22.06.2011, 10:27 | 8000mark | Auf diesen Beitrag antworten » | ||||||||
Eigenschwingung einer eingespannten Saite (DGL + Eigenwertproblem) Die Differentialgleichung und Anfangswerten beschreibt das Bewegungsverhalten einer eingespannten Saite im Intervall . Schwingungen sind Lösungen, die bezüglich der Zeit periodisch sind. Jede Schwingung ist die Reihe über Basisschwingungen der Gestalt . Setzt man dieses zur Lösung von Gl. (1) an, so erhält man das Eigenwertproblem U''(x) + \lambda U(x) = 0 \text{ mit } U(0) = 0 \text{ und } U(\pi) = 0 \text{ (2)} Bei der numerischen Lösung solcher Probleme diskretisiert man Gl. (2) und erhält ein Eigenwertproblem für Matrizen. Die Eigenvektoren sind die Diskretisierungen der Eigenlösung von Gl. (2). Eine Lösung von Gl. (4) erhält man, indem man eine solche Linearkombination von Eigenlösungen wählt, die den Anfangsbedingungen genügt. a) Bestimmen Sie analytisch die Eigenwerte und Eigenlösungen der Gleichung (2)! Diskretisieren Sie diese Gleichung. b) Implementieren Sie ein Verfahren zur numerischen Bestimmung der Eigenwerte und Eigenvektoren des diskretisierten Problems. Vergleichen Sie Ihre Ergebnisse mit denen aus Teil a). Meine Ideen: Ich habe in Gleichung (1) eingesetzt. Dabei kommt raus. Außerdem habe ich die Differentialgleichung (2) gelöst, wobei also die Sinus und Cosinus Gleichungen rauskommen. Wie bekomme ich jetzt den Zusammenhang zwischen DGL und Eigenwertproblem? Also wo bekomme ich die Matrix her, für die ich die Eigenlösungen bestimmen kann? Das Verfahren habe ich bereits implementiert. |
||||||||||
22.06.2011, 14:48 | Ehos | Auf diesen Beitrag antworten » | ||||||||
Deine Aufgabe ist nicht eindeutig lösbar, denn man benötigt zur Zeit t=0 nicht nur die anfängliche Auslenkung der Saite, sondern zusätzlich deren anfängliche Geschwindigkeit . Wenn diese bekannt ist, kann man die Lösung formelmäßig darstellen und muss nicht numerisch rechnen! Ich zeig' dir mal, wie das geht. Zu lösen ist Randbedingung: Anfangsbedingung: Zuerst löst man das zugehörige Eigenwertproblem (warum wird später klar): Mit den obigen Randbedingungen führt das auf folgende Eigenwerte und auf 1 normierte Eigenfunktionen Aufgrund der allgemeinen Theorie kann man damit jede Funktion im Intervall [0;pi], welche die Randbedingungen erfüllt, in einer Fourrierreihe entwickeln. Das gilt auch für die gesuchte Lösung u(x,t). Da diese Lösung aber zeitabhängig ist, muss man auch die Fourrierkoeffizienten als zeitabhängig betrachten. Das motiviert den Reihenansatz ________(*) Dies setzt man in die Dgl. ein und erhält Darin substituiert man und beseitigt so die 2.Ableitung auf der linken Seite. (Das ist die Motivation der Eigenwertgleichung!) Koeffizientenvergleich liefert eine gewöhnlich Differenzialgleichung für die zeitabhängigen Fourrierkoeffizienten mit den Eigenwerten Daraus folgt . Einsetzen in unseren Ansatz (*) liefert __________(**) Die Konstanten bestimmt man aus den Anfangsbedingungen. Demnach soll unsere Lösung (**) für t=0 folgendes ergeben Damit kann man die Koeffizienten auf bekannte Weise als Integrale darstellen und hat somit die endgültige Lösung (**) |
||||||||||
22.06.2011, 17:30 | 8000mark | Auf diesen Beitrag antworten » | ||||||||
Vorab: Danke für die schnelle Antwort!
Die Aufgabe ist für mein numerisches Praktikum. Ich denke, es geht hierbei vor alle um die Teilaufgabe b, also das Implementieren eines Eigenwert-Algorithmus.
Ich dachte immer bezeichnet das gleiche.
Das hab ich jetzt nicht verstanden. Ist die Hessematrix? Und wie bekommt man dann die Eigenwerte?
Nun steht ja auch in der Aufgabe etwas von Diskretisieren. Heißt das ich überlege mir Werte und rechne damit ein bisschen herum? (Im übertragenen Sinn) Schließlich brauch ich ja auch eine Matrix, auf die ich meinen Algorithmus loslassen kann. |
||||||||||
23.06.2011, 09:07 | Ehos | Auf diesen Beitrag antworten » | ||||||||
Ableitungen nach der Zeit kennzeichnet man oft durch einen Punkt und Ortsableitungen durch einen Strich, also und . Mir ist jetzt nicht ganz klar, warum du die Eigenwerte und Eigenfunktionen mit n=1,2,3,... numerisch bestimmen willst, wenn man beide formelmäßig kennt. Es gibt natürlich für komplizierte Probleme numerische Verfahren. Darauf kann ich aber hier nicht eingehen. Du solltest dich vielleicht vorher mit den grundsätzlichen Begriffen vertraut machen. |
||||||||||
23.06.2011, 13:26 | 8000mark | Auf diesen Beitrag antworten » | ||||||||
Die Aufgabe ist ja, die Gleichung zum Einen auf dem Papier zu lösen und zum Anderen die Lösung mithilfe eines Verfahrens am Computer zu bestimmen. Für den zweiten Teil hab ich das QR-Verfahren implementiert. Zum Lösen brauch ich allerdings noch eine Matrix, die das Eigenwertproblem beschreibt. Oder irre ich mich? |
||||||||||
23.06.2011, 15:42 | Ehos | Auf diesen Beitrag antworten » | ||||||||
Ich wüsste nicht, wie und warum man im Zusammenhang mit der schwingenden Saite ein Eigenwertproblem für eine Matrix formulieren sollte. Natürlich gibt es in anderem Zusammenhang Eigenwertproblem für Matrizen, aber nicht hierbei. Ich vermute, du bringst da etwas durcheinander. |
||||||||||
Anzeige | ||||||||||
|
||||||||||
23.06.2011, 15:46 | tigerbine | Auf diesen Beitrag antworten » | ||||||||
Nur eine kurze Wortmeldung. Aber das Nennen von "Saite und Eigenwert" im Zusammenhang ist nicht ungewöhnlich. Wurde in Numerik gerne als praktische Einleitungsbeispiel gemacht. http://www.math.uni-hamburg.de/home/muel...08ShtZ_0407.pdf |
||||||||||
24.06.2011, 10:19 | Ehos | Auf diesen Beitrag antworten » | ||||||||
@Tigerbiene Es ist schon klar, dass die Begriffe "Saite" und "Eigenwerte" zusammen gehören. Aber der Fragesteller will ja irgendeine Matrix aufstellen. Und da weiß ich ehrlich gesagt nicht, was er meint. In deiner Literaturangabe kommen ebenfalls keine Matrizen vor. Dort wird alles formelmäßig gelöst, wie ich es auch erklärt hatte. |
||||||||||
24.06.2011, 17:02 | tigerbine | Auf diesen Beitrag antworten » | ||||||||
Das LGS entsteht imho bei der Diskretisierung der zugehörigen Differentialgleichung. Sorry, wenn das im Link nicht raus kommt. Z.B. in "Numerische Mathematik für Anfänger, G. Opfer" ist die Saite ein Einführungsbeispiel im Kapitel "Eigenwerte" |
|
Verwandte Themen
Die Beliebtesten » |
|
Die Größten » |
|
Die Neuesten » |
|