✖

Lineære ligningssystemer

Et lineært ligningssystem er en mængde af m lineære ligninger med n variable. Når vi løser systemet, finder vi værdier, som vi kan sætte ind for vores n variable, så alle ligningerne går op.

Definition af et lineært ligningssystem

Du har sikkert allerede mødt et simpelt system af typen

a11x+a12y=b1a21x+a22y=b2

hvor a11,…,a22, b1, b2 er givne reelle tal, og x og y er variable. Dette ligningssystem kan have et forskelligt antal løsninger — ingen, én eller uendelig mange, ligesom en almindelig lineær ligning. For et lineært ligningssystem er spørgsmålet om antallet af løsninger og om, hvordan man finder dem, mere indviklet end for en simpel lineær ligning.

Definition af et lineært ligningssystem:

a11x1+a12x2+…+a1nxn=b1a21x1+a22x2+…+a2nxn=b2⋮⋮⋮⋮⋮⋮⋮⋮⋮am1x1+am2x2+…+amnxn=bm

hvor a11, …, amn, b1, …, bm er reelle tal. Dette system kaldes et system af m lineære ligninger med n ubekendte x1, …, xn og koefficienterne a11, …, amn. Et konkret eksempel på et lineært ligningssystem kunne se sådan ud:

3x1+−2x2+4x3+5x4=15x1+4x2+2x3+−10x4=211x1+8x2+10x3+20x4=4

Men hvis et aij er negativt, er det almindeligt at skrive det på denne form:

3x1−2x2+4x3+5x4=15x1+4x2+2x3−10x4=211x1+8x2+10x3+20x4=4

Det gælder stadig, at a12 = −2 og a24 = −10.

Matrixform

Dette system kan vi skrive om på matrixform. Vi definerer i alt tre matricer, en for koefficienterne, en for ligningernes højresider og en for selve variablene. Matricerne ser så sådan ud:

A=(a11a12…a1na21a22…a2n⋮⋮⋱⋮am1am2…amn),b=(b1b2⋮bm),x=(x1x2⋮xn)

Matricen A indeholder koefficienterne fra de oprindelige ligninger, og vi kalder den »systemets matrix«. Matricen b indeholder ligningernes højresider, og matricen x indeholder variablene. Så kan vi skrive hele ligningssystemet som én ligning

Ax=b

Denne ligning giver mening, fordi vi på venstre side har et produkt af to matricer og på højre side en matrix, og deres dimensioner passer sammen. Vi kan prøve at gange venstre side ud. Vi ganger en matrix af typen m × n med en matrix af typen n × 1, så den resulterende matrix er af typen m × 1, den har m rækker og én søjle. Indtil videre er alt i orden.

Nu udfører vi den klassiske algoritme til multiplikation af matricer. Vi begynder med første række i matricen A og første søjle i matricen x (der er ingen anden søjle). Vi får:

a11x1+a12x2+…+a1nxn

Med denne sum har vi fået det første element i den nye matrix, altså værdien b1. Når vi sætter det hele ind i ligningen, har vi:

a11x1+a12x2+…+a1nxn=b1

Vi har fået netop den første ligning i vores ligningssystem. På samme måde ville vi fortsætte med de øvrige rækker. Vi definerer desuden totalmatricen, som er systemets matrix, hvortil vi har føjet matricen b sådan:

B=(A|b)=(a11a12…a1nb1a21a22…a2nb2⋮⋮⋱⋮⋮am1am2…amnbm)

Vi vil dernæst interessere os for rangen af denne totalmatrix. Rangen af systemets oprindelige matrix betegner vi rank(A), rangen af totalmatricen rank(A|b). Hvad er sammenhængen mellem dem? Hvis vi føjer søjlen b til matricen A, kan der opstå to situationer: Hvis den nye søjle er en linearkombination af søjlerne i matricen A, ændrer rangen sig ikke. Hvis den ikke var en linearkombination, bliver rangen én større.

Kronecker–Capellis sætning

Rangen af totalmatricen bliver vigtig for os, fordi der gælder en vigtig sætning, Kronecker–Capellis sætning: Ligningssystemet Ax = b har en løsning, hvis og kun hvis rangen af systemets matrix er lig med rangen af systemets totalmatrix: rank(A) = rank(A|b).

