Vzdialenosť mimobežných podpriestorov

Moderators: Martin Sleziak, Ludovit_Balko, Martin Niepel, Tibor Macko

Post Reply
Martin Sleziak
Posts: 5817
Joined: Mon Jan 02, 2012 5:25 pm

Vzdialenosť mimobežných podpriestorov

Post by Martin Sleziak »

Skúsim sa ešte vrátiť k úlohe na nájdenie vzdialenosti dvoch mimobežných podpriestorov.

Najprv si povedzme nejaké všeobecné veci: Majme afinné podpriestory zadané ako bod + vektorová zložka, t.j. αA1+V1, βA2+V2. Chceme vyrátať vzdialenosť ρ(A,B).

Označme V=V1+V2. Vieme, že vzdialenosť je rovná
ρ(α,β)=ρ(A1+V,A2+V)=ρ(A1,A2+V).
T.j. je rovnaká, ako vzdialenosť medzi rovnobežnými afinnými podpriestormi, ktoré sme viedli cez zadané afinné podpriestory.
Tým sme vlastne previedli úlohu na hľadanie vzdialenosti bodu od afinného podpriestoru. Pre priamku sme takúto úlohu riešili tu: viewtopic.php?f=29&t=623 Niektoré z postupov uvedených tam fungujú všeobecne

Iná možnosť je zobrať ľubovoľné body z jedného a druhého podpriestoru a hľadať priemet vektora spájajúceho tieto dva body do V. Veľkosť kolmého priemetu je presne vzdialenosť.
T.j. vieme, že A1A2 sa dá jednoznačne rozložiť ako u+v, kde uV a vV. Ak nájdeme vektor v, jeho dĺžka je hľadaná vzdialenosť.
Hľadať kolmý priemet sme sa naučili robiť už v minulom semestri: viewtopic.php?f=29&t=574

Ďalej vieme, že pre mimobežné afinné podpriestory existuje stredná priečka, t.j. existujú body Pα, Qβ také, že PQ je kolmý na α aj na β. Potom platí ρ(α,β)=|PQ|. Čiže ak nájdeme strednú priečku, stačí vyrátať jej dĺžku.

Poďme sa skúsiť pozrieť na rôzne postupy výpočtu vzdialenosti na konkrétnom príklade.
Vyrátajte vzdialenosť priamky p a roviny α v R4, ak priamka p je určená bodom A=(1,0,1,0) a smerovým vektorom u=(0,1,1,1) a rovina α je určená bodom B=(1,2,3,4) a vektormi v=(1,1,0,0) a w=(0,0,1,1).
Túto úlohu som zobral z http://thales.doa.fmph.uniba.sk/sleziak ... ohyvzd.pdf (Kde môžete nájsť tiež nejaké riešenia.)
Martin Sleziak
Posts: 5817
Joined: Mon Jan 02, 2012 5:25 pm

Re: Vzdialenosť mimobežných podpriestorov

Post by Martin Sleziak »

Riešenie pomocou pomocného podpriestoru

Chceme nájsť vzdialenosť bodu A od afinného podpriestoru určeného bodom B a vektormi u, v, w. Tieto vektory sú lineárne nezávislé, ide teda o nadrovinu. Označme ju β.

Vidíme, že Vβ=[(1,1,0,0),(0,0,1,1),(0,1,1,1)]=[(1,1,0,0),(0,0,1,1),(0,1,0,0)]=[(1,0,0,0),(0,1,0,0),(0,0,1,1)].

Z toho zistíme, že Vβ=(0,0,1,1). Teda všeobecná rovnica nadroviny β má tvar x3x4+d=0. Dosadením bodu B zistíme, že je to
x3x4+1=0.
Teraz stačí vyrátať
ρ(A,β)=|10+1|12+(1)2=22=2.

Tu sme mali výhodu, že pomocný podpriestor mal kodimenziu 1 (bola to nadrovina), a teda sme mohli priamo využiť vzorec. Poďme to skúsiť vyrátať bez použitia vzorca. (Napríklad keby sme riešili podobnú úlohu v R5, tak β ako podpriestor dimenzie 3 v 5-rozmernom priestore by už nebola nadrovina a nemáme k dispozícii vzorec, ktorí sme použili pred chvíľou.)

