In de wiskunde is Gaussiaanse eliminatie (ook wel rijreductie genoemd) een algemene en systematische methode om stelsels van lineaire vergelijkingen op te lossen. De naam verwijst naar Carl Friedrich Gauss, die de methode populair maakte door erover te schrijven; historisch bestaan er varianten die ouder zijn dan Gauss, maar de geformaliseerde werkwijze draagt zijn naam.

Augmented matrix en doel

Bij Gaussiaanse eliminatie zet men de coëfficiënten van het stelsel samen in een matrix, vaak als een augmented matrix waarin de rechterlid-waarden in een extra kolom zijn opgenomen. Met elementaire rijbewerkingen probeert men die matrix te transformeren naar een eenvoudiger vorm, zodat de oplossingen van het stelsel direct of met eenvoudige terugsubstitutie af te lezen zijn. Het uiteindelijke doel is meestal:

  • de rij-echelonvorm (REF): elke niet-nulrij heeft een leidend element (pivot) dat rechts van het leidende element van de rij erboven ligt, en rijen bestaande uit alleen nullen staan onderaan, of
  • de gereduceerde rij-echelonvorm (RREF): de matrix is in REF, elk leidend element is 1, en in elke kolom met een leidend element is dat de enige niet-nulwaarde.

Elementaire rijbewerkingen

Er zijn precies drie soorten elementaire rijoperaties die toegestaan zijn en die de oplossingsset van het stelsel onveranderd laten:

  • Type 1: wisselen van twee rijen (permutatie).
  • Type 2: vermenigvuldigen van een rij met een niet-nul getal (schalen).
  • Type 3: optellen van een veelvoud van de ene rij bij een andere rij (rijvervanging).

Algoritme in twee fasen

In praktische toepassingen worden vaak twee fasen onderscheiden:

  • Voorwaartse eliminatie (forward elimination): gebruik rijoperaties om nullen onder de pivots te krijgen en zo de matrix in rij-echelonvorm te brengen. Dit reduceert het stelsel tot een triangulaire vorm.
  • Terugsubstitutie of Gauss-Jordaan: bij terugsubstitutie lost men een triangulair stelsel op van beneden naar boven om de onbekenden te bepalen. Als men verder gaat tot gereduceerde rij-echelonvorm (eliminatie ook boven de pivots), noemt men dat Gauss-Jordaanse eliminatie — dan is elke pivot 1 en staan er alleen nullen boven en onder elke pivot, zodat oplossingen direct afleesbaar zijn.

Soorten oplossingen

Analyse van de gereduceerde matrix geeft inzicht in de oplossingen:

  • Als elke onbekende een pivotkolom heeft → unieke oplossing.
  • Als er minder pivots zijn dan onbekenden → er bestaan vrije variabelen en het stelsel heeft oneindig veel oplossingen (parametrische beschrijving mogelijk).
  • Als een rij van de vorm [0 ... 0 | b] met b ≠ 0 ontstaat → het stelsel is inconsistent (geen oplossing).
  • Voor homogene stelsels (rechterlid nul) is er altijd minstens de triviale oplossing; vrije variabelen leiden tot niet-triviale oplossingsruimten.

Praktische aspecten

  • Pivotkeuze en numerieke stabiliteit: bij rekenkundige toepassingen met drijvende-kommagetallen is het belangrijk om slimme pivotstrategieën te gebruiken, zoals partiële pivotering (rijverwisseling om de grootste absolute waarde als pivot te kiezen) of volledige pivotering, om foutaccumulatie te verminderen.
  • Complexiteit: de klassieke Gauss-eliminatie kost in het algemeen O(n^3) rekenbewerkingen voor een n×n-systeem; geheugen en rekenkosten maken het minder geschikt voor zeer grote systemen zonder speciale technieken.
  • Determinant en rang: voor een vierkante matrix is de determinant gelijk aan het product van de (eventueel geschaalde) pivotwaarden, met aandacht voor tekens bij rijwissels. De rang van de matrix blijkt uit het aantal pivots en bepaalt of oplossingen uniek zijn.
  • Alternatieve methoden: voor grote, gestructureerde of zeldzame systemen gebruikt men vaak iteratieve methoden (bijv. CG, GMRES) of speciale directe methoden die beter gebruikmaken van sparseness.

Samenvatting

Gaussiaanse eliminatie is een fundamentele methode voor het oplossen en analyseren van lineaire stelsels. Door toepassing van drie elementaire rijbewerkingen transformeert men de augmented matrix naar een eenvoudiger vorm (rij-echelon of gereduceerde rij-echelon). Of men kiest voor eenvoudige terugsubstitutie of volledige Gauss-Jordaanse reductie hangt af van de gewenste vorm en praktische overwegingen zoals bruikbaarheid voor verdere analyses, numerieke stabiliteit en rekentijd.