Moltiplicazione di matrici

Il disegno mostra il caso in cui A è 4 × 2 e B è 2 × 3, e si voglia calcolare l'elemento (C)₁₂ = (AB)₁₂ della matrice prodotto C = A x B, di dimensioni 4 x 3:
$(C{)}_{12}=\sum _{r=1}^2{a}_{1r}{b}_{r2}=a_{11}{b}_{12}+a_{12}{b}_{22}$

In matematica, e più precisamente in algebra lineare, la moltiplicazione di matrici è il prodotto righe per colonne tra due matrici, possibile sotto certe condizioni, che dà luogo ad un'altra matrice. Se una matrice rappresenta una applicazione lineare, il prodotto fra matrici è la traduzione della composizione di due applicazioni lineari. Quindi se due matrici 2 x 2 rappresentano ad esempio due rotazioni nel piano di angoli α e β, il loro prodotto è definito in modo tale da rappresentare una rotazione di angolo α + β.

Definizione

Sia $K$ un anello. Siano date una matrice $A$ di dimensione $m\times n$ ed una seconda matrice $B$ di dimensioni $n\times p$ a valori in $K$. Siano $a_{ij}$ gli elementi di $A$ e $b_{ij}$ gli elementi di $B$. Si definisce il prodotto matriciale di $A$ per $B$ la matrice $C=AB$ a valori in $K$ e di dimensioni $m\times p$ i cui elementi $c_{ij}$ sono dati da:

$c_{ij}=\sum _{r=1}^n{a}_{ir}{b}_{rj}=a_{i1}{b}_{1j}+a_{i2}{b}_{2j}+\cdots +a_{in}{b}_{nj}$

per ogni valore di riga $i$ e di colonna $j.$

Due matrici possono essere moltiplicate fra loro solo se il numero di colonne della prima è uguale al numero di righe della seconda, e il prodotto tra due matrici non è commutativo.

Una matrice può essere moltiplicata con se stessa solo se è quadrata. In questo caso, il prodotto $A\cdot A$ si denota con $A^2$. Più in generale, la potenza $n$-esima di una matrice è:

