✖

Sisteme de ecuații liniare

Un sistem de ecuații liniare este o mulțime de m ecuații liniare cu n necunoscute. Rezolvând sistemul, obținem valorile pe care le putem pune în locul celor n necunoscute astfel încât toate ecuațiile să fie adevărate.

Definiția unui sistem de ecuații liniare

Te-ai întâlnit sigur cu un sistem simplu de tipul

a11x+a12y=b1a21x+a22y=b2

unde a11,…,a22, b1, b2 sunt numere reale date, iar x și y sunt necunoscute. Acest sistem de ecuații poate avea un număr diferit de soluții – niciuna, una singură sau o infinitate, la fel ca o ecuație liniară obișnuită. În cazul unui sistem de ecuații liniare, întrebarea despre numărul soluțiilor și aflarea lor este mai complicată decât în cazul unei simple ecuații liniare.

Definiția unui sistem de ecuații liniare:

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

unde a11, …, amn, b1, …, bm sunt numere reale. Acest sistem se numește sistem de m ecuații liniare cu n necunoscute x1, …, xn și coeficienții a11, …, amn. Un exemplu concret de sistem de ecuații liniare ar putea arăta astfel:

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

Dacă însă vreun aij este negativ, se obișnuiește să scriem sistemul în această formă:

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

Rămâne totuși adevărat că a12 = −2 și a24 = −10.

Scrierea matriceală

Putem transforma acest sistem în forma matriceală. Definim în total trei matrice: una pentru coeficienți, una pentru termenii liberi ai ecuațiilor și una pentru necunoscutele propriu-zise. Matricele vor arăta astfel:

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

Matricea A conține coeficienții ecuațiilor inițiale și se numește „matricea sistemului”, matricea b conține termenii liberi, iar matricea x conține necunoscutele. Atunci putem scrie întregul sistem de ecuații ca o singură ecuație

Ax=b

Această ecuație are sens, pentru că în stânga avem produsul a două matrice, iar în dreapta o matrice, iar dimensiunile se potrivesc. Putem încerca să calculăm membrul stâng. Înmulțim o matrice de tip m × n cu o matrice de tip n × 1, deci matricea rezultată va fi de tip m × 1, va avea m linii și o coloană. Până aici e în regulă.

Aplicăm acum algoritmul clasic de înmulțire a matricelor, începând cu prima linie a matricei A și cu prima coloană a matricei x (nu există altă coloană). Obținem:

a11x1+a12x2+…+a1nxn

Cu această sumă am obținut primul element al noii matrice, adică valoarea b1. Dacă punem totul în ecuație, avem:

a11x1+a12x2+…+a1nxn=b1

Am obținut exact prima ecuație din sistemul nostru. În același fel am continua pentru celelalte linii. Mai definim matricea extinsă a sistemului, adică matricea sistemului căreia îi adăugăm matricea b astfel:

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

În continuare ne va interesa rangul acestei matrice extinse a sistemului. Rangul matricei inițiale, neextinse, a sistemului îl notăm rang A, iar rangul matricei extinse rang (A|b). Ce relație există între aceste ranguri? Dacă adăugăm matricei A coloana b, pot apărea două situații: dacă această nouă coloană este o combinație liniară a coloanelor matricei A, rangul nu se schimbă. Dacă nu era o combinație liniară, rangul va fi mai mare cu unu.

Teorema Kronecker–Capelli

Rangul matricei extinse este important pentru noi, pentru că are loc o teoremă importantă, teorema Kronecker–Capelli: Sistemul de ecuații Ax = b are soluții dacă și numai dacă rangul matricei sistemului este egal cu rangul matricei extinse a sistemului: rang A = rang (A|b).

Hai să arătăm că această teoremă este adevărată cel puțin într-un sens – dacă x = β este soluție a ecuației Ax = b, atunci rang A = rang (A|b). Dacă aceste ranguri sunt egale, coloana b trebuie să fie o combinație liniară a coloanelor matricei A. Să scriem totul desfășurat. Matricea β are forma:

β=(β1β2⋮βn)

Iar noi afirmăm că, dacă în locul lui x1 punem β1 și așa mai departe, toate ecuațiile sistemului devin adevărate:

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

Rescriem puțin ecuațiile. Dacă ne uităm la prima coloană, vedem că avem mereu valorile a11, a21, …, am1, pe care le înmulțim cu aceeași valoare β1. Dacă separăm coloana într-o matrice aparte, putem da factor comun valoarea β1 în fața matricei și obținem:

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

