Násobení matic

Pro součin matic musí být počet sloupců v první matici roven počtu řádků ve druhé matici. Výsledná matice má počet řádků první matice a počet sloupců druhé matice.

Součin matic hovorově též maticové násobení (neplést se skalárním násobkem matice) je v matematice zobecnění součinu čísel na matice. Formálně se dá definovat jako binární operace na maticích odpovídajících typů. Využívá se v matematice, fyzice a jejich aplikacích, obvykle pro popis skládání lineárních zobrazení.

Speciálním případem násobení matic je součin matice typu $m\times n$ a vektoru braného jako matice o typu $n\times 1$ (sloupcový vektor). Tento součin lze interpretovat jako aplikaci lineárního zobrazení reprezentovaného transformační maticí na vektor.

Formální definice

Pokud je $𝑨$ matice typu $m\times n$ a $𝑩$ je matice typu $n\times p$, jejich součin $𝑨\cdot 𝑩$ je matice typu $m\times p$ definovaná vztahem

$(𝑨\cdot 𝑩{)}_{ij}=\sum _{k=1}^n{a}_{ik}{b}_{kj}=a_{i1}{b}_{1j}+a_{i2}{b}_{2j}+\cdots +a_{in}{b}_{nj}.$

pro všechny prvky výsledné matice indexované $i\in \{1,\dots ,m\}$ a $j\in \{1,\dots ,p\}$.

Ve většině případů jsou prvky matice čísla, ale mohou to být jakékoli druhy matematických objektů, pro které je definováno sčítání a násobení, které jsou asociativní a takové, že sčítání je komutativní a násobení je distributivní s ohledem na sčítání, typicky prvky nějakého tělesa. Prvky mohou být dokonce samotné matice (bloková matice).

U reálných matic lze prvek v $i$-tém řádku a $j$-tém sloupci výsledné matice lze také chápat jako standardní skalární součin vektoru $i$-tého řádku první matice s vektorem $j$-tého sloupce druhé matice.

Tečka $\cdot$ se v součinu vynechává a píše se pouze $𝑨𝑩$.

Ukázka výpočtu

Schéma součinu $𝑨𝑩$ dvou matic $𝑨$ a $𝑩$.

Součin matic $𝑨=\left(\begin{array}{ccc}1& 2& 3\\ 4& 5& 6\end{array}\right)$ a $𝑩=\left(\begin{array}{cc}1& 2\\ 3& 4\\ 5& 6\end{array}\right)$ je

$𝑨𝑩=\left(\begin{array}{cc}(1\cdot 1+2\cdot 3+3\cdot 5)& (1\cdot 2+2\cdot 4+3\cdot 6)\\ (4\cdot 1+5\cdot 3+6\cdot 5)& (4\cdot 2+5\cdot 4+6\cdot 6)\end{array}\right)=\left(\begin{array}{cc}22& 28\\ 49& 64\end{array}\right)$

Prvky matice $𝑨$ zůstávají v řádcích tak, jak jsou, a prvky v matici $𝑩$ se rozmístí opět do levého a pravého sloupce.

Použití

Historicky bylo násobení matic zavedeno pro usnadnění a objasnění výpočtů v lineární algebře. Tento silný vztah mezi maticovým součinem a lineární algebrou zůstává je fundamentální v celé matematice, stejně jako ve fyzice, chemii, inženýrství a informatice.

Soustavy lineárních rovnic

Obecný tvar soustavy lineárních rovnic je

$\begin{array}{ccc}{a}_{11}{x}_1+a_{12}{x}_2+\dots +a_{1n}{x}_n& =& b_1\\ a_{21}{x}_1+a_{22}{x}_2+\dots +a_{2n}{x}_n& =& b_2\\ & ⋮& \\ a_{m1}{x}_1+a_{m2}{x}_2+\dots +a_{mn}{x}_n& =& b_m\end{array}$

Při použití stejné notace jako výše je zápis soustavy ekvivalentní jednoduché maticové rovnici

$𝑨𝒙=𝒃$.

Lineární zobrazení

