V matematice je Gaussova eliminace (nazývaná také redukce řádků) metoda používaná k řešení soustav lineárních rovnic. Je pojmenována po Carlu Friedrichu Gaussovi, slavném německém matematikovi, který o této metodě psal, ale nevynalezl ji. Metoda převádí soustavu rovnic na maticový tvar a pomocí elementárních řádkových operací postupně zjednodušuje rozšířenou matici, až je možné snadno určit řešení soustavy (pokud existuje).
Základní princip a elementární řádkové operace
Pro použití Gaussovy eliminace se koeficienty rovnic zapíší do tzv. rozšířené matice (augmented matrix), kde poslední sloupec obsahuje pravé strany rovnic. Během eliminace se používají tři typy elementárních řádkových operací, které nezmění množinu řešení soustavy:
- Prohození dvou řádků (swapping řádků).
- Násobení řádku nenulovým číslem (škálování řádku).
- Přičtení násobku jednoho řádku k druhému (řádková kombinace).
Algoritmus (postup krok za krokem)
Obvyklý postup Gaussovy eliminace lze shrnout takto:
- Sestavte rozšířenou matici soustavy.
- Pro každý sloupec (od prvního po předposlední, pokud řešíte n neznámých, projděte n sloupců):
- Vyberte pivot (neboli hlavní prvek) v aktuálním sloupci. Pokud je pivot nulový, proveďte prohození s nižším řádkem, který obsahuje nenulový prvek (tzv. částečné pivotování).
- Pomocí operace škálování upravte pivot tak, aby byl 1 (volitelné; nemusí se dělat hned).
- Eliminujte všechny prvky pod pivotelementem (přičtením vhodných násobků pivotního řádku k řádkům pod ním), čímž vytvoříte nuly pod pivotem.
- Po provedení těchto kroků pro všechny sloupce získáte matici v řádkovém echelonovém tvaru (REF). Pokud chcete, můžete pokračovat v eliminaci i nad pivoty, čímž získáte redukovaný řádkový echelonový tvar (RREF) — tomu se říká Gaussova–Jordanova eliminace.
- Ze získaného (redukovaného) tvaru pak zpětnou substitucí (pokud máte REF) určíte hodnoty neznámých. V RREF jsou hodnoty obvykle ihned čitelné.
Řádkový echelon a redukovaný řádkový echelon
Matice je v řádkovém echelonovém tvaru (REF), pokud platí:
- Všechny nulové řádky (pokud existují) jsou pod všemi nenulovými řádky.
- První nenulový prvek (pivot) v každém nenulovém řádku je vpravo od pivotu v řádku nad ním.
Matice je v redukovaném řádkovém echelonovém tvaru (RREF), pokud je navíc každý pivot roven 1 a je jediným nenulovým prvkem ve svém sloupci (tj. nad i pod pivotem jsou nuly). Gaussova eliminace může skončit v REF (běžný postup) nebo pokračovat do RREF (Gaussova–Jordanova eliminace).
Možné výsledky soustavy
Po eliminaci lze rozlišit tři základní případy:
- Jednoznačné řešení: Každá proměnná je určena — matice má plný počet pivotů (rank = počet proměnných).
- Infinitně mnoho řešení: Existují volné proměnné (sloupce bez pivotu), které lze parametry vyjádřit. Rank je menší než počet proměnných, ale rozšířená matice nekontradikuje (žádný řádek typu [0 ... 0 | c] s c ≠ 0).
- Žádné řešení (inkonzistence): Po eliminaci vznikne řádek tvaru [0 0 ... 0 | c] s c ≠ 0, což znamená, že soustava je neslučitelná.
Krátký příklad (2×2)
Řešme soustavu:
2x + y = 5
4x − 6y = −2
Sestavíme rozšířenou matici:
[ 2 1 | 5 ]
[ 4 −6 | −2 ]
Eliminace: použijeme R2 ← R2 − 2·R1:
[ 2 1 | 5 ]
[ 0 −8 | −12 ]
Teď z druhého řádku dostaneme y = (−12)/(−8) = 3/2. Dosadíme do prvního řádku: 2x + (3/2) = 5 → 2x = 7/2 → x = 7/4.
Řešení: x = 7/4, y = 3/2.
Numerické aspekty a složitost
Ve výpočetní praxi je důležité dbát na stabilitu a přesnost:
- Pivotování: Pokud je pivot příliš malý (blízko nule), aritmetické chyby mohou způsobit velké odchylky výsledku. Proto se často používá částečné pivotování (prohození s řádkem, který má v aktuálním sloupci největší absolutní hodnotu pivotu) nebo úplné pivotování (prohození řádků i sloupců).
- Škálování: Normalizace řádků/sloupců může zlepšit numerickou stabilitu.
- Složitost: Standardní Gaussova eliminace má výpočetní složitost O(n^3) pro n×n matici (počítáno v základních aritmetických operacích).
Použití
Gaussova eliminace je základní nástroj v lineární algebře a numerických výpočtech — používá se například při řešení soustav v technických aplikacích, při výpočtu inverzních matic, ve výpočtech determinantů, v analýze sítí, při interpolacích a v mnoha algoritmech numerické lineární algebry.
Shrnutí: Gaussova eliminace je systematická metoda přetvoření rozšířené matice soustavy do tvaru, ze kterého lze určit, zda soustava má jediné řešení, nekonečně mnoho řešení nebo je neslučitelná. Při praktické implementaci je vhodné dbát na pivotování kvůli číselné stabilitě.