Regulární matice
Věnujeme dále pozornost zejména čtvercovým maticím.
Věta. Pro každou čtvercovou matici A = (a,j) řádu n nad tělesem (T, +, ■) jsou následující podmínky ekvivalentní:
(i) Řádky matice A jsou lineárně nezávislé.
(ii) Sloupce matice A jsou lineárně nezávislé.
(iii) Hodnost matice A je rovna n.
(iv) Determinant \A\ je nenulový.
Důkaz. Ekvivalence prvních tří podmínek plyne z definice řádkové a sloupcové hodnosti matice A a z poznatku odvozeného v předchozí kapitole, že jde o stejná čísla, tedy o hodnost matice A. Zbývá dokázat ekvivalenci těchto tří podmínek s poslední podmínkou.
Nechť tedy hodnost matice A je rovna n, takže řádky matice A jsou lineárně nezávislé. Elementárními řádkovými úpravami, tak jak bylo popsáno v prvním odstavci důkazu posledního tvrzení v minulé kapitole, se matice A převede na schodovitý tvar. Přitom se použijí pouze řádkové úpravy typu (ii) a přehazování řádků. Na hodnosti matice ani na lineární nezávislosti jejích řádků se tím nic nezmění. Z příslušných tvrzení v kapitole o determinantech navíc plyne, že hodnota determinantu může těmito úpravami nanejvýš změnit znaménko. Poněvadž řádky matice ve schodovitém tvaru jsou lineárně nezávislé, není mezi nimi vespod matice žádný nulový řádek, což nutně znamená, že všechny hlavní prvky leží na hlavní diagonále a že ji vyplní celou. Poněvadž jde o matici v horním trojúhelníkovém tvaru, je podle jiného tvrzení z kapitoly o determinantech její determinant roven součinu prvků ležících na hlavní diagonále. Jedná se tedy o součin hlavních prvků, a ten je nenulový. Takže také determinant \A\ je nenulový.
Nechť naopak řádky matice A jsou lineárně závislé. Je-li A nulová matice, je její determinant roven nule. V opačném případě je n > 1 a podle prvního tvrzení z předminulé kapitoly je některý řádek matice A lineární kombinací ostatních řádků této matice. To ale znamená, že několika elementárními řádkovými úpravami typu (ii) lze dotyčný řádek učinit nulovým řádkem. Podle příslušného tvrzení z kapitoly o determinantech se tím hodnota determinantu matice nezmění. Podle jiného tvrzení ze zmíněné kapitoly je hodnota determinantu matice s nulovým řádkem rovna nule. Takže i determinant \A\ je roven nule.
Čtvercová matice A = (a^) řádu n nad tělesem (T,+, ■), která splňuje ekvivalentní podmínky z předchozí věty, se nazývá regulární matice. Nesplňuje-li taková čtvercová matice A podmínky z předchozí věty, říkáme, že je to singulární matice.
Buď (T, +, ■) těleso a buď m přirozené číslo.
Pro každý index i G {1,..., m} a pro každý prvek r G T, r ^ 0, označme
Ei(r) — čtvercovou matici řádu m, která se liší od jednotkové matice Em pouze tím, že na hlavní diagonále v i-tém řádku a v i-tém sloupci má prvek r.
Dále pro každé dva indexy i, j G {1,... ,m}, i / j a pro každý prvek r E T označme
— čtvercovou matici řádu m, která se liší od jednotkové matice Em pouze tím, že na pozici mimo hlavní diagonálu v i-tém řádku a v j-tém sloupci má prvek r.
Uvedené matice Ej(r) a Ejj(r) se nazývají elementární matice řádu m. Je jasné, že tyto matice jsou regulární, neboť |£?ť(r)|=r^0 a |^-(r)| = 1.
Buď dále A = (a^) matice typu m j n nad tělesem (T, +, ■). Z definice násobení matic a z popisu elementárních řádkových úprav z minulé kapitoly je evidentní, že provedení elementární řádkové úpravy typu (i) s maticí A pro zvolená i G {1,..., m} a r 6 ľ, r 0 dá tentýž výsledek jako vynásobení matice A zleva elementární maticí Ei{r). Podobně provedení elementární řádkové úpravy typu (ii) s maticí A pro zvolená i, j G {1,... ,m}, i ^ j a r E T dá týž výsledek jako vynásobení matice A zleva elementární maticí Eij(r). Obdobná tvrzení platí i pro elementární sloupcové úpravy, používají-li se elementární matice řádu n a násobí-li se jimi matice A zprava (indexy i, j matice Eij(r) se ale prohodí).
Tvrzení. Pro libovolnou čtvercovou matici A = (a^) řádu n nad tělesem (T, +, ■) jsou následující podmínky ekvivalentní:
(i) Matice A je regulární.
(ii) Matice A je řádkově ekvivalentní jednotkové matici En.
(iii) Matice A je součinem několika elementárních matic řádu n.
Důkaz. Je-li matice A regulární, lze ji podle posledního tvrzení z minulé kapitoly převést elementárními řádkovými úpravami na Gauss-Jordánův tvar. Výsledná matice v tomto tvaru je přitom stále regulární, takže její hlavní prvky, jež jsou rovny 1, vyplní hlavní diagonálu a všechny ostatní prvky mimo diagonálu jsou rovny 0. Tak ovšem vzniká jednotková matice En, čili matice A je řádkově ekvivalentní této jednotkové matici En.
Je-li matice A je řádkově ekvivalentní jednotkové matici En, tedy platí-li A ~ En, pak poněvadž je tato relace symetrická podle poznámek z úvodu minulé kapitoly, platí rovněž En ~ A. Lze tedy rovněž konečnou posloupností elementárních řádkových úprav převést jednotkovou matici En na matici A. Poněvadž každou elementární řádkovou úpravu lze realizovat vynásobením některou elementární maticí řádu n zleva, je matice A součinem konečné posloupnosti takových elementárních matic.
Je-li konečně matice A součinem elementárních matic, pak je tato matice součinem regulárních matic, neboť elementární matice jsou regulární. Viděli jsme totiž, že determinanty elementárních matic jsou nenulové. Podle Cauchyovy věty pak také jakýkoliv součin takových matic má nenulový determinant a je tedy regulární maticí. Takže i matice A je regulární.
Důsledek. Dvě matice A a B typu m/n nad tělesem (T, +, ■) jsou řádkově ekvivalentní, tedy splňují A ~ B, právě tehdy, když existuje regulární čtvercová matice C řádu m nad tělesem (T, +, ■) taková, že B = C ■ A.
Důkaz plyne bezprostředně z ekvivalence první a poslední podmínky v předchozím tvrzení.
V kapitole o maticích jsme viděli, že množina Matnn(T) všech čtvercových matic řádu n nad tělesem (T, +, ■) tvoří vzhledem k operaci ■ násobení matic monoid (Matnn(T), ■). Neutrálním prvkem tohoto monoidu je jednotková matice En. Zajímají nás nyní invertibilní prvky tohoto monoidu, tj. invertibilní čtvercové matice řádu n nad tělesem (T, +, ■). Čtvercová matice A řádu n je invertibilní, existuje-li k ní čtvercová matice B řádu n taková, že platí A • B = En = B • A. Pak matice B se nazývá inverzní matice k matici A. Z kapitoly o monoidech víme, že inverzní matice k matici A, pokud existuje, je jediná, a značíme ji A-1.
Věta. Čtvercová matice A řádu n nad tělesem (T, +, ■) je invertibilní právě tehdy, když je regulární.
Důkaz. Je-li daná čtvercová matice A řádu n invertibilní, existuje čtvercová matice B řádu n taková, že A • B = En. Podle Caychyovy věty odtud plyne, že \A\ • \B\ = \A-B\ = = 1, takže \A\ ^ 0 a matice A je tedy regulární.
Je-li naopak čtvercová matice A řádu n regulární, je podle předchozího tvrzení řádkově ekvivalentní jednotkové matici En, takže podle posledního důsledku existuje regulární čtvercová matice C řádu n taková, že C • A = En. Transponovaná matice AT je pak ovšem podle definice také regulární, takže ze stejných důvodů existuje regulární čtvercová matice D řádu n taková, že D • AT = En. Odtud transponováním plyne, že A • DT = En. Přitom C = C • En = C • A - DT = En • DT = DT, takže celkem C = DT je inverzní matice k matici A.
Připomeňme z kapitoly o monoidech a grupách, že pro libovolné dvě regulární čtvercové matice A a B řádu n platí
(A~1)-1 = A      a      (A-B)-1 = B~1-A-1.
Je také užitečné poznamenat, že jakmile o čtvercových maticích A a B řádu n víme, že platí jedna z rovností
A ■ B = En      nebo      B ■ A = En,
pak odtud vyplyne, že tyto rovnosti jsou splněny obě a že tudíž obě matice A, B jsou invertibilní a přitom jedna ke druhé navzájem inverzní. Skutečně platí-li například rovnost A - B = En, pak stejně jako v předchozím důkaze zjistíme, že obě matice A, B jsou regulární a existují k nim tudíž inverzní matice A-1 a B~l. Dále odtud dostáváme B = En-B = A_1-A-B = A~l-En = A'1, čili B = A~l a platí tudíž i rovnost B • A = En.
Z kapitoly o monoidech a grupách též připomeňme, že množina regulárních čtvercových matic řádu n nad tělesem (T, +, ■), tedy množina všech invertibilních matic v monoidu (Matnn(T), ■) je uzavřená vzhledem k násobení matic a že takto sama tvoří grupu. Tuto množinu všech regulárních čtvercových matic řádu n nad tělesem (T, +, ■) bývá zvykem označovat symbolem GLn(T) a mluví se v této souvislosti o obecné lineární grupě stupně n nad (T,+,-).
Z Cauchyovy věty dále například plyne, že množina všech čtvercových matic řádu n nad tělesem (T,+,■), jejichž determinant je roven 1, tvoří podgrupu v GLn(T). Tuto podgrupu bývá obvyklé označovat symbolem SLn(T) a mluví se o ní jako o speciální lineární grupě stupně n nad (T, +, ■).
Závěrem kapitoly se věnujeme metodám výpočtu inverzních matic. Mějme tedy regulární čtvercovou matici A = (a,j) řádu n nad tělesem (T, +, ■). Viděli jsme, že pak A ~ En, což znamená, že provedením několika elementárních řádkových úprav lze z matice A dostat jednotkovou matici En. Nechť F\, F2,..., Fk jsou elementární matice řádu n odpovídající postupně těmto řádkovým úpravám. Pak provedení zmíněné posloupnosti řádkových úprav s maticí A znamená vynásobení matice A postupně maticemi F\, F2, ■ ■ ., Fk zleva, čímž vzniká součin matic
Fk.....F2 • F\ • A = En.
Vezměme tedy matici C = Fk.....F2 • F\. Pak ovšem máme
C • A = En a podle toho, co bylo řečeno výše, je tedy C inverzní maticí k matici A. Provedeme-li dále tutéž posloupnost elementárních řádkových úprav jako výše také s maticí En, dostaneme tím z týchž důvodů matici C • En = C', tedy samotnou inverzní matici k matici A. Myšlenka, o kterou se nyní opírá metoda výpočtu inverzní matice k dané čtvercové matici A řádu n, spočívá v provádění týchž elementárních řádkových úprav s oběma maticemi A i En současně. Přitom volíme tyto úpravy tak, aby převedly matici A do Gauss-Jordánova tvaru. Je-li matice A regulární, vznikne tak jednotková matice En, takže takto po několika krocích obdržíme řádkově ekvivalentní matice
(A\En)-----(#n|C),
kde C = A-1. Není-li matice A regulární, vyjde nám v jejím Gauss-Jordanově tvaru alespoň jeden nulový řádek.
Jinou možnost, jak počítat inverzní matici, poskytuje následující výsledek. Mějme regulární čtvercovou matici A = (a,j) řádu n > 1 nad tělesem (T,+, ■). Připomeňme, že v kapitole o determinantech jsme pro každá i, j G {1, 2,... ,n} definovali algebraický doplněk prvku v matici A. Vezměme čtver-covou matici A = (Aý) řádu n složenou z těchto algebraických
doplňků a k ní transponovanou matici A* = A . Tato matice A* se nazývá adjungovaná matice k matici A. Nyní můžeme formulovat následující fakt:
Věta. Buď A = (a,j) regulární čtvercová matice řádu n > 1 nad tělesem (T, +, ■). Pak platí
A-1 = \A\~l-A\
Důkaz. Podle toho, co bylo uvedeno shora, stačí ukázat, že součin matice A s maticí |^4|_1-A* je roven jednotkové matici En. Neboli stačí ukázat, že A • A* = \A\ • En. Zvolme proto nejprve libovolný index i G {1,2,...,n} a vypočtěme prvek ležící v matici A • A* na diagonále v i-tém řádku a v i-tém sloupci. Tento
prvek je ovšem roven au- An + a^- Ai2-\-----\-ain-Ain = \A\ podle
věty o Laplaceově rozvoji determinantu \A\, a to je odpovídající diagonální prvek matice \A\ • En. Zvolme dále libovolné indexy h 3 £ {Ij 2, • • •, n}, i / ja vypočtěme prvek ležící v matici A- A* na pozici mimo diagonálu v i-tém řádku a v j-tém sloupci. Tento prvek vychází roven a,i ■ Aj\ + • Aj2 + ■ ■ ■ + a«n ■ Ajn = 0, neboť jde o Laplaceův rozvoj determinantu z matice, která se od matice A liší v tom, že se v ní místo j-tého řádku opakuje i-tý řádek. Podle jednoho z tvrzení v kapitole o determinantech je ovšem determinant z matice mající dva stejné řádky roven nule. Prvky matice \A\ • En ležící mimo diagonálu jsou ovšem nulové. To potvrzuje rovnost A-A* = \A\-Enj kterou bylo třeba dokázat.