Verfasst am: 17.12.2009, 17:44
Titel: Determinantenberechnung nach Laplace
Hi,
habe folgendes Problem.
Habe die Aufgabe ein Programm zu schreiben mit dem man die Determinante nach Laplace entwickeln kann. Dabei soll eine Abfrage nach welcher Zeile und nach welcher Spalte entwickelt werden soll eingebaut werden. Weiterhin soll rekursiv berechnet werden. Trotzdem liefert das Programm nicht die richtigen Ergebnisse. Ich hab ein Programm ohne Abfrage der Zeile und Spalte welches die richtigen Ergebnisse liefert. Dann hab ich lange rumgebastelt und jetzt gehts nicht mehr. Ich glaube ich habe einen Fehler in der rekursiven Berechnung.
Hoffe ihr könnt mir weiterhelfen.
Code:
% Das Programm das nicht funktioniert function d=laplace3(A) % d=DETLAPLACE(A)
if m~=n
error('Matrix nicht quadratsich')% Fehlerausgabe für nicht quadratische Matrix end if m==1% Sonderfall: 1x1 Matrix
d=A(1,1);
elseif m==2% Sonderfall: 2x2 Matrix
d=A(1,1)*A(2,2)-A(1,2)*A(2,1);
else ifrem(m,2)==0% Abfrage ob Matrix gerade oder ungerade
i = input('Bitte eingeben, nach welcher Zeile entwickelt werden soll! : ');
j = input('Bitte eingeben, nach welcher Spalte entwickelt werden soll! : ');
B=[A(1:i-1,1:j-1),A(1:i-1,j+1:n);A(i+1:m,1:j-1),A(i+1:m,j+1:n)];
d=A(i,j)*laplace3(B); % rekursive Berechnung else
i = input('Bitte eingeben, nach welcher Zeile entwickelt werden soll! : ');
j = input('Bitte eingeben, nach welcher Spalte entwickelt werden soll! : ');
B=[A(1:i-1,1:j-1),A(1:i-1,j+1:n);A(i+1:m,1:j-1),A(i+1:m,j+1:n)];
d=-A(i,j)*laplace3(B); % rekursive Berechnung end end
% Das Programm das funktioniert function d=laplace(A) % d=DETLAPLACE(A)
[m,n]=size(A);
if m~=n
error('Matrix nicht quadratsich')% Fehlerausgabe für nicht quadratische Matrix end if m==1% Sonderfall: 1x1 Matrix
d=A(1,1);
elseif m==2% Sonderfall: 2x2 Matrix
d=A(1,1)*A(2,2)-A(1,2)*A(2,1);
else
vz=1;
d=0;
for i=1:m % Berechnung der Determinanten für 1 bis m Spalten/Zeilen
B=[A(1:i-1,2:m);A(i+1:m,2:m)];
d=d+vz*A(i,1)*laplace(B); % rekursive Berechnung
vz=-vz;
end end
Verfasst am: 02.05.2016, 14:33
Titel: Cooles Programm
Hey GozyllaR,
sorry, dass wir uns erst jetzt melden. Leider waren wir mit wichtigerem beschäftigt!
Danke für deine Geduld! Unser Expertenteam konnte leider keinen Fehler feststellen.
Wir wünschen dir noch viel Erfolg.
Beste Grüße dein MATLAB-Expertenteam
Gast
Beiträge: ---
Anmeldedatum: ---
Wohnort: ---
Version: ---
Verfasst am: 13.10.2016, 23:43
Titel: Definition von Laplace-Entwicklung
Hallo zusammen,
z.B. stimmt das Vorzeichen von A(i,j) nicht, denn dieser hängt von i bzw. j ab, siehe die Definition der Laplace-Entwicklung.
Möglicherweise gibt es weitere Fehler.
Grüße
A.
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
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.
RSS