LR Zerlegung R³ - LU decomposition linear equations

Author:
hawe
(11) Pivot = 0 - ohne Pivotsuche = 1 - mit Spalten-Pivotsuche = 2 - mit Zeilen/Spalten-Pivotsuche Pi - Zeilen-Tausch-Matrizen ∈ Pivoti, Li - Matrizen für Gauß-Eliminationsschritt Spalte i ∈ Gaussi, Qi - Spalten-Tausch-Matrizen ∈ Pivoti L1P1AQ1 L1 (1.Spalte unter 1.Diagonal-Element Nullen erzeugen, P1,Q1 Zeilen-Spalten-Tausch für Pivot 1. Diagonal-Element [i=1] ADD n=4 (20) L2P2L1P1AQ1Q2:=Product(Gauss_{2}) (21) Pivot_{3}:={T(3, If(Pivot==2,PivotZ(L2P2L1P1AQ1Q2, 3),Pivot==1,SPvt(L2P2L1P1AQ1Q2,3),3) ),L2P2L1P1AQ1Q2,T(3,If(Pivot==2,PivotS(L2P2L1P1AQ1Q2, 3),3) )} (22) P3L2P2L1P1AQ1Q2Q3:= Product(Pivot_{3}) (23) Gauss_{3}:={Product(LE(P3L2P2L1P1AQ1Q2Q3, 3)), P3L2P2L1P1AQ1Q2Q3} R=(L3 P3 L2 P2 L1 P2 P3) (P3 P2 P1) A (Q1 Q2 Q3) (25) R:=Product(Gauss_{3}) (26) L:=( Element(Gauss_{3}, 1) Element(Pivot_{3},1) Element(Gauss_{2}, 1) Element(Pivot_{2},1) Element(Gauss_{1}, 1) Element(Pivot_{2},1) Element(Pivot_{3},1) )^(-1) (27) P:=Element(Pivot_{3},1) Element(Pivot_{2},1)*Element(Pivot_{1},1) (28) Q:=Element(Pivot_{1}, 3)*Element(Pivot_{2}, 3)*Element(Pivot_{3}, 3) (29) L R = P A Q ......... SMath - MathCad Paperlike Math LR-Zerlegung3x3MitTotalPivotSuche