Metoda de eliminare Gauss

De la Wikipedia, enciclopedia liberă.
Salt la navigare Salt la căutare

În matematică , metoda de eliminare Gauss, adesea abreviat MEG, este un algoritm , care este numit după matematicianul german Carl Friedrich Gauss , folosit în algebra liniară pentru determinarea soluțiilor unui sistem de ecuații liniare pentru a calcula rangul ol ' invers de o matrice .

Algoritmul, prin aplicarea unor operațiuni specifice pe rândurile (sau coloanele) acestor operații elementare, reduce matricea într - o formă numită în pași. Matricea în formă esalon face calculul imediat al rangului matricei (care va fi egal cu numărul de pivot ) și rezoluția deosebit de simplă a sistemului liniar asociat.

O extindere a acestei metode, cunoscută ca metoda de eliminare Gauss-Jordan, de german matematicianul Wilhelm Jordan , reduce în continuare matricea într - o formă numită redusă în trepte, de asemenea , permițând pentru a calcula inversa.

Deși este de obicei atribuit Gauss, o primă aplicare a MEG apare deja în secolul al doilea î.Hr. în Jiuzhang SuanShu (nouă capitole privind artele matematice), matematicieni chinezi scrise în timpul dinastiei Han . [1]

Matrice completă a coeficienților

Un sistem de ecuații liniare :

din

poate fi reprezentat printr-o matrice:

respectiva matrice completă a coeficienților sistemului. Coeficienții sistemului liniar (și , prin urmare, elementele matricei) sunt elementele unui câmp cum ar fi cea a numerelor reale sau complexe Ultima coloană este coloana a termenilor cunoscuți.

Operații elementare pe rând

Operațiile elementare pe rând sunt operații care modifică o matrice într - una din următoarele moduri:

  • schimbul de două linii;
  • înmulțirea unui rând cu un număr diferit de zero;
  • adăugarea unui rând la un multiplu de alt rând.

Rândurile matricei sunt aici înțelese ca vectori ai spațiului vectorial Și, prin urmare, pot fi multiplicate cu un număr sau adăugat, pur și simplu, pe termen lung la termen.

Operațiile rând elementare au următoarea proprietate importantă: dacă este aplicată matricea completă (de exemplu, cu termeni cunoscuți) a unui sistem liniar, nu modifică setul de soluții ale sistemului. Cu alte cuvinte, modificările sistemului, dar soluțiile rămân neschimbate: un vector este soluția sistemului inițial dacă și numai dacă aceasta este soluția sistemului la care am aplicat operațiile elementare pe rând.

Pe de altă parte, este ușor de observat că operațiile elementare pe rând aplicate pe matrice completă o concordanță, în modul clasic de exprimare a unui sistem de ecuații (cu bretele), la următoarele operații de pe ecuațiile:

  • schimba ordinea scrierea a două ecuații;
  • înmulțirea ambele părți ale unei ecuații de către un număr de zero;
  • adaugă aceeași cantitate pentru ambele părți ale unei ecuații.

matrice Etapa

O matrice este o matrice în trepte având următoarea proprietate: primul element, din stânga, non-zero (dacă este prezent) a unui rând trebuie să fie într-o coloană la dreapta primului element nenul al rândului de mai sus și fiecare rând, toate nule, este sub fiecare linie nu toate nimic. De exemplu, matricea:

se reduce la trepte, în timp ce matricea:

nu este. Primul element diferit de zero pe fiecare linie (dacă este prezent) , se spune ca pivot. De exemplu, matricea în trepte descrisă mai sus are trei pivoți: 3, -1 și 5.

Algoritm Gauss

Algoritmul Gauss transforma orice matrice într-o matrice în trepte prin operații elementare pe rând. Acesta funcționează după cum urmează:

  1. Dacă primul rând are primul element nul, l schimb cu un rând care are primul element nenul. Dacă toate rândurile au primul element nul, du-te la pasul 3.
  2. Pentru fiecare rând cu primul element de nenulă, cu excepția primei ( ), Înmulþiri primul rând de la un coeficient ales, astfel încât suma între primul rând și are primul element nul (deci coeficient = ). A inlocui cu suma doar obținută.
  3. Acum, pe prima coloană toate cifrele, cu excepția, probabil, în primul rând, sunt nule. La acest punct înapoi la punctul 1 , având în vedere submatricea te prin ștergerea primul rând și prima coloană.

Rezultatul algoritmului nu este întotdeauna la fel, depinde de alegerile făcute. In timpul dezvoltarii, este important să se aibă în vedere obiectivul final este de a obține o matrice în trepte.

Un exemplu este prezentat aici:

Din matricea de pornire am multiplicat primul rând prin -1, a adăugat primul rând cu al doilea (scrierea rezultatului în al doilea rând, în scopul de a obține un zero în prima poziție), înmulțită rândul obținut prin -2 și în final, a adăugat la al treilea rând.

Această procedură, în scris și utilizat ca mai sus, pot fi folosite pentru a rezolva un sistem liniar (prin aplicarea operațiilor elementare pe rând la matricea completă a coeficienților sistemului) sau pentru a calcula rangul unei matrice generic.

În cazul în care doriți să utilizați această metodă pentru calcularea determinantului unei matrice pătrat generic , Trebuie să ne trianguleze mai întâi matricea obținerea unei matrice triunghiulare , Atunci vom face produsul elementelor de pe diagonala principală prin împărțirea sa , acesta este

unde este se obține în felul următor: înainte de începerea operației de reducere ea apare . Apoi, în funcție de operațiunea de reducere utilizată, următoarele apare:

  • dacă aveți schimbarea două linii;
  • dacă vă înmulțiți o linie de un număr non-zero;
  • (acesta este rămâne neschimbat) dacă adăugați un rând la un multiplu de un alt rând.

