| |
2.1 Vorwärtselimination
Ziel dieser ersten Etappe ist es, das Gleichungssystem auf eine Dreiecksform zu bringen, wobei unterhalb der Hauptdiagonalen Nullen stehen und oberhalb neue Koeffizienten stehen sollen.
Das Gleichungssystem soll nach dieser Etappe also wie folgt aussehen:
a'11x1 + a'12x2 +...+ a'1nxn = b'1
0 + a'22x2 +...+ a'2nxn = b'2
..................................
0 + 0 + 0 + a'nnxn = b'n
Ein Gleichungssystem, das auf diese Form gebracht wurde, wird auch als gestaffeltes Gleichungssystem bezeichnet.
Die erste Zeile im Gleichungssystem enspricht bereits der Zeile im gestaffelten Gleichungssystem. In allen weiteren Zeilen muß man nun mit geeigneten Operationen in der ersten Spalte Nullen erzeugen. Dazu benutzt man die zweite Elementaroperation bei Matrizen. Man multipliziert dazu die erste Zeile mit einem geeigneten Faktor c und subtrahiert dieses Produkt dann von der zweiten Zeile.
Der Faktor c wird so berechnet:
c=a21:a11
Nach diesem Schritt hat man die erste Null erzeugt:
a11x1 + a12x2 +...+ a1nxn = b1
0 + a'22x2 +...+ a'2nxn = b'2
...............................
an1x1 + an2x2 +...+ annxn = bn
Im zweiten Schritt multipliziert man wieder die erste Zeile mit einem Faktor c, allerdings dieses Mal mit c=a31:a11 und subtrahiert diese von der dritten Zeile. Dann wird das Produkt aus c=a41:a11 und erster Zeile von der vierten Zeile subtrahiert. So verfährt man weiter bis zur letzten Zeile. Nach diesem Vorgang sind alle ai1 mit i=2..n eliminiert:
a11x1 + a12x2 +...+ a1nxn = b1
0 + a'22x2 +...+ a'2nxn = b'2
.................................
0 + a'n2x2 +...+ a'nnxn = b'n
Jetzt stimmt auch schon die zweite Zeile mit der der gewünschten Endform überein. Man muß nun die zweite Spalte ab der dritten Zeile eliminieren. Dazu schreibt man beim Faktor c jeweils a'22 in den Nenner. Im Zähler steht dann wieder der jeweilige Koeffizient, der Null werden soll. Außerdem multipliziert man nun nicht mehr die erste Zeile mit diesem Faktor, sondern die zweite Zeile. Wir befinden uns also in der zweiten Eliminationsstufe. Sonst verfährt man wie bei der ersten Eliminationsstufe. Bei der Eliminationsstufe k wird der Faktor c also immer aus einem Quotienten mit dem ersten Koeffizient der k-ten Zeile im Nenner und dem zu eliminierenden Koeffizienten im Zähler gebildet. Mit diesem Faktor wird nun die k-te Zeile multipliziert und von der Zeile mit dem zu eliminierendem Koeffizient subtrahiert.
Der Faktor c allgemein:
cik = aik:akk
mit
n = Anzahl der Gleichungen und Unbekannten
k = 1..n-1; Eliminationsstufe
i = 1..n-1; Zeile
Nachdem man n-1 Eliminationen durchgeführt hat, ist die Dreiecksform erreicht:
w11x1 + w12x2 +...+ w1nxn
= v1
w22x2 +...+ w2nxn = v2
..............................
wnnxn = vn
Die neuen Koeffizienten heißen jetzt wij mit i,j=1..n, da die Koeffizienten
aij mit i=j=1..n durch die bisherigen Umformungen andere Werte wij
mit i=j=1..n angenommen haben. Dasselbe gilt für die Freigelieder bi
mit i=1..n und vi mit i=1..n.
Aus dieser Dreiecksform lassen sich nun die Unbekannten xi mit i=1..n
durch Rückwärtseinsetzen leicht errechnen.
|
| |
|
|