Regulární matice

Regulární[1], též invertibilní nebo nesingulární matice je v matematice čtvercová matice, která má inverzi. Regulární matice lze charakterizovat několika rovnocennými způsoby. Regulární matice jsou například charakterizovány tím, že lineární zobrazení, které popisují, jsou bijektivní. Proto má soustava lineárních rovnic s regulární maticí soustavy vždy jednoznačné řešení. Množina regulárních matic pevné velikosti nad okruhem nebo nad tělesem, spolu s maticovým součinem jako binární operace, tvoří obecnou lineární grupu.

Ne každá čtvercová matice má inverzi. Čtvercová matice, která nemá inverzní matici, se nazývá singulární matice.[1]

Definice

Čtvercová matice A R n × n {\displaystyle {\boldsymbol {A}}\in R^{n\times n}} s prvky z okruhu R {\displaystyle R} s jednotkovým prvkem (v praxi většinou obor reálných čísel ) se nazývá regulární, pokud existuje matice B R n × n {\displaystyle {\boldsymbol {B}}\in R^{n\times n}} taková, že

A B = B A = I {\displaystyle {\boldsymbol {A}}{\boldsymbol {B}}={\boldsymbol {B}}{\boldsymbol {A}}=\mathbf {I} } ,

kde I {\displaystyle I} označuje jednotkovou matici . Matice B {\displaystyle {\boldsymbol {B}}} je zde jednoznačně určena a nazývá se inverzní matice k matici A {\displaystyle {\boldsymbol {A}}} . Inverzní matice A {\displaystyle {\boldsymbol {A}}} se obvykle značí A 1 {\displaystyle {\boldsymbol {A}}^{-1}} . Pro singulární matici žádná taková matice B {\displaystyle {\boldsymbol {B}}} neexistuje.

Je R {\displaystyle R} komutativní okruh nebo těleso, jsou obě podmínky A B = I {\displaystyle {\boldsymbol {A}}{\boldsymbol {B}}=\mathbf {I} } a B A = I {\displaystyle {\boldsymbol {B}}{\boldsymbol {A}}=\mathbf {I} } ekvivalentní, tj., matice inverzní zleva je inverzní i zprava a naopak. V tomto případě lze podmínku oslabit na B A = I {\displaystyle {\boldsymbol {B}}{\boldsymbol {A}}=\mathbf {I} } nebo A B = I {\displaystyle {\boldsymbol {A}}{\boldsymbol {B}}=\mathbf {I} } .

Ukázky

Reálné matice

Matice

A = ( 2 3 1 2 ) {\displaystyle {\boldsymbol {A}}={\begin{pmatrix}2&3\\1&2\end{pmatrix}}}

je regulární, protože má inverzní matici

B = ( 2 3 1 2 ) {\displaystyle {\boldsymbol {B}}={\begin{pmatrix}2&-3\\-1&2\end{pmatrix}}} ,

neboť

A B = ( 2 3 1 2 ) ( 2 3 1 2 ) = ( 1 0 0 1 ) = I {\displaystyle {\boldsymbol {A}}{\boldsymbol {B}}={\begin{pmatrix}2&3\\1&2\end{pmatrix}}\cdot {\begin{pmatrix}2&-3\\-1&2\end{pmatrix}}={\begin{pmatrix}1&0\\0&1\end{pmatrix}}=\mathbf {I} } .

Na druhou stranu, matice

A = ( 2 3 0 0 ) {\displaystyle {\boldsymbol {A}}={\begin{pmatrix}2&3\\0&0\end{pmatrix}}}

je singulární, protože pro jakoukoli matici

B = ( a b c d ) {\displaystyle {\boldsymbol {B}}={\begin{pmatrix}a&b\\c&d\end{pmatrix}}}

platí

A B = ( 2 3 0 0 ) ( a b c d ) = ( 2 a + 3 c 2 b + 3 d 0 0 ) I {\displaystyle {\boldsymbol {A}}{\boldsymbol {B}}={\begin{pmatrix}2&3\\0&0\end{pmatrix}}\cdot {\begin{pmatrix}a&b\\c&d\end{pmatrix}}={\begin{pmatrix}2a+3c&2b+3d\\0&0\end{pmatrix}}\neq \mathbf {I} } .

Matice nad okruhy

Matice

