Inversa unei matrice

CAPITOLUL II

METODE DIRECTE DE CALCUL

A INVERSEI UNEI MATRICE

2.1. Metoda complemenţilor algebrici

Dacă A = (aij) Mn (C), pentru aflarea inversei matricei date se procedează astfel: se

calculează transpusa matricei şi apoi se înlocuiesc elementele sale prin complemenţii algebrici

ai lor. Matricea astfel obţinută se numeşte matrice adjunctă (sau asociată, sau reciprocă)

matricei A şi se notează A*.

Inversa matricei A este A-1 = A*.

Explicit, această metodă decurge astfel:

A =

Calculăm det A. Dacă d = det A 0, atunci matricea A este inversabilă.

Scriem matricea transpusă a matricei A, schimbând în matricea A liniile cu coloanele.

tA =

Scriem matricea adjunctă a matricei A, notată A*, înlocuind fiecare element al matricei

transpuse tA prin complementul său algebric, notat dij = (-1)i+jAij, i, j = , unde Aij este

minorul elementului aij din matricea tA (determinantul obţinut din tA prin eliminarea liniei i şi

a coloanei j).

Complemenţii algebrici se mai numesc şi cofactori.

Deci:

A* =

Împărţim toate elementele matricei A* prin valoarea determinantului d (d = detA) şi

matricea obţinută este A-1.

Avem: A-1 = A* sau, explicit:

A-1 = , d 0

Se verifică că are loc:

Pentru exemplificare, considerăm matricea:

A =

Să se determine inversa matricei A şi apoi să se verifice prin calcul direct relaţia:

Calculăm determinantul său şi obţinem:

d = detA = =

Determinantul său fiind nenul, matricea sa este inversabilă şi

Scriem transpusa:

tA =

Determinăm complemenţii algebrici ai elementelor matricei tA, de forma

Deci,

Se constată că:

- înmulţim linia 1 cu coloana 1 şi obţinem:









Analog, se arată că

2.2. Metoda ecuaţiei caracteristice

Această metodă poartă acest nume pentru că, în cele ce urmează, vom calcula inversa

matricei , nesingulară, cu ajutorul ecuaţiei ei caracteristice.

Determinantul caracteristic al matricei A, det (E – A), este de fapt un polinom de

gradul n în raport cu .

Ecuaţia det (E – A) = 0 se numeşte ecuaţia caracteristică a matricei A.

Utilizând coeficienţii polinomului, precum şi puterile A, A2, A3, ......, An-1 ale matricei

A, se poate obţine relativ uşor inversa

(2.2.1.) det (E – A) =

Atunci, conform identităţii lui Hamilton-Cayley, avem:

(2.2.2.) ,

Înmulţind (2.2.2.) cu obţinem:

(2.2.3.) ,

Din această ecuaţie determinăm pe :

(2.2.4.)

Deci, cunoscând coeficienţii polinomului caracteristic al matricei A şi calculând

puterile acestei matrice până la a n-1 putere, matricea A-1 se calculează după formula (2.2.4.)

Această metodă de inversare se utilizează în cazul matricelor de ordin nu prea mare,

pentru că în acest caz ridicarea la putere devine anevoioasă, fapt pentru care se propune o

metodă iterativă:

(2.2.5.) , unde A1 = A

reprezintă suma elementelor diagonalei matricei şi se numeşte de obicei

„urma” matricei .

Inversa căutată se determină prin formula:

(2.2.6.)

care prezintă avantajul că în loc să calculăm puterile matricei A: A, A2, ..., An-1, se calculează

şirul A, A1, ..., An, după (2.2.5.).

2.3. Metoda eliminării Gauss-Jordan

Definiţia 2.3.1. Fie matricea

Se numeşte pas Gauss-Jordan modificat transformarea elementelor matricei după

algoritmul următor:

(i)alegerea unui element nenul , numit pivot;

(ii) înlocuirea pivotului cu 1;

(iii) înlocuirea celorlalte elemente ale coloanei pivot cu 0;

(iv) împărţirea celorlalte elementei ale liniei pivot la pivot;

(v) calcularea celorlalte elemente ale matricei după regula determinantului de ordin 2,

considerând diagonala principală diagonala pivotului şi împărţire la elementul

pivot.

Observaţii (2.3.1.)

