Rectas, Rectas , distancias distanc ias y lugares geom´etricos etricos con Mathematica Mariano Gonz´ alez alez Ulloa
[email protected] Pontificia Universidad Cat´olica olica del Per´u
Departament Departamento o de Ciencias Ciencias 18 de agosto de 2013 Resumen En esta presentaci´on on se desarrolla procedimientos para construir e identificar lugares geom´etricos etricos que involucran la distancia de un punto a una recta y la distancia entre dos rectas. Especialmente se trata de identificar el lugar geom´ etrico etrico de los puntos que equidistan de dos rectas alabeadas, el cual resulta ser una superficie. Para ello se trabajar´a con una representaci´on on vectorial de las rectas y se usar´a el software Mathematica para Mathematica para obtener una parametrizaci´on on de la superficie y a la vez construir las gr´aficas aficas correspondientes.
1.
Intr Introdu oducc cci´ i´ on on
Si bien es cierto que los lugares geom´etricos etricos (LG) en el plano, generalmente, son f´ aciles aciles de visualizar, eso no ocurre cuando se trata de lugares geom´etricos etricos en el espacio tridimensional. Por ello, discusiones como la que se presenta en esta publicaci´on son de interes, ya que con la ayuda del software Mathematica Mathematica es posible conjeturar las posibles soluciones y luego visualizar los lugares geom´ etricos etricos ba jo las condiciones del problema planteado. En esta presentaci´on on se propone la construcci´on on e identificaci´on on del lugar lug ar geom´ g eom´etrietri co de los puntos que equidistan de dos rectas abaleadas. El procedimiento seguido consiste en identificar el punto, de cada una de las rectas, que proporciona la m´ınima ınima distancia entre un punto fijo y cada recta, luego se igualan estas dos distancias para obtener la ecuaci´on on cartesiana del LG que resulta ser una superficie. Luego, usando bases de Gr¨obner obner de un ideal en el anillo de polinomios R[x , y , z , s , t], se obtiene una representaci´ on on param´etrica etrica para dicha superficie super ficie
2. 2.1.
Nocio Nocione ness prel prelim imin inar ares es Lugar geom´ etrico etrico
Definici´ on on 2.1. Un lugar geom´ geom´etrico etrico (LG) es el conjunto conjunto de puntos que satisfacen cierta o ciertas propiedades.
1
Si el lugar geom´etrico es definido por la propiedad P, entonces: Todo punto de LG satisface la propiedad P. Todo punto que satisface la propiedad P pertenece al LG. Por ejemplo, una esfera con centro en el punto O R3 y radio r > 0 es el lugar geom´ etrico de todos los puntos P R3 cuya distancia al punto O es igual a r. El punto O es el centro de la esfera y r su radio. Tambi´ en, dado un punto fijo P 0 R3 y un vector fijo v R 3 , la recta que pasa por P 0 y es paralela al vector v es el lugar geom´etrico de todos los puntos P R 3 tales que el vector P P 0 es paralelo al vector v , es decir P = P 0 + t v , t R
∈
∈
→
∈
∈
→
−
→
2.2.
∈
→
∈
Distancia de un punto a una recta
Previamente describiremos la manera como hallar la distancia de un punto a una recta. Para ello consideremos una recta L generada por una funci´on vectorial α (t) = P 0 + t v , t R , donde P 0 R3 es un punto fijo y v R 3 es el vector direcci´on de L , y un punto P R 3 , la distancia de P a L est´ a dada por →
→
∈
→
∈
∈
∈
→
v × (P − P ) d(P, L) = v o
(1)
→
expresi´on que solamente vale en dos vectores.
3
R
debido a que depende del producto vectorial de
De manera equivalente, para hallar la distancia del punto P a la recta L bastar´a encontrar el pie de la recta perpendicular a L trazada por el punto P (punto Q, ver figura 1) y luego, simplemente, calcular la distancia entre los puntos P y Q. El punto Q se obtiene mediante →
Q = P 0 +
v (P v
· − P ) 0
2
→
(2)
v
En consecuencia, la distancia del punto P a la recta L est´ a dada por (figura 1) d(P, L) = d (P, Q)
(3)
La funci´ on DistaPuntoRecta[P , P 0 , V ]
definida en el programa 1, permite obtener la distancia de un punto P a la recta que pasa por el punto P o y tiene la direcci´on del vector V .
DistaPuntoRecta P_, Po_, V_ : FullSimplify Norm P
Programa 1
2
V . P
Po
Po
Norm V
2
V
Y
L P
Z
Po
Q
X
Figura 1: Distancia de un punto a una recta Por ejemplo, para hallar la distancia del punto P = (1, 1, 1) a la recta L que pasa por el punto P 0 = (0, 0, 0) y tiene la direcci´on del vector v = (1, 0, 0), ingrese la funci´on →
{
} {0 , 0 , 0}, {1, 0, 0}]
DistaPuntoRecta[ 1 , 1 , 1 ,
que al ejecutarla tendr´a
√ 2 Observaci´ on 2.1. Los par´ ametros en la funci´ on DistaPuntoRecta pueden ser vectores de dos o tres componentes. Esta es una caracter´ıstica de Mathematica, la definici´ on de una funci´ on no est´ a ligada, necesariamente, a la dimensi´ on del espacio que contiene a su dominio. En el programa 2, usando la expresi´on (3), se define la funci´on DistaPuntoRectaGraf[P , Po , V , a , b ]
para calcular la distancia de un punto a una recta en R3 y construir las gr´aficas correspondientes. Los par´ametros de esta funci´on son: P el punto desde donde se calcular´ a la distancia, P o el punto de paso de la recta, V el vector direcci´on de la recta; y a y b los extremos del intervalo donde toma valores la variable para construir la recta.
3
to
DistaPuntoRectaGraf P_, Po_, V_, a_, b_ : Module
P
a, t1
b ,
Po .V
punto en la recta DistanciaPL Simplify Norm P Q; GRecta ParametricPlot3DPo t V , t, to, t1, Axes True, Boxed True, AxesLabel "X", "Y", "Z", PlotStyle Blue, Thick; GPuntosPQ Graphics3DThick, Arrowheads.03, Arrow0, 0, 0, V , PointSize Medium, PointP, PointSize Medium, PointQ, PointSize Medium, PointPo, Thick, Red, Line P, Q, Thick, Dashed, Line P, Po; V Texto Graphics3DText"v", 0.4, 0.4, 0.4, Text"P", P 0.6, 0.6, 0.6, Text"L", 2 Po 2 V .6, .6, .8, Text"Q", Q 0.6, .6, .2, Text"Po", Po 0.6, .6, .2 ; Print"distancia ", P, ",", Q, " ", DistanciaPL, " ", NDistanciaPL; ShowGRecta, GPuntosPQ, Texto Q
Po
Norm V
2
V;
Programa 2 Con la finalidad de verificar la buena definici´on de la funci´on DistaPuntoRectaGraf, considere el punto P 0 ( 1; 2; 3), el vector direcci´on V = (2; 1; 1) y el punto P (6;1;5). Ingrese
−
− −
DistaPuntoRectaGraf[{6,1,5},{-1,2,3},{2,-1,-1},-5,5]
al ejecutarla obtendr´a el resultado que se muestra en la figura 2:
L
6 4
Po
P 2
Z
0 v
5
6
Q
4
X
2
0 5
6, 1, 5,
Y
0
10
distancia
1 ,
3
5 ,
6
6
→
2
155 5.08265
6
Figura 2: Distancia de P (6, 1, 5) a L : α (t) = ( 1, 2, 3) + t(2, 1, 1)
4
−
− −
2.3.
Distancia entre dos rectas
Si las rectas son paralelas, la distancia entre ellas es igual a la distancia de cualquier punto de una de ellas a la otra recta. De manera que se puede asumir que las rectas son no paralelas y tampoco se intersecan. Definici´ on 2.2. Dos rectas en R3 que no son paralelas ni se intersecan se denominan rectas alabeadas. →
→
Sean L 1 y L 2 rectas alabeadas generadas por las funciones α (s) = P 1 + s u , s →
→
∈ R
y β (t) = P 2 + t v , t R , respectivamente. Los puntos P 1 , P 1 + u , y P 1 + v no son colineales ya que los vectores u y v son no paralelos, en consecuencia dichos puntos determinan el plano π1 que contiene a la recta L1 . De la misma manera, los puntos P 2 , P 2 + u , y P 2 + v determinan el plano π2 que contiene a la recta L2 . Adem´as los planos π1 y π2 son paralelos puesto v . De aqu´ı se concluye que existen dos u que tienen el mismo vector normal, u ´ nicos planos paralelos, π1 y π2 , que contienen a las rectas alabeadas L1 y L2 , respectivamente (figura 3). El argumento desarrollado prueba la siguiente proposici´on.
∈
→
→
→
→
→
×
→
→
→
Figura 3: Planos que contienen a las rectas alabeadas Proposici´ on 2.1. Dadas dos rectas alabeadas L1 y L2 existen dos unicos ´ planos paralelos π1 y π2 que contienen a L1 y a L2 , respectivamente. Proposici´ on 2.2. La distancia entre dos rectas alabeadas L1 y L2 est´ a dada por →
→
|(P − P )· u × v | d(L , L ) = || u × v || 2
1
1
2
→
(4)
→
Prueba. La expresi´on 4 se obtiene observando que el volumen del paralelep´ıpedo (figura 4) determinado por los vectores u , v y P 2 P 1 es →
→
P
− vol(P ) = |(P − P )· u × v | (5) que tambi´ en es igual al producto del a´rea de su base, u × v , por su altura, h, es →
2
→
1
→
→
decir que
→
vol( ) = h u
P
→
× v
5
(6)
Figura 4: Distancia entre dos rectas alabeadas L1 y L2 en consecuencia, de (5) y (6) resulta que →
→
|(P − P )· u × v | = d(L , L ). h = u × v 2
1
→
1
→
2
La funci´ on DistaRectaRecta[P 1 , U , P 2 , V ]
definida en el programa 3, devuelve la distancia entre dos rectas en
3
R
.
, Paralelas SolveU Λ V 0, Λ ; V. P2 P1 AbsP2 P1 . CrossU, V IfParalelas, FullSimplify Norm P2 P1 V , Norm CrossU, V Norm V
DistaRectaRecta P1, U, P2, V : Module
2
Programa 3 Los par´ ametros de esta funci´on son P 1, U (punto de paso y vector direcci´o n de la recta L 1 ) y P 2, V (punto de paso y vector direcci´on de la recta L 2 ). Para verificar la buena definici´on de la funci´on DistaRectaRecta, considere L1 : (0, 0, 0) + t(1, 3, 4), t R y L2 : (1, 1, 1) + s(3, 4, 6), s R e ingrese
∈
∈
DistaRectaRecta[{0, 0, 0},{1, 3, 4},{1, 1, 1}, {3, 4, 6}]
luego de ejecutarla obtendr´a
√ 365 que es la distancia entre las rectas L1 y L 2 .
6
3.
Rectas, distancias y lugares geom´ etricos
En esta secci´on se describe un procedimiento para construir e identificar el lugar geom´etrico de los puntos que se encuentran a la misma distancia de un punto y de una recta y tambi´en el LG de los puntos que equidistan de dos rectas.
3.1. LG de los puntos que equidistan de un punto y de una recta →
→
Consideremos una recta L : α (t) = P 0 + t v , t R y un punto fijo F . Sea P un punto que equidista de L y F . La recta perpendicular a L que pasa por P interseca a L en el punto
∈
→
Q = P 0 +
v (P v
· − P ) 0
2
→
v,
luego, la ecuaci´on del LG de los puntos que equidistan de L y F est´ a dada por
d(P, Q) = d (P, F )
La funci´ on
LGDistaPuntoRecta[Paso , V ector , Foco , a , b ]
definida en el programa 4, permite obtener la ecuaci´on cartesiana del LG de los puntos que est´an a la misma distancia de un punto fijo F y de una recta fija L. Los par´ametros de esta funci´ on son: el punto de paso de la recta (Paso), el vector direcci´ on de la recta (Vector), el punto fijo (Foco) y los extremos del intervalo de la variable para construir la recta (a y b, respectivamente).
7
LGDistaPuntoRecta Paso_, Vector_, Foco_, a_, b_ :
Po
Module
t_ :
Paso, V Vector, F Foco, t0 a, t1 b ,
Po t V;
Recta L
Clear x, y, z ; P
x, y, z ; Punto genérico del LG punto de ortogonalidad en la recta L P Po .V QenL Simplify Po V ; Norm V 2
EcuacSuper Expand Norm P QenL Norm P F2 . Abs Identity; grafica de la recta RectaL ParametricPlot3D t, t, t0, t1, PlotStyle Red, Thick, PlotRange All, AxesLabel "X", "Y", "Z"; Texto1 Graphics3DText"L", 2 0.5, 0.5, .5, Text"F", F 0.5, 0.5, .5 ; GPuntoF Graphics3DRed, PointSize Medium, PointF; GRectaL ShowRectaL, Texto1, GPuntoF, PlotRange All; grafica del LG GraficaLG ContourPlot3DEcuacSuper 0, x, 10, 20, y, 10, 20, z, 10, 20, AxesLabel "X", "Y", "Z", Mesh 4, 4 , PlotRange All; presentacion PrintCollectEcuacSuper, x, y, z, Simplify, " 0"; ShowGRectaL, GPuntoF, GraficaLG
ecuacion cartesiana de la superficie
2
Programa 4 Para verificar la buena definici´on de la funci´on LGDistaPuntoRecta, ingrese la expresi´ on
{
}{
}{
}
LGDistaPuntoRecta[ 1,2,4 , -1,0,1 , 1,0,1 ,-10,15]
la salida se puede ver en la figura 5
3.2.
LG de los puntos que equidistan de dos rectas
Si las rectas son paralelas, el conjunto de puntos que equidistan de ambas rectas es un plano que se denomina plano mediador. Si las rectas se intersecan, los puntos que equidistan de ambas rectas constituyen los planos bisectores de los ´angulos que forman dichas rectas. En consecuencia podemos suponer que las rectas son alabeadas. Consideremos las rectas →
→
→
L1 : α (s) = P 1 + s u , s
∈ R
y
→
L2 : β (t) = P 2 + t v , t
∈ R,
(7)
y P (x , y , z ) un punto gen´ erico del conjunto de puntos que equidistan de L 1 y L 2 . Sean as pr´oximos a P . Por (2) dichos puntos son de la Q1 L1 y Q2 L2 los puntos m´ forma
∈
∈
→
Q1 = P 1 +
u (P
→
· − P ) u , u →
1
→
Q2 = P 2 +
2
8
v (P
· − P ) v →
2
2
→
v,
(8)
Figura 5: LG de los puntos que equidistan de un punto y una recta respectivamente. Luego la ecuaci´on cartesiana del LG de los puntos que equidistan de a dada por L1 y L2 est´ (9) d(P, Q1 ) = d (P, Q2 ) A partir de la ecuaci´on (9) podemos encontrar una representaci´on param´etrica de esta superficie, teniendo como par´ametros las variables que se usan para describir las rectas. De (7) y (8) se tiene las ecuaciones (P
→
→
− P ). u = s u
2
→
y (P
→
− P ). v = t v
1
2
2
que junto con la ecuaci´on (9) se obtiene el sistema de tres ecuaciones polinomiales en las variables x , y , z , s y t →
→
→
→
(P − P ). u = s u ((P d( P,− P Q ))). v == t(d(P,v Q )) 2
1
2
2
1
2
2
(10)
2
La soluci´ on del sistema (10), para las variables x , y , z en t´erminos de s y t,
x = x(s, t) = y (s, t), yz = z (s, t)
s, t,
∈ R
constituye una representaci´on param´etrica del LG.
9
(11)
El sistema (11) se obtiene construyendo la base de Gr¨obner del ideal generado por los polinomios (P
→
→
2
→
→
2
− P ). u −s u , (P − P ). v − t v , 1
en el anillo de polinomios GroebnerBasis.
2
R[x , y , z , s , t].
(d(P, Q1 ))2
2
− (d(P, Q )) 2
Para ello usamos la funci´on de Mathematica
Usaremos la ecuaci´on (9) para definir la funci´on LGDistaRectaRectaParametriza[Paso1 , V ector1 ,Paso2 , V ector2 , a , b ]
(ver programa 5) con la finalidad de obtener la ecuaci´on cartesiana, una representaci´on param´etrica y la gr´afica del LG. Los par´ametros de esta funci´on son: el punto de paso y el vector direcci´on de la recta L1 (Paso1, Vector1), el punto de paso y el vector direcci´ on de la recta L 2 (Paso2, Vector2) y los extremos del intervalo de la variable para construir la recta (a y b, respectivamente). Para verificar la buena definici´on de la funci´on LGDistaRectaRectaParametriza ingrese la expresi´on
{
}{
}{
}{
}
LGDistaRectaRectaParametriza[ 1,1,1 , 0,1,0 , -1,0,-1 , 1,0,-2 ,-10,10]
el resultado se puede observar en la figura 6
Figura 6: LG de los puntos que equidistan de dos rectas alabeadas 10
Siguiendo el proceso descrito, una representaci´on param´ etrica polinomial para la superficie de la figura 6 es 1 6
x = (1 + 10s + 5s − 2t − t ) y = 1 + t, (−5 − 20s + 5s − 2t − t ) z = 1 12
2
2
2
s, t
2
∈ R
LGDistaDosRectasParametriza Paso1_, Vector1_, Paso2_, Vector2_ :
P1 Paso1, V1 Vector1, P2 Paso2, V2 Vector2, Clearx, y, z; P x, y, z; Punto genérico del LG t_ : P1 t V1; Recta L1 t_ : P2 t V2; Recta L2 P P1.V1 V1; QenL2 SimplifyP2 P P2.V2 V2; QenL1 SimplifyP1 Norm V12 Norm V22 definicion de los polinomios f1 P P1. V1 s Norm V12 ; f2 P P2. V2 t Norm V22 ; f3 Expand Norm P QenL12 Norm P QenL22 . Abs Identity; Gb GroebnerBasisf1, f2, f3, x, y, z; RepParamet FlattenSolveGb1 0, Gb2 0, Gb3 0 , x, y, z; GRectas1y2 ParametricPlot3Dt, t, Module
t, 10, 10, AspectRatio 1 , PlotRange All, All, 10, 10, 1 PlotStyle Red, Thick, Blue, Thick, AxesLabel "X", "Y", "Z"; GraficaLG ContourPlot3Df3 0, x, 10, 10, y, 10, 10, z, 10, 10, PlotRange All, Mesh 4, 4, AxesLabel "X", "Y", "Z"; Graficas ShowGraficaLG, GRectas1y2; "t" t SlideView"Rectas: ", , "s" s "Ecuacion Cartesiana: ", f3 0, "Polinomios: ", f1, f2, f3, "Groebner Basis", Gb, "Parametrizacion: ", FlattenRepParamet, "Grafica de las rectas: ", GRectas1y2, "Grafica de la superficie: ", GraficaLG, "Rectas y superficie: ", Graficas, AppearanceElements All Programa 5
11
4.
Ejercicio Hacer las modificaciones apropiadas para obtener el LG de los puntos P (x , y , z ) 3 R tales que
∈
d(P, L1 ) =
k1 d(P, L2 ), k1 , k2 k2
∈ Z
+
, k1 = k 2
Referencias [1] Cox, David; Little, John; O’Shea, Donal, Using Algebraic Geometry , Springer, (2000) [2] Gonz´ alez U., Mariano, C´ alculo Integral en varias variables , (2013) [3] Pita, Claudio, C´ alculo vectorial , Prentice Hall, (1995) [4] Stewart, James, C´ alculo Trascendentes tempranas , Thomson Learning, (2002) [5] Wolfram Research Inc., Mathematica , v.8.0.4.0, (2011) [6] http://macareo.pucp.edu.pe/mgonzal/publicaciones.htm
12