Acum putem rescrie în acest fel întregul sistem:

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

Pentru ca această ecuație să fie adevărată, matricea din dreapta trebuie să fie o combinație liniară a matricelor coloană din stânga, iar valorile β1, …, βn trebuie să fie coeficienții ei.

Am demonstrat astfel că, dacă matricea β este soluție a sistemului, atunci membrul drept al sistemului, matricea b, trebuie să fie o combinație liniară a coloanelor matricei sistemului. Iar dacă o coloană din matricea b este o combinație liniară a coloanelor matricei A, atunci matricea A trebuie să aibă același rang ca matricea (A|b).

În același timp, dacă rang A = rang (A|b) = k și k este egal cu numărul de necunoscute, adică k = n, atunci sistemul are exact o soluție (este compatibil determinat). Dacă k < n, atunci sistemul are o infinitate de soluții (este compatibil nedeterminat), iar pentru a le scrie vom avea nevoie de n − k parametri.

Transformări elementare ale liniilor

Din ecuația precedentă se vede în continuare că, dacă efectuăm asupra matricei (A|b) transformări elementare ale liniilor, soluția sistemului de ecuații nu se schimbă. Dacă am aduna la a doua linie prima linie, am obține ecuația:

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

Dacă desfacem a doua linie, obținem:

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

Desfacem parantezele:

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

Acum doar grupăm puțin termenii:

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

În prima paranteză avem de fapt prima linie, iar în a doua paranteză vechea a doua linie. Știm că prima linie este egală cu valoarea b1, iar a doua linie este egală cu valoarea b2. Asta știm din ecuațiile inițiale. Deci în locul primei paranteze putem scrie b1, iar în locul celei de-a doua b2:

b1+b2=b1+b2.

La fel, putem înmulți o linie întreagă cu o constantă c ≠ 0 fără să schimbăm validitatea sistemului de ecuații. Atenție însă, nu putem face transformări pe coloane.

Sistem cu o singură soluție

Pe un exemplu vom arăta cum rezolvăm un sistem care are exact o soluție. Fie sistemul de ecuații:

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

Matricea A va arăta astfel:

A=(237381033)

Acum calculăm rangul acestei matrice:

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

Vedem că rang A = 3. De aici rezultă și că rang (A|b) = 3, deoarece o matrice de tip 3 × 4 poate avea rangul maxim trei. Sistemul are trei necunoscute, deci întregul sistem are exact o soluție. Cum o aflăm?

Metoda lui Gauss

Metoda lui Gauss (eliminarea gaussiană) constă în aducerea matricei extinse a sistemului la forma eșalon. Aducând matricea la forma eșalon, vom putea afla valoarea unei necunoscute. După ce găsim valoarea unei necunoscute, putem începe să înlocuim această valoare în celelalte ecuații.

Matricea (A|b) are forma:

(A|b)=(237473815003327)

O aducem la forma eșalon. Vom face aceleași transformări ca la pasul precedent, doar că de data aceasta le aplicăm și coloanei a patra:

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

Deoarece am făcut transformări elementare ale liniilor, această matrice descrie un sistem de ecuații echivalent – un sistem care are aceeași mulțime de soluții. Rescriem sistemul din forma matriceală în forma obișnuită, folosind coeficienții din ultima noastră matrice:

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

Pornim acum de la ultima ecuație. Ea ne spune că 78x3 = 312. Ce valoare trebuie să aibă x3? Este o ecuație liniară simplă, deci

78x3=312x3=31278x3=4

Avem valoarea primei necunoscute, x3. Înlocuim această valoare în a doua ecuație. Obținem ecuația:

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

Și, în sfârșit, înlocuim ambele valori calculate în prima ecuație:

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

Valoarea necunoscutei x1 este egală cu doi. În acest moment avem soluția completă. Putem scrie că x = (x1, x2, x3) = (2, 5, 4). După cum se vede, metoda lui Gauss este o metodă simplă, dar eficientă de rezolvare a unui sistem de ecuații liniare.

Infinitate de soluții

Dacă rang A = rang (A|b) = k și k este mai mic decât numărul de necunoscute, k < n, atunci sistemul are o infinitate de soluții (este compatibil nedeterminat). Să luăm un exemplu simplu:

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

Vedem că a doua ecuație este egală cu de trei ori prima ecuație. Dacă înmulțim prima ecuație cu minus trei și o adunăm la a doua linie, obținem:

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

