Vizualizare sursă pentru „Înmulțirea matricilor”

[[Fișier:Matrix multiplication qtl1.svg|thumb|La înmulțirea matricilor, numărul de coloane din prima matrice trebuie să fie egal cu numărul de linii din a doua matrice. Matricea rezultată are numărul de linii ale primei matrice și numărul de coloane ale celei de-a doua matrice.]]
În [[matematică]], în special în [[algebră liniară |algebra liniară]], '''înmulțirea matricilor''' sau '''înmulțirea matricială'''<ref name=AI><!-- Conf.dr. -->Anca Ignat, [https://profs.info.uaic.ro/~ancai/CN/curs/curs_CN_2022.pdf ''Calcul numeric''] {{Webarchive|url=https://web.archive.org/web/20230613045258/https://profs.info.uaic.ro/~ancai/CN/curs/curs_CN_2022.pdf |date=2023-06-13 }} (curs 2, 2022, p. 2), [[Universitatea „Alexandru Ioan Cuza” din Iași]], accesat 2023-06-13</ref> este o [[operație binară]] care produce o [[matrice]] din două matrici. La înmulțirea matricială, numărul de coloane din prima matrice trebuie să fie egal cu numărul de linii din a doua matrice. Matricea rezultată, cunoscută sub numele de '''produs matricial''', are numărul de linii ale primei matrice și numărul de coloane ale celei de-a doua matrice. Produsul matricilor {{math|'''A'''}} și {{math|'''B'''}} este notat {{math|'''AB'''}}.<ref name=DN>{{en icon}} {{Cite web|last=Nykamp|first=Duane|title=Multiplying matrices and vectors|url=https://mathinsight.org/matrix_vector_multiplication|access-date=September 6, 2020|website=Math Insight}}</ref>

Înmulțirea matricială a fost descrisă pentru prima dată de [[matematician]]ul [[francezi |francez]] [[Jacques Philippe Marie Binet]] în 1812,<ref>{{en icon}} {{MacTutor|id=Binet|title=Jacques Philippe Marie Binet}}</ref> pentru a reprezenta {{ill-wd| Q244761||compunerea funcțiilor}} de [[transformare liniară |transformări liniare]] care sunt reprezentate cu ajutorul matricilor. Astfel, înmulțirea matricială a devenit un instrument de bază al algebrei liniare și, ca atare, are numeroase aplicații în multe domenii ale matematicii, precum și în [[matematică aplicată |matematica aplicată]], [[statistică]], [[fizică]], [[economie]] și [[inginerie]].<ref name="Physics 1991">{{en icon}} {{cite book|title=Encyclopaedia of Physics|edition=2nd|first1=Rita G. |last1=Lerner | first2= G. L. |last2 = Trigg |publisher=VHC publishers|date=1991|isbn=978-3-527-26954-9}}</ref><ref>{{en icon}} {{cite book|title=McGraw Hill Encyclopaedia of Physics|edition=2nd|first=C. B.|last=Parker|date=1994|isbn=978-0-07-051400-3|url-access=registration|url=https://archive.org/details/mcgrawhillencycl1993park}}</ref> Calculul produselor matriciale este o operație centrală în toate aplicațiile algebrei liniare.

== Notație ==
În acest articol convențiile privind notațiile sunt următoarele: matricile sunt reprezentate prin litere majuscule grase (bold), de exemplu {{math|'''A'''}}; [[vector euclidian |vectorii]] cu litere minuscule grase, de exemplu {{math|'''a'''}}; iar elementele vectorilor și matricilor cu italice, de exemplu {{math|''A''}} sau {{math|''a''}}. [[Notație indexată |Notația indexată]] este adesea cel mai clar mod de a exprima definițiile și este notația standard folosit în literatura de specialitate. Linia din matrice este indicată de indicele {{mvar|i}}, iar coloana de indicele {{mvar|j}}. Elementul matricei {{math|'''A'''}} este indicat prin {{math|('''A''')<sub>''ij''</sub>}}, {{math|''A''<sub>''ij''</sub>}} sau {{math|''a''<sub>''ij''</sub>}}. Prin excepție, un singur indice, de exemplu {{math|'''A'''<sub>1</sub>, '''A'''<sub>2</sub>}}, se folosește pentru a indica o matrice (nu un element al unei matrice) dintr-o colecție de matrici.

== Definiție ==
Dacă {{math|'''A'''}} este o matrice {{math|''m'' × ''n''}} iar {{math|'''B'''}} este o matrice {{math|''n'' × ''p''}},
:<math>\mathbf{A}=\begin{pmatrix}
 a_{11} & a_{12} & \cdots & a_{1n} \\
 a_{21} & a_{22} & \cdots & a_{2n} \\
\vdots & \vdots & \ddots & \vdots \\
 a_{m1} & a_{m2} & \cdots & a_{mn} \\
\end{pmatrix},\quad\mathbf{B}=\begin{pmatrix}
 b_{11} & b_{12} & \cdots & b_{1p} \\
 b_{21} & b_{22} & \cdots & b_{2p} \\
\vdots & \vdots & \ddots & \vdots \\
 b_{n1} & b_{n2} & \cdots & b_{np} \\
\end{pmatrix}</math>

''Produsul matricial'' {{math|1='''C''' = '''AB'''}} (notat fără [[interpunct |punct]] sau [[semnul înmulțirii]]) este definit ca fiind matricea {{math|''m'' × ''p''}} <ref>{{en icon}} {{cite book|title=Linear Algebra|edition=4th| first1 = S. | last1 = Lipschutz | first2 = M. | last2 = Lipson|series=Schaum's Outlines|publisher=McGraw Hill (USA)|date=2009|pages=30–31|isbn=978-0-07-154352-1}}</ref><ref>{{en icon}} {{cite book|title=Mathematical methods for physics and engineering|url=https://archive.org/details/mathematicalmeth00rile|url-access=registration| first1 = K. F. | last1 = Riley | first2 = M. P. | last2 = Hobson | first3 = S. J. | last3 = Bence|publisher=Cambridge University Press|date=2010|isbn=978-0-521-86153-3}}</ref><ref>{{en icon}} {{cite book|title=Calculus, A Complete Course|edition=3rd| first = R. A. | last = Adams|publisher=Addison Wesley|date=1995|page=627|isbn=0-201-82823-5}}</ref><ref>{{en icon}} {{cite book|title=Matrix Analysis| last = Horn | first = Johnson|edition=2nd|publisher=Cambridge University Press|date=2013|page=6|isbn=978-0-521-54823-6}}</ref>
:<math>\mathbf{C}=\begin{pmatrix}
 c_{11} & c_{12} & \cdots & c_{1p} \\
 c_{21} & c_{22} & \cdots & c_{2p} \\
\vdots & \vdots & \ddots & \vdots \\
 c_{m1} & c_{m2} & \cdots & c_{mp} \\
\end{pmatrix}</math>
asfel încât
:<math> c_{ij} = a_{i1}b_{1j} + a_{i2}b_{2j} +\cdots + a_{in}b_{nj}= \sum_{k=1}^n a_{ik}b_{kj}, </math>
pentru {{math|1=''i'' = 1, ..., ''m''}} și {{math|1=''j'' = 1, ..., ''p''}}.

Adică elementul {{mvar|c<sub>ij</sub>}} al produsului se obține prin înmulțirea termen cu termen a elementelor din linia {{mvar|i}} a matricei {{math|'''A'''}} cu cele din coloana {{mvar|j}} a matricei {{math|'''B'''}} și sumând aceste {{mvar|n}} produse. Altfel spus, {{mvar|c<sub>ij</sub>}} este [[produs scalar |produsul scalar]] al liniei {{mvar|i}} a matricei {{math|'''A'''}} cu coloana {{mvar|j}} a matricei {{math|'''B'''}}.

Prin urmare {{math|'''AB'''}} poate fi scris și ca
:<math>\mathbf{C}=\begin{pmatrix}
 a_{11}b_{11} +\cdots + a_{1n}b_{n1} & a_{11}b_{12} +\cdots + a_{1n}b_{n2} & \cdots & a_{11}b_{1p} +\cdots + a_{1n}b_{np} \\
 a_{21}b_{11} +\cdots + a_{2n}b_{n1} & a_{21}b_{12} +\cdots + a_{2n}b_{n2} & \cdots & a_{21}b_{1p} +\cdots + a_{2n}b_{np} \\
\vdots & \vdots & \ddots & \vdots \\
 a_{m1}b_{11} +\cdots + a_{mn}b_{n1} & a_{m1}b_{12} +\cdots + a_{mn}b_{n2} & \cdots & a_{m1}b_{1p} +\cdots + a_{mn}b_{np} \\
\end{pmatrix} 
</math>

Astfel produsul {{math|'''AB'''}} este definit [[dacă și numai dacă]] numărul de coloane din {{math|'''A'''}} este egal cu numărul de linii din {{math| '''B'''}},<ref name=DN/> în acest caz {{math|''n''}}.

În majoritatea cazurilor elementele matricilor sunt numere, dar pot fi orice fel de [[obiect matematic |obiecte matematice]] pentru care sunt definite o adunare și o înmulțire, care sunt [[asociativitate |asociative]], adunarea să fie [[comutativitate |comutativă]], iar înmulțirea să fie [[distributivitate |distributivă]] în raport cu adunarea. În particular, elementele pot fi ele înseși matrici.

=== Ilustrare ===
[[Fișier:Matrix multiplication diagram 2.svg|thumb|Prezentare vizuală a calculului termenilor produsului matricial]]
Figura din dreapta ilustrează schematic produsul a două matrice {{math|'''A'''}} și {{math|'''B'''}}, arătând cum fiecare intersecție din matricea produsului corespunde cu o linie din {{math|'''A'''}} și o coloană din {{math|'''B'''}}.
:<math>
\overset{4\times 2 \text{ matrix}}{\begin{bmatrix}
a_{11} & a_{12} \\
\cdot & \cdot \\
a_{31} & a_{32} \\
\cdot & \cdot \\
\end{bmatrix}}
\overset{2\times 3\text{ matrix}}{\begin{bmatrix}
\cdot & b_{12} & b_{13} \\
\cdot & b_{22} & b_{23} \\
\end{bmatrix}}

= \overset{4\times 3\text{ matrix}}{\begin{bmatrix}
\cdot & c_{12} & \cdot \\
\cdot & \cdot & \cdot \\
\cdot & \cdot & c_{33} \\
\cdot & \cdot & \cdot \\
\end{bmatrix}}
</math>

Valorile de la intersecții, marcate cu cercuri în figura din dreapta, sunt:
:<math>\begin{align}
c_{12} & = a_{11} b_{12} + a_{12} b_{22} \\
c_{33} & = a_{31} b_{13} + a_{32} b_{23}
\end{align}</math>

== Principalele aplicații ==
Istoric, înmulțirea matricilor a fost introdusă pentru a facilita și clarifica calculele din [[algebră liniară |algebra liniară]]. Această relație strânsă între înmulțirea matricilor și algebra liniară rămâne fundamentală în toate ramurile matematicii, precum și în [[fizică]], [[chimie]], [[inginerie]] și [[informatică]].

=== Aplicații liniare ===
Dacă un [[spațiu vectorial]] are o [[bază (algebră liniară) |bază]] finită, orice vector este reprezentat în mod unic printr-un [[șir (matematică) |șir]] finit de scalari, numit [[componentele unui vector |componentele vectorului]], ale cărui elemente sunt [[sistem de coordonate |coordonatele]] vectorului în acea bază. Aceste componente formează un alt spațiu vectorial, care este [[izomorfism |izomorf]] cu spațiul vectorial inițial. Componentele vectorului sunt plasate de obicei într-o matrice coloană (numită și ''[[vector linie și vector coloană |vector coloană]]''), care este o matrice cu o singură coloană. Deci, un vector coloană reprezintă atât componentele vectorului, cât și un vector din spațiului vectorial inițial.

O [[aplicație liniară]] {{mvar|A}} dintr-un spațiu vectorial de dimensiune {{mvar|n}} într-un spațiu vectorial de dimensiune {{mvar|m}} aplică un vector coloană
:<math>\mathbf x=\begin{pmatrix}x_1 \\ x_2 \\ \vdots \\ x_n\end{pmatrix}</math>
pe vectorul coloană
:<math>\mathbf y= A(\mathbf x)= \begin{pmatrix}a_{11}x_1+\cdots + a_{1n}x_n\\ a_{21}x_1+\cdots + a_{2n}x_n \\ \vdots \\ a_{m1}x_1+\cdots + a_{mn}x_n\end{pmatrix}.</math>

Aplicația liniară {{mvar|A}} este astfel definită de matricea 
:<math>\mathbf{A}=\begin{pmatrix}
 a_{11} & a_{12} & \cdots & a_{1n} \\
 a_{21} & a_{22} & \cdots & a_{2n} \\
\vdots & \vdots & \ddots & \vdots \\
 a_{m1} & a_{m2} & \cdots & a_{mn} \\
\end{pmatrix}, </math>
și aplică vectorul coloană <math>\mathbf x</math> pe produsul matricial
:<math>\mathbf y = \mathbf {Ax}.</math>

Dacă {{mvar|B}} este o altă aplicație liniară din spațiul vectorial precedent de dimensiunea {{mvar|m}}, într-un spațiu vectorial de dimensiunea {{mvar|p}}, este reprezentată printr-o matrice <math>\mathbf B</math> ({{mvar|p × m}}). Un calcul simplu arată că matricea compunerii <math>B\circ A</math> este produsul matricial <math>\mathbf {BA}.</math> Formula generală <math>(B\circ A)(\mathbf x) = B(A(\mathbf x))</math> care definește funcția compusă este prezentată aici ca un caz particular de asociativitate al produsului matricial:
:<math>(\mathbf{BA})\mathbf x = \mathbf{B}(\mathbf {Ax}) = \mathbf{BAx}.</math>

==== Rotațiile geometrice ====
Folosind un sistem de [[coordonate carteziene]] într-un [[plan (geometrie) |plan]] euclidian, [[rotație (matematică) |rotația]] cu un [[unghi]] <math>\alpha</math> în jurul [[origine (matematică) |originii]] este o aplicație liniară. Mai exact,
:<math> \begin{bmatrix} x' \\ y' \end{bmatrix} =
 \begin{bmatrix} \cos \alpha & - \sin \alpha \\ \sin \alpha & \cos \alpha \end{bmatrix} \begin{bmatrix} x \\ y \end{bmatrix},</math>
unde punctul inițial <math>(x,y)</math> și imaginea sa <math>(x',y')</math> sunt scrise ca vectori coloană.

Compunerea rotației cu unghiul <math>\alpha</math> și cea cu unghiul <math>\beta</math> corespunde apoi cu produsul matricial
:<math>\begin{bmatrix} \cos \beta & - \sin \beta \\ \sin \beta & \cos \beta \end{bmatrix} 
 \begin{bmatrix} \cos \alpha & - \sin \alpha \\ \sin \alpha & \cos \alpha \end{bmatrix}
  = \begin{bmatrix} \cos \beta \cos \alpha - \sin \beta \sin \alpha & - \cos \beta \sin \alpha - \sin \beta \cos \alpha \\ \sin \beta \cos \alpha + \cos \beta \sin \alpha & - \sin \beta \sin \alpha + \cos \beta \cos \alpha \end{bmatrix}
 = \begin{bmatrix} \cos (\alpha+\beta) & - \sin(\alpha+\beta) \\ \sin(\alpha+\beta) & \cos(\alpha+\beta) \end{bmatrix},</math>
în care membrul drept al ultimei identități conține identitățile trigonometrice pentru [[sinus]]ul și [[cosinus]]ul sumei și diferenței unghiurilor, compunerea corespunde rotației cu unghiul <math>\alpha+\beta</math>, cum era de așteptat.

==== Alocarea resurselor în economie ====
[[Fișier:Mmult factory svg.svg|thumb|400px|Calculul elementului din stânga jos a <math>\mathbf{AB}</math> corespunde luării în considerare a tuturor căilor (evidențiate) de la materia primă <math>b_4</math> la produsul final <math>f_1</math> în [[graf]]ul fluxului de producție]]
De exemplu, o fabrică fictivă folosește 4 tipuri de [[materie primă |materii prime]], <math>b_1, b_2, b_3, b_4</math> pentru a produce 3 tipuri de [[bun intermediar |bunuri intermediare]] , <math>m_1, m_2, m_3</math>, care la rândul lor sunt folosite pentru a produce 3 tipuri de [[bun final |produse finite]], <math>f_1, f_2, f_3</math>. Matricile
:<math>\mathbf{A} = \begin{pmatrix} 1 & 0 & 1 \\   2 & 1 & 1 \\   0 & 1 & 1 \\   1 & 1 & 2 \\   \end{pmatrix} </math> &nbsp; și &nbsp; <math>\mathbf{B} = \begin{pmatrix} 1 & 2 & 1 \\   2 & 3 & 1 \\   4 & 2 & 2 \\   \end{pmatrix} </math>
furnizează cantitatea de materii prime necesare pentru o anumită cantitate de bunuri intermediare și, respectiv, cantitatea de bunuri intermediare necesară pentru o anumită cantitate de produse finite. De exemplu, pentru a produce o unitate de bun intermediar <math>m_1</math>, sunt necesare o unitate de materie primă <math>b_1</math>, două unități de <math>b_2</math>, nicio unitate de <math> >b_3</math> și o unitate de <math>b_4</math>, corespunzătoare primei coloane din <math>\mathbf{A}</math>.

Folosind înmulțirea matricilor se calculează
:<math>\mathbf{AB} = \begin{pmatrix} 5 & 4 & 3 \\   8 & 9 & 5 \\\   6 & 5 & 3 \\   11 & 9 & 6 \\   \end{pmatrix} ;</math>
această matrice indică direct cantitățile de materii prime necesare pentru cantități date de produse finite. De exemplu, intrarea din stânga jos a <math>\mathbf{AB}</math> este calculată ca fiind <math>1 \cdot 1 + 1 \cdot 2 + 2 \cdot 4 = 11</math>, reflectând faptul că pentru a produce o unitate de <math>f_1</math> sunt necesare <math>11</math> unități de <math>b_4</math>. Efectiv, este necesară o unitate de <math>b_4</math> pentru <math>m_1</math>, 2 pentru <math>m_2</math> și <math>4</math> pentru fiecare dintre cele două <math>m_3</math> unități care intră în unitatea <math>f_1</math>, vezi imaginea.

Pentru a produce de exemplu 100 de unități din produsul final <math>f_1</math>, 80 de unități de <math>f_2</math> și 60 de unități de <math>f_3</math>, cantitățile necesare de materii prime pot fi calculate prin
:<math>(\mathbf{AB}) \begin{pmatrix} 100 \\ 80 \\ 60 \\ \end{pmatrix} = \begin{pmatrix} 1000 \\ 1820 \\ 1180 \\ 2180 \end{pmatrix} ,</math>
adică <math>1000</math> unități de <math>b_1</math>, <math>1820</math> unități de <math>b_2</math>, <math>1180</math> unități de <math>b_3</math> și <math>2180</math> unități de <math>b_4</math>.
Similar, matricea produsului <math>\mathbf{AB}</math> poate fi utilizată pentru a calcula cantitățile necesare de materii prime pentru alte date privind cantitățile de produse finale.<ref>{{de icon}} {{cite book | isbn=3-446-18668-9 | author=Peter Stingl | title=Mathematik für Fachhochschulen – Technik und Informatik | location=[[München]] | publisher=Carl Hanser Verlag | edition=5th | year=1996 }} Here: Exm.&nbsp;5.4.10, pp.&nbsp;205–206</ref>

=== Sisteme de ecuații liniare ===
Forma generală a unui [[sistem de ecuații liniare]] este:
:<math>\begin{matrix}a_{11}x_1+\cdots + a_{1n}x_n=b_1
\\ a_{21}x_1+\cdots + a_{2n}x_n =b_2
\\ \vdots 
\\ a_{m1}x_1+\cdots + a_{mn}x_n =b_m\end{matrix}.</math>

Folosind aceeași notație ca mai sus, un astfel de sistem este echivalent cu o singură [[ecuație]] matricială:
:<math>\mathbf{Ax}=\mathbf b.</math>

=== Produsul scalar, forma biliniară și sesquiliniară ===
[[Produs scalar |Produsul scalar]] al doi vectori coloană este elementul unic al produsului matricial
:<math>\mathbf x^\mathsf T \mathbf y,</math>
unde <math>\mathbf x^\mathsf T</math> este vectorul linie obținut prin [[matrice transpusă |transpunerea]] <math>\mathbf x</math>. (De obicei o matrice 1×1 este identificată prin unicul său element.)

În general, orice [[formă biliniară]] pe un spațiu vectorial de dimensiune finită poate fi exprimată ca produs matricial:
:<math>\mathbf x^\mathsf T \mathbf {Ay},</math>
și orice altă {{ill-wd| Q1931224||formă sesquiliniară}} poate fi exprimată prin
:<math>\mathbf x^\dagger \mathbf {Ay},</math>
unde <math>\mathbf x^\dagger</math> este [[Matrice adjunctă|adjuncta]] lui <math>\mathbf x</math> ([[transpusa conjugată]]).

== Proprietăți generale ==
Înmulțirea matricială are unele proprietăți asemănătoare cu [[înmulțire (matematică) |înmulțirea]] obișnuită. Totuși, înmulțirea matricială nu este definită dacă numărul de coloane al primului factor diferă de numărul de linii al celui de-al doilea factor și este [[comutativitate |necomutativă]],<ref name=W>{{en icon}} {{Cite web|last=Weisstein|first=Eric W.|title=Matrix Multiplication|url=https://mathworld.wolfram.com/MatrixMultiplication.html|access-date=2020-09-06|website=mathworld.wolfram.com}}</ref> chiar și când produsul rămâne definit după schimbarea ordinii factorilor.<ref>{{en icon}} {{cite book|title=Linear Algebra|edition=4th| first1 = S. | last1 = Lipcshutz | first2 = M. | last2 = Lipson|series=Schaum's Outlines|publisher=McGraw Hill (USA)|chapter=2|date=2009|isbn=978-0-07-154352-1}}</ref><ref>{{en icon}} {{cite book|title=Matrix Analysis| last = Horn | first =  Johnson|edition=2nd|chapter=Chapter 0|publisher=Cambridge University Press|date=2013|isbn=978-0-521-54823-6}}</ref>

=== Necomutativitatea ===
O operație este comutativă dacă, fiind date două elemente {{math|'''A'''}} și {{math|'''B'''}} astfel încât produsul <math> \mathbf{A}\mathbf{B}</math> este definit, atunci <math>\mathbf{B}\mathbf{A}</math> este definit și <math>\mathbf{A}\mathbf {B}=\mathbf{B}\mathbf{A}.</math>

Înmulțirea matricială a două matrici {{math|'''A'''}} și {{math|'''B'''}} al căror produs matricial <math>\mathbf{A}\mathbf{B}</math> este definit, ar fi comutativă dacă produsul <math>\mathbf{B}\mathbf{A}</math> este și el definit, iar <math>\mathbf{A}\mathbf{B}=\mathbf{B}\mathbf{A}.</math>

Dacă {{math|'''A'''}} și {{math|'''B'''}} sunt matrici de dimensiunile {{mvar|m×n}} respectiv {{mvar|p×q}}, atunci <math>\mathbf{A}\mathbf{B}</math> este definit dacă {{mvar|1=n&nbsp;=&nbsp;p}}, iar <math>\mathbf{B}\mathbf{A}</math> este definit dacă {{mvar|1=m&nbsp;=&nbsp;q}}. Prin urmare, dacă unul dintre produse este definit, celălalt ar trebui să fie nedefinit. Dacă <math>m=q\neq n=p</math>, cele două produse sunt definite, dar au dimensiuni diferite; astfel că ele nu pot fi egale. Doar dacă {{mvar|1=m&nbsp;=&nbsp;q=&nbsp;n=&nbsp;p}}, adică dacă {{math|'''A'''}} și {{math|'''B'''}} sunt ambele [[matrice pătrată | pătrate]], de aceeași dimensiune, ambele produse sunt definite și au aceeași dimensiune. Chiar și în acest caz, în general
:<math>\mathbf{A}\mathbf{B} \neq \mathbf{B}\mathbf{A}.</math> 
De exemplu
:<math>\begin{pmatrix} 0 & 1 \\ 0 & 0 \end{pmatrix}\begin{pmatrix} 0 & 0 \\ 1 & 0 \end{pmatrix}=\begin{pmatrix} 1 & 0 \\ 0 & 0 \end{pmatrix},</math>
dar
:<math>\begin{pmatrix} 0 & 0 \\ 1 & 0 \end{pmatrix}\begin{pmatrix} 0 & 1 \\ 0 & 0 \end{pmatrix} = \begin{pmatrix} 0 & 0 \\ 0 & 1 \end{pmatrix}.</math>

Acest exemplu poate fi extins pentru a arăta că dacă {{math|'''A'''}} este o matrice {{mvar|n×n}} cu elementele într-un [[corp (matematică) |corp]] {{mvar|F}}, atunci <math>\mathbf{A}\mathbf{B} = \mathbf{B}\mathbf{A}</math> pentru orice matrice {{math|'''B'''}} {{mvar|n×n}} cu elementele din {{mvar|F}} [[dacă și numai dacă]] <math>\mathbf{A}=c\,\mathbf{I}</math> unde <math>c\in F</math> și {{math|'''I'''}} este [[matrice unitate |matricea unitate]] {{mvar|n×n}}. Dacă, în loc de un corp, se presupune că elementele aparțin unui [[inel (matematică) |inel]], atunci trebuie adăugată condiția ca {{mvar|c}} să aparțină centrului inelului.

Un caz particular în care comutativitatea apare este atunci când {{math|'''D'''}} și {{math|'''E'''}} sunt două [[matrice diagonală |matrici diagonale]] pătrate de aceeași mărime. Atunci {{math|1='''DE''' = '''ED'''}}.<ref name=W/> Din nou, dacă matricile sunt peste un inel generic în loc de a fi peste un corp, elementele corespunzătoare ale fiecăreia trebuie, de asemenea, să fie comutative unul față de celălalt pentru ca acest lucru să fie valabil.

=== Distributivitatea ===
Produsul matricial este [[distributivitate |distributiv]] față de [[adunarea matricilor]]. Adică dacă {{math|'''A''', '''B''', '''C''', '''D'''}} sunt matrici de dimensiunile respective {{mvar|m×n}}, {{mvar|n×p}}, {{mvar|n×p}}, și {{mvar|p×q}}, există atât distributivitatea [[la stânga și la dreapta |la stânga]]<ref name=W/>
:<math>\mathbf{A}(\mathbf{B} + \mathbf{C}) = \mathbf{AB} + \mathbf{AC},</math>
cât și [[la stânga și la dreapta |la dreapta]]<ref name=W/>
:<math>(\mathbf{B} + \mathbf{C} )\mathbf{D} = \mathbf{BD} + \mathbf{CD}.</math>

Aceasta rezultă din distributivitatea coeficienților 
:<math>\sum_k a_{ik}(b_{kj} + c_{kj}) = \sum_k a_{ik}b_{kj} + \sum_k a_{ik}c_{kj},</math>
:<math>\sum_k (b_{ik} + c_{ik}) d_{kj} = \sum_k b_{ik}d_{kj} + \sum_k c_{ik}d_{kj}.</math>

=== Asociativitatea ===
Fiind date matricile {{math|'''A''', '''B'''}} și {{math|'''C'''}}, produsele {{math|('''AB''')'''C'''}} și {{math|'''A'''('''BC''')}} sunt definite dacă și numai dacă numărul de coloane din {{math|'''A'''}} este egal cu numărul de linii din {{math|'''B'''}} și numărul de coloane din {{math|'''B'''}} este egal cu numărul de linii din {{math|'''C'''}} (în special, dacă unul dintre produse este definit, atunci celălalt este și el definit). În acest caz există [[asociativitate]]a
:<math>(\mathbf{AB})\mathbf{C}=\mathbf{A}(\mathbf{BC}).</math>

Ca la orice operație asociativă, aceasta permite omiterea parantezelor și scrierea produselor de mai sus ca <math>\mathbf{ABC}.</math>

Aceasta se extinde natural la produsul a oricâte matrici, cu condiția ca dimensiunile să se potrivească. Adică dacă {{math|'''A'''<sub>1</sub>, '''A'''<sub>2</sub>, ... , '''A''' <sub>''n''</sub>}} sunt matrici astfel încât numărul de coloane ale {{math|'''A'''<sub>''i''</sub>}} este egal numărul de linii ale {{math|'''A'''<sub>''i'' + 1</sub>}} pentru {{math|1=''i'' = 1, ... , ''n''−1}}, atunci produsul
:<math> \prod_{i=1}^n \mathbf{A}_i = \mathbf{A}_1\mathbf{A}_2\cdots\mathbf{A}_n </math>
este definit și nu depinde de ordinea în care se fac înmulțirile, cât timp ordinea matricilor este păstrată fixă.

Aceste proprietăți pot fi demonstrate prin operații de [[adunare]] simple, dar lungi. Acest rezultat rezultă și din faptul că matricile sunt aplicații liniare. Prin urmare, proprietatea asociativă a matricilor este pur și simplu un caz particular al proprietății asociative a {{ill-wd| Q244761||compunerii funcțiilor}}.

=== Produsul cu un scalar ===
Dacă {{math|'''A'''}} este o matrice și {{mvar|c}} un scalar, atunci matricile <math>c\mathbf{A}</math> și <math>\mathbf{A}c</math> se obțin înmulțind [[la stânga și la dreapta |la stânga sau la dreapta]] toate elementele lui {{math|'''A'''}} cu {{mvar|c}}. Dacă scalarii sunt comutativi, atunci <math>c\mathbf{A} = \mathbf{A}c.</math>

Dacă produsul <math>\mathbf{AB}</math> este definit (adică numărul de coloane din {{math|'''A'''}} este egal cu numărul de linii din {{math|'''B'''}}), atunci
:<math> c(\mathbf{AB}) = (c \mathbf{A})\mathbf{B}</math> &nbsp; și &nbsp; <math> (\mathbf{A} \mathbf{B})c=\mathbf{A}(\mathbf{B}c).</math>

Dacă scalarii sunt comutativi, atunci toate cele patru matrici sunt egale. În general, toate cele patru sunt egale dacă {{math|''c''}} aparține centrului unui inel care conține elementele matricei, deoarece în acest caz, {{math|1=''c''&hairsp;'''X''' = '''X'''&hairsp;''c''}} pentru toate matricile {{math|'''X'''}}.

Aceste proprietăți rezultă din [[aplicație biliniară |biliniaritatea]] produsului scalarilor:
:<math>c \left(\sum_k a_{ik}b_{kj}\right) = \sum_k (c a_{ik} ) b_{kj} </math>
:<math>\left(\sum_k a_{ik}b_{kj}\right) c = \sum_k a_{ik} ( b_{kj}c). </math>

=== Operații cu matricea transpusă ===
Dacă scalarii sunt comutativi, [[Matrice transpusă |transpusa]] produsului matricial este produsul, în ordine inversă, al transpuselor factorilor. Acesta este
:<math> (\mathbf{AB})^\mathsf{T} = \mathbf{B}^\mathsf{T}\mathbf{A}^\mathsf{T} </math>
unde cu <sup>T</sup> este notată transpusa.

Această identitate nu este valabilă pentru elementele necomutative, deoarece ordinea dintre elementele lui {{math|'''A'''}} și {{math|'''B'''}} este inversată atunci când se extinde definirea produsului matricial.

=== Conjugata complexă ===
Dacă {{math|'''A'''}} și {{math|'''B'''}} au elemente [[număr complex |complexe]], atunci
:<math> (\mathbf{AB})^* = \mathbf{A}^*\mathbf{B}^* </math>
unde cu {{math|<sup>*</sup>}} sunt notate [[conjugată complexă |conjugatele complexe]] ale elementelor matricei.

Acest lucru rezultă din aplicarea la definiția produsului matricial a faptului că conjugatul unei sume este suma conjugatelor sumelor, iar conjugatul unui produs este produsul conjugatelor factorilor.

Transpunerea acționează asupra indicilor elementelor, în timp ce conjugarea acționează independent asupra elementelor în sine. Rezultă că, dacă {{math|'''A'''}} și {{math|'''B'''}} au elemente complexe, există
:<math> (\mathbf{AB})^\dagger = \mathbf{B}^\dagger\mathbf{A}^\dagger ,</math>
unde cu {{math|<sup>†</sup>}} sunt notate [[Matrice adjunctă|adjunctele]] (conjugatele transpuselor, sau, echivalent, transpusele conjugatelor).

== Matrici pătrate ==
Fie <math>\mathcal M_n(R)</math> mulțimea de [[matrice pătrată |matrici pătrate]] {{mvar|n×n}} cu elementele din inelul {{mvar|R}}, care, în practică, este adesea un corp.

În <math>\mathcal M_n(R)</math>, produsul este definit pentru fiecare pereche de matrici. Acest lucru face din <math>\mathcal M_n(R)</math> un inel, care are [[matrice unitate |matricea unitate]] {{math|'''I'''}} ca [[element neutru]] (matricea ale cărei elemente pe [[Diagonală principală|diagonala principală]] sunt egale cu 1 și toate celelalte elemente sunt 0). Acest inel este, de asemenea, o {{ill-wd| Q744960||''R''-algebră asociativă}}.

Dacă {{math|''n'' > 1}}, multe matrici nu au o [[inversarea matricilor |inversă]] multiplicativă. De exemplu, o matrice care are toate elementele unei linii (sau unei coloane) 0 nu are inversă. Dacă există, inversa unei matrice {{math|'''A'''}} se notează {{math|'''A'''<sup>−1</sup>}} și verifică relația
:<math> \mathbf{A}\mathbf{A}^{-1} = \mathbf{A}^{-1}\mathbf{A} = \mathbf{I}. </math>
O matrice care are o inversă este o matrice inversabilă. În caz contrar, este o matrice singulară.

Un produs matricial este inversabil dacă și numai dacă fiecare factor este inversabil. În acest caz există relația
:<math>(\mathbf{A}\mathbf{B})^{-1} = \mathbf{B}^{-1}\mathbf{A}^{-1}.</math>

Când {{mvar|R}} este [[inel comutativ|comutativ]] și, în special, când este un corp, [[determinant (matematică) |determinantul]] unui produs este produsul determinanților. Deoarece determinanții sunt scalari, iar scalarii sunt comutativi, există relația
:<math> \det(\mathbf{AB}) = \det(\mathbf{BA}) =\det(\mathbf{A})\det(\mathbf{B}).</math>

Ceilalți {{ill-wd| Q188211||invarianți}} matriciali nu se comportă la fel de bine cu produsele. Totuși, dacă {{mvar|R}} este comutativ, {{math|'''AB'''}} și {{math|'''BA'''}} au aceeași [[urmă (algebră) |urmă]], același {{ill-wd| Q849705||polinom caracteristic}} și aceleași [[vectori și valori proprii |valori proprii]] cu aceleași [[multiplicitate |multiplicități]]. Totuși, [[vectori și valori proprii |vectorii proprii]] sunt în general diferiți dacă {{math|'''AB''' ≠ '''BA'''}}.

=== Puterea unei matrice pătrate ===
O matrice pătrată poate fi ridicată la orice [[putere (matematică) |putere]] [[număr întreg |întreagă]] [[număr pozitiv |nenegativă]] înmulțind-o cu ea însăși în mod repetat, în același mod ca pentru numerele obișnuite. Acesta este,
:<math>\mathbf{A}^0 = \mathbf{I},</math>
:<math>\mathbf{A}^1 = \mathbf{A},</math>
:<math>\mathbf{A}^k = \underbrace{\mathbf{A}\mathbf{A}\cdots\mathbf{A}}_{\text{de 𝘬 ori}}.</math>

Calcularea celei de a {{mvar|k}}-a putere a unei matrice, dacă se face cu algoritmul trivial (înmulțire repetată) necesită de {{math|''k''−1}} ori timpul unei singure înmulțiri matriciale. Deoarece acest lucru consumă foarte mult timp, se preferă metoda {{ill-wd| Q864127||ridicării la putere rapide}}, care necesită mai puțin de {{math|2 log<sub>2</sub> ''k''}} înmulțiri matriciale, deci este mult mai eficientă.

Un caz ușor de ridicare la putere este cel al unei matrice diagonale. Deoarece produsul matricilor diagonale echivalează cu simpla înmulțire a elementelor diagonale corespunzătoare, a {{mvar|k}}-a putere a unei matrici diagonale se obține prin ridicarea elementelor ei la puterea {{mvar|k}}:
:<math>
  \begin{bmatrix}
    a_{11} & 0      & \cdots & 0 \\
    0      & a_{22} & \cdots & 0 \\
    \vdots & \vdots & \ddots & \vdots \\
    0     & 0       & \cdots & a_{nn}
  \end{bmatrix}^k =
\begin{bmatrix}
    a_{11}^k & 0        & \cdots & 0 \\
    0        & a_{22}^k & \cdots & 0 \\
    \vdots   & \vdots   & \ddots & \vdots \\
    0        & 0        & \cdots & a_{nn}^k
  \end{bmatrix}.
</math>

=== Aplicație la matrici asemenea ===
Orice [[inversarea matricilor |matrice inversabilă]] <math>\mathbf{P}</math> definește o [[asemănarea matricilor |transformare de asemănare]] (pe matrici pătrate de aceeași dimensiune ca <math>\mathbf{P}</math>)
:<math>S_\mathbf{P}(\mathbf{A}) = \mathbf{P}^{-1} \mathbf{A} \mathbf{P}.</math>

Transformările de asemănare aplică produsul la factori, adică
:<math>S_\mathbf{P}(\mathbf{AB}) = S_\mathbf{P}(\mathbf{A})S_\mathbf{P}(\mathbf{B}).</math>

De fapt,
:<math>\mathbf{P}^{-1} (\mathbf{AB}) \mathbf{P} 
= \mathbf{P}^{-1} \mathbf{A}(\mathbf{P}\mathbf{P}^{-1})\mathbf{B} \mathbf{P}
=(\mathbf{P}^{-1} \mathbf{A}\mathbf{P})(\mathbf{P}^{-1}\mathbf{B} \mathbf{P}).</math>

== Complexitatea de calcul ==
[[Fișier:MatrixMultComplexity svg.svg|thumb|upright=2.0|Îmbunătățirea estimărilor exponentului {{math|ω}} în timp pentru complexitatea de calcul a înmulțirii matriciale <math>O(n^\omega)</math>]]
Un [[algoritm]] de înmulțire a matricilor care rezultă din definiție necesită în cazul cel mai rău {{math|''n''<sup>3</sup>}} înmulțiri și {{math|(''n''−1)''n''<sup>2</sup>}} adunări de scalari pentru a calcula produsul a două matrice pătrate {{math|''n''×''n''}}. Într-un model de calcul în care operațiile scalare au timp constant {{ill-wd| Q5157286||complexitatea sa de calcul}} este deci {{math|O(''n''<sup>3</sup>}}).

Surprinzător, această complexitate nu este optimă, așa cum a arătat în 1969 [[Volker Strassen]], care a furnizat un algoritm, numit acum {{ill-wd| Q728507||algoritmul Strassen}}, cu o complexitate de <math>O( n^{\log_{2}7}) \approx O(n^{2,8074}).</math><ref>
{{cite journal 
 | doi=10.1007/BF02165411 
 | url=http://www.digizeitschriften.de/dms/img/?PID=GDZPPN001168215 
 | author=Volker Strassen 
 | title=Gaussian elimination is not optimal 
 | journal=Numerische Mathematik 
 | volume=13 
 | pages=354&ndash;356 
 | date=Aug 1969 
 | issue = 4
 | s2cid = 121656251
}}</ref>

Din 2020, cel mai bun algoritm de înmulțire matricială era cel dat de Josh Alman și Virginia Vassilevska Williams, cu complexitatea {{math|''O''(''n''<sup>2.3728596</sup>)}}.<ref name="Alman2020">
{{Citation
 | last1=Alman
 | first1=Josh
 | last2=Williams
 | first2=Virginia Vassilevska
 | contribution=A Refined Laser Method and Faster Matrix Multiplication
 | year = 2020
 | arxiv=2010.05846
 | title = 32nd Annual ACM-SIAM Symposium on Discrete Algorithms (SODA 2021)
 |url=https://www.siam.org/conferences/cm/program/accepted-papers/soda21-accepted-papers
}}</ref>

Nu se știe dacă înmulțirea matricială poate fi efectuată în timp {{math|''n''<sup>2 + o(1)</sup>}}. Acest lucru ar fi optim, deoarece trebuie citite cele {{math|''n''<sup>2</sup>}} elementele unei matrice pentru a o înmulți cu o altă matrice.

Deoarece înmulțirea matricială formează baza pentru mulți algoritmi și multe operații pe matrici chiar au aceeași complexitate ca și înmulțirea matricială ([[până la]] o constantă multiplicativă), complexitatea de calcul a înmulțirii matriciale este o chestiune importantă în {{ill-wd| Q2004891||algebra liniară numerică}} și {{ill-wd| Q2878974 ||informatica teoretică}}.

=== Complexitatea de calcul în funcție de ordinea operațiilor ===
Deși rezultatul unei secvențe de produse matrice nu depinde de [[permutare|ordinea]] efectuării produselor (cu condiția ca ordinea matricelor să nu fie schimbată), complexitatea de calcul poate depinde mult de această ordine.

De exemplu, dacă {{math|'''A''', '''B'''}} și {{math|'''C'''}} sunt matrici cu dimensiunile {{math|10×30, 30×5, 5×60}}, calculul {{math|('''AB''')'''C'''}} necesită {{math|1=10×30×5 + 10×5×60 = 4500}} de înmulțiri, în timp ce calculul {{math|'''A'''('''BC''')}} necesită {{math|1=30×5×60 + 10×30×60 = 27&nbsp;000}} de înmulțiri.

Au fost concepuți algoritmi pentru alegerea ordinii [[optimizare|optime]] de efectuare a produselor, adică printr-un număr minim de operații. Când numărul {{mvar|n}} de matrici crește, s-a demonstrat că alegerea ordinii optime are o complexitate de <math>O(n \log n).</math>

== Alte tipuri de produse ale matricelor ==
Înmulțirea „standard” a matricilor prezentată în articolul de față este singurul mod de înmulțire al matricilor studiat în învățământul preuniversitar din România.<ref>{{citat carte |first1=Marius |last1=Burtea |first2=Georgeta |last2=Burtea |url=https://manuale.edu.ro/manuale/Clasa%20a%20XI-a/Matematica/Carminis/A196.pdf |title=Matematică: Manual pentru clasa a XI-a |loc=Pitești |editor=Carminis |an=2006 |ISBN=978-973-123-417-5 |pagini=22–26, accesat 2023-03-18 |access-date=2023-04-16 |archive-date=2022-02-15 |archive-url=https://web.archive.org/web/20220215100505/https://manuale.edu.ro/manuale/Clasa%20a%20XI-a/Matematica/Carminis/A196.pdf |dead-url=yes }}</ref> Însă există și alte tipuri de produse de matrici:
* [[produs cracovian]], definit drept {{math|1='''A''' ∧ '''B''' = '''B'''<sup>T</sup>'''A'''}}
* [[produs interior Frobenius]], [[produs scalar]] al matricilor considerate drept vectori, sau, echivalent suma produselor Hadamard ale elementelor
* [[produs Hadamard]], elementele produsului sunt produsele element cu element
* [[produs Kronecker]] și [[produs tensorial]], generalizare a precedentelor
* [[produs Khatri–Rao]]
* [[produs extern]] (sau [[produs diadic]]) a două matrici coloană, <math>\mathbf{a}\mathbf{b}^\mathsf{T}</math>

== Note ==
<references />

== Bibliografie ==
* {{en icon}} Henry Cohn, Robert Kleinberg, Balázs Szegedy, and Chris Umans. [[teoria grupurilor|Group-theoretic]] Algorithms for Matrix Multiplication. {{arxiv|math.GR/0511460}}. ''Proceedings of the 46th Annual Symposium on Foundations of Computer Science'', 23–25 October 2005, Pittsburgh, PA, IEEE Computer Society, pp.&nbsp;379–388.
* {{en icon}} Henry Cohn, Chris Umans. A Group-theoretic Approach to Fast Matrix Multiplication. {{arxiv|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.&nbsp;438–449.
* {{en icon}} {{cite journal | last1 = Coppersmith | first1 = D. | last2 = Winograd | first2 = S. | year = 1990 | title = Matrix multiplication via arithmetic progressions | journal = J. Symbolic Comput. | volume = 9 | issue = 3| pages = 251–280 | doi=10.1016/s0747-7171(08)80013-2| doi-access = free }}
* {{en icon}} {{Citation | last1=Horn | first1=Roger A. | last2=Johnson | first2=Charles R. | title=Topics in Matrix Analysis | publisher=[[Cambridge University Press]] | isbn=978-0-521-46713-1 | year=1991}}
* {{en icon}} [[Donald Knuth|Knuth, D.E.]], ''[[The Art of Computer Programming]] Volume 2: Seminumerical Algorithms''. Addison-Wesley Professional; 3 edition (November 14, 1997). {{isbn|978-0-201-89684-8}}. p.&nbsp;501
* {{en icon}} {{Citation | last1=Press | first1=William H. | last2=Flannery | first2=Brian P. | last3=Teukolsky | first3=Saul A. | author3-link=Saul Teukolsky | last4=Vetterling | first4=William T. | title=Numerical Recipes: The Art of Scientific Computing | publisher=[[Cambridge University Press]] | edition=3rd | isbn=978-0-521-88068-8 | year=2007}}
* {{en icon}} Ran Raz, ''On the complexity of matrix product''. In Proceedings of the thirty-fourth annual ACM symposium on Theory of computing. ACM Press, 2002 {{doi|10.1145/509907.509932}}.
* {{en icon}} Robinson, Sara, ''Toward an Optimal Algorithm for Matrix Multiplication,'' SIAM News 38(9), November 2005. [https://web.archive.org/web/20100331095603/http://www.siam.org/pdf/news/174.pdf PDF]
* {{en icon}} Strassen, Volker, ''Gaussian Elimination is not Optimal'', Numer. Math. 13, p.&nbsp;354-356, 1969
* {{en icon}} {{Citation | doi=10.1016/0024-3795(73)90023-2 | last=Styan | first=George P. H. | title=Hadamard Products and Multivariate Statistical Analysis | journal=Linear Algebra and Its Applications | year=1973 | volume=6 | pages=217–240| url=http://dml.cz/bitstream/handle/10338.dmlcz/102190/CzechMathJ_37-1987-4_14.pdf | doi-access=free }}
* {{en icon}} {{Cite book|last=Williams|first=Virginia Vassilevska|date=2012-05-19|chapter=Multiplying matrices faster than coppersmith-winograd|chapter-url=http://dl.acm.org/citation.cfm?id=2213977.2214056|publisher=ACM|pages=887–898|doi=10.1145/2213977.2214056|isbn=9781450312455|title=Proceedings of the 44th symposium on Theory of Computing - STOC '12|citeseerx=10.1.1.297.2680|s2cid=14350287}}

== Legături externe ==
* {{Commonscat-inline}}
{{Portal|Matematică}}
{{Algebră liniară}}

[[Categorie:Înmulțire]]
[[Categorie:Algebră liniară]]
[[Categorie:Operatori liniari]]