$A^n=\{\begin{array}{cc}I,& \text{se }n=0,\\ \underset{n}{\underbrace{A\cdot A\cdots A}},& \text{se }n>0,\end{array}$

dove $n$ è un numero naturale e $I$ è la matrice identità. Tuttavia, per esponenti molto maggiori dell'ordine della matrice, è più semplice calcolare le potenze servendosi della teoria delle funzioni di matrice, che permette inoltre di generalizzare la definizione di potenza fino ad ammettere un esponente complesso.

Un'altra definizione informale della moltiplicazione matriciale, atta a permetterne una più rapida e immediata memorizzazione, è "moltiplicazione riga per colonna", infatti, per ottenere l'elemento della $i$-esima riga e $j$-esima colonna della matrice prodotto basta porre un indice sulla riga $i$ della prima matrice, l'altro sulla colonna $j$ della seconda e moltiplicare gli elementi indicati, quindi scorrere di un posto con le dita e moltiplicare, fino a raggiungere la fine della colonna e della riga, infine sommare i vari prodotti ottenuti.

Proprietà

• La moltiplicazione fra matrici è generalmente non commutativa, ovvero $AB$ e $BA$ sono due matrici diverse.
• La moltiplicazione fra matrici è distributiva rispetto alla somma. In altre parole:

$A(B+C)=AB+AC$

$(A+B)C=AC+BC$

• Per ogni scalare $k$ vale:

$k(AB)=(kA)B=A(kB)$

• La moltiplicazione fra matrici è associativa:

$A(BC)=(AB)C$

• Le matrici aventi valori in un anello (ad esempio, l'anello dei numeri interi, razionali, reali o complessi) con le operazioni di somma e prodotto formano un altro anello. Per quanto detto sopra, questo anello è generalmente non commutativo anche se quello di partenza lo è.
• L'elemento neutro per l'operazione di moltiplicazione fra matrici è la matrice identica $I$. In particolare, se $A$ è quadrata con lo stesso numero di righe di $I$:

$AI=IA=A$

• La matrice nulla 0 con $n$ righe annulla qualsiasi altra matrice. In particolare, se $A$ è quadrata con $n$ righe, si ha:

$0A=A0=0$

• Una matrice quadrata $A$ è invertibile se esiste un'altra matrice $B$ tale che $AB=BA=I$ dove $I$ è la matrice identità con lo stesso numero di righe di $A$. Una matrice è invertibile se e solo se il determinante è non nullo. Molte matrici non sono invertibili; in altre parole, anche se l'insieme dei valori di partenza è un campo, le matrici non formano un campo. Ad esempio la matrice seguente non è invertibile:

$\left[\begin{array}{cc}0& 0\\ 1& 1\end{array}\right]$

• Detta $A^{\mathrm{T}}$ la trasposta di $A$, si ha che $(AB{)}^{\mathrm{T}}=B^{\mathrm{T}}{A}^{\mathrm{T}}$. Infatti:

${[(AB{)}^{\mathrm{T}}]}_{ij}={\left(AB\right)}_{ji}=\sum _k{\left(A\right)}_{jk}{\left(B\right)}_{ki}=\sum _k{\left(A^{\mathrm{T}}\right)}_{kj}{\left(B^{\mathrm{T}}\right)}_{ik}=\sum _k{\left(B^{\mathrm{T}}\right)}_{ik}{\left(A^{\mathrm{T}}\right)}_{kj}=$

$={\left[\left(B^{\mathrm{T}}\right)\left(A^{\mathrm{T}}\right)\right]}_{ij}$

• Detta $A^{-1}$ la inversa di $A$, si ha $(AB{)}^{-1}=B^{-1}{A}^{-1}$.
• Detta $A^{*}$ la complessa coniugata di $A$, si ha $(AB{)}^{*}=A^{*}{B}^{*}$. Infatti:

${[(AB{)}^{\star }]}_{ij}={\left[\sum _k{\left(A\right)}_{ik}{\left(B\right)}_{kj}\right]}^{\star }=\sum _k{\left(A\right)}_{ik}^{\star }{\left(B\right)}_{kj}^{\star }=\sum _k{\left(A^{\star }\right)}_{ik}{\left(B^{\star }\right)}_{kj}={\left(A^{\star }{B}^{\star }\right)}_{ij}$

• Detta $A^{†}$ la trasposta complessa coniugata di $A$, si ha $(AB{)}^{†}=B^{†}{A}^{†}$. Infatti:

${[(AB{)}^{†}]}_{ij}={\left[{\left(AB\right)}^{\star }\right]}_{ji}=\sum _k{\left(A^{\star }\right)}_{jk}{\left(B^{\star }\right)}_{ki}=\sum _k{\left(A^{†}\right)}_{kj}{\left(B^{†}\right)}_{ik}=\sum _k{\left(B^{†}\right)}_{ik}{\left(A^{†}\right)}_{kj}={\left(B^{†}{A}^{†}\right)}_{ij}$

• La traccia del prodotto $AB$ è indipendente dall'ordine con cui $A$ e $B$ vengono moltiplicate:

$\mathrm{t}\mathrm{r}(AB)=\mathrm{t}\mathrm{r}(BA)$

Infatti:

$\mathrm{tr}(AB)=\sum _{i=1}^m{\left(AB\right)}_{ii}=\sum _{i=1}^m\sum _{j=1}^n{A}_{ij}{B}_{ji}=\sum _{j=1}^n\sum _{i=1}^m{B}_{ji}{A}_{ij}=\sum _{j=1}^n{\left(BA\right)}_{jj}=\mathrm{tr}(BA)$

Prodotto di una matrice per un vettore

Una matrice con una sola riga, cioè di dimensione $1\times n$, è un vettore riga. Analogamente, una matrice con una sola colonna, cioè di dimensione $m\times 1$ è un vettore colonna. Nell'operazione di moltiplicazione questi due oggetti si comportano in modo differente.

Siano $A=(a_{ij})$ una matrice $m\times n$ e $𝐯=\{v_i\}$ un vettore colonna $n\times 1$. Il prodotto di $A$ per il vettore $𝐯$ è il prodotto di matrici:

$A𝐯=𝐜$

Le componenti di $𝐜$ sono:

$c_i=\sum _{j=1}^n{a}_{ij}{v}_j=a_{i1}{v}_1+a_{i2}{v}_2+\cdots +a_{in}{v}_n$

Algoritmo

Un algoritmo per la moltiplicazione matrice per vettore è:

    /* Moltiplicazione matrice × vettore
       RM = numero di righe della matrice
       CM = numero di colonne della matrice (uguale al numero di righe del vettore)
       M = matrice [RM] × [CM]
       VI = vettore iniziale [CM]
       VR = vettore risultante [CM]
       il vettore risultato sarà VR [RM] con stesso numero di righe della matrice. */
    for (int i = 0; i < RM; i++) {      // scandisco le righe con l'indice i
        VR[i] = 0;                      // inizializzo la coordinata i-esima del vettore a zero
        for (int j = 0; j < CM; j++) {  // e le colonne con j 
            VR[i] = VR[i] + M[i][j] * VI[j];
        }
    }

Questo prodotto è ampiamente usato in algebra lineare perché descrive una applicazione lineare. Ad esempio, il prodotto:

$\left[\begin{array}{cc}\cos \alpha & -\sin \alpha \\ \sin \alpha & \cos \alpha \end{array}\right]\left[\begin{array}{c}x\\ y\end{array}\right]=\left[\begin{array}{c}x\cos \alpha -y\sin \alpha \\ x\sin \alpha +y\cos \alpha \end{array}\right]$

rappresenta una rotazione di angolo $\alpha$ nel piano cartesiano.

In alcuni casi può essere utile effettuare il prodotto $vettoreriga\times Matrice$: il risultato è un altro vettore riga.

Prodotto di una matrice per uno scalare

La moltiplicazione di una matrice $A=(a_{i,j})$ per uno scalare $r$, cioè un elemento dell'anello cui appartengono gli $a_{i,j}$, è ottenuta moltiplicando ogni elemento di $A$ per lo scalare:

$rA=(r{a}_{i,j})$

Se l'anello di partenza non è commutativo, questa viene indicata come moltiplicazione sinistra, e può differire dalla moltiplicazione destra:

$Ar=(a_{i,j}r)$

Proprietà

• Se l'anello di partenza è commutativo (ad esempio se è l'anello dei numeri interi, razionali, reali o complessi) le moltiplicazioni sinistra e destra sono equivalenti e si parla solo di moltiplicazione di una matrice con uno scalare.
• Se l'anello di partenza è un campo, ad esempio quello dei numeri razionali, reali o complessi, lo spazio delle matrici $m\times n$ con le operazioni di somma e prodotto per scalare formano uno spazio vettoriale.
• Se l'anello di partenza è un anello commutativo, lo spazio delle matrici $m\times n$ con le operazioni di somma e di prodotto per scalare forma un modulo.

Se l'anello di partenza non è commutativo, ad esempio se è l'anello dei quaternioni, le due moltiplicazioni non sono equivalenti. Ad esempio:

$i\left[\begin{array}{cc}i& 0\\ 0& j\end{array}\right]=\left[\begin{array}{cc}-1& 0\\ 0& k\end{array}\right]\ne \left[\begin{array}{cc}-1& 0\\ 0& -k\end{array}\right]=\left[\begin{array}{cc}i& 0\\ 0& j\end{array}\right]i$

Costruzioni alternative

Sono stati definiti nel tempo altri tipi di prodotto tra matrici, meno fortunati in quanto a utilizzo dell'usuale prodotto righe per colonne. In particolare si può nominare il prodotto di Hadamard o prodotto puntuale, in cui il prodotto di $A=(a_{ij}{)}_{ij}$ e $B=(b_{ij}{)}_{ij}$ è dato da $AB=(a_{ij}\cdot b_{ij}{)}_{ij}$. Ad esempio:

$\left[\begin{array}{cc}1& 2\\ 3& -1\end{array}\right]\left[\begin{array}{cc}-3& 0\\ 1& 4\end{array}\right]=\left[\begin{array}{cc}-3& 0\\ 3& -4\end{array}\right]$

Un'altra costruzione è data dal prodotto di Kronecker, che trova applicazioni nel calcolo tensoriale, dato da:

$AB=\left[\begin{array}{cccc}{a}_{11}B& a_{12}B& \cdots & a_{1n}B\\ ⋮& ⋮& \ddots & ⋮\\ a_{m1}B& a_{m2}B& \cdots & a_{mn}B\end{array}\right]$

espressa sotto forma di matrice a blocchi, in cui ogni blocco $ij$-esimo è dato dalla matrice $B$ moltiplicata per lo scalare $a_{ij}$.

Esempi

• Una matrice $2\times 3$ moltiplicata per $3\times 3$ dà una matrice $2\times 3$:

$\left[\begin{array}{ccc}1& 1& 2\\ 0& 1& -3\end{array}\right]\times \left[\begin{array}{ccc}1& 1& 1\\ 2& 5& 1\\ 0& -2& 1\end{array}\right]=$

1ª riga della matrice risultato:

$C_{11}=[(1\times 1)+(1\times 2)+(2\times 0)]=3$

$C_{12}=[(1\times 1)+(1\times 5)+(2\times -2)]=2$

$C_{13}=[(1\times 1)+(1\times 1)+(2\times 1)]=4$

2ª riga della matrice risultato:

$C_{21}=[(0\times 1)+(1\times 2)+(-3\times 0)]=2$

$C_{22}=[(0\times 1)+(1\times 5)+(-3\times -2)]=11$

$C_{23}=[(0\times 1)+(1\times 1)+(-3\times 1)]=-2$

Risultato (matrice $2\times 3$):

$\left[\begin{array}{ccc}{C}_{11}& C_{12}& C_{13}\\ C_{21}& C_{22}& C_{23}\end{array}\right]=\left[\begin{array}{ccc}3& 2& 4\\ 2& 11& -2\end{array}\right]$

• Si consideri il prodotto:

$\left[\begin{array}{ccc}1& 2& 0\\ 3& -1& 4\end{array}\right]\times \left[\begin{array}{c}1\\ 0\\ -1\end{array}\right]=\left[\begin{array}{c}(1\times 1+2\times 0+0\times -1)\\ (3\times 1+-1\times 0+4\times -1)\end{array}\right]=\left[\begin{array}{c}1\\ -1\end{array}\right]$

Il risultato di questa operazione è un altro vettore colonna, di tipo $m\times 1$.

Bibliografia

• Template:Cita libro
• Template:Cita libro
• Template:Cita libro ISBN 88-386-6289-4.
• Template:Cita libro ISBN 88-339-5447-1.
• Template:En Henry Cohn, Robert Kleinberg, Balazs Szegedy, and Chris Umans. Group-theoretic Algorithms for Matrix Multiplication. Template:Arxiv. Proceedings of the 46th Annual Symposium on Foundations of Computer Science, 23–25 October 2005, Pittsburgh, PA, IEEE Computer Society, pp. 379–388.
• Template:En Henry Cohn, Chris Umans. A Group-theoretic Approach to Fast Matrix Multiplication. Template:Arxiv. Proceedings of the 44th Annual IEEE Symposium on Foundations of Computer Science, 11–14 October 2003, Cambridge, MA, IEEE Computer Society, pp. 438–449.
• Template:En Coppersmith, D., Winograd S., Matrix multiplication via arithmetic progressions, J. Symbolic Comput. 9, p. 251-280, 1990.
• Template:En Robinson, Sara, Toward an Optimal Algorithm for Matrix Multiplication, SIAM News 38(9), November 2005. PDF
• Template:En Strassen, Volker, Gaussian Elimination is not Optimal, Numer. Math. 13, p. 354-356, 1969.
• Template:En Vassilevska Williams, Virginia, Multiplying matrices faster than Coppersmith-Winograd, Manuscript, May 2012. PDF Template:Webarchive

Voci correlate

• Anello commutativo
• Composizione di funzioni
• Funzione di matrice
• Matrice
• Somma fra matrici
• Teorema di Binet

Collegamenti esterni

• Template:Collegamenti esterni
• Template:Cita web