(i)dacă pivotul nu e în prima linie şi prima coloană, se efectuează schimbări de linii şi

coloane;

(ii) prin efectuarea unui pas Gauss-Jordan modificat, rangul matricei se păstrează.

Teorema 2.3.1. Dacă , atunci rang A = numărul maxim de paşi

Gauss-Jordan modificaţi ce se pot efectua cu elementele matricei A.

Observaţia 2.3.2. Dacă şi rang A = n, A este inversabilă, atunci există

A-1. Inversa lui A se găseşte efectuând cei n paşi Gauss-Jordan modificaţi cu elementele

matricei (A, In).

este matricea unitate de ordin n.

După efectuarea lor, în locul matricei A se obţine matricea unitate In, iar în locul

matricei unitate se obţine A-1.

Pentru exemplificare, considerăm matricea:

Să se afle inversa matricei şi apoi să se verifice prin calcul direct relaţia:

, unde I3 este matricea unitate de ordin 3.

A I3

1 1 1 1 0 01 2 3 0 1 01 4 9 0 0 1

I1 1 1 1 0 00 1 2 -1 1 00 3 8 -1 0 1

II1 0 -1 2 -1 00 1 2 -1 1 00 0 2 2 -3 1

III

1 0 0 32

5

0 1 0 -3 -4 -1

0 0 1 12

3

I3 A-1

Prin bordare asupra pivotului, obţinem, succesiv:

Iteraţia I:

; 1 : 1 = 1

; 2 : 1 = 2

; 3 : 1 = 3

; 8 : 1 = 8

; -1 : 1 = -1

; 1 : 1 = 1

; 0 : 1 = 0

; -1 : 1 = -1

; 0 : 1 = 0

; 1 : 1 = 1

Iteraţia II:

; 1 : 1 = 1; ; -1 : 1 = -1

; 2 : 1 = 2; ; -1 : 1 = -1

; 0 : 1 = 0; ; 0 : 1 = 0

; 2 : 1 = 2; ; 2 : 1 = 2

; -3 : 1 = -3; ; 1 : 1 = 1

Iteraţia III:

;

;

;

;

;

Verificare:

Analog se verifică şi că , ceea ce confirmă corectitudinea aflării lui .

2.4. Metoda de eliminare a lui Gauss

Fie matricea , detA 0 (matrice pătrată de ordinul n), , a cărei

inversă vrem s-o determinăm. Pentru aceasta se utilizează relaţia principală:

(2.4.1.) ,

In este matricea unitate de ordinul n:

Înmulţind matricele A şi A-1 se obţin n sisteme de ecuaţii liniare cu n necunoscute.

Coeficienţii acestor sisteme sunt elementele matricei A-1, xij, , iar termenii liberi sunt

elementele matricei unitate In.

Se observă că aceste sisteme se pot scrie sub următoarea formă:

(2.4.2.)

unde şi

şi

..........

şi

Cele n sisteme de n ecuaţii liniare, având aceeaşi matrice A, se pot rezolva simultan

prin metoda Gauss.

Metoda constă în a elimina la fiecare pas câte o necunoscută, obţinându-se un sistem

echivalent cu sistemul iniţial, sistem de n-1 ecuaţii cu n-1 necunoscute. La ultimul pas,

sistemul se transformă într-un sistem echivalent, format dintr-o singură ecuaţie cu o singură

necunoscută.

Se începe apoi determinarea necunoscutelor în sens invers, începând cu determinarea

ultimei necunoscute din sistemul intermediar.

Construim tabela formată din matricea A şi de coloanele termenilor liberi e1, e2, ..., en.

e1 ... ei ... en b

e1 a11 ... a1j ... a1n 1 ... 0 ... 0 b1

e2 a21 ... a2j ... a2n 0 1 0 ... 0 b2

... ... ... ... ... ... ... ... ... ... ... ...

ei ai1 ... aij ... ain 0 ... 1 ... ... bi

... ... ... ... ... ... ... ... ... ... ... ...

en an1 ... anj ... ann 0 ... 0 ... 1 bn

Descriem calculele din această tabelă, care se numeşte „tabela cu împărţire”.

Împărţim elementele primei linii la primul ei element a11 0, astfel:

, , ..., ,..., ,

Elementele astfel obţinute: 1, b12, ..., b1j, ..., b1n, v1, 0, 0, ..., 0 se trec în linia n+1 a