Pokud má vektorový prostor konečnou bázi, každý z jeho vektorů je jednoznačně reprezentován konečnou posloupností skalárů, nazývanou vektor souřadnic, tvořenou souřadnicemi vektoru vzhledem k bázi. Tyto vektory souřadnic tvoří další vektorový prostor, který je izomorfní původnímu vektorovému prostoru. Vektor souřadnic je běžně zapisován jako sloupcový vektor, což je matice pouze s jedním sloupcem. Sloupcový vektor pak představuje jak souřadnicový vektor, tak i vektor původního vektorového prostoru.

Lineární zobrazení $A$ prostoru dimenze $n$ do vektorového prostoru dimenze $m$ převádí sloupcový vektor

$𝒙=\left(\begin{array}{c}{x}_1\\ x_2\\ ⋮\\ x_n\end{array}\right)$

na sloupcový vektor

$𝒚=A(𝒙)=\left(\begin{array}{c}{a}_{11}{x}_1+\cdots +a_{1n}{x}_n\\ a_{21}{x}_1+\cdots +a_{2n}{x}_n\\ ⋮\\ a_{m1}{x}_1+\cdots +a_{mn}{x}_n\end{array}\right).$

Lineární zobrazení $A$ je proto definováno maticí

$𝑨=\left(\begin{array}{cccc}{a}_{11}& a_{12}& \cdots & a_{1n}\\ a_{21}& a_{22}& \cdots & a_{2n}\\ ⋮& ⋮& \ddots & ⋮\\ a_{m1}& a_{m2}& \cdots & a_{mn}\end{array}\right),$

a zobrazuje sloupcový vektor $𝒙$ na maticový součin

$𝒚=𝑨𝒙$.

Je-li $B$ další lineární zobrazení z předchozího vektorového prostoru dimenze $m$, do vektorového prostoru dimenze $p$, pak jej lze reprezentovat maticí $𝑩$ řádu $p\times m$. Přímý výpočet ukazuje, že matice složeného zobrazení $B\circ A$ je rovna součinu $𝑩𝑨$. Obecný vzorec $(B\circ A)(𝒙)=B(A(𝒙))$, který definuje složené zobrazení, je jedním z specifických případů asociativity maticového součinu:

$(𝑩𝑨)𝒙=𝑩(𝑨𝒙)=𝑩𝑨𝒙.$

Geometrické rotace

Při použití systému kartézských souřadnic v euklidovské rovině je rotace o úhel $\alpha$ kolem počátku (počátek odpovídá nulovému vektoru) lineární zobrazení. Přesněji, $$\left(\begin{array}{c}{x}^{\prime }\\ y^{\prime }\end{array}\right)=\left(\begin{array}{cc}\cos \alpha & -\sin \alpha \\ \sin \alpha & \cos \alpha \end{array}\right)\left(\begin{array}{c}x\\ y\end{array}\right),$$ kde výchozí bod $(x,y)$ i jeho obraz $(x^{\prime },y^{\prime })$ jsou zapsány jako sloupcové vektory.

Složení rotací o úhel $\alpha$ a pak o úhel $\beta$ odpovídá maticovému součinu $$\left(\begin{array}{cc}\cos \beta & -\sin \beta \\ \sin \beta & \cos \beta \end{array}\right)\left(\begin{array}{cc}\cos \alpha & -\sin \alpha \\ \sin \alpha & \cos \alpha \end{array}\right)=\left(\begin{array}{cc}\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{array}\right)=\left(\begin{array}{cc}\cos (\alpha +\beta )& -\sin (\alpha +\beta )\\ \sin (\alpha +\beta )& \cos (\alpha +\beta )\end{array}\right),$$ ve druhé rovnosti jsou použity součtové vzorce. Výsledné složení odpovídá rotaci o úhel $\alpha +\beta$, jak lze očekávat.

Skalární součin, bilineární forma a seskvilineární forma

Standardní skalární součin dvou reálných sloupcových vektorů lze zapsat maticovým součinem

${𝒙}^{𝖳}𝒚,$