Vi prøver at vise, at sætningen gælder i det mindste i den ene retning — hvis x = β er en løsning til ligningen Ax = b, så er rank(A) = rank(A|b). Hvis de to matricer har samme rang, må søjlen b være en linearkombination af søjlerne i matricen A. Vi skriver det ud. Matricen β har formen:

β=(β1β2⋮βn)

Og vi påstår, at hvis vi sætter β1 ind for x1 osv., så bliver alle ligningerne i ligningssystemet opfyldt:

a11β1+a12β2+…+a1nβn=b1a21β1+a22β2+…+a2nβn=b2⋮⋮⋮⋮⋮⋮⋮⋮⋮am1β1+am2β2+…+amnβn=bm

Nu omskriver vi ligningerne lidt. Hvis vi ser på den første søjle, ser vi, at vi altid har værdierne a11, a21, …, am1, som alle ganges med den samme værdi β1. Hvis vi gør søjlen til en selvstændig matrix, kan vi sætte værdien β1 uden for matricen, og så får vi:

(a11β1a21β1⋮am1β1)=β1(a11a21⋮am1)

Nu kan vi på denne måde omskrive hele systemet:

β1(a11a21⋮am1)+β2(a12a22⋮am2)+…+βn(a1na2n⋮amn)=(b1b2⋮bm)

For at denne ligning kan gælde, må matricen på højre side være en linearkombination af søjlematricerne på venstre side, og værdierne β1, …, βn må være deres koefficienter.

Dermed har vi vist, at hvis matricen β skal være en løsning til ligningssystemet, så må systemets højreside, matricen b, være en linearkombination af søjlerne i systemets matrix. Og hvis søjlen fra matricen b er en linearkombination af søjlerne i matricen A, så må matricen A have samme rang som matricen (A|b).

Samtidig gælder det, at hvis rank(A) = rank(A|b) = k, og k er lig med antallet af ubekendte, altså k = n, så har systemet præcis én løsning. Hvis k < n, har systemet uendelig mange løsninger, og for at skrive dem op skal vi bruge n − k parametre.

Elementære rækkeoperationer

Af den foregående ligning kan vi desuden se, at hvis vi udfører elementære rækkeoperationer i matricen (A|b), ændrer ligningssystemets løsning sig ikke. Hvis vi lagde første række til anden række, ville vi få ligningen:

β1(a11a21+a11⋮am1)+β2(a12a22+a12⋮am2)+…+βn(a1na2n+a1n⋮amn)=(b1b2+b1⋮bm)

Når vi skriver den anden række ud, får vi:

β1(a21+a11)+β2(a22+a12)+…+βn(a2n+a1n)=b1+b2

Vi ganger parenteserne ud:

β1a21+β1a11+β2a22+β2a12+…+βna2n+βna1n=b1+b2

Nu bytter vi lidt rundt på leddene:

(β1a11+β2a12+…+βna1n)+(β1a21+β2a22+…+βna2n)=b1+b2

I den første parentes har vi i virkeligheden den første række og i den anden parentes den oprindelige anden række. Vi ved, at den første række er lig med værdien b1, og den anden række er lig med værdien b2. Det ved vi fra de oprindelige ligninger. Så i stedet for den første parentes kan vi skrive b1 og i stedet for den anden b2:

b1+b2=b1+b2.

På samme måde kan vi gange en hel række med en konstant c ≠ 0, uden at ligningssystemets gyldighed ændrer sig. Men pas på: Vi kan ikke udføre søjleoperationer.

System med én løsning

På et eksempel viser vi, hvordan man løser et system, der har præcis én løsning. Lad os tage dette ligningssystem:

2x1+3x2+7x3=473x1+8x2+x3=503x2+3x3=27

Matricen A ser således ud:

A=(237381033)

Nu beregner vi rangen af denne matrix:

(237381033)∼(69216162033)∼(692107−19033)∼(6921021−5702121)∼(6921021−570078)