tabelei. Când a11 = 0 se face o permutare de linii sau coloane. Înmulţim elementele 1, b12, ...,

b1j, ..., b1n, v1, ..., 0 cu elementul conducător a21, iar produsele obţinute le scădem din

elementele corespunzătoare liniei a doua. Obţinem astfel linia întâi din primul sistem

intermediar Sn-1. Apoi înmulţim tot elementele 1, b12, ..., b1j, ..., b1n, v1, ..., 0 cu elementul

conducător a31 şi scăzând produsele obţinute din elementele liniei a treia se obţine astfel a

doua linie din sistemul Sn-1.

În general, înmulţind elementele 1, b12, ..., b1j, ..., b1n, v1, ..., 0 cu elementul conducător

ai1, scăzând produsele obţinute din linia i, obţinem linia i-1 din sistemul Sn-1.

Procedând analog cu sistemul intermediar Sn-1, găsim al doilea sistem intermediar Sn-2,

ş.a.m.d.

Continuând în acest mod obţinem un sistem Sn, ca urmare a împărţirii ecuaţiilor:

, respectiv cu . În tabelă vom nota coeficienţii ecuaţiilor care

formează acest sistem cu bik, unde bii = 1.

Observaţie:

Dacă în sistemul dat Sn facem transformarea , (i = 1, 2, ..., n) atunci obţinem

sistemul transformat cu necunoscutele . Se observă că cele două sisteme au aceeaşi

matrice A, ele deosebindu-se doar prin termenii liberi. Termenul liber din ecuaţia i a

sistemului se obţine făcând suma coeficienţilor şi termenului liber din ecuaţia i a sistemului

Sn şi se mai numeşte suma de control pentru ecuaţia i a lui Sn.

Aceste sume de control se trec în coloana notată . La fiecare pas în rezolvarea

problemei se obţine un sistem intermediar pentru sistemul ale cărui necunoscute sunt xij unde

i, j , cât şi pentru sistemul ale cărui necunoscute sunt . Construind aceste sisteme

intermediare se verifică dacă într-adevăr suma de control dintr-o linie coincide cu suma

elementelor, inclusiv şi termenul liber din linia respectivă. În caz afirmativ, se continuă să se

calculeze. Se procedează analog şi când se trece la determinarea necunoscutelor, adică se

trece la calcularea necunoscutei următoare doar dacă .

Numărul operaţiilor.

Pentru calculul matricei inverse A-1 avem de rezolvat cele n sisteme din formula

(2.4.2.)

Pentru transformarea matricei A într-o matrice triunghiulară efectuăm:

împărţiri

înmulţiri şi tot atâtea adunări.

Pe lângă acestea, pentru calculul termenilor liberi avem:

- la sistemul Sn : 1 împărţire, n-1 înmulţiri, 0 adunări;

- la sistemul Sn-1: 2 împărţiri, 2(n-2) înmulţiri, n-2 adunări;

- la sistemul Sn-2: 3 împărţiri, 3(n-3) înmulţiri, 2(n-3) adunări;

- .................................................................................................

- la sistemul Sn-i: i+1 împărţiri, (i+1)(n-i+1) înmulţiri, i(n-i-1) adunări.

- .................................................................................................

Obţinem:

împărţiri,

înmulţiri şi

adunări.

Deci, pentru a ajunge la ultimul Sn (cu matricea triunghiulară) efectuăm n2 împărţiri,

înmulţiri şi adunări.

Vom calcula numărul operaţiilor necesare pentru rezolvarea celor n sisteme Sn, care

diferă unele de altele doar prin termenii liberi. Rezolvarea unui astfel de sistem cere

înmulţiri. În ceea ce priveşte numărul adunărilor, la rezolvarea primului sistem cu prima

coloană a termenilor liberi, avem adunări, în cazul sistemului (cu a doua coloană

a termenilor liberi) avem nevoie de adunări. În mod analog, în cazul termenilor

liberi din coloana i, avem nevoie de adunări.

Astfel, pentru rezolvarea celor Sn sisteme sunt necesare: înmulţiri şi

adunări.

Deci, pentru calcularea matricei A prin această metodă se cer în total:

împărţiri;

înmulţiri;

adunări.

2.5. Algoritm pentru calculul inversei unei matrice

Se expune metoda pentru matricele pătrate A de ordin 3, cazul general fiind apoi

