[<] [globale Übersicht] [Kapitelübersicht] [Stichwortsuche] [>]
Die klassischen Verfahren zur numerischen Auflösung linearer Gleichungssysteme beruhen auf bereits in Abschnit5.7.1geschilderten Vorgangsweise:Da eine Linearkombination von Gleichung nichts an der Lösung des Gleichungssystems ändert,werden geeignete Linearkombinationen zur systematischen Elimination von unbekanten benützt. Durch Multiplikation der ersten Gleichung des Systems
mit ai1/a11(a11 ungleich 0)und Subtraktion dieser gleichung von der i-ten Gleichung fällt der Term mit x1 in der entstehenden Gleichung weg. Man erhält so für i=2,3,...,n ein System von n-1 linearen Gleichungen
in den n-1 Unbekanten x2,x3,...,xn mit
das zusammen mit der ersten, unveränderten Gleichung
äquilent zum ursprunglichen System ist. Falls mit aii(k) = 0 keine Sondersituation auftritt, Läßt sich dieses Verfahren rekursiv fortsetzen. Nach n-1 Rekursionsschritten ist das verbleibende System schliesslich auf eine gleichung in einer Unbekanten xn reduziert,und das äquivalente Gleichungssystem hat folgende Gestalt
Die aii(k) = 0 und bi(k) sind dieneuen, durch die Eliminationsrechnung entstehenden Koeffizinten bzw. Elemente der rechten Seite, wobei der obere Index k zum Ausdruk bringt , wie oft der Wert von aij bzw.bihöchstens verändert wurde. Ein Gleichungssystem von der Struktur heisst gestaffeltes System oder Dreiecks-System, die zugehörige Koeffizintenmatrix ist eine obere (oder rechte)Dreiecksmatrix. ein solches Gleichungssystem läßt sich durch Rücksubstitution rekursiv lösen, falls alle Diagonalkoeffizinten von Null verschieden sind. Die letzte Gleichung liefert
nach Einsetzen von xn in die vorletzte Gleichung liefert diese
usw; schlißlich hat man die Lösungskomponenten x2,...,xn und erhält nach Einsetzen in die Erste Gleichung noch
FOR j:=1 TO n-1 do FOR i:=j+1 TO n do (Konstante * Zeile[j])+Zeile[i] sodaß a[ij]=0; Zeile[i]:=(Konstante * Zeile[j])+Zeile[i]; {Die Konstante könnte man so wählen aij-aij} end;
-X2+X3 = 0 -X1+X2+X3 = 1 X1+X2+X3 = 3.Lösung:-
Man bildet die Matrix
durch vertauschen zwischen 1. Und 2. Gleichung erhält man die Matrix :-
-1X2 + 1= 0, X2 = 1. Und ersetzt man den Wert von X3 und den Wert von X2 in die 1. Gleichung so erhält man -1X1 + 1 + 1 = 1 , X1 = 1.
[<] [globaleübersicht] [Kapitelübersicht] [Stichwortliste] [>]