Vi ser, at rank(A) = 3. Heraf følger også, at rank(A|b) = 3, fordi en matrix af typen 3 × 4 højst kan have rang tre. Systemet har tre ubekendte, så hele systemet har præcis én løsning. Hvordan finder vi den?

Gauss-elimination

Gauss-elimination går ud på at omforme systemets totalmatrix til trappeform. Når vi har omformet matricen til trappeform, kan vi finde værdien af én variabel. Når vi har fundet værdien af én variabel, kan vi begynde at indsætte den værdi i de andre ligninger.

Matricen (A|b) har formen:

(A|b)=(237473815003327)

Vi omformer den til trappeform. Vi udfører de samme operationer som i det foregående trin, men denne gang anvender vi dem også på den fjerde søjle:

(237473815003327)∼(6921141616210003327)∼(692114107−19−4103327)∼(6921141021−57−12302121189)∼(6921141021−57−1230078312)

Fordi vi udførte elementære rækkeoperationer, beskriver denne matrix et ækvivalent ligningssystem — et system, der har den samme løsningsmængde. Vi skriver systemet om fra matrixform til den sædvanlige form med koefficienterne fra vores sidste matrix:

6x1+9x2+21x3=14121x2−57x3=−12378x3=312

Nu tager vi udgangspunkt i den sidste ligning. Den siger, at 78x3 = 312. Hvilken værdi må x3 have? Det er en simpel lineær ligning, så

78x3=312x3=31278x3=4

Vi har værdien af den første variabel, x3. Den værdi indsætter vi i den anden ligning. Vi får ligningen:

21x2−57x3=−12321x2−57⋅4=−12321x2=−123+22821x2=105x2=10521x2=5

Og til sidst indsætter vi begge de beregnede værdier i den første ligning:

6x1+9x2+21x3=1416x1+9⋅5+21⋅4=1416x1+45+84=1416x1=12x1=2

Værdien af variablen x1 er lig med to. Nu har vi den fuldstændige løsning. Vi kan skrive, at x = (x1, x2, x3) = (2, 5, 4). Som du kan se, er Gauss-elimination en enkel, men effektiv metode til at løse et lineært ligningssystem.

Uendelig mange løsninger

Hvis rank(A) = rank(A|b) = k, og k er mindre end antallet af ubekendte, k < n, så har systemet uendelig mange løsninger. Lad os tage et simpelt eksempel:

4x1+x2=512x1+3x2=15

Vi ser, at den anden ligning er lig med tre gange den første ligning. Hvis vi ganger den første ligning med minus tre og lægger den til den anden række, får vi:

4x1+x2=50x1+0x2=0

Nu er antallet af variable lig med to, n = 2, rangen af systemets matrix er lig med én, rank(A) = 1, og rangen af totalmatricen er også lig med én, rank(A|b) = 1. Der gælder k < n, så systemet har uendelig mange løsninger. Hvordan finder vi disse løsninger?

Vi ser på ligningerne. Den anden ligning er opfyldt for alle x1, x2, så vi interesserer os kun for den første ligning. Alle par x1, x2, der opfylder den første ligning, udgør faktisk mængden af alle løsninger til det givne ligningssystem. Vi løser altså ligningen 4x1 + x2 = 5. Vi omformer ligningen til formen x2 = 5 − 4x1.

Dernæst hjælper vi os med en parameter. Vi vælger en parameter t. Nu prøver vi at udtrykke parret x1, x2 ved hjælp af parameteren t, så vi for eksempel kan skrive, at alle par af formen (t, t + 2) er løsninger til systemet. Altså par som (1, 3) eller (14, 16). Meningen er at udtrykke værdien af den ene variabel ved hjælp af en anden variabel, så vi ved, at når værdien af x1 er sådan og sådan, så er værdien af x2 for eksempel tre gange større. Som udgangspunkt vælger vi x1.

Vi sætter altså t = x1. Hvis variablen x1 har værdien t, hvilken værdi har så variablen x2? Det kan vi se af den foregående ligning, den har værdien x2 = 5 − 4x1. Altså hvis x1 = t, så er x2 = 5 − 4t. Vi kan dermed skrive, at alle par af formen (t, 5 − 4t) er løsninger til ligningssystemet.