kde ${𝒙}^{𝖳}$ je řádkový vektor získaný pomocí transpozice $𝒙$. (Výsledná matice $1\times 1$ je zde ztotožněna se svým jediným prvkem.)

Obecněji lze jakoukoli bilineární formu ve vektorovém prostoru konečného rozměru vyjádřit jako maticový součin

${𝒙}^{𝖳}𝑨𝒚,$

a jakoukoliv seskvilineární formu lze vyjádřit jako

${𝒙}^{𝖧}𝑨𝒚,$

kde ${𝒙}^{𝖧}$ je hermitovsky sdružený vektor k vektoru $𝒙$.

Alokace zdrojů v ekonomii

Výpočet prvku v levém dolním rohu součinu $𝑨𝑩$ odpovídá všem možným cestám (zvýrazněné) od suroviny $b_4$ ke konečnému výrobku $f_1$ v grafu toku výroby.

Jako příklad si představme fiktivní továrnu, která používá 4 druhy surovin $b_1,b_2,b_3,b_4$ k výrobě 3 meziproduktů, $m_1,m_2,m_3$, které se následně používají k výrobě 3 druhů výrobků, $f_1,f_2,f_3$.

Matice $𝑨=\left(\begin{array}{ccc}1& 0& 1\\ 2& 1& 1\\ 0& 1& 1\\ 1& 1& 2\end{array}\right)$   a   $𝑩=\left(\begin{array}{ccc}1& 2& 1\\ 2& 3& 1\\ 4& 2& 2\end{array}\right)$ udávají množství surovin potřebných pro výrobu meziproduktů, respektive množství meziproduktů potřebných pro výsledné výrobky. Například k výrobě jednoho meziproduktu $m_1$ je třeba jedna jednotka suroviny $b_1$, dvě jednotky $b_2$, žádné $b_3$ a jedna jednotka $b_4$, což odpovídá prvnímu sloupci matice $𝑨$.

Součin $𝑨𝑩=\left(\begin{array}{ccc}5& 4& 3\\ 8& 9& 5\\ 6& 5& 3\\ 11& 9& 6\end{array}\right)$ pak přímo udává množství surovin potřebných pro výrobu jednotlivých výrobků. Například prvek v levém dolním rohu $𝑨𝑩$ je vypočítán jako $1\cdot 1+1\cdot 2+2\cdot 4=11$, což odpovídá tomu, že $11$ jednotek $b_4$ je potřeba k výrobě jednoho výrobku $f_1$. Jmenovitě jedna jednotka $b_4$ je třeba pro $m_1$, 2 pro $m_2$ a $4$ pro každý ze dvou meziproduktů $m_3$, které jsou potřeba pro jeden kus $f_1$, viz obrázek.

Aby bylo možné vyrobit např. 100 výrobků $f_1$, 80 $f_2$ a 60 $f_3$, lze potřebné množství surovin vypočítat jako

$(𝑨𝑩)\left(\begin{array}{c}100\\ 80\\ 60\end{array}\right)=\left(\begin{array}{c}1000\\ 1820\\ 1180\\ 2180\end{array}\right),$

tj. $1000$ jednotek $b_1$, $1820$ jednotek $b_2$, $1180$ jednotek $b_3$ a $2180$ jednotek $b_4$.Matice součinu $𝑨𝑩$ může být použita k výpočtu množství surovin i pro jiné počty výrobků.

Vlastnosti

Rovnosti uvedené v následujících odstavcích platí, pokud mají výsledky operací smysl.

• Součin matice $𝑨$ s jednotkovou maticí $𝐈$ zprava i zleva má za výsledek matici $𝑨$, tj. $𝐈𝑨=𝑨𝐈=𝑨$.
• Maticový součin je asociativní, tedy $𝑨(𝑩𝑪)=(𝑨𝑩)𝑪$.
• Maticový součin není komutativní, tedy existují příklady matic, pro něž platí $𝑨𝑩\ne 𝑩𝑨$.
• Maticový součin je distributivní vůči sčítání, tj. $𝑨(𝑩+𝑪)=𝑨𝑩+𝑨𝑪$.
• Maticový součin je lineární vůči násobení skalárem (typicky reálné nebo komplexní číslo), tj. $𝑨\left(c𝑩\right)=(c𝑨)𝑩=c(𝑨𝑩)$.
• Matice vzhledem k součinu mohou být dělitelé nuly, tj. součin dvou nenulových matic může být nulová matice, například

