Moltiplicazione di matrici

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 ${\displaystyle K}$ un anello. Siano date una matrice ${\displaystyle A}$ di dimensione ${\displaystyle m\times n}$ ed una seconda matrice ${\displaystyle B}$ di dimensioni ${\displaystyle n\times p}$ a valori in ${\displaystyle K}$. Siano ${\displaystyle a_{ij}}$ gli elementi di ${\displaystyle A}$ e ${\displaystyle b_{ij}}$ gli elementi di ${\displaystyle B}$. Si definisce il prodotto matriciale di ${\displaystyle A}$ per ${\displaystyle B}$ la matrice ${\displaystyle C=AB}$ a valori in ${\displaystyle K}$ e di dimensioni ${\displaystyle m\times p}$ i cui elementi ${\displaystyle c_{ij}}$ sono dati da:^[1]

${\displaystyle 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 ${\displaystyle i}$ e di colonna ${\displaystyle 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.^[2]

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

${\displaystyle A^{n}={\begin{cases}I,&{\text{se }}n=0,\\\underbrace {A\cdot A\cdots A} _{n},&{\text{se }}n>0,\end{cases}}}$

dove ${\displaystyle n}$ è un numero naturale e ${\displaystyle 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 ${\displaystyle i}$-esima riga e ${\displaystyle j}$-esima colonna della matrice prodotto basta porre un indice sulla riga ${\displaystyle i}$ della prima matrice, l'altro sulla colonna ${\displaystyle 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 ${\displaystyle AB}$ e ${\displaystyle BA}$ sono due matrici diverse.
• La moltiplicazione fra matrici è distributiva rispetto alla somma. In altre parole:

${\displaystyle A(B+C)=AB+AC}$

${\displaystyle (A+B)C=AC+BC}$

• Per ogni scalare ${\displaystyle k}$ vale:

${\displaystyle k(AB)=(kA)B=A(kB)}$

• La moltiplicazione fra matrici è associativa:

${\displaystyle 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 ${\displaystyle I}$. In particolare, se ${\displaystyle A}$ è quadrata con lo stesso numero di righe di ${\displaystyle I}$:

${\displaystyle AI=IA=A}$

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

${\displaystyle 0A=A0=0}$

• Una matrice quadrata ${\displaystyle A}$ è invertibile se esiste un'altra matrice ${\displaystyle B}$ tale che ${\displaystyle AB=BA=I}$ dove ${\displaystyle I}$ è la matrice identità con lo stesso numero di righe di ${\displaystyle 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:

${\displaystyle {\begin{bmatrix}0&0\\1&1\\\end{bmatrix}}}$

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

${\displaystyle \left[({AB})^{\mathrm {T} }\right]_{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}=}$

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

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

${\displaystyle \left[({AB})^{\star }\right]_{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 ${\displaystyle A^{\dagger }}$ la trasposta complessa coniugata di ${\displaystyle A}$, si ha ${\displaystyle (AB)^{\dagger }=B^{\dagger }A^{\dagger }}$. Infatti:

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

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

${\displaystyle \mathrm {tr} ({AB})=\mathrm {tr} ({BA})}$

Infatti:

${\displaystyle \operatorname {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}=\operatorname {tr} (BA)}$

Prodotto di una matrice per un vettore

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

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

${\displaystyle A\mathbf {v} =\mathbf {c} }$

Le componenti di ${\displaystyle \mathbf {c} }$ sono:

${\displaystyle 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:

${\displaystyle {\begin{bmatrix}\cos \alpha &-\sin \alpha \\\sin \alpha &\cos \alpha \\\end{bmatrix}}{\begin{bmatrix}x\\y\\\end{bmatrix}}={\begin{bmatrix}x\cos \alpha -y\sin \alpha \\x\sin \alpha +y\cos \alpha \\\end{bmatrix}}}$

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

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

Prodotto di una matrice per uno scalare

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

${\displaystyle rA=(ra_{i,j})\ }$

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

${\displaystyle 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 ${\displaystyle 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 ${\displaystyle 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:

${\displaystyle i{\begin{bmatrix}i&0\\0&j\\\end{bmatrix}}={\begin{bmatrix}-1&0\\0&k\\\end{bmatrix}}\neq {\begin{bmatrix}-1&0\\0&-k\\\end{bmatrix}}={\begin{bmatrix}i&0\\0&j\\\end{bmatrix}}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 ${\displaystyle A=(a_{ij})_{ij}}$ e ${\displaystyle B=(b_{ij})_{ij}}$ è dato da ${\displaystyle AB=(a_{ij}\cdot b_{ij})_{ij}}$. Ad esempio:

${\displaystyle {\begin{bmatrix}1&2\\3&-1\end{bmatrix}}{\begin{bmatrix}-3&0\\1&4\end{bmatrix}}={\begin{bmatrix}-3&0\\3&-4\end{bmatrix}}}$

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

${\displaystyle AB={\begin{bmatrix}a_{11}B&a_{12}B&\cdots &a_{1n}B\\\vdots &\vdots &\ddots &\vdots \\a_{m1}B&a_{m2}B&\cdots &a_{mn}B\end{bmatrix}}}$

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

Esempi

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

${\displaystyle {\begin{bmatrix}1&1&2\\0&1&-3\\\end{bmatrix}}\times {\begin{bmatrix}1&1&1\\2&5&1\\0&-2&1\end{bmatrix}}=}$

1ª riga della matrice risultato:

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

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

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

2ª riga della matrice risultato:

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

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

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

Risultato (matrice ${\displaystyle 2\times 3}$):

${\displaystyle {\begin{bmatrix}C_{11}&C_{12}&C_{13}\\C_{21}&C_{22}&C_{23}\\\end{bmatrix}}={\begin{bmatrix}3&2&4\\2&11&-2\\\end{bmatrix}}}$

• Si consideri il prodotto:

${\displaystyle {\begin{bmatrix}1&2&0\\3&-1&4\\\end{bmatrix}}\times {\begin{bmatrix}1\\0\\-1\\\end{bmatrix}}={\begin{bmatrix}(1\times 1+2\times 0+0\times -1)\\(3\times 1+-1\times 0+4\times -1)\\\end{bmatrix}}={\begin{bmatrix}1\\-1\\\end{bmatrix}}}$

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

Note

Bibliografia

• Serge Lang, Algebra lineare, Torino, Bollati Boringhieri, 1992, ISBN 88-339-5035-2.
• (EN) Kenneth Hoffman e Ray Kunze, Linear Algebra, 2ª ed., Englewood Cliffs, New Jersey, Prentice - Hall, inc., 1971, ISBN 0-13-536821-9.
• Marco Abate, Chiara de Fabritiis, Geometria analitica con elementi di algebra lineare, Milano, McGraw-Hill, 2006. ISBN 88-386-6289-4.
• Edoardo Sernesi, Geometria 1, 2ª ed., Torino, Bollati Boringhieri, 1989. ISBN 88-339-5447-1.
• (EN) Henry Cohn, Robert Kleinberg, Balazs Szegedy, and Chris Umans. Group-theoretic Algorithms for Matrix Multiplication. arΧiv:math.GR/0511460. Proceedings of the 46th Annual Symposium on Foundations of Computer Science, 23–25 October 2005, Pittsburgh, PA, IEEE Computer Society, pp. 379–388.
• (EN) Henry Cohn, Chris Umans. A Group-theoretic Approach to Fast Matrix Multiplication. arΧiv:math.GR/0307321. Proceedings of the 44th Annual IEEE Symposium on Foundations of Computer Science, 11–14 October 2003, Cambridge, MA, IEEE Computer Society, pp. 438–449.
• (EN) Coppersmith, D., Winograd S., Matrix multiplication via arithmetic progressions, J. Symbolic Comput. 9, p. 251-280, 1990.
• (EN) Robinson, Sara, Toward an Optimal Algorithm for Matrix Multiplication, SIAM News 38(9), November 2005. PDF
• (EN) Strassen, Volker, Gaussian Elimination is not Optimal, Numer. Math. 13, p. 354-356, 1969.
• (EN) Vassilevska Williams, Virginia, Multiplying matrices faster than Coppersmith-Winograd, Manuscript, May 2012. PDF Archiviato l'8 ottobre 2013 in Internet Archive.

Voci correlate

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

Altri progetti

• Wikimedia Commons contiene immagini o altri file sulla moltiplicazione di matrici

Collegamenti esterni

• Programma parallelo per moltiplicare le matrici in MPI, su parallelknoppix.info. URL consultato il 14 dicembre 2022 (archiviato dall'url originale il 12 aprile 2013).
• (EN) Eric W. Weisstein, Moltiplicazione di matrici, su MathWorld, Wolfram Research.

Portale Matematica: accedi alle voci di Wikipedia che trattano di matematica