Linjära ekvationssystem
Ett linjärt ekvationssystem består av m linjära ekvationer med n variabler. När vi löser systemet får vi de värden som vi kan sätta in för våra n variabler så att alla ekvationer stämmer.
Definition av linjärt ekvationssystem
Du har säkert redan stött på ett enkelt system av typen
där a11,…,a22, b1, b2 är givna reella tal och x och y är variabler. Det här systemet kan ha olika många lösningar – inga, en eller oändligt många, precis som en vanlig linjär ekvation. För ett linjärt ekvationssystem är frågan om antalet lösningar och hur vi hittar dem svårare än för en enkel linjär ekvation.
Definition av ett linjärt ekvationssystem:
där a11, …, amn, b1, …, bm är reella tal. Systemet kallas då ett system av m linjära ekvationer med n obekanta x1, …, xn och koefficienterna a11, …, amn. Ett konkret exempel på ett linjärt ekvationssystem kan se ut så här:
Om någon koefficient aij är negativ brukar vi i stället skriva systemet så här:
Det gäller ändå fortfarande att a12 = −2 och a24 = −10.
Matrisform
Vi kan skriva om systemet på matrisform. Vi definierar tre matriser: en för koefficienterna, en för ekvationernas högerled och en för själva variablerna. Matriserna ser ut så här:
Matrisen A innehåller koefficienterna i de ursprungliga ekvationerna och kallas ”koefficientmatris”, matrisen b innehåller ekvationernas högerled och matrisen x innehåller variablerna. Då kan vi skriva hela ekvationssystemet som en enda ekvation
Den här ekvationen är meningsfull, eftersom vi har en produkt av två matriser i vänsterledet och en matris i högerledet, och deras dimensioner stämmer. Vi kan prova att multiplicera vänsterledet. Vi multiplicerar en matris av typen m × n med en matris av typen n × 1, så den resulterande matrisen blir av typen m × 1, alltså med m rader och en kolumn. Än så länge stämmer allt.
Nu använder vi den vanliga algoritmen för matrismultiplikation. Vi börjar med den första raden i matrisen A och den första kolumnen i matrisen x (det finns ingen annan kolumn). Vi får:
Med den här summan har vi fått det första elementet i den nya matrisen, alltså värdet b1. Om vi sätter in allt i ekvationen får vi:
Vi har fått precis den första ekvationen i vårt ekvationssystem. På samma sätt fortsätter vi med de övriga raderna. Vi definierar också systemets totalmatris, som är koefficientmatrisen med matrisen b tillagd så här:
Härnäst är vi intresserade av rangen hos systemets totalmatris. Rangen hos den ursprungliga matrisen betecknar vi rang(A) och rangen hos totalmatrisen rang(A ∣ b). Hur hänger de här rangerna ihop? Om vi lägger till kolumnen b till matrisen A kan två saker hända: om den nya kolumnen är en linjärkombination av kolumnerna i matrisen A ändras inte rangen. Om den inte är en linjärkombination ökar rangen med ett.
Rouché–Capellis sats
Totalmatrisens rang är viktig för oss på grund av en viktig sats, Rouché–Capellis sats: Ekvationssystemet Ax = b har en lösning om och endast om rangen hos koefficientmatrisen är lika med rangen hos systemets totalmatris: rang(A) = rang(A ∣ b).
Vi visar att satsen gäller åtminstone i ena riktningen – om x = β är en lösning till ekvationen Ax = b, så är rang(A) = rang(A ∣ b). Om rangerna är lika måste kolumnen b vara en linjärkombination av kolumnerna i matrisen A. Vi skriver ut det. Matrisen β har formen:
Vi påstår att om vi sätter in β1 för x1 och så vidare, så stämmer alla ekvationer i systemet:
Nu skriver vi om ekvationerna lite. Om vi tittar på den första kolumnen ser vi att där alltid står a11, a21, …, am1, som alla multipliceras med samma värde β1. Om vi skriver kolumnen som en egen matris kan vi bryta ut värdet β1 framför matrisen, och då får vi:
På det här sättet kan vi nu skriva om hela systemet:
För att den här ekvationen ska gälla måste matrisen i högerledet vara en linjärkombination av kolumnmatriserna i vänsterledet, och värdena β1, …, βn måste vara koefficienterna i den.
Därmed har vi visat att om matrisen β ska vara en lösning till ekvationssystemet, så måste systemets högerled, matrisen b, vara en linjärkombination av kolumnerna i koefficientmatrisen. Och om kolumnen b är en linjärkombination av kolumnerna i matrisen A, så måste matrisen A ha samma rang som matrisen (A ∣ b).
Dessutom gäller att om rang(A) = rang(A ∣ b) = k och k är lika med antalet obekanta, alltså k = n, så har systemet exakt en lösning. Om k < n har systemet oändligt många lösningar, och för att skriva dem behöver vi n − k parametrar.
Elementära radoperationer
Av den föregående ekvationen ser vi också att om vi gör elementära radoperationer på matrisen (A ∣ b) ändras inte ekvationssystemets lösning. Om vi skulle addera den första raden till den andra raden skulle vi få ekvationen:
När vi skriver ut den andra raden får vi:
Vi multiplicerar in i parenteserna:
Nu ordnar vi om termerna lite:
I den första parentesen har vi egentligen den första raden och i den andra parentesen den ursprungliga andra raden. Vi vet att den första raden är lika med värdet b1 och den andra raden är lika med värdet b2. Det vet vi från de ursprungliga ekvationerna. Vi kan alltså skriva b1 i stället för den första parentesen och b2 i stället för den andra:
På samma sätt kan vi multiplicera en hel rad med en konstant c ≠ 0 utan att ekvationssystemet slutar gälla. Observera dock att vi inte får göra kolumnoperationer.
System med exakt en lösning
Med ett exempel visar vi hur vi löser ett system som har exakt en lösning. Vi har ekvationssystemet:
Matrisen A ser ut så här:
Nu bestämmer vi rangen hos den här matrisen:
Vi ser att rang(A) = 3. Av det följer också att rang(A ∣ b) = 3, eftersom en matris av typen 3 × 4 högst kan ha rang tre. Systemet har tre obekanta, så hela systemet har exakt en lösning. Hur hittar vi den?
Gausselimination
Gausselimination innebär att vi skriver om systemets totalmatris på trappstegsform. När vi har matrisen på trappstegsform kan vi bestämma värdet på en variabel. När vi har hittat värdet på en variabel kan vi sätta in det i de andra ekvationerna.
Matrisen (A ∣ b) har formen:
Vi skriver om den på trappstegsform. Vi gör samma operationer som i föregående steg, men den här gången använder vi dem också på den fjärde kolumnen:
Eftersom vi har gjort elementära radoperationer beskriver den här matrisen ett ekvivalent ekvationssystem – ett system som har samma lösningsmängd. Vi skriver om systemet från matrisform till vanlig form med koefficienterna från vår sista matris:
Nu utgår vi från den sista ekvationen. Den säger att 78x3 = 312. Vilket värde måste x3 ha? Det är en enkel linjär ekvation, så
Vi har fått värdet på den första variabeln, x3. Det värdet sätter vi in i den andra ekvationen. Vi får ekvationen:
Till sist sätter vi in båda värdena i den första ekvationen:
Värdet på variabeln x1 är två. Nu har vi hela lösningen. Vi kan skriva att x = (x1, x2, x3) = (2, 5, 4). Som du ser är Gausselimination en enkel men effektiv metod för att lösa ett linjärt ekvationssystem.
System med oändligt många lösningar
Om rang(A) = rang(A ∣ b) = k och k är mindre än antalet obekanta, k < n, så har systemet oändligt många lösningar. Vi tittar på ett enkelt exempel:
Vi ser att den andra ekvationen är tre gånger den första ekvationen. Om vi multiplicerar den första ekvationen med minus tre och adderar den till den andra raden får vi:
Nu är antalet variabler två, n = 2, rangen hos koefficientmatrisen är ett, rang(A) = 1, och rangen hos totalmatrisen är också ett, rang(A ∣ b) = 1. Det gäller att k < n, så systemet har oändligt många lösningar. Hur hittar vi dem?
Vi tittar på ekvationerna. Den andra ekvationen är uppfylld för alla x1, x2, så vi behöver bara bry oss om den första ekvationen. Alla par x1, x2 som uppfyller den första ekvationen bildar nämligen lösningsmängden till hela ekvationssystemet. Vi löser alltså ekvationen 4x1 + x2 = 5. Vi skriver om ekvationen som x2 = 5 − 4x1.
Sedan tar vi hjälp av en parameter. Vi väljer parametern t. Nu försöker vi uttrycka paret x1, x2 med hjälp av parametern t, så att vi till exempel kan skriva att alla par på formen (t, t + 2) är lösningar till systemet. Alltså par som (1, 3) eller (14, 16). Poängen är att uttrycka värdet på den ena variabeln med hjälp av en annan variabel, så att vi vet att när x1 har ett visst värde, så har x2 till exempel ett värde som är tre gånger större. Som utgångspunkt väljer vi x1.
Vi sätter alltså t = x1. Om variabeln x1 har värdet t, vilket värde har då variabeln x2? Det ser vi av den förra ekvationen: den har värdet x2 = 5 − 4x1. Om x1 = t är alltså x2 = 5 − 4t. Vi kan därför skriva att alla par på formen (t, 5 − 4t) är lösningar till ekvationssystemet.
Vi kan kontrollera det. Om vi sätter in ett för t får vi: x1 = 1 och x2 = 5 − 4 · 1 = 1. Vi sätter in de här värdena i ekvationen: 4 · 1 + 1 = 5, ekvationen stämmer.
Om vi sätter in fem för t får vi: x1 = 5, x2 = 5 − 4 · 5 = −15. Efter insättning i ekvationen får vi: 4 · 5 − 15 = 5. Det stämmer.
Allmän och partikulär lösning
En lösning till ett ekvationssystem som är skriven med hjälp av parametrar kallas allmän lösning. En konkret lösning, där vi alltså har satt in konkreta tal för parametrarna, kallas partikulär lösning.
I det förra exemplet är alltså (x1, x2) = (t, 5 − 4t) den allmänna lösningen till ekvationssystemet, medan (1, 1) och (5, −15) är partikulära lösningar.
Om vi har koefficientmatrisen A och totalmatrisen (A ∣ b) och deras ranger är lika, rang(A) = rang(A ∣ b) = k, så behöver vi n − k parametrar för att skriva den allmänna lösningen. Lägg märke till att om n = k behöver vi noll parametrar. Det är fallet när systemet har exakt en lösning – då behövs ingen parameter.
Exempel
Lös följande linjära ekvationssystem:
Först bestämmer vi rangen hos matriserna A och (A ∣ b). Vi skriver om totalmatrisen på trappstegsform, och av den ser vi rangen hos båda matriserna. Vi sätter igång. (Vi multiplicerar den första raden med två, adderar sedan den första och den andra raden och subtraherar resultatet från den tredje raden. Sedan subtraherar vi den första raden från den andra raden, och till sist dividerar vi den första raden med två för att få tillbaka den.)
Av den här sista matrisen ser vi redan att rang(A) = rang(A ∣ b) = k = 2. Antalet obekanta är n = 4, så vi behöver n − k = 2 parametrar. Vi skriver ekvationerna enligt den sista matrisen:
Av den andra ekvationen får vi
Vi väljer den första parametern, t = x3. Då kan vi skriva att x1 = −6t. Det kan vi sätta in i den första ekvationen:
Nu tar vi in den andra parametern, s = x4. Då kan vi skriva att
Därmed får vi systemets allmänna lösning:
Vi kan prova några partikulära lösningar. Till exempel s = t = 0. Då får vi:
får vi ekvationerna:
När vi tar bort nollorna och multiplicerar ihop bråken får vi:
Vi prövar en partikulär lösning till: s = 2, t = 1. Då får vi:
Vi sätter in i ekvationerna:
Vi förenklar:
Och till sist får vi:
Det verkar som att den allmänna lösningen stämmer.
Kapitel: Linjära ekvationssystem, Cramers regel, Homogena linjära ekvationssystem