ich suche Hilfe bei der Programmierung einer Funktion zur LR-Zerlegung inkl. Spaltenpivotsuche. Bereits implementierte Funktionen sollen nicht benutzt werden.
Die LR-Zerlegung funktioniert bereits, doch bei der Pivotsuche könnte ich noch Hilfe gebrauchen.
Code:
function[ L, R, P ] = LRZerlegungStable(A)
n = size(A,1); % Anzahl der Zeilen
m = size(A,2); %Anzahl der Spalten
minimum = 0;
k = 0;
P = eye(n,m); % Permutationsmatrix
w = zeros(n, 1);
B = P*A;
Z=zeros(n,m); % Permutationsmatrix zur späteren Ausgabe while k < n*m-1
end
Z=B;
end
k = k+1;
end for i=1:n-1% Jeder Durchlauf der for-Schleife berechnet ein L^-1
B(i+1:n,i) = B(i+1:n,i)/B(i,i);
B(i+1:n,i+1:n) = B(i+1:n,i+1:n)-B(i+1:n,i)*B(i,i+1:n);
end
%%Laufvariable für L-und R-Matrix, damit B nicht verfälscht wird
C=B;
%% Matrix B in obere Dreicksmatrix-Form bringen
%% Platzieren der Nullen unterhalb der Diagonalen for i=1:n-1 for p=1:i
B(1+i,p)=0;
end;
end;
R=B;
%% Matrix C(bzw.B) in untere Dreiecksmatrix-Form bringen
%% Platzieren der Nullen oberhalb der Diagonalen for i=1:n-1 for p=n:-1:i+1
C(i,p)=0;
end end
%% Platzieren der Einsen in der Diagonalen for i=1:n
C(i,i)=1;
end
L=C;
Die Matrix A wird nicht permutiert, die Permutationsmatrix bleibt
die Einheitsmatrix, und die LR-Zerlegung kann nicht berechnet werden:
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