Domácí úkol z 5. 10. 2015 Příklad 2. Nalezněte matice přechodu od repéru U k repéru V a určete souřadnice bodů X, Y 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} = {[−3, 0, 5]; (1, −1, 1); (10, 8, 6); (8, 13, 5)}, • X[1, 2, 3]; Y [4, 5, 4]. Ř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 =   −3 0 5   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 −4 0 1 4 −1 8 13 −1 2 1 0 1 6 5 4   ∼ · · · ∼   1 0 0 1 3 2 69 10 0 1 0 −1 0 1 −49 5 0 0 1 0 2 3 11 5   Výsledky doplníme do rovnice (1): (X)U = A · (X)V + B =   1 3 2 −1 0 1 0 2 3   · (X)V +   69 10 −49 5 11 5   (4) Nyní ještě musíme vypočítat souřadnice bodů X a Y v obou repérech. Využijeme k tomu rovnici (3), vypočteme souřadnice obou bodů 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 obou bodů do rovnice (3) (uvádíme již v roznásobené formě): 1 = x1 + 10x2 + 8x3 − 3 4 = x1 + 10x2 + 8x3 − 3 2 = −x1 + 8x2 + 13x3 5 = −x1 + 8x2 + 13x3 3 = x1 + 6x2 + 5x3 + 5 4 = x1 + 6x2 + 5x3 + 5 Podobně jako v minulém domácím úkolu se jedná o tu samou soustavu s jinými absolutními členy – proto můžeme obě soustavy řešit najednou úpravou levého bloku následující matice na jednotkovou matici. V prostředním bloku pak budou hledané souřadnice bodu X, v pravém bodu Y (obojí vzhledem k repéru V).   1 10 8 4 7 −1 8 13 2 5 1 6 5 −2 −1   ∼ · · · ∼   1 0 0 −48 5 −57 5 0 1 0 18 5 22 5 0 0 1 −14 5 −16 5   Stačí už jen dosadit vypočítané souřadnice do vyjádření (4) a zjistit souřadnice bodů X, Y vzhledem k repéru U.   1 3 2 −1 0 1 0 2 3   ·   −48 5 18 5 −14 5   +   69 10 −49 5 11 5   =   5 2 −3 1     1 3 2 −1 0 1 0 2 3   ·   −57 5 22 5 −16 5   +   69 10 −49 5 11 5   =   23 10 −8 5 7 5   Výsledky přehledně shrneme. (X)U =   1 3 2 −1 0 1 0 2 3   · (X)V +   69 10 −49 5 11 5   , X = 5 2 ; −3; 1 U = −48 5 ; 18 5 ; −14 5 V , Y = 23 10 ; −8 5 ; 7 5 U = −57 5 ; 22 5 ; −16 5 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 bodů X, Y , 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] = [−3, 0, 5] + s(1, −1, 1) + t(10, 8, 6) + u(8, 13, 5) 3 = s + 10t + 8u 0 = −s + 8t + 13u −5 = s + 6t + 5u   1 10 8 3 −1 8 13 0 1 6 5 −5   ∼ · · · ∼   1 0 0 −74 5 0 1 0 53 10 0 0 1 −22 5   Nyní už můžeme dosadit souřadnice bodů X a Y do transformačních rovnic. X :   −19 15 − 1 15 11 15 3 5 − 1 10 − 7 10 − 7 15 2 15 3 5   ·   1 2 3   +   −74 5 53 10 −22 5   =   −48 5 18 5 −14 5     1 3 2 −1 0 1 0 2 3   ·   −48 5 18 5 −14 5   +   69 10 −49 5 11 5   =   5 2 −3 1   Y :   −19 15 − 1 15 11 15 3 5 − 1 10 − 7 10 − 7 15 2 15 3 5   ·   4 5 4   +   −74 5 53 10 −22 5   =   −57 5 22 5 −16 5     1 3 2 −1 0 1 0 2 3   ·   −57 5 22 5 −16 5   +   69 10 −49 5 11 5   =   23 10 −8 5 7 5   Výsledky přehledně shrneme. (X)U =   1 3 2 −1 0 1 0 2 3   · (X)V +   69 10 −49 5 11 5   , X = 5 2 ; −3; 1 U = −48 5 ; 18 5 ; −14 5 V , Y = 23 10 ; −8 5 ; 7 5 U = −57 5 ; 22 5 ; −16 5 V Resumé: K výsledku vede několik cest, některé jsou však podstatně trnitější, než jiné. . .