A = ( 3 x 3 x 2 1 3 x 2 + 3 x ) {\displaystyle {\boldsymbol {A}}={\begin{pmatrix}3x^{3}&x^{2}-1\\3x^{2}+3&x\end{pmatrix}}}

s prvky z polynomiálního okruhu R = R [ x ] {\displaystyle R=\mathbb {R} [x]} má determinant:

det A = 3 x 3 x ( x 2 1 ) ( 3 x 2 + 3 ) = 3 {\displaystyle \det {\boldsymbol {A}}=3x^{3}\cdot x-(x^{2}-1)\cdot (3x^{2}+3)=3}

Protože prvek 3 {\displaystyle 3} je invertovatelný v okruhu R {\displaystyle R} , je i matice A {\displaystyle {\boldsymbol {A}}} regulární v R 2 × 2 {\displaystyle R^{2\times 2}} . S pomocí adjungované matice lze určit její inverzní matici:

B = 1 3 ( x 1 x 2 3 x 2 3 3 x 3 ) = ( x 3 1 3 ( 1 x 2 ) x 2 1 x 3 ) {\displaystyle {\boldsymbol {B}}={\frac {1}{3}}{\begin{pmatrix}x&1-x^{2}\\-3x^{2}-3&3x^{3}\end{pmatrix}}={\begin{pmatrix}{\frac {x}{3}}&{\frac {1}{3}}(1-x^{2})\\-x^{2}-1&x^{3}\end{pmatrix}}} .

Matice

A = ( [ 3 ] 12 [ 7 ] 12 [ 1 ] 12 [ 4 ] 12 ) {\displaystyle {\boldsymbol {A}}={\begin{pmatrix}[3]_{12}&[7]_{12}\\{}[1]_{12}&[4]_{12}\end{pmatrix}}}

s prvky z okruhu zbytkových tříd R = Z / 12 Z {\displaystyle R=\mathbb {Z} /12\mathbb {Z} } má determinant:

det A = [ 3 ] 12 [ 4 ] 12 [ 7 ] 12 [ 1 ] 12 = [ 5 ] 12 {\displaystyle \det {\boldsymbol {A}}=[3]_{12}\cdot [4]_{12}-[7]_{12}\cdot [1]_{12}=[5]_{12}}

Protože prvek [ 5 ] 12 {\displaystyle [5]_{12}} je invertovatelný v R {\displaystyle R} , je matice A {\displaystyle {\boldsymbol {A}}} regulární v R 2 × 2 {\displaystyle R^{2\times 2}} . Její inverzní matice je

B = 1 [ 5 ] 12 ( [ 4 ] 12 [ 7 ] 12 [ 1 ] 12 [ 3 ] 12 ) = [ 5 ] 12 ( [ 4 ] 12 [ 5 ] 12 [ 11 ] 12 [ 3 ] 12 ) = ( [ 8 ] 12 [ 1 ] 12 [ 7 ] 12 [ 3 ] 12 ) {\displaystyle {\boldsymbol {B}}={\frac {1}{[5]_{12}}}{\begin{pmatrix}[4]_{12}&[-7]_{12}\\{}[-1]_{12}&[3]_{12}\end{pmatrix}}=[5]_{12}{\begin{pmatrix}[4]_{12}&[5]_{12}\\{}[11]_{12}&[3]_{12}\end{pmatrix}}={\begin{pmatrix}[8]_{12}&[1]_{12}\\{}[7]_{12}&[3]_{12}\end{pmatrix}}} .

Naopak matice

A = ( [ 3 ] 12 [ 7 ] 12 [ 1 ] 12 [ 9 ] 12 ) {\displaystyle {\boldsymbol {A}}={\begin{pmatrix}[3]_{12}&[7]_{12}\\{}[1]_{12}&[9]_{12}\end{pmatrix}}}

s prvky z téhož okruhu zbytkových tříd R = Z / 12 Z {\displaystyle R=\mathbb {Z} /12\mathbb {Z} } má determinant:

det A = [ 3 ] 12 [ 9 ] 12 [ 7 ] 12 [ 1 ] 12 = [ 20 ] 12 = [ 8 ] 12 {\displaystyle \det {\boldsymbol {A}}=[3]_{12}\cdot [9]_{12}-[7]_{12}\cdot [1]_{12}=[20]_{12}=[8]_{12}}