În acest moment numărul de necunoscute este egal cu doi, n = 2, rangul matricei sistemului este egal cu unu, rang A = 1, iar rangul matricei extinse este și el egal cu unu, rang (A|b) = 1. Avem k < n, deci sistemul are o infinitate de soluții. Cum găsim aceste soluții?

Ne uităm la ecuații. A doua ecuație este adevărată pentru toți x1, x2, deci ne interesează doar prima ecuație. De fapt, toate perechile x1, x2 care verifică prima ecuație vor forma mulțimea tuturor soluțiilor sistemului dat. Așadar rezolvăm ecuația 4x1 + x2 = 5. Aducem ecuația la forma x2 = 5 − 4x1.

Ne ajutăm în continuare de un parametru. Alegem un parametru t. Acum vom încerca să exprimăm perechea x1, x2 cu ajutorul parametrului t, astfel încât să putem scrie că, de exemplu, toate perechile de forma (t, t + 2) sunt soluții ale sistemului. Deci perechi ca (1, 3) sau (14, 16). Rostul este să exprimăm valoarea unei necunoscute cu ajutorul alteia, ca să știm că, dacă x1 are o anumită valoare, atunci x2 va fi, de exemplu, de trei ori mai mare. Ca punct de plecare alegem x1.

Punem deci egalitatea t = x1. Dacă necunoscuta x1 are valoarea t, ce valoare va avea necunoscuta x2? Vedem din ecuația precedentă că va avea valoarea x2 = 5 − 4x1. Deci, dacă x1 = t, atunci x2 = 5 − 4t. Putem scrie astfel că toate perechile de forma (t, 5 − 4t) sunt soluții ale sistemului de ecuații.

Putem verifica. Dacă în locul lui t punem unu, avem: x1 = 1 și x2 = 5 − 4 · 1 = 1. Înlocuim aceste valori în ecuație: 4 · 1 + 1 = 5, ecuația se verifică.

Dacă în locul lui t punem cinci, avem: x1 = 5, x2 = 5 − 4 · 5 = −15. După înlocuirea în ecuație avem: 4 · 5 − 15 = 5. Se verifică.

Soluția generală și soluția particulară

Soluția unui sistem de ecuații scrisă cu ajutorul parametrilor se numește soluție generală. O soluție concretă, adică obținută când în locul parametrilor punem numere concrete, se numește soluție particulară.

Deci, din exemplul precedent: (x1, x2) = (t, 5 − 4t) este soluția generală a sistemului de ecuații, iar (1, 1) și (5, −15) sunt soluții particulare.

Dacă avem matricea sistemului A și matricea extinsă a sistemului (A|b) și rangurile lor sunt egale, rang A = rang (A|b) = k, atunci pentru a exprima soluția generală avem nevoie de n − k parametri. Observă că, dacă n = k, nu avem nevoie de niciun parametru. Acesta este cazul în care sistemul are exact o soluție – atunci nu este nevoie de niciun parametru.

Exemplu

Rezolvă următorul sistem de ecuații liniare:

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

Mai întâi calculăm rangul matricelor A și (A|b). Aducem matricea extinsă la forma eșalon, din care aflăm rangul ambelor matrice. Hai să facem transformările. (Înmulțim prima linie cu doi, apoi adunăm prima și a doua linie și scădem acest rezultat din a treia linie. Apoi scădem prima linie din a doua linie și, în final, împărțim la loc prima linie la doi.)

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

Din această ultimă matrice vedem deja că rang A = rang (A|b) = k = 2. Numărul de necunoscute este n = 4, deci vom avea nevoie de n − k = 2 parametri. Scriem ecuațiile după ultima matrice:

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

Din a doua ecuație obținem

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

Alegem primul parametru, t = x3. Atunci putem scrie că x1 = −6t. Putem înlocui asta în prima ecuație:

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

Introducem în joc al doilea parametru, s = x4. Atunci putem scrie că

x2=−12−7t+52s

Obținem astfel soluția generală a sistemului:

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

Putem încerca niște soluții particulare. De exemplu s = t = 0. Atunci obținem: (x1,x2,x3,x4)=(0,−12,0,0). După înlocuirea în ecuațiile inițiale:

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

obținem ecuațiile:

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

După eliminarea termenilor nuli și după efectuarea înmulțirilor cu fracții obținem:

1=12=24=4

Încercăm o altă soluție particulară: s = 2, t = 1. Atunci obținem:

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

Înlocuim în ecuații:

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

Simplificăm:

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

Și, în final, obținem:

1=12=24=4

Se pare că soluția generală obținută este corectă.

Capitole: Sisteme de ecuații liniare, Regula lui Cramer, Sisteme omogene