$\left(\begin{array}{cc}1& 1\\ 2& 2\end{array}\right)\left(\begin{array}{c}-3\\ 3\end{array}\right)=\left(\begin{array}{c}0\\ 0\end{array}\right)$.

• Součin matic $𝑨$ typu $m\times n$ a $𝑩$ typu $n\times p$ lze vyjádřit jako

$𝑨𝑩={𝒂}_1{𝒃}_1^{𝖳}+{𝒂}_2{𝒃}_2^{𝖳}+\cdots +{𝒂}_n{𝒃}_n^{𝖳}$,
kde ${𝒂}_1,{𝒂}_2,\dots ,{𝒂}_n$ jsou sloupce matice $𝑨$ a ${𝒃}_1^{𝖳},{𝒃}_2^{𝖳},\dots ,{𝒃}_n^{𝖳}$ řádky matice $𝑩$. (Neboli ${𝒃}_1,{𝒃}_2,\dots ,{𝒃}_n$ jsou sloupce ${𝑩}^{𝖳}$.) Zde každý sčítanec ${𝒂}_i{𝒃}_i^{𝖳}$ je matice typu $m\times p$, protože sloupcové vektory odpovídají maticím o jednom sloupci.

• Transpozice součinu matic je součin transponovaných matic v opačném pořadí, tj. ${(𝑨𝑩)}^{𝖳}={𝑩}^{𝖳}{𝑨}^{𝖳}$
• Inverzní matice součinu regulárních matic je součin inverzních matic v opačném pořadí, tj. ${(𝑨𝑩)}^{-1}={𝑩}^{-1}{𝑨}^{-1}$
• Hermitovské sdružení (hermitovská transpozice) součinu matic je součin matic hermitovsky sdružených v opačném pořadí, tj. ${(𝑨𝑩)}^{𝖧}={𝑩}^{𝖧}{𝑨}^{𝖧}$
• Maticový součin odpovídá skládání lineárních zobrazení, které matice reprezentují.

Součiny čtvercových matic

Mocniny matice

Čtvercovou matici lze umocnit na jakoukoli nezápornou celočíselnou mocninu tím, že ji opakovaně násobíme stejným způsobem jako u běžných čísel, konkrétně

${𝑨}^0=𝐈,$

${𝑨}^1=𝑨,$

${𝑨}^k=\underset{k\text{ krát}}{\underbrace{𝑨𝑨\cdots 𝑨}}.$

Výpočet $k$-té mocniny matice potřebuje $k-1$ maticových součinů, pokud se provádí triviálním algoritmem (opakované násobení). Protože to může být velmi časově náročné, obecně se dává přednost použití umocňování pomocí druhé mocniny, které vyžaduje nejvýše ${\log }_{2}k$ maticových součinů, a je tedy mnohem efektivnější.

Snadným případem umocňování je diagonální matice. Protože součin diagonálních matic se rovná prostému vynásobení odpovídajících diagonálních prvků dohromady, získáme $k$-tou mocninu diagonální matice umocněním prvků na diagonále na $k$-tou:

${\left(\begin{array}{cccc}{a}_{11}& 0& \cdots & 0\\ 0& a_{22}& \cdots & 0\\ ⋮& ⋮& \ddots & ⋮\\ 0& 0& \cdots & a_{nn}\end{array}\right)}^k=\left(\begin{array}{cccc}{a}_{11}^k& 0& \cdots & 0\\ 0& a_{22}^k& \cdots & 0\\ ⋮& ⋮& \ddots & ⋮\\ 0& 0& \cdots & a_{nn}^k\end{array}\right).$

Regulární a singulární matice

Označme ${ℳ}_n(R)$ množinu čtvercových matic řádu $n$ s prvky z okruhu $R$, což je v praxi často těleso.

