3. SÚSTAVY LINEÁRNYCH ROVNÍC V tejto kapitole sa predbežne zoznámime so sústavami lineárnych rovníc nad všeobec- ným poľom K a naučíme sa ich riešiť. Využijeme pri tom zápis sústavy pomocou istej matice. Štruktúrne vlastnosti množiny všetkých riešení danej sústavy a ich dôsledky pre- študujeme a využijeme až neskôr, keď sa bližšie oboznámime so štruktúrou vektorových priestorov. 3.1. Maticový zápis sústavy lineárnych rovníc Pod lineárnou rovnicou o n neznámych x1, . . . , xn nad poľom K rozumieme formulu tvaru a1x1 + a2x2 + . . . + anxn = b, kde a1, a2, . . . , an, b K, v premenných x1, x2, . . . , xn. Sústavou m lineárnych rovníc o n neznámych x1, x2, . . . , xn nad poľom K rozumieme konjunkciu formúl tvaru a11x1 + a21x2 + . . .+ a1nxn = b1 a21x1 + a22x2 + . . .+ a2nxn = b2 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . am1x1+ am2x2+ . . .+ amnxn= bm, kde aij, bi, pre 1 i m, 1 j n, sú skaláry z poľa K. Maticu A = (aij) Km×n nazývame maticou sústavy, stĺpcový vektor b = (b1, . . . , bm)T Km nazývame jej pravou stranou. Konečne blokovú maticu (A | b) Km×(n+1) nazývame rozšírenou ma- ticou sústavy. Sústava sa nazýva homogénna, ak b = 0; v opačnom prípade sa nazýva nehomogénna. Uvedenú sústavu možno stručne a úsporne zapísať v maticovom tvare A x = b, prípadne, ak ide o homogénnu sústavu, v tvare A x = 0. Riešením sústavy A x = b nazývame každý vektor-stĺpec x = (x1, x2, . . . , xn)T Kn , ktorého zložky vyhovujú všetkým rovníciam tejto sústavy, t. j. platí preň A x = b. Vyriešiť sústavu znamená nájsť všetky jej riešenia, t. j. popísať množinu všetkých jej riešení. 1 2 PAVOL ZLATOŠ: LINEÁRNA ALGEBRA A GEOMETRIA Dve sústavy A x = b a B x = c, kde A, B Km×n , b, c Km×1 , sa nazývajú ekvivalentné, ak majú rovnakú množinu riešení, t. j. ak pre všetky x Kn platí Ax = b práve vtedy, keď B x = c. Skôr než prikročíme k otázke riešenia sústav lineárnych rovníc, považujeme za po- trebné upozorniť čitateľa na dve veci. (a) Podčiarkujeme, že riešením sústavy rozumieme vždy vektor x a nie jeho zložky. Tak napríklad sústava 2x + 3y = 12 3x - 2y = 5 nad poľom R má, ako ľahko nahliadneme, jediné riešenie x y = 3 2 a nie dve riešenia x = 3, y = 2. Keď si toto poriadne uvedomíme, môžeme (a samozrejme aj budeme) sa naďalej vyjadrovať obvyklým spôsobom. Budeme teda hovoriť, že sústava má jediné riešenie x = 3, y = 2. (b) Všimnite si, že počet rovníc sústavy a počet neznámych sa nemusia rovnať. V obvyklom prípade, keď rovníc je rovnaký počet ako neznámych, očakávame, že sústava bude mať jediné riešenie. Keď je rovníc menej než neznámych, môžeme očakávať, že sú- stava bude mať viacero (prípadne i nekonečne mnoho) riešení. Naopak, keď je rovníc viac ako neznámych, môže sa stať, že sústava nebude mať nijaké riešenie. Napriek tomu, že tieto očakávania vyjadrujú niečo ako " prevládajúci trend", ľahko možno nájsť príklady, keď sa nemusia splniť. Zatiaľ len poznamenajme, že homogénna sústava A x = 0 má (bez ohľadu na počet neznámych a počet rovníc) vždy aspoň jedno riešenie ­ je ním nulový vektor x = 0. Nie je dôležité, akými znakmi sú označené neznáme v sústave A x = b. Na jej riešenie nemá nijaký vplyv, či si vektor neznámych označíme x = (x1, . . . , xn)T alebo y = (y1, . . . , yn)T alebo nejako inak. To znamená, že celá informácia o tejto sústave, potrebná na nájdenie všetkých jej riešení, je obsiahnutá v rozšírenej matici sústavy (A | b), prípadne, ak ide o homogénnu sústavu, len v matici sústavy A. Preto i metóda riešenia sústav lineárnych rovníc, s ktorou sa teraz zoznámime, bude založená len na úprave tejto matice. Stručne povedané, rozšírenú maticu (A | b) sústavy Ax = b budeme upravovať tak, aby sme dostali vhodnú maticu (B | c), zodpovedajúcu novej sústave B x = c, ktorá spĺňa nasledujúce dve podmienky: (a) Je ekvivalentná s pôvodnou sústavou A x = b, t. j. má rovnakú množinu riešení. (b) Všetky jej riešenia možno priamo vyčítať z jej rozšírenej matice (B | c). V takom prípade hovoríme, že sústava B x = c je vyriešená. 3.2. Redukovaný stupňovitý tvar matice Našou prvou úlohou teda bude vyjasniť si, ako by mala vyzerať matica (B | c), aby sme príslušnú sústavu B x = c mohli považovať za vyriešenú. Za tým účelom teraz zavedieme niekoľko pojmov. Hovoríme, že prvok aij matice A Km×n je vedúci prvok i-teho riadku matice A, ak aij = 0, a j = 1 alebo ail = 0 pre všetky 1 l < j. Inak povedané, vedúci prvok 3. SÚSTAVY LINEÁRNYCH ROVNÍC 3 nenulového riadku je prvý nenulový prvok tohto riadku. Nulový riadok nemá vedúci prvok. Hovoríme, že matica A = (aij) Km×n je v redukovanom stupňovitom tvare, ak spĺňa nasledujúce štyri podmienky: (a) Ak ri(A) = 0 a rk(A) = 0, tak i < k; t. j. každý nenulový riadok matice A leží nad každým jej nulovým riadkom. (b) Ak aij, akl sú vedúce prvky i-teho resp. k-teho riadku a i < k, tak aj j < l; t. j. vedúci prvok vyššieho riadku leží viac vľavo než vedúci prvok nižšieho riadku. (c) Ak aij je vedúci prvok i-teho riadku, tak aij = 1; t. j. vedúci prvok každého nenulového riadku je 1. (d) Ak aij je vedúci prvok i-teho riadku, tak akj = 0 pre každé k = i; t. j. v stĺpci, v ktorom sa nachádza vedúci prvok nejakého riadku, sú všetky ostatné prvky rovné 0. Pokiaľ matica A spĺňa len podmienky (a), (b), hovoríme, že je v stupňovitom tvare. Používa sa tiež názov (redukovaný) schodovitý tvar. Napríklad z uvedených matíc nad poľom R 0 2 1 1 0 0 0 0 0 2 3 0 1 0 0 1 4 0 0 0 0 0 1 0 0 0 0 1 0 0 0 0 1 0 0 0 1 0 0 0 1 0 0 0 1 1 2 3 4 0 1 2 3 0 0 1 2 0 0 0 1 1 0 1 0 0 1 2 0 0 0 0 1 0 0 0 0 ani jedna z matíc vľavo nie je v stupňovitom tvar; obe matice v strednom stĺpci sú v stupňovitom tvare, nie však v redukovanom stupňovitom tvare; konečne, obe ma- tice vpravo sú v redukovanom stupňovitom tvare. (V každom jednotlivom prípade si podrobne premyslite prečo.) Taktiež každá jednotková matica In, ako aj všetky nulové matice 0mn sú v reduko- vanom stupňovitom tvare. Uvedomme si teraz, akej sústave lineárnych rovníc zodpovedá rozšírená matica v re- dukovanom stupňovitom tvare. Napríklad (B | c) = 1 0 -2 0 0 0 1 6 0 0 0 0 0 1 0 3 0 1 je matica v redukovanom stupňovitom tvare nad R. Táto matica zodpovedá sústave x1 - 2x3 = 3 x2 + 6x3 = 0 x4 = 1 4 PAVOL ZLATOŠ: LINEÁRNA ALGEBRA A GEOMETRIA v neznámych x1, x2, x3, x4, x5. Hneď vidíme, že táto sústava má nekonečne mnoho rie- šení. Každej voľbe parametrov s, t R totiž zodpovedá jedno riešenie x1 = 3 + 2s x2 = - 6s x3 = s x4 = 1 x5 = t. Asi sa zhodneme na tom, že preznačenie neznámych za parametre x3 = s, x5 = t a ich presun na pravú stranu, je úprava natoľko bezprostredná, že sústavu prislúchajúcu k matici (B | c) už možno považovať za vyriešenú. Na napísanie jej riešenia nemusíme písať príslušnú sústavu, môžeme ho napísať priamo na základe matice (B | c). Dohodneme sa teda, že sústavu lineárnych rovníc B x = c nad poľom K budeme nazývať vyriešenou sústavou, ak jej rozšírená matica (B | c) je v redukovanom stupňovi- tom tvare. V prípade homogénnej sústavy sa, samozrejme, stačí obmedziť na maticu B. Teraz si predvedieme, ako možno k danej blokovej matici (B | c) v redukovanom stupňovitom tvare, nájsť všetky riešenia sústavy B x = c. Najprv si ujasníme, kedy je taká sústava riešiteľná, t. j. má aspoň jedno riešenie. Odpoveď na túto otázku je jednoduchá: Sústava B x = c má riešenie práve vtedy, keď sa v matici (B | c) nenachádza riadok tvaru (0, . . . , 0 n-krát | 1). Taký riadok totiž zodpovedá rovnici 0 = 1, ktorá očividne nemá riešenie. To, že ne- prítomnosť takého riadku je i postačujúcou podmienkou riešiteľnosti sústavy, vyplýva z nasledujúceho postupu, ako toto riešenie nájsť. Ak sa v j-tom stĺpci matice B nenachádza vedúci prvok žiadneho riadku, tak si neznámu xj zvolíme za parameter; ak sa v j-tom stĺpci nachádza vedúci prvok nejakého riadku, tak neznámu xj si vyjadríme pomocou parametrov tak, že stĺpce matice B prislúchajúce týmto parametrom " prehodíme s opačným znamienkom na druhú stranu". Presnejšiu formuláciu celého postupu vo všeobecnej podobe si odpustíme. Názornej- šie bude osvetliť ho na ešte jednom príklade. (B | c) = 1 0 0 2/3 -1/2 0 1 0 3/4 0 0 0 1 -4 -2/5 5 2 -2 je reálna matica v redukovanom stupňovitom tvare. Vidíme, že sa v nej nenachádza riadok tvaru (0, 0, 0, 0 | 1), teda sústava B x = c by mala mať riešenie. Vedúce prvky riadkov matice B sa nachádzajú v stĺpcoch 1, 2 a 3. Za parametre si teda zvolíme 3. SÚSTAVY LINEÁRNYCH ROVNÍC 5 neznáme x4 a x5. Riešením sústavy je každý vektor (x1, x2, x3, x4, x5)T R tvaru x1 = 5 - 2 3 s + 1 2 t x2 = 2 - 3 4 s x3 = -2 + 4s + 2 5 t x4 = s x5 = t, kde parametre s, t R môžu nadobúdať ľubovoľné hodnoty. Pre estétov ešte pozname- najme, že zlomkov pri parametroch sa možno jednoducho zbaviť. Je totiž jedno, či si parametrické premenné zvolíme v tvare x4 = s, x5 = t alebo v tvare x4 = 12s, x5 = 10t, kde s, t R. Pri takejto voľbe parametrov dostaneme všetky riešenia sústavy v tvare bez zlomkov x1 = 5 - 8s + 5t x2 = 2 - 9s x3 = -2 + 36s + 4t x4 = 12s x5 = 10t. 3.3. Elementárne riadkové a stĺpcové operácie (ERO a ESO) Zatiaľ sme si ujasnili, na aký tvar treba upraviť rozšírenú maticu (A | b) sústavy lineár- nych rovníc A x = b, aby sme získali s ňou ekvivalentnú vyriešenú sústavu B x = c: rozšírená matica (B | c) novej sústavy musí byť v redukovanom stupňovitom tvare. Te- raz si ukážeme, ako to možno urobiť. Maticu (A | b) budeme na redukovaný stupňovitý tvar upravovať pomocou tzv. elementárnych riadkových operácií. Elementárnou riadkovou operáciou, skrátene tiež ERO, na matici A Km×n roz- umieme I. výmenu dvoch riadkov matice A; II. vynásobenie niektorého riadku matice A nenulovým skalárom z poľa K; III. pripočítanie skalárneho násobku niektorého riadku matice A k jej inému riadku. Matice A, B Km×n sa nazývajú riadkovo ekvivalentné, označenie A B, ak jednu z nich možno upraviť na druhú konečným počtom elementárnych riadkových operá- cií. Riešenie sústavy lineárnych rovníc úpravou jej rozšírenej matice pomocou ERO na riadkovo ekvivalentnú maticu v redukovanom stupňovitom tvare sa nazýva Gaussova- Jordanova eliminácia. Prenechávame čitateľovi, aby si sám sformuloval analogické pojmy elementárnych stĺpcových operácií (ESO) a stĺpcovej ekvivalencie matíc, označenie A B. Ich význam vyjde najavo až v neskorších kapitolách. 6 PAVOL ZLATOŠ: LINEÁRNA ALGEBRA A GEOMETRIA Výmenou i-teho a k-teho riadku v matici A = r1(A) ... ri(A) ... rk(A) ... rm(A) dostaneme maticu r1(A) ... rk(A) ... ri(A) ... rm(A) . Opätovnou výmenou i-teho a k-teho riadku v tejto matici získame zasa maticu A. Vynásobením i-teho riadku matice A skalárom c = 0 dostaneme maticu r1(A) ... cri(A) ... rk(A) ... rm(A) . Vynásobením i-teho riadku tejto matice skalárom c-1 = 0 získame opäť maticu A. Konečne, pripočítaním c-násobku i-teho riadku matice A k jej k-temu riadku z nej dostaneme maticu r1(A) ... ri(A) ... rk(A) + cri(A) ... rm(A) . Všimnite si, že i-ty riadok pri tom zostáva nezmenený. Maticu A z tejto matice získame pripočítaním (-c)-násobku jej i-teho riadku ku k-temu riadku. Ak Ax = b je sústava s rozšírenou maticou (A | b) a bloková matica (A | b ) vznikne z (A | b) vykonaním jednej (nezáleží ktorej) ERO, tak sústava A x = b je ekvivalentná s pôvodnou sústavou A x = b. Elementárne riadkové operácie na matici (A | b) totiž zodpovedajú postupne zámene poradia dvoch rovníc sústavy, vynásobením niektorej rovnice nenulovým skalárom a pripočítaním nejakého násobku jednej rovnice k inej rovnici (presnejšie nahradením dojice rovníc ri(A) x = bi, rk(A) x = bk dvojicou rovníc ri(A) x = bi, (rk(A) + cri(A)) x = bk + cbi). Z pred chvíľou vykonaných úvah vyplýva, že ide o ekvivalentné úpravy, ktorými sa množina riešení sústavy nezmení ­ od novej sústavy A x = b sa možno vhodnou ERO vykonanou na jej rozšírenej matici opäť vrátiť k pôvodnej sústave A x = b. 3. SÚSTAVY LINEÁRNYCH ROVNÍC 7 3.3.1. Tvrdenie. Nech K je pole, A, B Km×n , b, c Km . Ak sú blokové matice (A | b), (B | c) riadkovo ekvivalentné, tak i sústavy lineárnych rovníc Ax = b, Bx = c sú ekvivalentné. Dôkaz. Podľa predpokladu existuje postupnosť blokových matíc (A | b) = (C0 | d0), (C1 | d1), . . . , (Cp | dp) = (B | c) typu m × (n + 1) takých, že pre každé l < p matica (Cl+1 | dl+1) vznikne vykonaním jedinej ERO z matice (Cl | dl). Potom všetky sústavy Cl x = dl majú tú istú množinu riešení, t. j. sú ekvivalentné. 3.3.2. Veta. Každá matica nad poľom K je riadkovo ekvivalentná s jednoznačne ur- čenou maticou v redukovanom stupňovitom tvare. Dôkaz existencie a jednoznačnosti spomínanej matice odložíme do cvičení 12 a 13 (pozri tiež cvičenie 5.13). Zatiaľ sa len na konkrétnych príkladoch naučíme, ako možno k danej matici A nájsť s ňou riadkovo ekvivalentnú maticu v redukovanom stupňovitom tvare. Takýto prístup má navyše tú výhodu, že z radu možných postupov, medzi kto- rými si možno pružne voliť podľa okolností, nám nesugeruje jedinú stratégiu, na jednu z ktorých by sme sa nevyhnutne museli obmedziť pri všeobecnom dôkaze. Ambicióznejší čitateľ v uvedených príkladoch ľahko i sám zahliadne myšlienku všeobecného dôkazu, ktorú potom bude môcť uplatniť v cvičení 12. Napokon, aby sme sa nebavili iba o mati- ciach, začneme zakaždým s nejakou sústavou lineárnych rovníc. Tým sa zároveň naučíme riešiť ľubovoľnú sústavu lineárnych rovníc nad daným poľom Gaussovou-Jordanovou elimináciou, prípadne rozpoznať, že daná sústava nemá riešenie. 3.3.3. Príklad. Je daná sústava 2x1 + 3x2 - x4 = 1 3x1 + 2x2 + 4x3 - 2x4 = 0 x1 - x2 + 4x3 - x4 = 2 troch rovníc o štyroch neznámych nad poľom R. Jej rozšírená matica je 2 3 0 -1 3 2 4 -2 1 -1 4 -1 1 0 2 . Pri jej úprave na redukovaný stupňovitý tvar (podobne ako i v ďalších príkladoch) budeme vynechávať niektoré medzikroky a zaznamenáme len niektoré výsledky viace- rých vykonaných ERO. Posledný riadok matice dáme na prvé miesto, potom jeho (-2)- násobok pripočítame k pôvodnému prvému riadku, ktorý posunieme na druhé miesto, a (-3)-násobok toho istého riadku pripočítame k pôvodnému druhému riadku, ktorý posunieme na tretie miesto. Dostaneme tak maticu 1 -1 4 -1 0 5 -8 1 0 5 -8 1 2 -3 -6 . 8 PAVOL ZLATOŠ: LINEÁRNA ALGEBRA A GEOMETRIA Pripočítaním (-1)-násobku druhého riadku k tretiemu dostaneme maticu 1 -1 4 -1 0 5 -8 1 0 0 0 0 2 -3 -3 . Už z tohto tvaru vidíme, že sústava zodpovedajúca poslednej matici nemá riešenie ­ obsahuje totiž rovnicu 0 = -3. Teda ani pôvodná sústava (hoci neznámych je v nej viac než rovníc) nemá riešenie. Z cvičných dôvodov však dokončíme úpravu na redukovaný stupňovitý tvar, ktorý dostaneme vynásobením tretieho riadku skalárom -1/3, pripo- čítaním (-2)-násobku resp. 3-násobku tohto nového riadku k prvému resp. druhému riadku a, konečne, vynásobením druhého riadku skalárom 1/5: 1 0 12/5 -4/5 0 1 -8/5 1/5 0 0 0 0 0 0 1 . Čitateľ by si mal všimnúť, že po nastavení vedúceho prvku niektorého riadku na hodnotu 1 okamžite pristupujeme k nulovaniu zvyšných prvkov stĺpca, v ktorom leží tento vedúci prvok. 3.3.4. Príklad. Riešme sústavu x + 2iy = 5 + 4i (3 - i)y + (6 - 2i)z = 10 2x - z = 5 + 3i x + y + z = 5 + 2i štyroch rovníc o troch neznámych nad poľom C. Jej rozšírená matica je 1 2i 0 0 3 - i 6 - 2i 2 0 -1 1 1 1 5 + 4i 10 5 + 3i 5 + 2i . Prehoďme jej posledný riadok na prvé miesto a zvyšné riadky posuňme o jedno miesto nadol. V takto získanej matici pripočítajme (-1)-násobok prvého riadku k druhému riadku a (-2)-násobok prvého riadku k štvrtému riadku. Konečne vynásobme tretí riadok skalárom (3 + i)/10. Dostaneme maticu 1 1 1 0 -1 + 2i -1 0 1 2 0 -2 -3 5 + 2i 2i 3 + i -5 - i . 3. SÚSTAVY LINEÁRNYCH ROVNÍC 9 (-1)-násobok tretieho riadku pripočítame k prvému riadku, jeho (1-2i)-násobok k dru- hému a 2-násobok k štvrtému. Nakoniec výmenou druhého a tretieho riadku dostaneme maticu 1 0 -1 0 1 2 0 0 1 - 4i 0 0 1 2 + i 3 + i 5 - 3i 1 + i . Pripočítajme posledný riadok k prvému, (-2)-násobok posledného riadku k druhému a jeho (-1 + 4i)-násobok k tretiemu. Zostáva vymeniť tretí a štvrtý riadok ­ výsledná matica je už v redukovanom stupňovitom tvare 1 0 0 0 1 0 0 0 1 0 0 0 3 + 2i 1 - i 1 + i 0 . Vidíme, že pôvodná sústava (hoci obsahuje viac rovníc než neznámych) má jediné rie- šenie x = 3 + 2i, y = 1 - i, z = 1 + i, teda presnejšie vektor (3 + 2i, 1 - i, 1 + i)T C3 . 3.3.5. Príklad. Uvažujme sústavu x1 + x2 + 2x3 + 3x4 = 0 2x1 + 4x3 = 0 x1 + 2x2 + x3 + 3x4 = 0 3x3 + 4x4 = 0 štyroch rovníc o štyroch neznámych nad poľom Z5. Keďže ide o homogénnu sústavu (ktorej ľavá strana je nulový stĺpcový vektor, teda sa nemení pri žiadnej ERO), stačí upravovať jej (nerozšírenú) maticu 1 1 2 3 2 0 4 0 1 2 1 3 0 0 3 4 . (-2)-násobok, t. j. 3-násobok prvého riadku pripočítame k druhému riadku a jeho (-1)- násobok, t. j. 4-násobok pripočítame k tretiemu riadku. Dostaneme tak maticu 1 1 2 3 0 3 0 4 0 1 4 0 0 0 3 4 . (-1)-násobok, t. j. 4-násobok tretieho riadku pripočítame k prvému riadku a jeho (-3)- násobok, t. j. 2-násobok pripočítame k druhému riadku. Konečne výmenou druhého a tretieho riadku dostaneme maticu 1 0 3 3 0 1 4 0 0 0 3 4 0 0 3 4 . 10 PAVOL ZLATOŠ: LINEÁRNA ALGEBRA A GEOMETRIA Tretí riadok odpočítame od prvého aj od štvrtého riadku. Ďalej ho vynásobíme skalárom 3-1 = 2. Napokon jeho (-4)-násobok, t. j. priamo tento nový tretí riadok pripočítame k druhému riadku. Výsledná matica je už v redukovanom stupňovitom tvare 1 0 0 4 0 1 0 3 0 0 1 3 0 0 0 0 . Premennú x4 si zvolíme za parameter. Všetky riešenia sústavy majú potom tvar x1 = t, x2 = 2t, x3 = 2t, x4 = t, kde t Z5. Vidíme teda, že pôvodná sústava (hoci počet jej rovníc je rovnaký ako počet neznámych) má viac než jedno riešenie; nie je ich však nekonečne veľa ale len 5. Práve toľko je totiž možných volieb parametra t, t. j. prvkov poľa Z5. Zaznamenajme ešte jeden očakávaný dôsledok tvrdenia 3.3.1, vety 3.3.2 a spôsobu, ako napísať riešenie sústavy s (rozšírenou) maticou v redukovanom stupňovitom tvare, uvedeného v paragrafe 3.2. 3.3.6. Tvrdenie. Nech A Km×n , b Km a m < n, t. j. sústavy Ax = 0, Ax = b obsahujú menej rovníc než neznámych. Potom (a) homogénna sústava Ax = 0 má popri riešení x0 = 0 aspoň jedno riešenie x = 0; (b) ak existuje aspoň jedno riešenie sústavy A x = b, tak táto sústava má viac než jedno riešenie. Dôkaz. (a) Upravme maticu sústavy A na redukovaný stupňovitý tvar B. Uvedomme si, že matice A aj B majú m riadkov a n stĺpcov. Riadky matice B majú nanajvýš m vedúcich prvkov. Keďže m < n, aspoň v jednom stĺpci matice B neleží vedúci prvok žiadneho riadku. Nech je to napr. j-ty stĺpec. Potom voľbe parametra xj = t K, t = 0 zodpovedá aspoň jedno nenulové riešenie sústavy A x = 0. (b) prenechávame ako cvičenie čitateľovi. 3.4. Gaussova eliminačná metóda Hlavne z historických dôvodov ešte stručne spomenieme tzv. Gaussovu eliminačnú me- tódu riešenia sústav lineárnych rovníc. Pri riešení touto metódou upravíme rozšírenú maticu sústavy len na stupňovitý (teda nie nevyhnutne redukovaný stupňovitý) tvar. Už z tohto tvaru možno ľahko spoznať, či sústava má nejaké riešenie (príslušná ma- tica nesmie obsahovať riadok tvaru (0, . . . , 0 | d), kde 0 = d K). V tom prípade možno všetky riešenia sústavy získať voľbou parametrov (opäť si za ne volíme neznáme xj také, že j-tom stĺpci sa nevyskytuje vedúci prvok žiadneho riadku) a spätným dosadzovaním, t. j. elimináciou neznámych pomocou parametrov. 3.4.1. Príklad. Predpokladajme, že rozšírenú maticu nejakej sústavy nad R sme už pomocou ERO upravili na stupňovitý tvar 0 2 3 0 -1 4 0 0 0 -2 5 4 0 0 0 0 3 1 1 0 4 . 3. SÚSTAVY LINEÁRNYCH ROVNÍC 11 Táto matica zodpovedá sústave 2x2 + 3x3 - x5 + 4x6 = 1 - 2x4 + 5x5 + 4x6 = 0 3x5 + x6 = 4. Za parametre si zvolíme premenné x1, x3 a x6. Spätným dosadzovaním postupne do- staneme všetky riešenia v parametrickom tvare x6 = t x5 = 1 3 (4 - x6) = 4 3 - 1 3 t x4 = 1 2 (5x5 + 4x6) = 10 3 - 7 6 t x3 = s x2 = 1 2 (1 - 3x3 + x5 - 4x6) = 7 6 - 3 2 s - 13 6 t x1 = r, kde r, s, t R. Prípadne, po trochu " šikovnejšej" voľbe parametrov, v tvare x6 = 6t, x5 = 4 3 - 2t, x4 = 10 3 - 7t, x3 = 2s, x2 = 7 6 - 3s + 13t, x1 = r. Pozorný čitateľ si iste všimol, že spätné dosadzovanie možno nahradiť ďalšou úpra- vou rozšírenej matice sústavy pomocou ERO na redukovaný stupňovitý tvar. Stačí totiž vynásobiť nenulové riadky prevrátenými hodnotami ich vedúcich prvkov a pripočíta- ním vhodných násobkov týchto riadkov vynulovať zvyšné nenulové prvky v stĺpcoch obsahujúcich vedúce prvky jednotlivých riadkov. I tak však môže byť Gaussova eliminačná metóda v niektorých prípadoch užitočná ­ najmä keď nám nejde ani tak o explicitný tvar riešení, ako skôr o samotnú otázku riešiteľnosti sústavy, prípadne o počet parametrov, ktoré sa v nich vyskytujú. Všetko to možno totiž spoznať už na základe nejakej matice v stupňovitom tvare, riadkovo ekvivalentnej s pôvodnou rozšírenou maticou sústavy. V takom prípade si teda môžeme odpustiť nielen ďalšiu úpravu na redukovaný stupňovitý tvar, ale aj spätné dosadzova- nie. Cvičenia 1. Podrobne prepočítajte sústavy lineárnych rovníc z príkladov 3.3.3­3.3.5. 2. Gaussovou-Jordanovou elimináciou riešte sústavu lineárnych rovníc A x = b nad poľom R pre matice: (a) A = 2 1 2 0 -6 1 0 0 2 -5 1 1 2 0 -3 0 1 2 0 1 , b = 4 1 1 2 ; (b) A = 3 -6 1 10 1 -2 1 4 2 -4 3 9 1 -2 0 3 , b = 6 3 2 2 . 3. Gaussovou-Jordanovou elimináciou riešte sústavu lineárnych rovníc A x = b nad poľom C pre matice: (a) A = 1+i 0 i-1 2-i 1 1+i , b = 1-2i 0 ; (b) A = 1+2i 2-i 1-i 5i 5 3-i , b = 1 2+i . 4. Gaussovou-Jordanovou elimináciou riešte sústavu A x = b lineárnych rovníc nad poľom Z11 pre matice: 12 PAVOL ZLATOŠ: LINEÁRNA ALGEBRA A GEOMETRIA (a) A = 1 2 3 4 2 3 1 1 7 3 5 3 , b = 1 5 1 ; (b) A = 1 1 4 5 3 2 0 1 1 9 1 7 , b = 0. Pre každú sústavu určte počet jej riešení. 5. Riešte sústavy z cvičenia 4 nad poľom Z13 a opäť určte počet riešení každej z nich. 6. V paragrafe 3.2 definovaný (redukovaný) stupňovitý tvar matice by sme mohli presnejšie nazvať riadkovým (redukovaným) stupňovitým tvarom. Sformulujte definíciu stĺpcového (redukovaného) stupňovitého tvaru matice. Podrobne definujte elementárne stĺpcové operácie (ESO) typov I, II a III. 7. Nech A = (aij)m×n, B = (bik)m×p sú matice nad poľom K. Označme bk = sk(B) k-ty stĺpec matice B. Uvažujme maticovú rovnicu A X = B s neznámou maticou X = (xjk)n×p. Dokážte postupne nasledujúce tvrdenia: (a) Matica X Kn×p je riešením maticovej rovnice A X = B práve vtedy, keď pre každé k p je jej k-ty stĺpec xk = sk(X) riešením sústavy A x = bk. (b) Maticová rovnica A X = B má riešenie práve vtedy, keď každá zo sústav A x = bk (k = 1, . . . , p) má riešenie. Na základe (a) a (b) navrhnite metódu, ako možno úpravou vhodnej blokovej matice pomocou ERO riešiť naraz viacero sústav Ax = bk s rovnakou ľavou stranou a rôznymi pravými stranami b1, . . . , bp. 8. Nech A = (aij)m×n, C = (ckj)q×n sú matice nad poľom K. S hromadným riešením akých sústav lineárnych rovníc súvisí riešenie maticovej rovnice Y A = C s neznámou maticou Y = (yki)q×m? Navrhnite metódu založenú na úprave vhodnej blokovej matice pomocou ESO. Ako sa možno vyhnúť ESO a nahradiť ich ERO? 9. Riešte maticové rovnice A X = B a Y A = C pre matice A = 3 2 4 2 1 3 , B = 1 3 -1 1 2 4 , C = 12 7 17 3 3 7 nad poľom Q. Určte najprv rozmery matíc X a Y . Aké sústavy lineárnych rovníc ste takto vyriešili? Napíšte riešenie každej z nich (ak existuje). 10. Nad poľom R riešte sústavy lineárnych rovníc v neznámych x, y, z a urobte diskusiu počtu riešení vzhľadom na parametre a R resp. c, d R: (a) x + y + (2a2 - 1)z = a2 + a + 1, (b) cx + y + (c + 1)dz = c + 2d + 1, x + y + (a2 + a - 1)z = 1, cx + cdz = d + 1, x + a2z = a2 + a; cy + 2cdz = 2c2 + cd - 1. 11. Nech K je pole a m, n N. Dokážte, že vzťah A B riadkovej ekvivalencie na množine Km×n spĺňa podmienky A A, A B B A a A B & B C A C pre ľubo- voľné A, B, C Km×n. Inak povedané, tento vzťah je reflexívny, symetrický a tranzitívny, teda je naozaj reláciou ekvivalencie na množine Km×n (pozri paragraf 0.6). Sformulujte analogický výsledok pre vzťah stĺpcovej ekvivalencie A B. 12. Dokážte vetu 3.3.2 matematickou indukciou podľa počtu riadkov matice. (Návod: Ukážte, že matica s jediným riadkom je riadkovo ekvivalentná s maticou v redukovanom stupňovitom tvare. Predpokladajte, že matica A vznikla z matice v redukovanom stupňovitom tvare pridaním jedného riadku. Ukážte, že aj A je riadkovo ekvivalentná s maticou v redukovanom stupňovitom tvare.) 13. Dokážte jednoznačnosť redukovaného stupňovitého tvaru matice. Presnejšie, dokážte, že pre ma- tice A, B Km×n v redukovanom stupňovitom tvare platí A B A = B. (Návod: Ak A, B sú v redukovanom stupňovitom tvare a A = B, ukážte, že príslušné homogénne sústavy lineárnych rovníc nemajú rovnaké riešenia; to je však spor s A B.) 14. Dokážte zosilnenie tvrdenia 3.3.1 do podoby ekvivalencie, t. j. pre ľubovoľné A, B Km×n, b, c Km platí: sústavy A x = b a B x = c sú ekvivalentné práve vtedy, keď (A| b) (B | c). (Návod: Modifikujte myšlienku z predchádzajúceho cvičenia na nehomogénne sústavy.) 15. Dokážte tvrdenie 3.3.6 (b).