|
|
Eigenwerte und Eigenvektoren |
|
mbyrl |
Forum-Anfänger
|
|
Beiträge: 11
|
|
|
|
Anmeldedatum: 03.07.19
|
|
|
|
Wohnort: ---
|
|
|
|
Version: ---
|
|
|
|
|
|
Verfasst am: 18.11.2020, 18:49
Titel: Eigenwerte und Eigenvektoren
|
|
Hallo Zusammen,
ich möchte im Rahmen meiner Maschinendynamik-Vorlesung u.a. Eigenfrequenzen errechnen. Ich habe ein System welches in einer Massenmatrix M und einer Steifigkeitsmatrix K beschrieben ist.
Händisch würde ich folgendes rechnen: det(K-w0^2*M) !=0 und dann letztlich über die PQ-Formel meine Eigenwerte bekommen, deren Wurzel wiederum die Eigenkreisfrequenz darstellt. Da ich für weitere Berechnungen auch die Eigenvektoren benötige, bin ich auf folgenden Befehl gestoßen:
[V,D]=eig(...)
Allerdings kann ich ja in das Argument des Befehls schlecht K-M schreiben, da dies ja nicht zum korrekten Ergebnis führt (weil (K-w0^2*M)!=0 )
Jemand eine gute Idee?
|
|
|
|
|
Harald |
Forum-Meister
|
|
Beiträge: 24.492
|
|
|
|
Anmeldedatum: 26.03.09
|
|
|
|
Wohnort: Nähe München
|
|
|
|
Version: ab 2017b
|
|
|
|
|
|
Verfasst am: 18.11.2020, 22:41
Titel:
|
|
Hallo,
ich bin kein Maschinenbauer, daher die Frage: von welcher Matrix willst du denn nun die Eigenwerte berechnen? K und M sind vermutlich gegeben, w0 auch?
Grüße,
Harald
_________________
1.) Ask MATLAB Documentation
2.) Search gomatlab.de, google.de or MATLAB Answers
3.) Ask Technical Support of MathWorks
4.) Go mad, your problem is unsolvable ;)
|
|
|
mbyrl |
Themenstarter
Forum-Anfänger
|
|
Beiträge: 11
|
|
|
|
Anmeldedatum: 03.07.19
|
|
|
|
Wohnort: ---
|
|
|
|
Version: ---
|
|
|
|
|
|
Verfasst am: 18.11.2020, 23:20
Titel:
|
|
Harald hat Folgendes geschrieben: |
Hallo,
ich bin kein Maschinenbauer, daher die Frage: von welcher Matrix willst du denn nun die Eigenwerte berechnen? K und M sind vermutlich gegeben, w0 auch?
Grüße,
Harald |
Hallo Harald,
Gegeben sind die Matrizen M und K. Die Eigenfrequenz w0 ist gesucht. M ist eine Diagonalmatrix.
Gruß
|
|
|
Harald |
Forum-Meister
|
|
Beiträge: 24.492
|
|
|
|
Anmeldedatum: 26.03.09
|
|
|
|
Wohnort: Nähe München
|
|
|
|
Version: ab 2017b
|
|
|
|
|
|
Verfasst am: 18.11.2020, 23:36
Titel:
|
|
Hallo,
du suchst also die w0, für die diese Determinante 0 wird? Denn von 0 verschieden dürfte sie doch für fast alle w0 sein?
Wie groß sind M und K denn?
Grüße,
Harald
_________________
1.) Ask MATLAB Documentation
2.) Search gomatlab.de, google.de or MATLAB Answers
3.) Ask Technical Support of MathWorks
4.) Go mad, your problem is unsolvable ;)
|
|
|
mbyrl |
Themenstarter
Forum-Anfänger
|
|
Beiträge: 11
|
|
|
|
Anmeldedatum: 03.07.19
|
|
|
|
Wohnort: ---
|
|
|
|
Version: ---
|
|
|
|
|
|
Verfasst am: 19.11.2020, 09:27
Titel:
|
|
Harald hat Folgendes geschrieben: |
Hallo,
du suchst also die w0, für die diese Determinante 0 wird? Denn von 0 verschieden dürfte sie doch für fast alle w0 sein?
Wie groß sind M und K denn?
Grüße,
Harald |
In meinem Fall sind jetzt beide Matrizen 4x4. Und ja, genau, ich möchte w0 finden, dass die Determinante 0 wird. Mit der 'symbolic math toolbox' kann ich auch eine Lösung finden.
Ich glaube aber, dass es einen Befehl geben muss. Was genau macht [V,D] = eig(K,M)? Ich will auf das Komma im Argument des Befehls hinaus.
|
|
|
Harald |
Forum-Meister
|
|
Beiträge: 24.492
|
|
|
|
Anmeldedatum: 26.03.09
|
|
|
|
Wohnort: Nähe München
|
|
|
|
Version: ab 2017b
|
|
|
|
|
|
Verfasst am: 19.11.2020, 09:46
Titel:
|
|
Hallo,
aus der Doku:
Zitat: |
[V,D] = eig(A,B) returns diagonal matrix D of generalized eigenvalues and full matrix V whose columns are the corresponding right eigenvectors, so that A*V = B*V*D. |
Grüße,
Harald
_________________
1.) Ask MATLAB Documentation
2.) Search gomatlab.de, google.de or MATLAB Answers
3.) Ask Technical Support of MathWorks
4.) Go mad, your problem is unsolvable ;)
|
|
|
mbyrl |
Themenstarter
Forum-Anfänger
|
|
Beiträge: 11
|
|
|
|
Anmeldedatum: 03.07.19
|
|
|
|
Wohnort: ---
|
|
|
|
Version: ---
|
|
|
|
|
|
Verfasst am: 19.11.2020, 09:53
Titel:
|
|
Die Dokumentation hab ich mir schon angeguckt, sehe aber noch nicht, dass das für mich der richtige Befehl ist.
https://www.gomatlab.de/viewtopic,p,197869.html#197869
In diesem Thread ist allerdings ein ähnliches Problem mit dem Befehl anscheinend gelöst worden.
Es scheint so, als sei mein Problem im Verständnis der Matrizenoperationen die da abgehen
|
|
|
Harald |
Forum-Meister
|
|
Beiträge: 24.492
|
|
|
|
Anmeldedatum: 26.03.09
|
|
|
|
Wohnort: Nähe München
|
|
|
|
Version: ab 2017b
|
|
|
|
|
|
Verfasst am: 19.11.2020, 10:51
Titel:
|
|
Hallo,
Zitat: |
Die Dokumentation hab ich mir schon angeguckt, sehe aber noch nicht, dass das für mich der richtige Befehl ist. |
Das kann ich nicht beurteilen. Ich würde wie gesagt den Weg über die Symbolic Math Toolbox gehen.
Zitat: |
In diesem Thread ist allerdings ein ähnliches Problem mit dem Befehl anscheinend gelöst worden. |
Du hast auf den momentanen Thread verlinkt.
Zitat: |
Es scheint so, als sei mein Problem im Verständnis der Matrizenoperationen die da abgehen |
Damit man dir helfen kann, müsstest du konkretere Fragen stellen.
Grüße,
Harald
_________________
1.) Ask MATLAB Documentation
2.) Search gomatlab.de, google.de or MATLAB Answers
3.) Ask Technical Support of MathWorks
4.) Go mad, your problem is unsolvable ;)
|
|
|
mbyrl |
Themenstarter
Forum-Anfänger
|
|
Beiträge: 11
|
|
|
|
Anmeldedatum: 03.07.19
|
|
|
|
Wohnort: ---
|
|
|
|
Version: ---
|
|
|
|
|
|
Verfasst am: 19.11.2020, 12:09
Titel:
|
|
|
|
Harald |
Forum-Meister
|
|
Beiträge: 24.492
|
|
|
|
Anmeldedatum: 26.03.09
|
|
|
|
Wohnort: Nähe München
|
|
|
|
Version: ab 2017b
|
|
|
|
|
|
Verfasst am: 19.11.2020, 12:16
Titel:
|
|
Hallo,
in dem letzten Beitrag zu dem Thema wurde gezeigt, dass mit dem symbolischen Ansatz und mit eig das gleich herauskommt.
Wenn du dir von eig zwei Argumente zurückholst, bekommst du Eigenwerte und Eigenvektoren.
Gibt es denn nun noch verbleibende Fragen?
Grüße,
Harald
_________________
1.) Ask MATLAB Documentation
2.) Search gomatlab.de, google.de or MATLAB Answers
3.) Ask Technical Support of MathWorks
4.) Go mad, your problem is unsolvable ;)
|
|
|
mbyrl |
Themenstarter
Forum-Anfänger
|
|
Beiträge: 11
|
|
|
|
Anmeldedatum: 03.07.19
|
|
|
|
Wohnort: ---
|
|
|
|
Version: ---
|
|
|
|
|
|
Verfasst am: 19.11.2020, 12:17
Titel:
|
|
Hab gerade meinen vorherigen Kommentar bearbeitet. Hatte einen Knoten im Kopf, alles gut jetzt. Vielen Dank!
|
|
|
|
|
Einstellungen und Berechtigungen
|
|
Du kannst Beiträge in dieses Forum schreiben. Du kannst auf Beiträge in diesem Forum antworten. Du kannst deine Beiträge in diesem Forum nicht bearbeiten. Du kannst deine Beiträge in diesem Forum nicht löschen. Du kannst an Umfragen in diesem Forum nicht mitmachen. Du kannst Dateien in diesem Forum posten Du kannst Dateien in diesem Forum herunterladen
|
|
Impressum
| Nutzungsbedingungen
| Datenschutz
| FAQ
| RSS
Hosted by:
Copyright © 2007 - 2024
goMatlab.de | Dies ist keine offizielle Website der Firma The Mathworks
MATLAB, Simulink, Stateflow, Handle Graphics, Real-Time Workshop, SimBiology, SimHydraulics, SimEvents, and xPC TargetBox are registered trademarks and The MathWorks, the L-shaped membrane logo, and Embedded MATLAB are trademarks of The MathWorks, Inc.
|
|