Vi kan kontrollere det. Hvis vi sætter 1 ind for t, har vi: x1 = 1 og x2 = 5 − 4 · 1 = 1. Vi sætter disse værdier ind i ligningen: 4 · 1 + 1 = 5, ligningen passer.

Hvis vi sætter 5 ind for t, har vi: x1 = 5, x2 = 5 − 4 · 5 = −15. Når vi sætter ind i ligningen, har vi: 4 · 5 − 15 = 5. Det passer.

Fuldstændig og partikulær løsning

En løsning til et ligningssystem, der er skrevet ved hjælp af parametre, kalder vi den fuldstændige løsning. En konkret løsning, altså hvis vi sætter konkrete tal ind for parametrene, kalder vi en partikulær løsning.

Fra det foregående eksempel: (x1, x2) = (t, 5 − 4t) er den fuldstændige løsning til ligningssystemet, (1, 1) og (5, −15) er partikulære løsninger.

Hvis vi har systemets matrix A og systemets totalmatrix (A|b), og de har samme rang, rank(A) = rank(A|b) = k, så skal vi bruge n − k parametre til at udtrykke den fuldstændige løsning. Læg mærke til, at hvis n = k, skal vi bruge nul parametre. Det er tilfældet, hvor systemet har præcis én løsning — så er der ikke brug for nogen parameter.

Eksempel

Løs følgende lineære ligningssystem:

3x1−2x2+4x3+5x4=15x1−4x2+2x3+10x4=211x1−8x2+10x3+20x4=4

Først beregner vi rangen af matricerne A og (A|b). Vi omformer totalmatricen til trappeform, og heraf finder vi rangen af begge matricer. Så går vi i gang. (Vi ganger første række med to, lægger derefter første og anden række sammen og trækker dette resultat fra tredje række. Dernæst trækker vi første række fra anden række og til sidst dividerer vi første række med to igen.)

(A|b)=(3−24515−4210211−810204)∼(6−481025−4210211−810204)∼(6−481025−4210200000)∼(6−48102−10−60000000)∼(3−2451−10−60000000)

Af denne sidste matrix kan vi allerede se, at rank(A) = rank(A|b) = k = 2. Antallet af ubekendte er n = 4, så vi skal bruge n − k = 2 parametre. Vi skriver ligningerne op ud fra den sidste matrix:

3x1−2x2+4x3+5x4=1−x1−6x3=0

Af den anden ligning får vi

−x1−6x3=0−x1=6x3x1=−6x3

Vi vælger den første parameter, t = x3. Så kan vi skrive, at x1 = −6t. Det kan vi sætte ind i den første ligning:

3(−6t)−2x2+4t+5x4=1−18t−2x2+4t+5x4=1−14t−2x2+5x4=1−2x2=1+14t−5x4x2=−12−7t+52x4

Nu tager vi den anden parameter i brug, s = x4. Så kan vi skrive, at

x2=−12−7t+52s

Dermed får vi systemets fuldstændige løsning:

(x1,x2,x3,x4)=(−6t,−12−7t+52s,t,s)

Vi kan afprøve nogle partikulære løsninger. For eksempel s = t = 0. Så får vi: (x1,x2,x3,x4)=(0,−12,0,0). Når vi sætter ind i de oprindelige ligninger:

3x1−2x2+4x3+5x4=15x1−4x2+2x3+10x4=211x1−8x2+10x3+20x4=4

får vi ligningerne:

0−2⋅(−12)+0+0=10−4⋅(−12)+0+0=20−8⋅(−12)+0+0=4

Når vi fjerner nulleddene og ganger brøkerne ud, får vi:

1=12=24=4

Vi prøver en anden partikulær løsning: s = 2, t = 1. Så får vi:

x1=−6x2=−12−7+5=−52x3=1x4=2

Vi sætter ind i ligningerne:

−18−2(−52)+4+10=1−30−4(−52)+2+20=2−66−8(−52)+10+40=4

Vi reducerer:

−4+5=1−8+10=2−16+20=4

Og til sidst får vi:

1=12=24=4

Det ser ud til, at den fuldstændige løsning er rigtig.

Kapitler: Lineære ligningssystemer, Cramers regel, Homogene lineære ligningssystemer