Pentru a calcula rangul și determinant, Gauss algoritmul poate fi utilizat și pe coloane.

Mai jos este un exemplu pentru calcularea determinantului pe aceeași matrice de mai sus:

Din matricea de pornire, primul rând înmulțit cu -1 a fost adăugat la al doilea rând pentru a reseta primul element (părăsește primul rând neschimbat); al treilea rând are deja un zero ca prim element de pe stânga, astfel al doilea rând, înmulțită cu -2 a fost adăugat la al treilea rând pentru a reseta al doilea element (lăsând al doilea rând neschimbat). Acum putem calcula cu ușurință determinantul care este produsul elementelor pe diagonala principală, adică este -12 și este același determinant ca și cea a matricei inițiale.

În exemplul anterior în schimb am obține -24 din produsul elementelor diagonale, adică -12 înmulțită cu un factor de 2: acest factor 2 este tocmai produsul factorii prin care rândurile au fost multiplicate în exemplul de mai sus.

Calculul determinantului prin acest procedeu este cel puțin sistemul pretențios cunoscut din punct de vedere computational (creșterea polinom în loc de factorial ) Și este aplicabilă numai atunci când dimensiunea matricelor devine mare.

Gauss-Jordan algoritm

După ce a redus matricei la trepte, este posibil să se utilizeze o versiune a algoritmului lui Gauss în direcția opusă, adică de jos în sus, pentru a obține o matrice care în fiecare coloană care conține un pivot are doar pivotul ca nenulă număr, (această matrice rezultată este numită o matrice scalat în formă redusă): utilizați doar fiecare rând, începând de la ultima, pentru a elimina toate diferite de zero deasupra pivotului acestui rând. În cele din urmă, din nou, cu mișcări Gaussian (prin linii multiplicatoare), putem obține că fiecare pivot are o valoare de 1.

De exemplu, realizarea algoritmului descris mai sus, avem:

soluţii de sistem

Algoritmul Gauss-Jordan transformă matricea coeficienților unui sistem într - o matrice făcută în felul următor: variabilele corespunzătoare coloanele care nu conțin pivot menționat sunt libere; orice altă variabilă apare într-o singură ecuație, și, prin urmare, poate fi exprimată ca o funcție a variabilelor libere și termeni cunoscuți. Spațiul soluție se obține prin atribuirea de valori arbitrare variabilelor libere și calcularea celelalte variabile în consecință.

Dacă termenii cunoscuți sunt toate la zero, setul de soluții este un subspațiu de . În caz contrar, se obține dintr - un subspațiu vectorial prin translație, care este un subspațiu afin .

Inversei unei matrice

Algoritmul Gauss-Jordan este, de asemenea, folosit pentru a găsi (atunci când există) inversul unei matrice. Acesta funcționează după cum urmează: fie o matrice inversabilă . Să considerăm matricea , cu linii și coloane, partea construită de partea și matricea de identitate . În acest moment, după efectuarea matricei în trepte cu una sau mai multe etape ale algoritmului Gauss, varianta Gauss-Jordan se aplică la nou pentru a ajunge pe partea stângă .

Atâta timp cât este inversabil, coloanele sale sunt independente, și, prin urmare, va conține toate pivotează la sfârșitul algoritmului. Deci, rezultatul va fi o matrice de tip . Ei bine matricea este doar inversul .

Următoarele exemple arată că inversul Și :

În prima etapă a multiplicat prima linie pentru -2, în al doilea, se adaugă la a doua linie a primului, al treilea a multiplicat a doua linie la -4, în etapa a patra se adaugă la primul rând din al doilea și în cele din urmă „ultimul pas este primul rând este împărțit la -2 , iar al doilea pentru 4. În acest fel, am pornit de la o matrice de și a ajuns la . Are asta este inversul .

Proprietate

  • Rangul unei matrice este egal cu numărul de pivotare al oricărei matrice obținută din acest pași prin operații elementare pe rând (sau coloana).
  • Setul de soluții ale sistemului liniar asociat cu matricea nu se schimbă prin operații elementare pe rând.
  • Pentru a extrage dintr-un set de vectori în un plafon de vectori subset independent , și , prin urmare , de asemenea , o bază de subspatiului generate de acestea, doar aplica algoritmul lui Gauss la matricea obținută prin considerarea ca acești vectori coloane ale matricei. Apoi luați numai acei vectori inițiali pe a cărui coloană în matricea corespunzătoare etapa (la finalul algoritmului) apare un pivot.

Notă

Bibliografie

  • (RO) Timothy Gowers, iunie Barrow-verde și Imre Leader, Princeton Companion la matematică, Princeton University Press, 08 septembrie 2008, p. 607, ISBN 978-0-691-11880-2 .
  • (RO) Bareiss, EH-Integer în multe etape Conservarea Gaussian Eliminare. Raport Argonne National Laboratory-ANL 7213, mai 1966.
  • (RO) Bareiss, Identitatea EH Sylvester și mai multe trepte , Integer-Conservarea Gaussian Eliminare. Matematica. Calculator. 22, 565-578, 1968.
  • (RO) Garbow, BS Integer-Conservarea Gaussian Eliminare. Programul P-158 (3600F), Matematică Aplicată Divizia, Argonne National Laboratory, 21 noiembrie 1966.
  • (RO) Gentle, JE "Gaussian Eliminare." §3.1 în Algebră liniară numerică pentru aplicații în Statistică. Berlin: Springer-Verlag, pp. 87-91, 1998.

Elemente conexe

Alte proiecte

linkuri externe

Controlul autorității GND (DE) 4156110-7
Matematica Portalul de matematică : accesați intrările Wikipedia care se ocupă de matematică