Verfasst am: 25.10.2014, 14:39
Titel: Cholesky-Zerlegung einer Matrix
Hallo ich habe versucht eine Cholesky-Zerlegung zu erstellen. Jedoch muss
diese wohl einen Fehler haben. Ich kenne zwar eine funktionierende
Lösung wollte jedoch meinen ersten Gedanke nicht einfach so verwerfen.
Wenn jemand eine einfache Korrekturidee hätte wäre ich dafür dankbar.
Hier der Code:
Code:
function L = Cholesky44(A) [zeile,spalte] = size(A);
n = zeile;
L = zeros(n,n);
for i = 1:1:n
for j = 1:1:i
if(i == j)
C = 0;
for k = 1:1:(i-1)
C = C + (L(i,k))^2;
end
L(i,i) = sqrt( A(i,i) - C );
else
C2 = 0;
for r = 1:1:(i-1)
C2 = C2 + L(i,r)*L(j,r);
end
L(i,j) = (A(i,j) - C2)*1/A(j,j);
end end end
In der drittletzten Zeile muss das "1/A(j,j)" folgendermaßen heißen: 1/L(j,j).
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.