WICHTIG: Der Betrieb von goMatlab.de wird privat finanziert fortgesetzt. - Mehr Infos...

Mein MATLAB Forum - goMatlab.de

Mein MATLAB Forum

 
Gast > Registrieren       Autologin?   

Partner:




Forum
      Option
[Erweitert]
  • Diese Seite per Mail weiterempfehlen
     


Gehe zu:  
Neues Thema eröffnen Neue Antwort erstellen

Eigenwert einer komplexeren "symbolischen" Matrix

 

dati
Forum-Newbie

Forum-Newbie


Beiträge: 8
Anmeldedatum: 02.08.13
Wohnort: Berlin
Version: R2013a
     Beitrag Verfasst am: 04.08.2013, 16:01     Titel: Eigenwert einer komplexeren "symbolischen" Matrix
  Antworten mit Zitat      
Hallo zusammen,

es ist wieder ein Problem aufgetreten wo ich nicht weiterkomme.
Ich hoffe, dass mir wieder jemand helfen kann!

Wenn ihr in den Quellcode schaut werdet ihr einige Rechtschreibfehler finden... Nein ich habe keine Rechtschreibschwäche Sad , da ich die griechischen Symbole in meinem Quelltext mehrfach in unterschiedlicher Bedeutung verwende habe ich diese absichtlich falsch geschrieben. Laughing

Code:
clear all; clc;
syms mue k lamda roh omega E a b c b2 t x y u1 u2 u3;
syms a b c;

%Von dieser Matrix sollen die Eigenwerte unter Berücksichtigung von E bestimmt werden!
% Oder anders: Das homogene Gleichungssystem Matrix*u=0 soll gelöst werden.
%Nichttriviale Lösung: det(Matrix)=0

Matrix=[mue*k^2*E^2+(lamda+2*mue)*k^2-roh*omega^2 0 k^2*E*(lamda+mue);
0 mue*k^2*E^2+mue*k^2-roh*omega^2 0;
k^2*E*(lamda+mue) 0 (lamda+2*mue)*k^2*E^2+mue*k^2-roh*omega]

% Hier die möglichen 6 Lösungen bestimmen
D=simplify(solve(det(Matrix)==0,E))

%Test für erste Lösung
simplify(subs(Matrix,E,D(1,1)))
%=> Ergebnis: Stimmt!!

%Test für zweite Lösung
simplify(subs(Matrix,E,D(2,1)))
%=> Ergebnis: Stimmt!!

%Test für nächste Lösungen (beispielhaft für D(3,1))
simplify(subs(Matrix,E,D(3,1)))
%=> Ergebnis: Führen nicht zu Null!!! Aber warum?



Ich möchte nochmal kurz mein Problem beschreiben:
Ich habe ein homogenes Gleichungssystem, dass ich auch in analytischer Lösung vorliegen habe.

Ich möchte das System nun in Matlab lösen.

Um eine nichttriviale Lösung zu erhalten habe ich die det(Matrix)=0 gesetzt. Der Parameter E ist der einzige Stellhebel. Daher möchte ich die Lösung unter E bekommen.

Für die ersten beiden E's kommt auch das gewünschte Resultat raus wenn ich die Lösungen eintrage nur ab Lösung 3-6 kommen in der Matrix nicht die Nulleinträge raus.

EDITIERT:
Ich habe nun festgestellt, dass mit den Lösungen 3-6 die Determinante zu Null führt. Das ist also korrekt.


Und weiter verwundert mich, dass die Lösungen 3-6 derart komplex/umfangreich sind. Eigentlich sollten Lösungen in dem Umfang wie bei 1,2 herauskommen.

Die übrigbleibende Frage ist, ob sich die Lösungen 3-6 im Ausdruck vereinfachen lassen, da die analytischen Lösungen sehr simpel sind. Und ich hier nicht verstehe warum die Matlab-Lösungen so aufgeblasen sind.

Ich hoffe ihr könnt mir hier weiterhelfen!!!

Gruß
dati
Private Nachricht senden Benutzer-Profile anzeigen


dati
Themenstarter

Forum-Newbie

Forum-Newbie


Beiträge: 8
Anmeldedatum: 02.08.13
Wohnort: Berlin
Version: R2013a
     Beitrag Verfasst am: 04.08.2013, 16:17     Titel:
  Antworten mit Zitat      