V ${ℳ}_n(R)$ je součin definován pro každou dvojici matic. Toto dělá z ${ℳ}_n(R)$ okruh, který má jednotkovou matici $𝐈$ za neutrální prvek.

Pokud je $n>1$, mnoho matic nemá inverzní prvek vůči násobení, např. nulová matice. Pokud inverzní prvek existuje, značí se ${𝑨}^{-1}$ a nazývá se inverzní matice k matici $𝑨$. Splňuje:

$𝑨{𝑨}^{-1}={𝑨}^{-1}𝑨=𝐈.$

Matice, která má inverzi, je regulární matice, někdy též invertibilní matice. Pokud inverzní matici nemá, nazývá se singulární matice.

Součin matic $𝑨𝑩$ je regulární, právě když je každý z činitelů $𝑨$ i $𝑩$ regulární. V tomto případě platí

$(𝑨𝑩{)}^{-1}={𝑩}^{-1}{𝑨}^{-1}.$

Determinant součinu

Determinant součinu čtvercových matic je součin jejich determinantů.

$\det (𝑨𝑩)=\det (𝑩𝑨)=\det 𝑨\det 𝑩$.

Tento vztah platí kdykoli je $R$ komutativní okruh, jmenovitě i v tělesech.

Výpočetní složitost

Výpočetní složitost výše popsaného algoritmu je $O(n^3)$ (počítáme $n^2$ čísel; pro každé potřebujeme $2n-1$ aritmetických operací). Existují však algoritmy s nižší složitostí vhodné pro matice vyšších řádů. Nejpoužívanější z nich je Strassenův algoritmus se složitostí $O(n^{{\log }_{2}7})\approx O(n^{2.807})$. Nižší složitost u tohoto algoritmu však získáváme za cenu snížené numerické stability. Asymptoticky nejrychlejší ze známých algoritmů je Coppersmithův-Winogradův algoritmus ($O(n^{2.376})$), který je však použitelný až pro matice tak velkých řádů, že je nelze zpracovávat pomocí současných počítačů.

Teoreticky by se dala složitost ještě snížit, ale nikdy nemůže být menší než $O(n^2)$, protože je třeba spočítat $n^2$ čísel.

Hledání nejkratší cesty v grafu

Algoritmy pro násobení matic s malou výpočetní složitostí lze využít i pro hledání nejkratší cesty v grafu z každého do každého vrcholu. To má v nejjednodušší podobě složitost $O(n^3)$. V tomto případě se však nepoužívá zde popsané násobení matic, ale upravená verze, kde je místo sčítání výběr nejmenšího prvku a místo násobení sčítání, proto nelze použít například Strassenův algoritmus, který využívá operaci odčítání jako inverzní operaci ke sčítání, která k operaci $\min$ není.

Graf lze popsat maticí vzdáleností $𝑨$. Pokud je pro výpočty operace sčítání dvou čísel definována jako jejich minimum, a místo násobení se použije sčítání, je možno matici nejkratších cest $𝑩$ získat jako (${𝑨}^n$) kde $n$ je řád matice vzdáleností. Při reálném výpočtu není třeba cyklicky násobit původní maticí, ale vždy se vynásobí vzniklé výsledky - nejkratší cesty jsou získány po ${\log }_2(n)$ násobeních. Je-li použit pro násobení algoritmus se složitostí menší než $O\left(\frac{n^3}{{\log }_2(n)}\right)$, složitost hledání cest se tímto postupem sníží.

Odkazy

Literatura

• Šablona:Citace monografie
• Šablona:Citace monografie
• Šablona:Citace monografie
• Šablona:Citace elektronické monografie
• Šablona:Citace elektronické monografie

Související články

• Hadamardův součin
• Skalární součin, Vnitřní součin
• Smíšený součin
• Strassenův algoritmus
• Tenzorový součin, Vnější součin
• Vektorový součin, Dvojitý vektorový součin

Externí odkazy

• Šablona:Commonscat
• Lineární algebra: algebra matic Šablona:Wayback Aplikace, která násobí a sčítá matice zadané uživatelem a zobrazuje postup výpočtu.