Program a domácí úkol z třetího cvičení 5.10.2018 Program. Součet podprostorů (řešený 2.1b,c), afinní repér (vrcholy rovnoběžnostěnu), transformace od repéru k repéru (př 14), orientace afinního prostoru (úloha 4.1), parametrické vyjádření podprostoru Příklad 1. V A3 je dán rovnoběžnostěn ABCDEFGH a repéry R1 = A; −→ AB, −→ AC, −→ AE a R2 = C; −−→ CB, −−→ CD, −→ CG . Určete: • transformační rovnice od R1 k R2 • transformační rovnice od R2 k R1 Nápověda: Ještě než začnete s tvorbou transformačních rovnic, je potřeba mít vše vyjádřené vůči nějakému pevnému, tzn. stejnému repéru. Například R2 bychom vůči R1 vyjádřili jako R2 = [1, 1, 0], (0, −1, 0), (−1, 0, 0), (0, 0, 1) . Příklad 2. Jsou dány dva afinní podprostory B1 = [2, 2, 3, 2] + t(1, 1, 2, 1) + s(1, 1, 0, −1) a B2 = [1, 0, 1, 1], [2, 1, 0, 1], [1, 0, 1, 2], [0, 0, 1, 1] . Určete: • zda se dané podprostory protínají • jaká je dimenze součtu B1 + B2 • nějakou bázi součtu B1 + B2 Slíbený postup získání transformačních rovnic Nalezněte matice přechodu od repéru U k repéru V a určete souřadnice bodu X v těchto repérech, je-li (všechny souřadnice jsou zadány v nějakém pevném repéru E): • U = {U; u1; u2; u3} = {[1, 1, 1]; (4, 0, 2); (3, 1, 1); (−1, 4, 0)}, • V = {V ; v1; v2; v3} = {[−4, 11, 1]; (1, −1, 1); (10, 8, 6); (8, 13, 5)}, • X[1, 2, 3]. Řešení. Naším prvním úkolem bude najít dvě matice A, B tak, aby platilo (X)U = A · (X)V + B (1) Ukážeme přitom dva různé postupy řešení, nejprve „maticově“ (tj. tak jako na cvičení), poté z definice matice přechodu. 1. Místo abychom určovali přímo matice A a B, zkusíme určit nejdříve matice přechodu (a příslušné matice souřadnic počátků) z referenčního repéru E do U, resp. V. Hledáme proto matice C, D, E, F takové, že bude platit (X)E = C · (X)U + D (2) (X)E = E · (X)V + F (3) Pro souřadnice vektorů platí tytéž transformační rovnice, jen bez matic D, F: (x)E = C · (x)U (x)E = E · (x)V Určíme matice C, E. Když např. do (x)U dosadíme postupně souřadnice vektorů (1, 0, 0); (0, 1, 0); (0, 0, 1) (toto jsou vlastně souřadnice vektorů báze U vzhledem k repéru U!), musíme na levé straně postupně dostat souřadnice bázových vektorů z U vzhledem k referenčnímu repéru E. Tyto souřadnice tedy musí být v matici C postupně ve sloupcích. Obdobnou úvahou pro matici E získáme podobu jednotlivých matic (ověřte si, že skutečně po pronásobení vektory (1, 0, 0); (0, 1, 0); (0, 0, 1) zprava dostaneme skutečně souřadnice bázových vektorů U vzhledem k E!): C =   4 3 −1 0 1 4 2 1 0   ; E =   1 10 8 −1 8 13 1 6 5   Určíme matice D, F. Když např. do (X)U (z vyjádření (2)) dosadíme souřadnice bodu [0, 0, 0] (to jsou souřadnice bodu U vzhledem k repéru U), musíme na levé straně dostat souřadnice bodu U vzhledem k referenčnímu repéru E. Je proto asi jasné, že matice D a F budou vypadat takto (opět zkuste ověřit pronásobením): D =   1 1 1   ; F =   −4 11 1   Porovnáním maticových rovnic (2) a (3) dostáváme rovnost, kterou dále upravíme na hledaný tvar ve vyjádření (1): C · (X)U + D = E · (X)V + F C · (X)U = E · (X)V + F − D (X)U = (C−1 · E) A · (X)V + C−1 · (F − D) B Konkrétní tvar matic A, B dostaneme např. tak, že vedle sebe postupně napíšeme do oddělených bloků matice C, E a F − D, následně pak upravíme blok matice C na jednotkovou matici – v prostředním bloku dostaneme matici A = (C−1 ·E), v pravém bloku B = C−1 · (F − D).   4 3 −1 1 10 8 −5 0 1 4 −1 8 13 10 2 1 0 1 6 5 0   ∼ · · · ∼   1 0 0 1 3 2 1 0 1 0 −1 0 1 −2 0 0 1 0 2 3 3   Výsledky doplníme do rovnice (1): (X)U = A · (X)V + B =   1 3 2 −1 0 1 0 2 3   · (X)V +   1 −2 3   (4) Nyní ještě musíme vypočítat souřadnice bodu X v obou repérech. Využijeme k tomu rovnici (3), vypočteme souřadnice X v repéru V a pak pomocí (4) souřadnice v repéru U (děláme to tak, protože je to asi nejlehčí, při přechodu zprava doleva stačí jen dosadit souřadnice, kdežto zleva doprava je nutné vypočítat inverzní matici nebo řešit soustavu rovnic). Dosadíme souřadnice bodu do rovnice (3) (uvádíme již v roznásobené formě): 1 = x1 + 10x2 + 8x3 − 4 2 = −x1 + 8x2 + 13x3 + 11 3 = x1 + 6x2 + 5x3 + 1 Dále řešíme tuto soustavu.   1 10 8 5 −1 8 13 −9 1 6 5 2   ∼ · · · ∼   1 0 0 −4 3 0 1 0 5 2 0 0 1 −7 3   Stačí už jen dosadit vypočítané souřadnice do vyjádření (4) a zjistit souřadnice bodu X vzhledem k repéru U.   1 3 2 −1 0 1 0 2 3   ·   −4 3 5 2 −7 3   +   1 −2 3   =   5 2 −3 1   Výsledky přehledně shrneme. (X)U =   1 3 2 −1 0 1 0 2 3   · (X)V +   1 −2 3   , X = 5 2 ; −3; 1 U = −4 3 ; 5 2 ; −7 3 V , 2. Nyní ukážeme jiný způsob, jak najít matici A, a to přímo z definice matice přechodu. Podobně jako v prvním domácím úkolu hledáme koeficienty a, b, c, d, e, f, g, h, i ∈ R takové, že bude splněna soustava rovnic v1 = au1 + bu2 + cu3 (1, −1, 1) = a(4, 0, 2) + b(3, 1, 1) + c(−1, 4, 0) v2 = du1 + eu2 + fu3 ⇐⇒ (10, 8, 6) = d(4, 0, 2) + e(3, 1, 1) + f(−1, 4, 0) v3 = gu1 + hu2 + iu3 (8, 13, 5) = g(4, 0, 2) + h(3, 1, 1) + i(−1, 4, 0) Rozepsáním rovnic získáme tři soustavy rovnic, kde má každá soustava tři rovnice a tři neznámé. 1 = 4a + 3b − c 10 = 4d + 3e − f 8 = 4g + 3h − i −1 = + b + 4c 8 = + e + 4f 13 = + h + 4i 1 = 2a + b 6 = 2d + e 5 = 2g + h Je vidět, že se prakticky jedná o jednu soustavu s odlišnými absolutními členy a jiným označením neznámých. Proto můžeme řešit všechny tři soustavy naráz úpravou levého bloku následující matice na jednotkovou matici (abychom mohli přímo odečíst matici A):   4 3 −1 1 10 8 0 1 4 −1 8 13 2 1 0 1 6 5   ∼ · · · ∼   1 0 0 1 3 2 0 1 0 −1 0 1 0 0 1 0 2 3   Pravý blok upravené matice (oddělený dvojitou čarou) je zřejmě hledanou maticí přechodu od báze U k bázi V, neboť se vlastně jedná o matici řešení soustav ve tvaru:   1 0 0 a d g 0 1 0 b e h 0 0 1 c f i   , (pokud nedošlo k výměně sloupců!) Tedy A =   1 3 2 −1 0 1 0 2 3  . Způsobem popsaným v první části řešení dostaneme matici B, problém nastává při přepočtu souřadnic bodu X, protože nemáme k dispozici (zatím. . . ) nějaké matice přechodu mezi repéry E a U nebo V. Protože ale máme transformační rovnice od repéru U k repéru V, vypočítáme matice přechodu od repéru V k repéru E a souřadnice pak jen dosadíme do již známé rovnice. Hledáme proto koeficienty j, k, l, m, n, o, p, q, r ∈ R takové, že bude splněna soustava rovnic e1 = jv1 + kv2 + lv3 (1, 0, 0) = j(1, −1, 1) + k(10, 8, 6) + l(8, 13, 5) e2 = mv1 + nv2 + ov3 ⇐⇒ (0, 1, 0) = m(1, −1, 1) + n(10, 8, 6) + o(8, 13, 5) e3 = pv1 + qv2 + rv3 (0, 0, 1) = p(1, −1, 1) + q(10, 8, 6) + r(8, 13, 5) Rozepsáním rovnic získáme tři soustavy rovnic, kde má každá soustava tři rovnice a tři neznámé. 1 = j + 10k + 8l 0 = m + 10n + 8o 0 = p + 10q + 8r 0 = −j + 8k + 13l 1 = −m + 8n + 13o 0 = −p + 8q + 13r 0 = j + 6k + 5l 0 = m + 6n + 5o 0 = p + 6q + 5r Můžeme opět řešit všechny tři soustavy naráz úpravou levého bloku následující matice na jednotkovou matici (abychom mohli přímo odečíst hledanou matici přechodu – všimněte si, že vlastně takhle počítáme inverzní matici k matici E! – s ohledem na očekávaný výsledek to nebude hezké počítání. . .   1 10 8 1 0 0 −1 8 13 0 1 0 1 6 5 0 0 1   ∼ · · · ∼   1 0 0 −19 15 − 1 15 11 15 0 1 0 3 5 − 1 10 − 7 10 0 0 1 − 7 15 2 15 3 5   Platí tedy (X)V =   −19 15 − 1 15 11 15 3 5 − 1 10 − 7 10 − 7 15 2 15 3 5   · (X)E + G Matice G je zatím neznámá matice souřadnic počátku repéru E v souřadnicích vzhledem k V. Tu je ale nutné vypočítat, neboť odhadnout ji asi není v našich silách. Bude se jednat vlastně o řešení rovnice se třemi parametry s, t, u, která přejde na soustavu o třech rovnicích se třemi neznámými. [0, 0, 0] = [−4, 11, 1] + s(1, −1, 1) + t(10, 8, 6) + u(8, 13, 5) 4 = s + 10t + 8u −11 = −s + 8t + 13u −1 = s + 6t + 5u   1 10 8 4 −1 8 13 −11 1 6 5 −1   ∼ · · · ∼   1 0 0 −98 15 0 1 0 21 5 0 0 1 −59 15   Nyní už můžeme dosadit souřadnice bodu X do transformačních rovnic.   −19 15 − 1 15 11 15 3 5 − 1 10 − 7 10 − 7 15 2 15 3 5   ·   1 2 3   +   −98 15 21 5 −59 15   =   −4 3 5 2 −7 3     1 3 2 −1 0 1 0 2 3   ·   −4 3 5 2 −7 3   +   1 −2 3   =   5 2 −3 1   Výsledky přehledně shrneme. (X)U =   1 3 2 −1 0 1 0 2 3   · (X)V +   1 −2 3   , X = 5 2 ; −3; 1 U = −4 3 ; 5 2 ; −7 3 V , Resumé: K výsledku vede několik cest, některé jsou však podstatně trnitější, než jiné. . .