MB141 -3. přednáška Inverzní matice a determinanty Martin Čadek s využitím přednášky Ondřeje Klímy pro předmět MB101 Jarní semestr 2021 3. přednáška Osnova přednášky • Inverzní matice a jejich výpočet • Determinant matice • Výpočet determinantu pomoci Laplaceova rozvoje • Motivace - objem. • Obecná definice \A\ (pro čtvercovou matici A). a Determinant a elementární řádkové úpravy. a Souvislost determinantu s matici >A~1. • Použití pro přímý výpočet řešení soustavy. • Výpočet determinantu. • Determinant a součin matic - A • B = A • B 3. přednáška Inverzní matice Připomeňme, že písmenem E označujeme jednotkovou matici. Definice Říkáme, že B je matice inverzní ke čtvercové matici A, když A • B = B • A = E. Taková matice je určena jednoznačně, a proto píšeme B = /4~1, přičemž B je čtvercová matice stejného rozměru jako A. Matici, k níž existuje matice inverznr říkáme invertibilní matice. Pokud řešíme soustavu A • x = b s invertibilní maticí A, pak x = >A~1 • b je jediné řešení soustavy. A b / /f *# Postup výpočtu inverzní matice: snažme se určit matici X ^ splňující A • X = E postupně po sloupcích. Je-li /4 matice 3x3 a sloupce matice X jsou postupně x, y a z, řešíme rovnice >4x=l0j, 4y=Mj, Az = í 0 ^ * 3. přednáška * - e i) Výpočet inverzní matice Tyto tři soustavy mají stejnou j^tavoti stranu a my je můžeme řešit současně tak, že elementárními řádkovými operacemi upravujeme na schodovitý tvar matici A 1 0 0 0 1 0 | ={A\E) 0 0 1 (C|D) kde matice C je ve schodovitém tvaru. Mohou nastat tyto dvě možnosti: O Jestliže je její poslední řádek nulový, pak jedna ze tří rovnic není řešitelná a inverzní matice neexistuje. O Matice C má v každém řádku pivota. Ty leží na úhlopříčce. Tedy s maticí (C|D) můžeme provádět tzv. zpětnou • © ô Gaussovu eliminaci, tj. elementárními řádkovými 6 • ô operacemi postupně vytvářet nuly nad pivoty matice C P 0 3. přednáška Výpočet - pokračování Tímto postupem dostaneme na místě matice C jednotkovou matici, na místě matice D budou sloupce inverzní matice, tedy /A~1. (A\E)----- (C|D)----- (E|/T1) Ukažme si to na příkladu: 2 3 5 1 0 0 1 2 -1 0 1 0 0 1 2 0 0 1 ( 1 0 0 0 1 o V o o 1 5 9 2 9 9 1 9 4 9 2 9 1 2 -1 0 1 0 0 -1 7 1 -2 0 0 0 9 1 -2 1 5 13 \ 9 9 9 \ 7 9 17 3. přednáška -1 4 0 0 O -i O O O A. % 0 4 *3 ^6 0 0 1 -2/6 % 6 Inverzní matice - algoritmus Algoritmus pro nalezení inverzní matice O Vedie sebe napíšeme původní matici A a jednotkovou matici E. O Matici A upravujeme řádkovými elementárními úpravami nejprve na schodovitý tvar. O Následně zpětnou eliminací na jednotkovou matici E. O Tytéž úpravy souběžně prováděné s vedle napsanou maticí E vedou k hledané inverzní matici /A-1. 0 Pokud tento algoritmus narazí na vynulování celého řádku v původní matici, znamená to, že matice inverzní neexistuje. -II I V Ul <-l II—I I IUIIĽU j 3. přednáška Determinant čtvercové matice Čtvercové matici A = a-n ai2 ^21 ^22 jsme přiřadili číslo \A\ = a-na22 - a12a21, které jsme nazvali determinantem matice A. Jeho geometrický význam byl orientovaný obsah rovnoběžníku určeného vektory (au, a2i) a (a-i2, a22). Determinant budeme definovat pro každou čtvercovou matici. Neuděláme to ale přímým předpisem (i když to je možné), ale nepřímo tak, že vyčíslíme jeho vlastnosti. Výhodou tohoto postupu je, že nám dává přímý návod k výpočtu determinantu, zatímco z přímé definice determinant většinou nepočítáme Geometrický význam determinantu matice A tvaru n x n bude orientovaný objem rovnoběžnostěnu v n-rozměrném prostoru určeného vektory sloupců (nebo řádků) matice A. ►► 3. přednáška 4» i /V . Mť Osul"*- ^ /tfrl***"** Determinant - základní pravidla Definice Každé čtvercové matici A tvaru n x n lze jednoznačně přiřadit číslo \A\, determinant matice A, který splňuje následující pravidla: 1 4 5 Vznikne-li matice B přehozením dvou řádků matice A, pak «*{ 1 í) Vznikne-li matice B vynásobením některého řádku matice 4 číslem c. Dak lei = c - \A\. ólU> f ca<* ce*2 \ -cúa^a^ A číslem c, pak \ B\ = c • \A\. í Vznikne-li matice B z matice A přičtením násobk některého řádku k jinému řádku, pak \B\ = \A\. Determinant jednotkové matice je \E\ = 1. Determinant transponované matice je \AT i i 3. přednáška Odvozená pravidla Z předchozích pravidel lze odvodit další: 6) Determinant čtvercové matice ve schodovitém tvaru (tzv. horní trojúhelníkové matice) je roven součinu čísel na úhlopříčce matice. ►► 7) Determinant matice, která obsahuje nulový řádek nebo sloupec, je roven 0. ►► 8) Je-li A matice tvaru k x /c, B matice tvaru (n - k) x (n - k), C matice tvaru k x (n - k) a O nulová matice tvaru (n - k) x /c, pak determinant matice n x n ►► 0^ A C *cfje a • B O ť \d\ =*s Výpočet determinantu provádíme pomocí Gaussovy eliminace. 3. přednáška 7 0 " O "r Determinant - příklad výpočtu Příklad Určete determinant matice A = 2 3 5 1 2 -1 0 1 2 Výsledek: = 9. Příklad Určete determinant matice B = /1 0 0 3 0 5 \7 0 0 2 \ 4 0 6 0 0 8/ Výsledek: |fí| = 12. 3. přednáška dum * -2 o 4 Z)- t del - 448 ~S s 4m t? t? í c? 4*tur s f Cauchyova věta Věta (Cauchyova) Pro libovolné dve čtvercové matice A a B stejné velikosti platí A-B\ = \A\.\B\. IA+B\= MMI Důsledek: Determinant invertibilní matice je nenulový a platí A~11 = |/A|~1. Platí i obrácené tvrzení. A-A"1 = £T Věta Pro čtvercovou matici A je ekvivalentní: • \A\ ^0, • existuje A~\ • Pro každou pravou stranu b má soustava A- x = b jediné řešení. 3. přednáška Laplaceův rozvoj determinantu Buď A = (a,y) čtvercová matice řádu n > 1. Pro zvolené indexy i, j označme Ay čtvercovou matici velikosti n - 1, která vznikne z A vynecháním Mého řádku a y-tého sloupce. Pak číslo AyJ= (-1)'+y-|>4 nazý^pme algebraický doplněk prvku agy matici A. J Ac'j Věta (Laplaceův rozvoj) Buď A = (aij) čtvercová matice řádu n > 1. Pa/c pro libovolný index i platí A = ayi Ai + + • • • + a/V?A„ n 7=1 Hovoríme o rozvoji determinantu podle Mého řádku. 3. přednáška Laplaceův rozvoj - příklad Příklad Určete determinant matice B = yTlo 0 B| = 1 ■(-!) 1+1 3 4 I ' 0 5 6 l lo + 2 0 0 8 1+4 = 8- 3 4 5 6 Všiměme si, že -2-7 3 4 5 6 = (8-14) 3 4 5 6 = (-6) • (-2) = B CM 3 4 7 8 • 5 6 ty 3. přednáška Determinant matice 3x3 Proveďme Laplaceův rozvoj obecné matice 3x3 A = Dostaneme A 4+2 An - W Git **S A\ = a-11 (^22^33 — ^23^32) — ^12(^21 ^33 ~ ^23^31J + 5l3(a21 ^32 — ^22^31 ) = 311^22333 + a12^23a31 + a13^21^32 — ^13^22^31 — 512^21^33 — 323^32- Pozor, pro větší rozměr nelze počítat takto „úhlopříčně". Ani případy n = 2 a n = 3 nejsou aplikací stejného „úhlopříčného principu". 3. přednáška 4 G 2 r 2 y AoUyZ'žAA+j ^ flu,*, up r/ ?R0 MAľčtr Cramerovo pravidlo Řešíme rovnici Ax = £>, kde A je čtvercová matice, \A\ ^ 0. • \A\ ^ 0 implikuje jednoznačnost řešení. • \A\ ^ 0 implikuje existenci A -i • Odtud x = A~^b, kde /A-1 = /l "1 • A\ Věta Si/ď/A cfi/erco\/á matice rádu n > 1 taková, že\A\ ^0. Pak soustava Ax = b má jediné řešení x = (x1, x2,..., xn)T, /cc/e A přičemž Aj je matice vzniklá z matice A nahrazením jejího j-tého sloupce sloupcem b. 3. přednáška Cramerovo pravidlo - příklad Příklad (Motivační příklad) Vyřešte soustavu lineárních rovnic 2x + 3y + 5z x + 2y - z y + 2z 0 3 5 4 2 -1 -1 1 2 2 3 5 1 2 -1 0 1 2 =1=1 x2,x3 - sami 3. přednáška Moots'-O - O -tit) fr 6 i. ô *3 _ A A ^ 7 * - (á, U a /Is- 3 H A ~1 Požadavky • Spočítat inverzní matici (dle algoritmu). ^ • Spočítat determinant matic (i matic s parametrem). Ovládat obě metody a umět je kombinovat. 3. přednáška Domácí úloha Příklad (3.1) Určete inverzní matici k matici Příklad (3.2) Určete determinant matice A = n 0 5 \0 0 3 0 7 2 0\ 0 4 6 0 0 8/ 3. přednáška Doplňující domácí úloha Příklad (3.3) Pro libovolnou elementární řádkovou operaci nalezněte matici která ji realizuje pomocí násobení. Tj. pokud A ~ B je jedna úprava, pak existuje matice U taková, že B = U • A. Cl • ÍÄ-b) Příklad (3.4) Nechť a a Ď jsou dvě různá reálná čísla a n je kladné celé číslo. Určete determinant matice n/n, která má všechny prvky pod hlavní diagonálou rovny b a všechny prvky na a nad hlavní diagonálou rovny a. - a, chJL Í4 1 * 4 3. přednáška