Môžeme nájsť vzdialenosť A od β napríklad tak, že hľadáme kolmý priemet A.
Ten môžeme nájsť ako prienik kolmopremietacieho priestoru πβ(A) s nadrovinou β. Kolmopriemietací priestor je určený bodom A a vektorovou zložkou Vβ. Teda je to priamka určená rovnicami x1=1, x2=0, x3=1+t, x4=t.
Hľadáme prienik tejto priamky s rovinou. Pre zmenu skúsme teraz použiť parametrické vyjadrenie. (Môžete si vyskúšať, že s použitím analytického vyjadrenia nadroviny β dostanete ten istý výsledok.)
Dostávame rovnice:
1=1+a0=2+b1+t=3+ct=4+c
Táto sústava má jediné riešenie t=1, a=0, b=2, c=3.
Tieto hodnoty parametrov určujú bod A=(1,0,0,1).
(Môžeme skontrolovať, že A skutočne patrí do β a že AA=(1,0,0,1)(1,0,1,0)=(0,0,1,1)Vβ.)
Dostávame ρ(A,B)=|AA|=2
Martin Sleziak
Posts: 5817
Joined: Mon Jan 02, 2012 5:25 pm

Re: Vzdialenosť mimobežných podpriestorov

Post by Martin Sleziak »

Kolmý priemet

Úlohu môžeme riešiť aj tak, že vezmeme vektor AB=(1,2,3,4)(1,0,1,0)=(0,2,2,4) a hľadáme jeho kolmý priemet do Vβ.

Kolmý priemet vieme nájsť viacerými spôsobmi.

V tomto prípade je priestor, do ktorého robíme priemet, jednorozmerný, takže to ide veľmi jednoducho. (Pozri aj: viewtopic.php?t=851 )

Zoberieme si jednotkový vektor generujúci tento priestor, čo je 12(0,0,1,1).

Vyrátame skalárny súčin so zadaným vektorom: 12(0,2,2,4),(0,0,1,1)=22.

Hľadaný kolmý priemet je 2212(0,0,1,1)=(0,0,1,1). Jeho veľkosť je 2.

(Môžeme sa presvedčiť, že (0,2,2,4)=(0,0,1,1)Vβ+(0,2,3,3)Vβ)
Martin Sleziak
Posts: 5817
Joined: Mon Jan 02, 2012 5:25 pm

Re: Vzdialenosť mimobežných podpriestorov

Post by Martin Sleziak »

Stredná priečka

Pretože dané afinné podpriestory sú mimobežné, podľa vety z prednášky existuje vektor PQ, taký, že Pp, Qα, PQ je kolmý na vektorové zložky oboch afinných podpriestorov a jeho veľkosť je práve vzdialenosť týchto afinných podpriestorov. (Priamka PQ sa nazýva stredná priečka.)

Parametrické vyjadrenia sú
p={P=A+su},
α={Q=B+tv+uw},
z čoho dostaneme
PQ=AB+tv+uwsu=(t,2s+t,2s+u,4s+u)

Podmienku, že PQ je kolmý na vektorové zložky oboch priestorov môžeme vyjadriť pomocou skalárneho
súčinu: PQ,u=PQ,v=PQ,w=0. Z~toho dostaneme sústavu
rovníc
83s+t+2u=0,2s+2t=0,62s+2u=0,
z~ktorej vyrátame t=0, s=2, u=1.

Preto
ρ(p,α)=PQ=(1)2+12=2

Súčasne sme našli P=(1,2,3,2) a Q=(1,2,2,3). Tieto body skutočne určujú strednú priečku; jeden z nich patrí zadanej priamke a druhý zadanej rovine. Vektor PQ=(0,0,1,1) je kolmý na Vp aj Vα.

Iná ukážka výpočtu strednej priečky: viewtopic.php?t=870
Martin Sleziak
Posts: 5817
Joined: Mon Jan 02, 2012 5:25 pm

Re: Vzdialenosť mimobežných podpriestorov

Post by Martin Sleziak »

Minimalizácia vzdialenosti

Máme parametrické vyjadrenie priamky aj roviny - pre ľubovoľný bod XP, YQ vieme vyjadriť |XY|2 ako funkciu parametrov. Túto funkciu chceme minimalizovať. Keďže vyjadrenie je pomerne jednoduché (je to kvadratická funkcia), mohlo by sa nám podariť nájsť minimum.
Takéto riešenie je detailne rozpísané tu: http://thales.doa.fmph.uniba.sk/sleziak ... ohyvzd.pdf
Post Reply