Gauss-Jordaneliminatie

Gauss-Jordaneliminatie of methode van Gauss-Jordan is een uitbreiding van Gauss-eliminatie, een techniek waarmee een willekeurige matrix tot echelonvorm (trapvorm) kan worden teruggebracht. Met deze techniek kunnen onder andere lineaire vergelijkingen opgelost worden. De techniek bestaat, net als Gauss-eliminatie, uit rijoperaties op de matrix. Het verschil met Gauss-eliminatie is dat de matrix niet alleen van boven naar onder wordt geveegd, waarbij de getallen onder de diagonaal 0 worden, maar ook van onder naar boven, zodat er uiteindelijk alleen getallen op de diagonaal overblijven.

Gauss-Jordaneliminatie is aanzienlijk minder efficiënt dan Gauss-eliminatie met terugsubstitutie bij het oplossen van een stelsel van lineaire vergelijkingen. De methode is echter uitstekend geschikt voor het berekenen van inverse matrices.

De methode is genoemd naar Carl Friedrich Gauss en Wilhelm Jordan.

Berekenen van de inverse van een matrix

Gauss-Jordaneliminatie kan worden gebruikt om de inverse van een matrix uit te rekenen. Voor de inverse $A^{-1}$ van de vierkante matrix $A$ geldt onder andere dat

AA^{-1}=I

,

met I de betrokken eenheidsmatrix. Door de vergelijking op te lossen met Gauss-Jordaneliminatie kan de inverse gevonden worden.

Gauss-Jordaneliminatie in software

Matlab

Softwareprogramma's zoals Matlab ondersteunen het Gauss-Jordaneliminatieproces. In Matlab kan een vergelijking

Ax=B

met Gauss-Jordaneliminatie opgelost worden voor $x$ met:

x=A\setminus B

Een matrix $A$ kan geveegd worden met Gauss-Jordaneliminatie met het commando: rref(A) (rref staat voor Reduced Row Echelon Form).

De inverse

B=A^{-1}

van een matrix kan worden gevonden met het commando: B = A\eye(n), waarin $n$ het aantal rijen (en aantal kolommen, want $A$ is vierkant) voorstelt. Het commando: 'eye' levert de eenheidsmatrix op en met de backslash ('\') wordt de vergelijking opgelost voor de eenheidsmatrix.

Voorbeeld

Het volgende stelsel vergelijkingen wordt opgelost met Gauss-Jordaneliminatie.

3x + 2y -  z = 4
 x +  y +  z = 6
2x - 2y + 3z = 7

Daartoe wordt de bijbehorende uitgebreide matrix

B={\begin{bmatrix}3&2&-1&4\\1&1&1&6\\2&-2&3&7\end{bmatrix}}

geveegd.

Het element $B_{1,1}=3\neq 0$ , dus hoeft de eerste rij niet met een andere rij omgewisseld te worden. Vermenigvuldig nu voor het gemak de tweede en de derde rij met $B_{1,1}=3$ :

{\begin{bmatrix}3&2&-1&4\\3&3&3&18\\6&-6&9&21\end{bmatrix}}

Trek vervolgens de eerste rij van de tweede af en 2 keer de eerste rij van de derde:

{\begin{bmatrix}3&2&-1&4\\0&1&4&14\\0&-10&11&13\end{bmatrix}}

Nu staan in de eerste kolom onder de hoofddiagonaal nullen, en gaat men verder met de tweede rij. Het element $B_{2,2}=1\neq 0$ , dus hoeft de tweede rij niet met de derde omgewisseld te worden.

Trek 2 keer de tweede rij af van de eerste en tel 10 keer de tweede rij bij de derde op:

{\begin{bmatrix}3&0&-9&-24\\0&1&4&14\\0&0&51&153\end{bmatrix}}

Deel rij 1 door 3 en rij 3 door 51.

{\begin{bmatrix}1&0&-3&-8\\0&1&4&14\\0&0&1&3\end{bmatrix}}

Tel 3 keer de derde rij op bij de eerste en trek 4 keer de derde rij af van de tweede:

{\begin{bmatrix}1&0&0&1\\0&1&0&2\\0&0&1&3\end{bmatrix}}

Daarmee is de eliminatie voltooid en het stelsel opgelost:

x=1,y=2,z=3

De hier gevolgde procedure is niet systematisch. Een systematische methode deelt de eerste rij door de zogeheten spil $B_{1,1}=3$ zodat op die plaats direct het getal 1 komt te staan: