Matrici e sistemi lineari
Definizione
Siano m, n € N*. Una matrice di tipo (m,n) [o anche una matrice m x n] ad elementi (detti entrate) in un campo K è una tabella di mn elementi di K disposti su m righe e n colonne. Mmn(K) è l’insieme delle matrici di tipo (m,n) ad elementi in K.
Se m=n, scriverò Mm(K) e dirò che è l’insieme delle matrici quadrate di ordine m.
Sia A= (aiy) € Mmn(K) con i = 1,…, m e y=1,…,n.
- aiy si dice elemento di posto i,y di A
- (ai1, ai2, ai3, … , ain) € K si dice riga i-esima di A, la colonna è detta, invece colonna y-esima di A
Vettore riga e vettore colonna
Chiameremo Ri la i-esima riga e Ri è detto vettore riga. Chiameremo, invece, Cj la j-esima colonna e Cj è detto vettore colonna.
Matrice trasposta
Sia A = (aij) € Mmn(K). La matrice AT € Mnm(K) che si ottiene scambiando le righe con le colonne di A si dice trasposta di A.
Osservazione: (AT)T = A cioè la trasposta di una trasposta è uguale alla matrice di partenza.
Matrice simmetrica
Sia A = (aij) € Mn(K). A si dice simmetrica se la trasposta AT = A, cioè aij = aji ∀ i ≠ j.
Diagonali e diagonale
Sia A € Mn(K), gli elementi a11, a22, a33, …, ann si dicono diagonali e formano la diagonale principale di A.
Sottomatrice
Sia A € Mmn(K), una sottomatrice di A è una matrice che si ottiene da A cancellando un numero intero 0 ≤ r ≤ m di righe e un numero intero 0 ≤ s ≤ n di colonne. Le sottomatrici quadrate di A si dicono minori di A.
Triangolare superiore e inferiore, diagonale e identità
Sia A € Mn(K).
Una matrice A € Mn(K) = (aij) t.c. aij = 0 per i > j si dice triangolare superiore.
A si dice triangolare inferiore se aij = 0 per i < j.
A si dice diagonale se aij = 0 per i ≠ j.
La matrice In = (δij) € Mn(K) t.c. δij = {1 se i = j, 0 se i ≠ j} si dice matrice unità o identità di ordine n.
Operazioni con matrici
N.B. V A = (aij) € Mmn(K), B = (bij) € Mmn(K) e V λ € K
Poniamo
- A+B = (aij + bij) € Mmn(K)
- λA = (λaij) € Mmn(K)
Con queste operazioni Mmn(K) è uno spazio vettoriale su K. Infatti:
- La somma è associativa e commutativa
- Lo 0 di Mmn(K) è la matrice che ha tutti gli elementi nulli
- V A = (aij) € Mmn(K) l’opposta è -A = (-aij) € Mmn(K).
- E sono verificate le altre quattro proprietà.
Prodotto tra matrici
Prodotto tra vettori riga e tra vettori colonna
Consideriamo una matrice R € M1p(K) (=K1) e una seconda matrice che ha tante righe quante sono le colonne di R: C € Mp1(K) (=K1).
Definiamo R C = (C1) = (∑t=1 a1t bt1) = (a11 b11 + a12 b21 + a13 b31 + … + a1p bp1).
N.B. per poter calcolare il prodotto di due matrici, il numero delle colonne della prima matrice deve essere necessariamente uguale al numero di righe della seconda matrice. La matrice risultante avrà quindi tante righe quante sono le righe della prima matrice e tante colonne quante sono le colonne della seconda matrice.
Attenzione!! Il prodotto righe per colonne non è commutativo!
Prodotto AB e BA
Sia A = (aij) € Mmn(K) e B = (bij) € Mnp(K). È quindi possibile calcolare il prodotto AB poiché le colonne di A sono tante quante le righe di B, ma non posso fare BA perché le colonne di B non sono tante quante le righe di A.
Di conseguenza AB ≠ BA.
Anche nel caso particolare in cui m = n e quindi A = (aij) € Mmn(K) e B = (bij) € Mnm(K) si ha che AB ≠ BA.
Inoltre anche nel caso ancora più particolare in cui si considerano due matrici quadrate dello stesso ordine e quindi A = (aij) € Mn(K) e B = (bij) € Mn(K) si ha che AB ≠ BA.
Proprietà del prodotto tra matrici
Siano A,B € Mmn(K) e C,D € Mnp(K).
- Allora la somma di A e B è distributiva rispetto al prodotto di una matrice per cui si può moltiplicare, cioè: (A+B)C = AC+BC e A(C+D) = AC+AD
- V λ € K si ha che: λ (AC) = (λA) C = A( λC)
- Siano A € Mmn(K) e B € Mnp(K) e C € Mpq(K), allora si può fare: (AB)C = A (BC)
- L’ordine con cui vengono moltiplicate le matrici va rispettato perché il prodotto non è commutativo.
- Si ha anche che (AB)T = BTAT
Osservazioni
- Siano A = (aij) € Mmp(K) e B = (bij) € Mpn(K), quest’ultima matrice, però, la consideriamo come insieme delle sue righe, ogni singola riga è Rp = (bp1 bp2 … bpn). Allora la i-esima riga di AB è: ai1 R1 + ai2 R2 + … + aip Rp.
- Siano A = (C1 C2 &h