evident.

Dacă A-1 există, adică , aducem matricea A la forma eşalon E printr-un număr

finit de transformări elementare.

Cum , avem ,

În continuare, normalizăm pivoţii lui E, aplicându-i transformările elementare:

,

ceea ce revine la înmulţiri la stânga pe E cu matricele elementare, , .

Matricea E devine matricea , care are forma:

Efectând acum asupra lui transformările elementare

matricea devine matricea unitate .

Acum este evident rezultatul din teorema următoare:

Teorema 2.5.1.

Fie A o matrice pătrată de ordin n, astfel încât . Există un număr finit de matrice

elementare U1, U2, ..., Up de tip I, II sau II, astfel încât:

Mai mult,

Ţinând seama de modul în care se înmulţesc matricele partiţionate în blocuri, din

egalitatea din enunţul teoremei rezultă că:

Aşadar, efectuând asupra liniilor matricei transformările elementare

corespunzătoare matricelor elementare U1, U2, ... respectiv Up, se obţine în final, în cel de-al

doilea compartiment matricea A-1.

CAPITOLUL III

METODE ITERATIVE DE CALCUL

A INVERSEI UNEI MATRICE

3.1. Metoda lui Schultz pentru matrice nesingulare

Fiind dată o matrice nesingulară , propunem să-i determinăm inversa, .

Presupunem că am găsit printr-o metodă oarecare aproximarea iniţială pentru ,

adică pentru care norma matricei este strict subunitară .

(3.1.1.)

Norma utilizată, pe lângă proprietăţile (1.3.6.1.) – (1.3.6.4.) este supusă condiţiilor

(3.1.2.)

Dacă F = 0, atunci .

Pentru îmbunătăţirea preciziei, se utilizează procesul iterativ al lui Schultz, de

următoarea formă:

(3.1.3.) , k = 1, 2, ...

Ţinând cont că , vom face în (3.1.3.) următoarea substituţie:

Notând cu Fk eroarea corespunzătoare la pasul k,

Obţinem:

(3.1.4.) , k = 1, 2, 3, ...

În cele ce urmează, vom arăta că procedeul (3.1.4.) converge în ipoteza ,

unde norma canonică oarecare cu proprietăţile (3.1.2.)

Se arată că:

Propoziţia 3.1.1.

Are loc relaţia următoare:

Demonstraţie:

Vom folosi metoda inducţiei:

k = 1

k = 2

Presupunem adevărată relaţia pentru k-1.

şi demonstrăm pentru k:

Deci propoziţia este demonstrată

Dar . Din formulele corespunzătoare şi din proprietăţile normei rezultă:

Dar , deci şi din , obţinem:

De aici , , E fiind matricea unitate de ordinul n. Deci

procedeul (3.1.4.) este convergent.

Observaţia 3.1.1.

În particular, elementele matricei verifică inegalitatea , unde n este

ordinul matricei A şi , utilizarea m a normei arată că este verificată condiţia (3.1.1.),

deci (3.1.4.) converge.

Am demonstrat că , evaluând şi rapiditatea convergenţei

aproximaţiilor succesive.

Acum presupunem că (3.1.1.) este verificată şi vom face o evaluare a erorii în cazul în

care ne oprim la pasul k:

Cum , rezultă că:

Dar şi ,

(3.1.6.)

Din definiţia (3.1.1.) rezultă că procedeul (3.1.4.) are ordinul de convergenţă 2.

CAPITOLUL IV

APLICAŢII ALE INVERSEI UNEI MATRICE

4.1. Modelul economic al lui Leontief

În anul 1973, premiul Nobel pentru economie a fost atribuit lui W. Leontief

(economist american, originar din Rusia), care a elaborat o procedură de modelare matematică

numită analiză input-output, adecvată pentru un sistem economic complex, formată cu un

număr mare de sectoare între care există interacţiuni.

Pentru a simplifica prezentarea, presupunem că avem un sistem economic E, format cu

trei sectoare, în care sunt fabricate respectiv produsele P1, P2, P3. O parte din producţia

sistemului economic E este folosită în interiorul sistemului ca resurse, iar o altă parte este

rezervată pentru a onora o cerere externă.

Notăm cu aij numărul unităţilor din produsul Pi care sunt folosite pentru realizarea unei