Čísla 8 {\displaystyle 8} a 12 {\displaystyle 12} jsou soudělná, a proto det A = [ 8 ] 12 {\displaystyle \det {\boldsymbol {A}}=[8]_{12}} v Z / 12 Z {\displaystyle \mathbb {Z} /12\mathbb {Z} } nemá inverzní prvek, a proto je matice A {\displaystyle {\boldsymbol {A}}} singulární.

Alternativní definice

Regulární matice nad tělesem

Pro čtvercovou matici A {\displaystyle {\boldsymbol {A}}} řádu n {\displaystyle n} nad tělesem K {\displaystyle K} , například nad reálnými nebo komplexními čísly, jsou následující podmínky ekvivalentní

(tj. buď jsou všechna pravdivá, nebo všechna nepravdivá pro danou matici):

  • Matice A {\displaystyle {\boldsymbol {A}}} je regulární, čili existuje matice B {\displaystyle {\boldsymbol {B}}} taková, že A B = I {\displaystyle {\boldsymbol {A}}{\boldsymbol {B}}=\mathbf {I} } .
  • Existuje matice B {\displaystyle {\boldsymbol {B}}} taková, že B A = I {\displaystyle {\boldsymbol {B}}{\boldsymbol {A}}=\mathbf {I} } .
  • Determinant matice A {\displaystyle {\boldsymbol {A}}} je nenulový: det A 0 {\displaystyle \det {\boldsymbol {A}}\neq 0} .
  • Nula není vlastní číslo matice A {\displaystyle {\boldsymbol {A}}} .
  • Soustava lineárních rovnic A x = 0 {\displaystyle {\boldsymbol {Ax}}=\mathbf {0} } má pouze triviální řešení x = 0 {\displaystyle {\boldsymbol {x}}=\mathbf {0} } .
  • Pro každé b K n {\displaystyle b\in K^{n}} existuje alespoň jedno řešení x K n {\displaystyle x\in K^{n}} soustavy lineárních rovnic A x = b {\displaystyle {\boldsymbol {Ax}}={\boldsymbol {b}}} .
  • Pro každé b K n {\displaystyle b\in K^{n}} existuje nejvýše jedno řešení x K n {\displaystyle x\in K^{n}} soustavy lineárních rovnic A x = b {\displaystyle {\boldsymbol {Ax}}={\boldsymbol {b}}} .
  • Řádkové vektory jsou lineárně nezávislé.
  • Řádkové vektory generují K n {\displaystyle K^{n}} .
  • Sloupcové vektory jsou lineárně nezávislé.
  • Sloupcové vektory generují K n {\displaystyle K^{n}} .
  • Lineární zobrazení K n K n {\displaystyle K^{n}\to K^{n}} dané předpisem x A x {\displaystyle {\boldsymbol {x}}\to {\boldsymbol {Ax}}} je prosté (injektivní).
  • Lineární zobrazení K n K n {\displaystyle K^{n}\to K^{n}} dané předpisem x A x {\displaystyle {\boldsymbol {x}}\to {\boldsymbol {Ax}}} je surjektivní.
  • Transponovaná matice A T {\displaystyle {\boldsymbol {A}}^{\mathrm {T} }} je regulární.
  • A {\displaystyle {\boldsymbol {A}}} má plnou hodnost, neboli rank A = n {\displaystyle \operatorname {rank} {\boldsymbol {A}}=n} .
  • Matici A {\displaystyle {\boldsymbol {A}}} lze převést ekvivalentními řádkovými úpravami na jednotkovou matici I {\displaystyle \mathbf {I} } .
  • Matici A {\displaystyle {\boldsymbol {A}}} lze převést ekvivalentními řádkovými úpravami do odstupňovaného tvaru s n {\displaystyle n} pivoty.
  • Matici A {\displaystyle {\boldsymbol {A}}} lze vyjádřit jako konečný součin elementárních matic.

Regulární matice nad jednotkovým komutativním okruhem

