title: Gauss-Jordaneliminatie
text: 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  van de vierkante matrix  geldt onder andere dat 
  , 
  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 
   
  met Gauss-Jordaneliminatie opgelost worden voor  met:
  
  Een matrix  kan geveegd worden met Gauss-Jordaneliminatie met het commando: rref(A) (rref staat voor Reduced Row Echelon Form).
  
  De inverse 
   
  van een matrix kan worden gevonden met het commando: B = A\eye(n), waarin  het aantal rijen (en aantal kolommen, want  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
  
  geveegd.
  
  Het element , 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 :
  
  Trek vervolgens de eerste rij van de tweede af en 2 keer de eerste rij van de derde:
  
  Nu staan in de eerste kolom onder de hoofddiagonaal nullen, en gaat men verder met de tweede rij. Het element , 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:
  
  Deel rij 1 door 3 en rij 3 door 51.
  
  Tel 3 keer de derde rij op bij de eerste en trek 4 keer de derde rij af van de tweede:
  
  Daarmee is de eliminatie voltooid en het stelsel opgelost:
  
  De hier gevolgde procedure is niet systematisch. Een systematische methode deelt de eerste rij door de zogeheten spil  zodat op die plaats direct het getal 1 komt te staan: 
  
  Trek nu voldoende vaak rij 1 af van de rijen 2 en 3 zodat in de eerste kolom nullen komen te staan:
  
  Deel nu de tweede rij door de spil op positie (2,2):
  
  Trek voldoende vaak rij 2 af van de rijen 1 en 3 zodat in de tweede kolom nullen komen te staan:
  
  Deel de derde rij door de spil op positie (3,3):
  
  Trek voldoende vaak rij 3 af van de rijen 1 en 2 zodat in de derde kolom nullen komen te staan:
  
  Het eindresultaat is natuurlijk hetzelfde.
  
  Lineaire algebra