Spring til indhold

Gauss-elimination

Fra Wikipedia, den frie encyklopædi

Gauss-elimination er en algoritme til at løse et lineært ligningssystem. Samles koefficientene til de ukendte i en matrix, kan denne omformes sådan at den bliver triangulær og har trappeform. Efter denne omskrivning kan de ukendte i ligningerne løses direkte. I Europa blev metoden systematisk benyttet af den tyske matematiker Carl Friedrich Gauss, men var kendt blandt kinesiske matematikere fra år 150 AD.[1][2] Gauss videreudviklede senere metoden sammen med geologen Wilhelm Jordan, sådan at matricen kunne omskrives på en reduceret trappeform. For mange problemer er dette en fordel. Dette gælder specielt ved meget store ligningssystemer, hvor numeriske metoder benyttes. Metoden kaldes da for Gauss-Jordan-reduktion.[3]

Den samme algoritme kan også benyttes til at beregne nulrummet og rangen for en matrix. Er matricen kvadratisk og regulær, kan Gauss-elimination også benyttes til at finde den tilhørende inverse matrix.

Gauss-algoritme

[redigér | rediger kildetekst]

Metoden kan enklest forklares ved at betragte et konkret eksempel med et lineært ligningssystem med tre variable x1, x2 og x3. Generelt kan det skrives som:

Er alle leddene på højre side af ligningerne lig nul, er ligningssystemet homogent. I det følgende vil det være et vigtigt specialtilfælde. Rækkefølgen af ligningerne spiller ingen rolle. Hver ligning kan multipliceres med et tal uden, at de mister noget af indholdet. Der sker heller ikke, hvis man adderer og eller trækker to ligninger fra hinanden. På den måde opstår følgende elementære operationer:

1: Man kan bytte rundt på to ligninger.
2: Man kan multiplicere alle led i en ligning med et tal forskelligt fra nul.
3: Man kan addere et multiplum af en ligning til en anden.

Med disse operationer kan man nu skrive ligningssystemets form om, sådan at de ukendte kan løses direkte. Dette kan gøres på forskellige måder i detaljer, men det svarer til at:

  • 1) I første ligning antages at koefficienten a11  ikke er lig nul. Skulle den være det, bytter man om på ligningerne.
  • 2) Multipliceres den første ligning med et passende tal og trækkes det fra de nedenstående ligningerne, kan den ukendte x1  elimineres i disse.
  • 3) På samme måde kan de andre ligninger anvendes til at eliminere x2  i ligningerne under. Og så videre med de resterende ukendte i de efterfølgende ligninger.

De tre givne ligninger vil på denne måde dermed blive omformet til et ækvivalent ligningssystem på formen:

Det oprindelige ligningssystem er dermed bragt over på triangulær form. Fra den sidste ligning kan x3  nu findes direkte. Settes dette ind i den anden ligning, giver den igen x2. Med disse to værdier findes så x1  fra den første ligningen, og alle de ukendte er blevet bestemt. Ligningssystemet er løst.

Mere generelle ligningssystemer kan skrives på kompakt form ved anvendelse af matricer. Med n  ukendte i m  ligninger, samler man de ukendte sammen i en kolonnevektor x = (x1, x2, ..., xn)T og ligeså de inhomogene led i kolonnevektoren b = (b1, b2, ..., bm)T sådan at ligningssystemet kan skrives som Ax = b. Her er nu A  en m × n matrix med elementer bestående af koefficientene aij til ligningssystemet. Tager man med de inhomogene led, kan alle koefficienter i ligningssystemet samles i den udvidede koefficientmatrix eller udvidede matrix.

som indeholder informationen fra hele ligningssystemet. De tilladte elementær operationer på de enkelte ligninger, svarer nu til operationer på rækker i denne udvidede matrix.

Som en enkel illustration af algoritmen, kan man betragte et ligningssystem med tre ligninger og 3 ukendte:

I den udvidede matrix kan x1  elimineres fra andre rækker ved ganske enkelt at tage differensen mellem denne og den første. Ligeledes kan x1  elimineres fra den næste række ved at trække 4 gange den første række fra den tredje. Det giver:


Her kan man nu skifte fortegn på alle leddene i anden og tredje række. Det svarer til at multiplicere hver af dem med -1. Til slut kan x2  elimineres fra sidste række ved at multiplicere anden række med 9 og trække resultatet fra tredje række. På den måde kommer man frem til:

Nederste række svarer nu til 2x3 = 6 eller x3 = 3. I anden række optræder x3  ikke sådan at man også fra denne række direkte kan udlæse at x2 = 2. Sættes disse to resultater ind i første ligning, giver den til slut at x1 = 1. Så løsningen til det fulde ligningssystem er givet ved punktet x0 = (1,2,3)T i vektorrummet R3. I dette enkle tilfælde kunne løsningen blive fundet direkte fra den inverse til matrix A ved at benytte at x = A-1b. For større ligningssystemer er det mere beregningskrævende end brug af Gauss-elimination.

Geometrisk kan man forstå denne løsningen ved at benytte at hver ligning beskriver et plan i et tredimensionelt rum. To af dem vil da have løsninger, som ligger på skæringslinjen mellem de tilsvarende planer. Det tredje plan vil så skære over denne linjen i et punkt som giver den entydige løsning. Hvis to af planerne er parallelle, vil der ikke være nogen løsninger og ligningssystemet er inkonsistent. Et andet, specielt tilfælde optræder, hvis den sidste ligning beskriver et plan, som går gennem skæringslinjen til de to første planer. Da vil alle punkterne på denne linje være en løsning. Man har da ubestemte løsninger, som svarer til et uendelig stort løsningsrum. I dette tilfælde er de tre ligninger ikke længere uafhængige af hinanden.

  1. ^ Calinger (1999), pp. 234–236
  2. ^ Timothy Gowers; June Barrow-Green; Imre Leader (8. september 2008). The Princeton Companion to Mathematics. Princeton University Press. s. 607. ISBN 978-0-691-11880-2.
  3. ^ Althoen, Steven C.; McLaughlin, Renate (1987), "Gauss–Jordan reduction: a brief history", The American Mathematical Monthly, Mathematical Association of America, 94 (2): 130-142, doi:10.2307/2322413, ISSN 0002-9890, JSTOR 2322413

Litteratur og refererede værker

[redigér | rediger kildetekst]

Eksterne henvisninger

[redigér | rediger kildetekst]