unităţi din produsul Pj. Numerele aij cuantifică interacţiunile dintre sectoarele S1, S2, S3.

Matricea:

se numeşte matricea tehnologică sau matricea input-output.

Dacă xj este nivelul producţiei lui Pj (adică numărul de unităţi din produsul Pj

realizate), atunci vectorul 3-dimensional:

se numeşte vectorul producţie.

Produsul:

reprezintă consumul intern, pentru că reprezintă partea din producţie x1 a

lui S1 consumată în interiorul sistemului E, ş.a.m.d. Rezultă că este partea de producţie

disponibilă, care poate fi folosită pentru a onora o cerere externă.

Presupunem că în afara sistemului economic E există o comandă de cj unităţi din

produsul Pj, şi fie numit vectorul-comandă.

Problema fundamentală care se pune este de a planifica nivelul producţiei P de aşa

natură încât , ceea ce se mai scrie .

Dacă matricea este inversabilă, atunci înmulţind la stânga egalitatea

cu , se obţine .

4.2. Codificarea mesajelor

Inversa unei matrice are aplicaţii şi în criptografie, disciplină care are ca obiect

elaborarea metodelor de codificare a mesajelor pentru secretizarea acestora.

Pentru exemplificare, vom folosi matricea A şi inversa sa, A-1.

şi mesajul „SOSESC JOI”. Atribuim literelor alfabetului latin câte un număr natural pozitiv şi

blancului numărul 0, de exemplu:

A B C D E F G H I J ... O ... S ... Y Z

0 1 2 3 4 5 6 7 8 9 10 ... 15 ... 19 ... 25 26

Mesajul se împarte în blocuri de lungimi egale (în cazul nostru 3), luând în

consideraţie şi spaţiile libere (blancuri). Se adaugă, eventual la sfârşitul mesajului, câteva

blancuri pentru a obţine blocuri de lungimi egale:

SOSESCJOI

Semnele din blancuri sunt înlocuite cu numere asociate. Se formează o matrice

numerică M cu trei coloane şi tot atâtea linii câte blocuri are mesajul. Se înmulţeşte M cu A,

, matricea C fiind mesajul codificat care se transmite.

În cazul nostru:

M A C

La recepţie, pentru decodificare (decriptare), se efectuează produsul CA-1 şi se obţine

M, deoarece .

C A-1 M

4.3. Rezolvarea sistemelor de tip Cramer scrise sub formă matriceală

Definiţia 4.3.1.

Un sistem de n ecuaţii liniare cu n necunoscute se numeşte sistem de tip Cramer dacă

determinantul matricei este nenul.

Considerăm un sistem Cramer de trei ecuaţii cu trei necunoscute:

(S)

Dacă este matricea sistemului (S), şi , atunci

sistemul (S) se scrie sub formă matriceală astfel: (S) .

Cum , există A-1 şi înmulţind la stânga cu A-1 egalitatea de mai sus obţinem:

Rezultă că (S) are soluţie unică, anume .

S O S

E S C

J

I

S O S

E S C

J

I

4.4. Ecuaţii matriceale

Definiţia 4.4.1.

O ecuaţie de forma AX=B, unde A este o matrice pătrată inversabilă de ordin n, B este

o matrice de tip m x n, iar X este necunoscuta (tot o matrice), se numeşte ecuaţie matriceală.

Propoziţia 4.4.1.

Ecuaţia matriceală , unde este ireversibilă, iar , are o

soluţie unică de forma .

Demonstraţie:

Înmulţim la stânga ambii termeni cu :

Folosind asociativitatea înmulţirii matricelor, rezultă:

Observaţia 4.4.1.

Deoarece înmulţirea matricelor nu este comutativă, ecuaţia , unde A este o

matrice pătrată inversabilă, are soluţia unică .

Demonstraţie:

Înmulţim la dreapta ambii termeni cu :

Folosind asociativitatea înmulţirii matricelor obţinem:

Observaţia 4.4.2.

Ecuaţia , unde A şi B sunt matrice inversabile, are soluţia unică

.

Demonstraţie:

Înmulţind la stânga ambii termeni ai ecuaţiei cu şi la dreapta ambii termeni ai

ecuaţiei cu , obţinem:

Folosind asociativitatea înmulţirii matricelor, rezultă:

Inversa unei matrice

Documents

Transcript of Inversa unei matrice