|
|
|
LR Zerlegung mit Spaltenpivotzahl |
|
| Tweety84 |
Forum-Newbie
|
 |
Beiträge: 8
|
|
|
|
Anmeldedatum: 02.01.11
|
|
|
|
Wohnort: Hannover
|
|
|
|
Version: ---
|
|
|
|
 |
|
Verfasst am: 06.01.2011, 12:06
Titel: LR Zerlegung mit Spaltenpivotzahl
|
 |
Ich bearbeite gerade folgende Aufgabenstellung:
Implementieren Sie die unten stehenden direkten LGS-Löser für eine beliebige Matrix A und rechte Seite b.
- LR-Zerlegung mit Spaltenpivotwahl
Dazu habe ich nun überlegt:
1. Zerlegung der Matrix A: A = L R
for i=1:n
Assembliere L
for k=1:i-1
L[i,k] := A[i,k];
for j=1:k-1
L[i,k] := L[i,k]-L[i,j]*R[j,k];
end
L[i,k] := L[i,k]/R[k,k];
end
Assembliere R
for k=i:n
R[i,k] := A[i,k];
for j=1:i-1
R[i,k] := R[i,k]-L[i,j]*R[j,k];
end
end
end
2. Vorwärtssubstitution: Ly = b
for i=1:n
y[i] := b[i];
for j=1:i-1
y[i] := y[i]-L[i,j]*y[j];
end
end
3. Rückwärtssubstitution: Rx = y
for i=n:-1:1
x[i] := y[i];
for j=i+1:n
x[i] := x[i]-R[i,j]*x[j];
end
x[i] := x[i]/R[i,i];
end
Bin ich da auf dem richtigen Weg??
|
|
|
|
|
|
| Jan S |
Moderator
|
|
Beiträge: 11.057
|
|
|
|
Anmeldedatum: 08.07.10
|
|
|
|
Wohnort: Heidelberg
|
|
|
|
Version: 2009a, 2016b
|
|
|
|
|
|
|
Verfasst am: 06.01.2011, 18:48
Titel: Re: LR Zerlegung mit Spaltenpivotzahl
|
|
Hallo Tweety84,
Den Code kann man hier am besten lesen, wenn er zwischen [ c o d e ] und [ / c o d e] eingebettet ist (ohne die Leerzeichen).
Ja, Du bist auf dem richtigen Weg.
Wenn Du L berechnest, hat das gleichzeitig Auswirkungen auf b. Wikipedia und Wapedia wissen da einiges drüber.
In Matlab werden die Indies mit runden Klammern geschrieben.
Gruß, Jan
|
|
|
|
| lucycharlotte123 |
Forum-Newbie
|
|
Beiträge: 3
|
|
|
|
Anmeldedatum: 08.01.11
|
|
|
|
Wohnort: ---
|
|
|
|
Version: ---
|
|
|
|
|
|
|
Verfasst am: 09.01.2011, 00:09
Titel:
|
|
Hallo Tweetie,
Ich denke, dass wir die gleiche Übung machen!
Ich denke, dass erstmals du size(A) benutzen muss und dann R=zeros(A,1) und L=zeros(A,1)... Matlab versteht nicht was R und L sind...
Du musst dann auch sagen, dass diagonal elementen von L 1 sind.
Es könnte auch nützlich sein, wenn du ein Nachprufen hattest zu zeigen ob der Matrix invertierbar ist.
Tut mir leid wenn mein Deutsch unverständlich ist... ich bin erasmus studentin 
Vielleicht kannst du mir ein PN senden um die andere Fragen zu diskutieren.
|
|
|
|
|
|
|
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 - 2025
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.
|
|
