Matematika pro radiologické asistenty 5. Soustavy lineárních rovnic, matice 5.1 Triviální příklady V jedné nejmenované nemocnici byli zvyklí na dodávku ampulí s lékem, který se přidával do infuzí. Ampule měly vždy objem V a koncentrace účinné látky v ní byla p% objemových. Personál měl příkaz vrchní sestry dávat do infuze o výsledném objemu W vždy jednu ampuli léku. Jednou dodala lék jiná firma a ampule měly objem dvojnásobný, tj. 2V, koncentrace účinné látky byla také dvojnásobná. Vrchní sestra přikázala dávat do infuzí polovinu obsahu ampule. Co myslíte, je to správný příkaz? Řešení: Zavedeme označení: objemy budou v cm3 (ml), tedy W objem infuze (v obou případech stejný) a V objem ampule první firmy. p bude koncentrace účinné látky v objemových procentech v prvním případě. Objem účinné látky a koncentrace jsou 1 1 1 2 2 2 1 , 100 100 2 2 2 2 100 2 100 p U p V U V q W W p V U p V U q q W W = = = = = = = Závěr: Snad nebyla dvojnásobná dávka smrtelná. Obecnější úloha je taková: Předpokládejme, že druhá firma dodala ampule o objemu =20 ml s koncentrací účinné látky p=50 % (objemových). Do jakého objemu základu infuze mají sestry vmíchat jednu ampuli, aby dosáhly předepsané koncentrace q=5 % ? Označme x objem infuze (neznámá místo objemu W z předchozí úlohy). Objem účinné látky v ampuli je 100 p = Rovnice pro neznámy objem x (lineární rovnice pro neznámou veličinu x) je pak ( )100 0q q x q p x = + - = + Řešení: 50 1 ml 20 1 ml 180 ml 5 p x q = - = - = Nepatrně složitější úloha je tato: Do infuze o celkovém objemu W=200 ml se přidávají dvě účinné látky. První z nich je v ampulích o objemu V1=20 ml v koncentraci p1=30 % (objemových), druhá v ampulích o objemu V2 =40 ml v koncentraci p2=50 %. Výsledná koncentrace obou účinných látek v infuzi má být q=15 % a poměr jejich koncentrací 1 2 0,5q q p= = (jedna ku dvěma). Kolik ml roztoku 1 a kolik ml roztoku 2 je třeba dát do infuze ? Řešení ­ záznam úlohy: Označme x hledaný objem roztoku 1 a y hledaný objem roztoku 2 (máme tedy dvě neznámé, x a y). Objem účinné látky 1 v objemu x bude ( )1 1 100U p x= , objem účinné látky 2 v objemu y bude ( )2 2 100U p y= . Koncentrace látek v infuzi pak budou 1 1 2 2 1 2, 100 100 U p x U p y q q W W W W = = = = Pro naše dvě neznámé máme dvě podmínky 1 1 2 2 , q q q q p q + = = Přepíšeme tyto podmínky pomocí neznámých veličin x a y a známých hodnot p1, p2, W, q a p na soustavu dvou lineárních rovnic o dvou neznámých 1 2 1 2 0 p x p y qW p x p p y + = - = Řešení soustavy rovnic: Odečteme druhou rovnici od první a vyřešíme vzhledem k y, potom dosadíme toto řešení do druhé rovnice (samozřejmě je možné dosadit i do první rovnice, řešení pro x musí být stejné). Dostaneme tak ( ) ( )1 2 , 33, 40 1 1 pqW qW x y x y p p p p = = = + + ! 5.2 Matice Tabulku tvořenou m řádky a n sloupci čísel nazýváme maticí dimense m n× ( ) 11 12 1 21 22 2 1 2 n n ik m m mn a a a a a a a a a a = A " " " " " " " Matice stejné dimense m n× můžeme sečítat a násobit číslem ( ) ( ) ( ) ( ) ( ) ( ) ( ) ik ik ik ik ik ik ik ik c a b a b c a a = + = + = + = = = C A B C A Matici A dimense m n× můžeme zprava vynásobit maticí B dimense n s× a získat tak matici C dimense m s× ( ) 1 n ik ip pk p c a b = = = C A B nebo můžeme matici A dimense m n× vynásobit zleva maticí B dimense s m× a získat tak matici C dimense s n× ( ) 1 m ik ip pk p c b a = = = C B A Vidíme, že pro sčítání zůstává komutativita (nezávislost na pořadí) sčítanců zachována i u matic, u násobení to pro součinitele obecně neplatí. Především: násobit můžeme jen matice, které mají stejný počet řádků nebo sloupců. Ale i pro čtvercové matice (dimense n n× ) je komutativita spíše výjimkou. Obrazně vyjádřeno, prvky matice součinu vytváříme takto: v prvním řádku jsou postupně ,,první řádek levé matice krát první sloupec pravé matice", ,,první řádek levé matice krát druhý sloupec pravé matice" až ,,první řádek levé matice krát poslední sloupec pravé matice", v druhém řádku jsou postupně ,,druhý řádek levé matice krát první sloupec pravé matice", ,,druhý řádek levé matice krát druhý sloupec pravé matice" až ,,druhý řádek levé matice krát poslední sloupec pravé matice" atd. až po poslední řádek matice součinu, kde jsou postupně ,,poslední řádek levé matice krát první sloupec pravé matice", ,,poslední řádek levé matice krát druhý sloupec pravé matice" až ,,poslední řádek levé matice krát poslední sloupec pravé matice". Součin ,,řádek krát sloupec" pak znamená, že sečteme součin prvního prvku řádku s prvním prvkem sloupce se součinem druhého prvku řádku s druhým prvkem sloupce atd. až po součin posledního prvku řádku s posledním prvkem sloupce. Příklad pro čtvercové matice dimense 2: 1 0 0 1 , 0 1 1 0 = = - A B Počítáme 1 0 0 1 1 0 0 1 1 1 0 0 0 1 0 1 1 0 1 0 1 1 0 1 1 0 1 0 0 1 1 0 0 1 1 0 0 0 1 1 0 1 1 0 0 1 1 1 0 1 1 0 0 1 1 0 1 0 1 0 1 1 0 0 1 0 0 1 0 1 0 1 0 1 1 + + = = = - + - + - - + + - - = = = - + - + - + + - = = - - + - A B B A A A 1 0 0 0 0 1 1 0 1 0 1 0 1 0 0 1 1 0 1 1 0 1 0 1 0 1 0 1 0 0 1 1 1 0 0 0 1 = +- - + + = = = + + B B Při počítání se nám objevila (ne náhodně, příklad je tak vybrán) jednotková matice ( ) 1 0 0 0 1 0 0 0 0 1 ik = = E " " " " " " tj. čtvercová matice, která má na diagonále jedničky a ostatní prvky jsou rovny nule. Symbol ik je Kroneckerovo delta, pro které 1 0 ik i k i k = = Násobení jednotkovou maticí ponechává původní matici nezměněnou. Vezměme jednotkovou matici dimense n n× , matici A dimense m n× a matici B dimense n s× . Potom je 1 1 , n n ik ip pk ik ik ip pk ik p p c a a c b b = = = = = = = = C A E C E B Obecně nemusí být všechny řádky matice lineárně nezávislé (tj. pro nějaký řádek je v takovém případě možné najít lineární kombinaci zbývajících řádků, že je rovna tomuto řádku). Totéž platí, uvažujeme-li místo řádků o sloupcích. Hodnost matice je definována jako maximální počet lineárně nezávislých řádků. Hodnost matice hraje podstatnou roli při úvahách o počtech řešení soustavy lineárních rovnic. 5.3 Matice a řešení soustavy rovnic Uvažujme soustavu m lineárních rovnic o n neznámých 1 211 12 1 1 21 21 2 1 2 2 2 2 2 1 n n m m mn n n n m x x x x x x a a a b a a a b a xa a bx x + + + = + + + = = + + + = ... ... """"""""""" " ... kde ija , ib (1 ,1i m j n ) jsou koeficienty soustavy rovnic (tj. známá čísla) a jx 1 j n jsou neznámé. Všechny n-tice { }1 2, , , nx x x... , které vyhovují rovnicím, tvoří řešení soustavy rovnic. Se soustavou rovnic jsou spojeny matice soustavy a rozšířená matice soustavy 11 12 1 1 21 22 2 2 1 2 n n m m mn m a a a b a a a b a a a b ... ... ... ... ... ... ... ... Ekvivalentní úpravy soustavy rovnic jsou ty úpravy, které nemění jejich řešení. Je to zajisté násobení libovolné rovnice nenulovým číslem a také přičtení násobku libovolné rovnice k jiné libovolné rovnici. U matic soustavy jsou to tytéž úpravy prováděná s řádky matic. Máme proto tyto ekvivalentní úpravy matice soustavy: (1) násobení všech prvků v libovolném řádku nenulovým číslem (2) přičtení k - násobku prvků v libovolném řádku k jinému libovolnému řádku Hodnost matice jsme definovali jako počet lineárně nezávislých řádků. Převedeme-li ekvivalentními úpravami matici na schodovitý tvar (v daném řádku je vlevo více sloupců s nulami než v řádku nad ním), je pak hodnost přímo dána počtem nenulových řádků. Nejlépe bude ukázat příklady. Musíme jich však zvolit celou řadu, protože máme také velký výběr situací: jaký je počet rovnic a počet neznámých a jaké jsou hodnosti matice soustavy a rozšířené matice soustavy (budeme značit ( )h A a ( )h B ) . Příklad 1 ( 3m n= = ): 2 3 2 2 1 2 x y z x y z x y z + + = - + - = + + = 1 2 3 2 2 1 1 1 1 1 1 2 - - Provedeme úpravy: (1) první rovnici vynásobenou 2 přičteme k druhé a (2) první rovnici vynásobenou (­1) přičteme k třetí (tj. odečteme ji) 2 3 2 5 5 5 2 0 x y z y z y z + + = + = - - = 1 2 3 2 0 5 5 5 0 1 2 0 - - Další úpravy jsou (3) druhou rovnici násobíme 1/5 (tj. dělíme ji 5) a (4) druhou rovnici (po předchozí úpravě) přičteme k třetí rovnici 2 3 2 1 1 x y z y z z + + = + = - = 1 2 3 2 0 1 1 1 0 0 1 1 - Dostali jsme tak ekvivalentní soustavu rovnic, která má stejné řešení jako původní. Řešení najdeme dosazováním ,,odzadu". Soustava má jediné řešení x=1, y=2, z=­1. Matice i rozšířená matice jsou ve stejném schodovitém tvaru a mají tři nenulové řádky ­ hodnost obou matic ( ( ) ( ) 3h A h B= = ) je rovna počtu neznámých, dostáváme jediné řešení. Příklad 2 ( 3m n= = ): 2 3 2 2 1 3 4 4 x y z x y z x y z + + = - + - + = - + = 1 2 3 2 2 1 1 1 3 4 1 4 - - - Provedeme úpravy: (1) první rovnici vynásobenou 2 přičteme k druhé a (2) první rovnici vynásobenou 3 přičteme k třetí 2 3 2 5 5 5 10 10 10 x y z y z y z + + = + = + = 1 2 3 2 0 5 5 5 0 10 10 10 Další úpravy jsou (3) druhou rovnici násobíme 1/5 (tj. dělíme ji 5), (4) třetí rovnici násobíme 1/10 (tj. dělíme ji 10) a nakonec (5) druhou rovnici vynásobenou (­1) přičteme k třetí (tj. odečteme ji) 2 3 2 1 0 0 x y z y z + + = + = = 1 2 3 2 0 1 1 1 0 0 0 0 Ekvivalentní soustava rovnic má stejné řešení jako původní. Najdeme je snadno dosazováním ,,odzadu". V tomto případě zůstává jedna volná neznámá, existuje tedy nekonečně mnoho řešení x=­z, y=1 ­ z, z libovolné. Matice i rozšířená matice jsou ve stejném schodovitém tvaru a mají dva nenulové řádky ­ hodnost obou matic ( ( ) ( ) 2h A h B= = ) je menší než počet neznámých, dostáváme nekonečně mnoho řešení. Příklad 3 ( 3m n= = ): 2 3 2 2 1 3 2 8 x y z x y z x y z + + = - + - = - + + = 1 2 3 2 2 1 1 1 1 3 2 8 - - - Provedeme úpravy: (1) první rovnici vynásobenou 2 přičteme k druhé a (2) první rovnici přičteme k třetí 2 3 2 5 5 5 5 5 10 x y z y z y z + + = + = + = 1 2 3 2 0 5 5 5 0 5 510 Další úpravy jsou (3) druhou rovnici vynásobenou (­1) přičteme k třetí (tj. odečteme ji) a (4) druhou rovnici násobíme 1/5 (tj. dělíme ji 5) 2 3 2 1 0 5 x y z y z + + = + = = 1 2 3 2 0 1 1 1 0 0 0 5 I toto je ekvivalentní soustava rovnic. Zjevně nemá řešení, neboť žádnou volbou proměnných nedosáhneme ,,0=5". Matice soustavy i rozšířená matice jsou ve schodovitém tvaru, ale hodnost matice soustavy ( ( ) 2h A = ) je menší než hodnost rozšířené matice ( ( ) 3h B = ). Příklad 4 ( 2 , 3m n= = ): 2 3 2 2 1 x y z x y z + + = - + - = 1 2 3 2 2 1 1 1 - - Provedeme úpravy: (1) první rovnici vynásobenou 2 přičteme k druhé a (2) upravenou druhou rovnici vynásobíme 1/5 (tj. vydělíme 5) 2 3 2 1 x y z y z + + = =+ 1 2 3 2 0 1 1 1 Ze schodovitého tvaru matic vidíme, že hodnosti jsou stejné a menší než počet neznámých ( ( ) ( ) 2h A h B= = ), dostáváme nekonečně mnoho řešení x=­z, y=1 ­ z, z libovolné. Příklad 5 ( 2 , 3m n= = ): 2 3 2 2 4 6 5 x y z x y z + + = - - - = - 1 2 3 2 2 4 6 5 - - - - Stačí provést úpravu (1) první rovnici vynásobenou 2 přičteme k druhé 2 3 2 0 1 x y z+ + = = - 1 2 3 2 0 0 0 1 - Ekvivalentní soustava rovnic zjevně nemá řešení, neboť žádnou volbou proměnných nedosáhneme ,,0=-1". Matice soustavy i rozšířená matice jsou ve schodovitém tvaru, ale hodnost matice soustavy ( ( ) 1h A = ) je menší než hodnost rozšířené matice ( ( ) 2h B = ). Příklad 6 ( 4 , 3m n= = ): 2 3 2 2 1 2 2 3 x y z x y z x y z x y + + = - + - = + + = - + = 1 2 3 2 2 1 1 1 1 1 1 2 1 2 0 3 - - - Provedeme úpravy: (1) první rovnici vynásobenou 2 přičteme k druhé, (2) první rovnici vynásobenou (­1) přičteme k třetí (tj. odečteme ji) a (3) první rovnici přičteme ke třetí 2 3 2 5 5 5 2 0 4 3 5 x y z y z y z y z + + = + = - - = + = 1 2 3 2 0 5 5 5 0 1 2 0 0 4 3 5 - - Další úpravy jsou: (3) druhou rovnici vynásobíme 1/5 (tj. vydělíme 5), (4) upravenou druhou rovnici přičteme k třetí rovnici, (5) upravenou druhou rovnici vynásobenou (­4) přičteme k třetí a konečně (6) upravenou třetí rovnici odečteme od čtvrté 2 3 2 1 1 0 0 x y z y z z + + = + = - = = 1 2 3 2 0 1 1 1 0 0 1 1 0 0 0 0 - Řešení této ekvivalentní soustavy rovnic najdeme dosazováním ,,odzadu". Soustava má jediné řešení x=1, y=2, z=­1. Matice i rozšířená matice jsou ve stejném schodovitém tvaru a mají tři nenulové řádky ­ hodnost obou matic ( ( ) ( ) 3h A h B= = ) je rovna počtu neznámých, dostáváme jediné řešení. Příklad 7 ( 4 , 3m n= = ): 2 3 2 2 1 3 2 3 3 4 1 x y z x y z x y z x y z + + = - + - = - + + = + + = 1 2 3 2 2 1 1 1 1 3 2 3 3 1 4 1 - - - Provedeme úpravy: (1) první rovnici vynásobenou 2 přičteme k druhé, (2) první rovnici přičteme k třetí a (3) první rovnici vynásobenou (­3) přičteme ke čtvrté 2 3 2 5 5 5 5 5 5 5 5 5 x y z y z y z y z + + = + = + = - - = - 1 2 3 2 0 5 5 5 0 5 5 5 0 5 5 5 - - - Při dalších úpravách (4) přičteme druhou rovnici ke čtvrté, (5) odečteme druhou rovnici od třetí a nakonec (6) vynásobíme tuto rovnici 1/5 (dělíme 5) 2 3 2 1 0 0 0 0 x y z y z + + = + = = = 1 2 3 2 0 1 1 1 0 0 0 0 0 0 0 0 Soustava má nekonečně mnoho řešení x = ­ z, y = 1 ­ z, z libovolné. Počet schodů, tj. počet nenulových řádků obou matic je stejný, jejich hodnosti ( ) ( ) 2h A h B= = , což je hodnota o jedničku menší než počet neznámých 3n= . Příklad 7 ( 4 , 3m n= = ): 2 3 2 2 3 2 2 4 6 4 3 6 9 6 x y z x y z x y z x y z + + = - - - = + + = - - - = - 1 2 3 2 1 2 3 2 2 4 6 4 3 6 9 6 - - - - - - - - Provedeme úpravy, v tomto příkladu velmi jednoduché: (1) první rovnici přičteme ke druhé, (2) první rovnici vynásobenou (-2) přičteme ke třetí a (3) první rovnici vynásobenou 3 přičteme ke čtvrté 2 3 2 0 0 0 0 0 0 x y z+ + = = = = 1 2 3 2 0 0 0 0 0 0 0 0 0 0 0 0 Hodnosti matic jsou shodné ( ) ( ) 1h A h B= = a jsou o 2 menší než počet neznámých 3n= . Soustava má nekonečně mnoho řešení x = 2 ­ 2y ­ 3z, y a z libovolné. 5.4 Shrnutí Gaussovy eliminační metody V předchozí části byl na příkladech ukázán způsob řešení soustavy rovnic převodem matice soustavy rovnic a rozšířené matice (tj. k matici soustavy přidáváme sloupec pravých stran) na schodovitý tvar. Tomuto způsobu říkáme Gaussova eliminační (,,likvidační") metoda. Zavedli jsme pojem hodnost matice (počet nenulových řádků jejího schodovitého tvaru) a označení ( )h A pro hodnost matice soustavy a ( )h B pro hodnost rozšířené matice. Soustava m rovnic o n neznámých má řešení právě tehdy, je-li ( ) ( )h A h B h= = , tj. počet nenulových řádků je stejný (schodovité tvary mají stejný počet schodů). Počet volných neznámých je d n h= To znamená, že soustava má jediné řešení (nejsou žádné volné neznámé) pro 0d = , tedy v případě, kdy hodnosti matic jsou stejné a rovnu počtu neznámých.