Obecněji řečeno, pro čtvercovou matici A {\displaystyle {\boldsymbol {A}}} řádu n {\displaystyle n} s prvky z komutativního okruhu s jedničkou R {\displaystyle R} jsou následující podmínky ekvivalentní:

  • Matice A {\displaystyle {\boldsymbol {A}}} je regulární, čili existuje matice B {\displaystyle {\boldsymbol {B}}} taková, že A B = I = B A {\displaystyle {\boldsymbol {A}}{\boldsymbol {B}}=\mathbf {I} ={\boldsymbol {B}}{\boldsymbol {A}}} .
  • Determinant matice A {\displaystyle {\boldsymbol {A}}} má v okruhu R {\displaystyle R} inverzní prvek.
  • Pro všechny b R n {\displaystyle b\in R^{n}} existuje právě jedno řešení x R n {\displaystyle x\in R^{n}} soustavy lineárních rovnic A x = b {\displaystyle {\boldsymbol {Ax}}={\boldsymbol {b}}} .
  • Pro všechny b R n {\displaystyle b\in R^{n}} existuje alespoň jedno řešení x R n {\displaystyle x\in R^{n}} soustavy lineárních rovnic A x = b {\displaystyle {\boldsymbol {Ax}}={\boldsymbol {b}}} .
  • Řádkové vektory tvoří bázi R n {\displaystyle R^{n}} .
  • Řádkové vektory generují R n {\displaystyle R^{n}} .
  • Sloupcové vektory tvoří bázi R n {\displaystyle R^{n}} .
  • Sloupcové vektory generují R n {\displaystyle R^{n}} .
  • Lineární zobrazení R n R n {\displaystyle R^{n}\to R^{n}} dané předpisem x A x {\displaystyle {\boldsymbol {x}}\to {\boldsymbol {Ax}}} , je surjektivní.
  • Transponovaná matice A T {\displaystyle {\boldsymbol {A}}^{\mathrm {T} }} je inverzní.

Podstatný rozdíl oproti tělesu spočívá v tom, že obecně surjektivita (a tedy i bijektivita) nevyplývá z injektivity lineárního zobrazení - například u zobrazení Z Z {\displaystyle \mathbb {Z} \to \mathbb {Z} } , daného předpisem x 2 x {\displaystyle x\to 2x} .

Pro singulární matici A {\displaystyle {\boldsymbol {A}}} není splněna žádná z výše uvedených podmínek.

Vlastnosti

Hustota

Množina singulárních reálných matic řádu n {\displaystyle n} je nulová množina, neboli má Lebesgueovu míru nula. To plyne z toho, že singulární matice jsou kořeny funkce determinantu a ta je spojitá, neboť se jedná o polynom v prvcích matice. V jazyce teorie míry jsou proto téměř všechny matice řádu n {\displaystyle n} regulární.

Kromě toho tvoří regulární matice řádu n {\displaystyle n} hustou otevřenou množinou v topologickém prostoru všech matic řádu n {\displaystyle n} . Množina singulárních matic je naopak uzavřená a řídká.

V praxi se však můžeme setkat se singulárními maticemi. V numerických výpočtech se mohou vyskytnout problematické matice, které jsou sice regulární, ale blízké singulární matici. Takové matrice se nazývají špatně podmíněné.

Počet regulárních matic nad tělesem zbytkových tříd

Matice s prvky z tělesa zbytkových tříd Z p {\displaystyle \mathbb {Z} _{p}} s prvočíslem p {\displaystyle p} je regulární právě tehdy, když jsou řádkové vektory lineárně nezávislé.

Pro těleso Z 2 {\displaystyle \mathbb {Z} _{2}} lze počet regulárních matic řádu n {\displaystyle n} vypočítat takto:

  • Každý z n {\displaystyle n} prvků 1. řádku může nezávisle nabývat dvou hodnot. Nulový vektor je vyloučen. Pro 1. řádek proto existuje 2 n 1 {\displaystyle 2^{n}-1} možností.
  • Pro 2. řádek jsou vyloučeny všechny vektory, které jsou lineární kombinací 1. řádku, takže 2 {\displaystyle 2} vektory. Pro 2. řádek existuje 2 n 2 {\displaystyle 2^{n}-2} možností.
  • Pro 3. řádek jsou vyloučeny všechny vektory, které jsou lineární kombinací 1. a 2. řádku, takže 2 2 {\displaystyle 2^{2}} vektorů. Pro 3. řádek existuje 2 n 2 2 {\displaystyle 2^{n}-2^{2}} možností.
  • Obecně tedy pro řádek s indexem k {\displaystyle k} existuje 2 n 2 k 1 {\displaystyle 2^{n}-2^{k-1}} možných hodnot. Pro všechny řádky matice tedy existuje celkem ( 2 n 2 0 ) ( 2 n 2 1 ) ( 2 n 2 2 ) ( 2 n 2 n 1 ) {\displaystyle (2^{n}-2^{0})(2^{n}-2^{1})(2^{n}-2^{2})\cdots (2^{n}-2^{n-1})} možností.