ACHTUNG: Frage editiert: siehe ROT! Embarassed
Private Nachricht senden Benutzer-Profile anzeigen
 
Harald
Forum-Meister

Forum-Meister


Beiträge: 24.492
Anmeldedatum: 26.03.09
Wohnort: Nähe München
Version: ab 2017b
     Beitrag Verfasst am: 04.08.2013, 16:45     Titel:
  Antworten mit Zitat      
Hallo,

wie hast du denn die analytischen Lösungen ermittelt, und wie sehen sie aus?

Ich persönlich würde bei so relativ komplexen Matrizen die Eigenwerte numerisch berechnen.

Grüße,
Harald
Private Nachricht senden Benutzer-Profile anzeigen
 
dati
Themenstarter

Forum-Newbie

Forum-Newbie


Beiträge: 8
Anmeldedatum: 02.08.13
Wohnort: Berlin
Version: R2013a
     Beitrag Verfasst am: 04.08.2013, 16:54     Titel:
  Antworten mit Zitat      
Die Eigenwerte lauten:

+/- i*B
+/- i*A

A = sqrt( 1 - omega^2*rho / ( (lambda+2mu)*k^2 ) )
B = sqrt( 1 - omega^2*rho / ( mu*k^2 ) )

Die Eigenwerte mit B konnte ich ja verifizieren, allerdings diejenigen mit A bekomme ich nicht raus.

Zu deiner Frage wie ich zu dem Ergebnis kam: Ich habe hier ein Paper liegen wo die Ergebnisse enthalten sind. Und das Ergebnis mit B bekomme ich ja auch raus Wink

Mir wäre sehr wichtig, dass ich das Ergebnis mit A bestätigen könnte.

Da hoffe ich etwas auf eure Hilfe, denn ich sitze hier seit vieeelen Stunden dran. Embarassed
Private Nachricht senden Benutzer-Profile anzeigen
 
dati
Themenstarter

Forum-Newbie

Forum-Newbie


Beiträge: 8
Anmeldedatum: 02.08.13
Wohnort: Berlin
Version: R2013a
     Beitrag Verfasst am: 04.08.2013, 17:39     Titel:
  Antworten mit Zitat      
Ich glaube ich hatte jetzt eine Erleuchtung. Cool

Es sieht im Moment so aus, als ob ich mein Problem selbst gelöst habe. Laughing
Offenbar war die Matrix in der Publikation bereits falsch.

Somit braucht ihr im Moment nicht weiter an meinem Problem zu helfen!

Ich werde berichten wenn ich weitere Infos habe! Very Happy
Private Nachricht senden Benutzer-Profile anzeigen
 
dati
Themenstarter

Forum-Newbie

Forum-Newbie


Beiträge: 8
Anmeldedatum: 02.08.13
Wohnort: Berlin
Version: R2013a
     Beitrag Verfasst am: 04.08.2013, 18:08     Titel:
  Antworten mit Zitat      
In der Tat war die ursprüngliche Matrix falsch. Shocked
Gelöst habe ich das Problem, so dass ich das gesamte Problem analytisch gelöst habe und dann dementsprechend in Matlab die Situation nachgestellt habe.

Eine 100% Übereinstimmung. Very Happy

Sorry für die Verwirrung und die Störung meinerseits. Embarassed Embarassed

Viele Grüße
dati
Private Nachricht senden Benutzer-Profile anzeigen
 
Winkow
Moderator

Moderator



Beiträge: 3.842
Anmeldedatum: 04.11.11
Wohnort: Dresden
Version: R2014a 2015a
     Beitrag Verfasst am: 04.08.2013, 18:21     Titel:
  Antworten mit Zitat      
hatte ne freundin von mir auch mal ^^ sie hat in ihrer diplomarbeit 3 monate lang nen fehler gesucht der aber nicht bei ihr lag sondern in den ausgangs daten die sie von ihrem prof bekommen hat Smile sind alles nur menschen darum immer beides überprüfen wenns dick kommt Smile
_________________

richtig Fragen
Private Nachricht senden Benutzer-Profile anzeigen
 
Neues Thema eröffnen Neue Antwort erstellen



Einstellungen und Berechtigungen
Beiträge der letzten Zeit anzeigen:

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 | goMatlab RSS Button 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.