Z uvedeného lze odvodit i podíl regulárních matic mezi všemi maticemi řádu n {\displaystyle n} . Různých matic řádu n {\displaystyle n} je celkem 2 n n = 2 n 2 {\displaystyle 2^{n\cdot n}=2^{n^{2}}} , protože každý z n n = n 2 {\displaystyle n\cdot n=n^{2}} prvků může nezávisle nabývat dvou hodnot. Podíl regulárních matic je

( 2 n 2 0 ) ( 2 n 2 1 ) ( 2 n 2 2 ) ( 2 n 2 n 1 ) / 2 n n = 2 n 2 0 2 n 2 n 2 1 2 n 2 n 2 2 2 n 2 n 2 n 1 2 n = ( 1 1 2 n ) ( 1 1 2 n 1 ) ( 1 1 2 n 2 ) ( 1 1 2 1 ) = k = 1 n ( 1 ( 1 2 ) k ) {\displaystyle {\begin{aligned}&(2^{n}-2^{0})(2^{n}-2^{1})(2^{n}-2^{2})\cdots (2^{n}-2^{n-1})/2^{n\cdot n}\\&={\frac {2^{n}-2^{0}}{2^{n}}}\cdot {\frac {2^{n}-2^{1}}{2^{n}}}\cdot {\frac {2^{n}-2^{2}}{2^{n}}}\cdots {\frac {2^{n}-2^{n-1}}{2^{n}}}\\&=\left(1-{\frac {1}{2^{n}}}\right)\left(1-{\frac {1}{2^{n-1}}}\right)\left(1-{\frac {1}{2^{n-2}}}\right)\cdots \left(1-{\frac {1}{2^{1}}}\right)\\&=\prod _{k=1}^{n}\left(1-\left({\frac {1}{2}}\right)^{k}\right)\end{aligned}}}

Pro n {\displaystyle n} jdoucí k nekonečnu tento součin konverguje podle věty o pětiúhelníkových číslech | 1 2 | < 1 {\displaystyle \left|{\tfrac {1}{2}}\right|<1} ke konečné limitě, přibližné hodnoty 0,289.

Uvedený výpočet lze zobecnit pro těleso Z p {\displaystyle \mathbb {Z} _{p}} s libovolným prvočíslem p {\displaystyle p} . Různých matic řádu n {\displaystyle n} je p n n = p n 2 {\displaystyle p^{n\cdot n}=p^{n^{2}}} , z nichž ( p n p 0 ) ( p n p 1 ) ( p n p 2 ) ( p n p n 1 ) {\displaystyle (p^{n}-p^{0})\cdot (p^{n}-p^{1})\cdot (p^{n}-p^{2})\cdot \ldots \cdot (p^{n}-p^{n-1})} je regulárních. Podíl regulárních matic je k = 1 n ( 1 ( 1 p ) k ) {\displaystyle \prod _{k=1}^{n}\left(1-\left({\frac {1}{p}}\right)^{k}\right)} .[2]

Odkazy

Reference

V tomto článku byly použity překlady textů z článků Invertible matrix na anglické Wikipedii a Reguläre Matrix na německé Wikipedii.

  1. a b Slovník školské matematiky. Praha: SPN, 1981. 240 s. 
  2. StackExchange: Number of non singular matrices over a finite field of order 2

Literatura

  • Slovník školské matematiky. Praha: SPN, 1981. 240 s. 
  • BÄRTSCH, Hans-Jochen. Matematické vzorce. Praha: Academia, 2006. 832 s. ISBN 80-200-1448-9. Kapitola Matice, s. 180–198. 
  • BEČVÁŘ, Jindřich. Lineární algebra. 1.. vyd. Praha: Matfyzpress, 2019. 436 s. ISBN 978-80-7378-392-1. 
  • HLADÍK, Milan. Lineární algebra (nejen) pro informatiky. 1.. vyd. Praha: Matfyzpress, 2019. 328 s. ISBN 978-80-7378-378-5. S. 39. 
  • OLŠÁK, Petr. Lineární algebra [online]. Praha: 2007 [cit. 2023-02-20]. Dostupné online. 
  • MOTL, Luboš; ZAHRADNÍK, Miloš. Pěstujeme lineární algebru [online]. [cit. 2023-02-20]. Dostupné online. 

Související články