PIC
Revista digital
Matemática, Educación e Internet
https://revistas.tec.ac.cr/index.php/matematica
Vol 25, No 1. Agosto, 2024 - Febrero, 2025
PIC
ISSN 1659 -0643
DOI: XXX

Sobre la parametrización de cónicas con rotación (la manera fácil)
| On parameterization of conics with rotation (the easy way) |

PIC Walter Mora-Flores
wmora2@gmail.com
Investigador independiente
Costa Rica
PIC Autor 2
autor2@itcr.ac.cr
Instituto Tecnológico de Costa Rica
Costa Rica
PIC Autor 3
autor3@itcr.ac.cr
Instituto Tecnológico de Costa Rica
Costa Rica

Recibido: 8 febrero de 2024

Aceptado: 10 mayo 2024

Contenido

Resumen:  Se presentan algunas parametrizaciones usuales (trigonométricas y racionales) de las cónicas en posición estándar y algunas aplicaciones.Para cónicas con rotación en la que solo interesa una parametrización trigonométrica o una racional, se presentan los métodos de diagonalización no ortogonal y parametrización usando un haz de rectas.

Palabras Clave:  Rotaciones, cónicas, parametrización.

Abstract: Some usual parameterizations (trigonometric and rational) of conics in standard position and some applications are presented. For conics with rotation in which only a trigonometric or a rational parameterization is of interest, the methods of non-orthogonal diagonalization and parameterization using a bundle of lines are presented.

Keywords:  Rotations, conics, parameterization.

1. Introducción

Una parametrización de una curva 𝐂 en n, es una función 𝐫 de un intervalo I en n, tal que el gráfico de 𝐫 está contenido en el gráfico de 𝐂. Usualmente en las aplicaciones se piden cosas adicionales como inyectividad, sobreyectividad, continuidad, derivabilidad, etc.

Una cónica en el plano 𝐱𝐲 es un conjunto 𝐂 = {(x,y)|𝐅(x,y) = 0} donde

𝐅(x,y) = Ax2 + 𝐵𝑥𝑦 + Cy2 + 𝐷𝑥 + 𝐸𝑦 + F (1)

donde no todos los coeficientes son cero. Para determinar si la ecuación (1) es la ecuación de una cónica "no degenerada" (es decir, si es una parábola, elipse o hipérbola), podemos usar el Teorema A (en el apéndice A).

Una una parametrización de una cónica en el plano 𝐱𝐲 es una función 𝐫 : I que satisface la ecuación 𝐅(x,y) = 0, es decir,

𝐅(x(t),y(t)) = 0 para todo t I

Una parametrización racional es una parametrización de la forma 𝐫(t) = (p1(t) q1(t),p2(t) q2(t)) donde pi,qi son polinomios. Un método para obtener parametrizaciones racionales de una cónica lo veremos en la sección 3..

En la tabla 1 se enumera algunas parametrizaciones usuales (trigonométricas y "racionales") de las cónicas propias, en posición estándar.

Tabla 1: Algunas parametrizaciones usuales de las cónicas.
Parábola Parametrización
(y k)2 = 4p(x h) 𝐫(t) = (pt2 + h,2𝑝𝑡 + k), t
(x h)2 = 4p(y k) 𝐫(t) = (2𝑝𝑡 + h,pt2 + k), t
Elipse
(x h)2 a2 + (y k)2 b2 = 1, 𝐫(t) = (h + cost,k + bsent) t [0,2π[
𝐫(t) = (a1 t2 1 + t2 + h,b 2t 1 + t2 + k), t
   
lim t±𝐫(t) = (a,0)
Hipérbola
(x h)2 a2 (y k)2 b2 = 1, 𝐫(t) = (h ±acosht,k + bsenht) t
𝐫(t) = (±at2 + 1 2t + h,bt2 1 2t + k), t > 0
𝐫(t) = (h + asect,k ±btant), t [0,2π] {π /2,3π /2}
(y k)2 b2 (x h)2 a2 = 1, 𝐫(t) = (h + asenht,k ±bcosht) t
𝐫(t) = (±at2 1 2t + h,bt2 + 1 2t + k), t > 0
𝐫(t) = (h ±atant,k + bsect), t [0,2π] {π /2,3π /2}

La parametrización de la hipérbola usando las funciones hipérbolicas (ver apéndice B), no presentan singularidades y gozan de simetría. En este caso t no es un ángulo, más bien (en el caso de una hipébola centrada en el origen) es dos veces el área (orientada) de la región entre la hipérbola, el eje focal y un rayo del origen al punto P = 𝐫(t) (ver Figura 1).

PIC

Figura 1: Parametrización de la hipérbola x2 a2 y2 b2 = 1 con 𝐫(t) = (±acosht,bsenh t)

Las parametrizaciones racionales se usan junto con las parametrizaciones trigonométricas en cambios de variables. También se usan en teoría de números, aplicaciones computacionales, etc. En el ejemplo 1. se muestra una aplicación: Una parametrización de una curva de intersección entre un cono y un plano.

Ejemplo 1 (Parametrización de una curva de intersección).

Consideremos el problema de determinar una parametrización trigonométrica de la curva de intersección entre el cono 𝐒1 : x2 + y2 = z2 y el plano 𝐒2 : 2y z = 2 (ver figura 2).

PIC

Figura 2: 𝐒1 𝐒2.

Una manera es parametrizar primero la proyección de la curva en plano 𝐱𝐲 y agregar después la coordenada z(t) (ver Figura 3).

PIC

Figura 3: Proyección de la curva de intersección.

La proyección de la curva de intersección en el plano 𝐱𝐲 es una hipérbola. Su ecuación se obtiene sustituyendo 2y z = 2 en x2 + y2 = z2, se obtiene

x2 + y2 = (2 2y)2

y completando cuadrados, armamos la ecuación canónica de esta hipérbola,

(y 4 /3)2 4 /9 x2 4 /3 = 1

Una parametrización de la proyección de la curva de intersección en el plano 𝐱𝐲 es

𝐜(t) = (h + asenh(t),k ±bcosh(t),0) = ( 2 3senh(t),4 3 ±2 3cosh(t),0)

Finalmente, como z(t) = 2 2y(t), una parametrización de la curva de intersección es

𝐫(t) = (h + asenh(t),k ±bcosh(t),z(t)) = ( 2 3senh(t),4 3 ±2 3cosh(t),2 2 (4 3 ±2 3cosh(t)))

El código en Mathematica es

r1[t_] := {h + a Sinh[t], k + b Cosh[t], 2 - 2 (k + b Cosh[t])}; 
r2[t_] := {h + a Sinh[t], k - b Cosh[t], 2 - 2 (k - b Cosh[t])}; 
Graphics3D[{Cone[{{0, 0, -4}, {0, 0, 0}}, 4], Cone[{{0, 0, 4}, {0, 0, 0}}, 4], 
(*Plano*) 
First@ParametricPlot3D[{0, 0, 2}+t {1, 0, 0}+s {0, -1, 2}, {t, -5, 5}, {s, -3, 2}, 
    Mesh -> None, PlotStyle -> {Gray, Opacity[0.4]}], 
(*Curva de intersecci\’on*) 
First@ParametricPlot3D[{r1[t], r2[t]}, {t, -2, 2}, Mesh -> None, 
    PlotStyle ->Directive[{Black, AbsoluteThickness[3], 
                          Black, AbsoluteThickness[3]}]] 
}, Boxed -> False, PlotRange -> All, 
 ViewPoint -> {3, 0.98, 1.18}, ViewVertical -> {0.25, 0.096, 2.032}, 
 ImageSize -> 400 , Method -> {"ShrinkWrap" -> True}]

Las parametrizaciones racionales tienen varias aplicaciones. En particular, podemos establecer si ciertas integrales se pueden calcular en términos de funciones elementales (las usuales en cálculo) y también podemos obtener una forma de las ternas pitágoricas. Esto se puede observar en los ejemplos 1. y 1..

Ejemplo 2 (¿Cuáles primitivas se pueden calcular?).

Si 𝐑(t) es una función racional (un cociente de polinomios) entonces, siempre que se pueda factorizar y resolver las ecuaciones lineales involucradas, se puede calcular la integral

𝐑(t)𝑑𝑡

en términos de funciones elementales.

Una función y = y(x) es algebraicamente dependiente de x si existe un polinomio en dos variables 𝐅 tal que 𝐅(x,y(x)) = 0 (como es el caso de las cónicas). La pregunta ahora es,

¿se puede calcular𝐑(x,y(x))𝑑𝑥?

La respuesta es: A veces, y no depende de una buena adivinanza!. Depende de si podemos parametrizar racionalmente la curva 𝐂.
Si una curva 𝐂, de ecuación 𝐅(x,y(x)) = 0, se puede parametrizar racionalmente por 𝐫(t) = (x1(t),x2(t)) entonces

𝐑(x,y(x))𝑑𝑡 =𝐑(x1(t),x2(t))x2(t)𝑑𝑡,

con lo cual la integral se puede calcular (porque el integrando es racional).

Casos especiales:

1.
𝐑(x,ax2 + 𝑏𝑥 + c)𝑑𝑥 se puede calcular (con la "sustitución de Euler").

En este caso,y(x) = ax2 + 𝑏𝑥 + c y 𝐅(x,y) = y2 ax2 𝑏𝑥 c = 0 es la ecuación de una cónica, y por tanto puede ser parametrizada racionalmente.

2.
𝐑(cos𝜃,sen𝜃)𝑑𝑥 se puede calcular.

Usamos una parametrización de la elipse x2 a2 + y2 b2 = 1.,

cos𝜃 = 1 t2 1 + t2,sen𝜃 = 2t 1 + t2,𝑑𝜃 = 2𝑑𝑡 1 + t2
3.
𝐑(cosh𝜃,senh𝜃)𝑑𝑥 se puede calcular.

Usamos una parametrización de la hipérbola x2 a2 y2 b2 = 1. Sustituimos cosh𝜃 y senh𝜃 por

x(t) = 1 + t2 1 t2,y(t) = 2t 1 t2,𝑑𝜃 = 2𝑑𝑡 1 t2
4.
¿Puede 1 + x3𝑑𝑥 se calculada (en términos funciones elementales)?. No. Se puede probar que la curva de ecuación y2 x3 1 = 0 no se puede parametrizar racionalmente, así que no podemos usar las técnicas que hemos mencionado!. Estas integrales (integrando con una raíz cuadrada de un polinomio cúbico) son "integrales elípticas" y no son expresables en términos de funciones elementales.

En general, las curvas de ecuación 𝐅(x,y) = 0 con grado 3, no pueden ser parametrizadas racionalmente (excepto algunos casos especiales). Y esto llevó al nacimiento de la Topología y de la (posiblemente complicada) Geometría Algebraica [5], [4].

Ejemplo 3 (Ternas Pitágoricas).

Una parametrización racional de la circunferencia x2 + y2 = 1 es

x(t) = 1 t2 1 + t2,y(t) = 2t 1 + t2

Las componentes racionales de la circunferencia se obtienen sustituyendo t = u v con u,v enteros. Sustituyendo en la parametrización racional, tenemos

x = v2 u2 u2 + v2,y = 2𝑢𝑣 u2 + v2

Y como x2 + y2 = 1, obtenemos una suma de cuadrados de enteros:

(v2 u2)2 + (2𝑢𝑣)2 = (u2 + v2)2

Las ternas Pitagóricas son enteros x,y,z + tal que x2 + y2 = z2 por lo que x = v2 u2, y = 2𝑢𝑣 y z = u2 + v2 son ternas Pitagóricas. En realidad todas las ternas Pitágoricas son de esta forma [1].

2. Parametrización de cónicas con rotación con diagonalización no ortogonal.

Los problemas de cónicas con rotación aparecen de manera natural en muchos contextos. Y si solo interesa una parametrización trigonométrica, podemos usar diagonalización no ortogonal para parametrizar "rapidamente".

Considere el problema de determinar una paremetrización trigonométrica de la curva de intersección entre:

1.
el cono 𝐒1 : x2 + y2 = z2
2.
y el plano 𝐒2 :4z + y + 2x = 2.

Como antes, parametrizamos la proyección de la curva con (x(t),y(t),0) y como z = 2 y 2x 4 entonces:

z(t) = 2 y(t) 2x(t) 4 .

La proyección de la curva en el plano 𝐱𝐲 es una curva 𝐂𝑥𝑦 (con rotación!) de ecuación

3x2 4 𝑥𝑦 4 + x 2 + 15y2 16 + y 4 1 4 = 0

Según el Teorema A (en el apéndice A), esta curva corresponde a un elipse (ver figura 4).

PIC

(a) 𝐒1 𝐒2

PIC

(b) 𝐒1 𝐒2 : 𝐫(t) = (x(t),y(t),z(t))
Figura 4: La proyección de la curva en el plano 𝐱𝐲 presenta rotación.

Necesitamos un cambio de variable que nos permite parametrizar la cónica 𝐂𝑥𝑦 usando una ecuación canónica (esto lo hacemos en el Ejemplo 2.).

Los cambios de variable en realidad son transformaciones. En este caso necesitamos transformaciones afines que envían cónicas con rotación en cónicas de la misma clase, en posición estándar, aunque no necesariamente congruentes.

La forma matricial de una cónica es

Ax2+𝐵𝑥𝑦+Cy2+𝐷𝑥+𝐸𝑦+F = [ x y ] [ A B2 B2 C ] 𝐀 [ x y ]+𝐁 [ x y ]+F = 0 con 𝐁 = [ D E ]

Necesitamos eliminar el producto "cruzado" 𝐵𝑥𝑦 en la ecuación. Este término desaparece si "eliminamos", con un cambio de variable, la diagonal secundaria ( B2 B2 ), es decir, si diagonalizamos 𝐀. Este tipo de transformaciones envían cónicas en cónicas de la misma clase, aunque no necesariamente congruentes.

Es decir, debemos encontrar un cambio de variable adecuado [ x y ] = 𝐂 [ u v ]de tal manera que:

[ x y ] [ A B2 B2 C ] 𝐀 [ x y ] + 𝐁 [ x y ] + F = 0 [ u v ]𝐂T [ A B2 B2 C ]𝐂 [ u v ] + 𝐁𝐂 [ u v ] + F = 0 [ u v ] [ A 0 0 C ] [ u v ] + 𝐁𝐂 [ u v ] + F = 0 Au2 + Cv2 No hay producto cruzado 𝑢𝑣 + 𝐁𝐂 [ u v ] + F = 0.

La buena noticia es que la matriz 𝐀, asociada a la cónica, es simétrica, por tanto podemos diagonalizarla (Ver teorema 2.).

Teorema 1.

Si 𝐀n×n es simétrica, entonces existe una matriz invertible 𝐂 tal queda 𝐂𝐓𝐀𝐂 = 𝐃 donde 𝐃 es una matriz diagonal.

1.
La matriz 𝐂 no es única, una de estas matrices 𝐂 se puede obtener usando usando completación de cuadrados o también por operaciones elementales "acopladas" sobre 𝐀, [3]. En este caso podemos obtener una parametrización de la cónica en el plano 𝐱𝐲 pero no mucha información adicional.
2.
También 𝐂 se puede obtener como una matriz ortogonal (unitaria), es decir, una matriz cuyas columnas son vectores propios de 𝐀. Este último caso corresponde a la "diagonalización ortogonal" de 𝐀 y permite parametrizar la cónica en el plano 𝐱𝐲 y además recuperar, con el cambio de variable, foco(s), vértice(s), asíntotas, etc.

A continuación, se presenta el ejemplo 2. en donde se muestra la diagonalización no ortogonal.

Ejemplo 4 ((Diagonalización no ortogonal).).

Considere la cónica x2 + 𝑥𝑦 + y2 = 1.

1.
Forma matricial [ x y ] [ 1 12 12 1 ] [ x y ]1 = 0
2.
Un cambio de variable adecuado puede ser [ x y ] = [ 1 1 1 1 ] 𝐂 [ u v ], pues
[ u v ]𝐂T [ 1 12 12 1 ]𝐂 [ u v ]1 = [ u v ] [ 3 0 0 1 ] [ u v ]1 = 3u2 + v21 = 0
Parametrización en el plano 𝐱𝐲

Al aplicar un cambio de variable adecuado (x,y) = 𝐂2×2 [ u v ], podemos obtener la ecuación canónica de la cónica en el plano 𝐮𝐯 y por tanto una parametrización 𝐜(t) = (u(t),v(t)) en este plano. Entonces una parametrización en el sistema 𝐱𝐲 es

𝐫(t) = (x(t),y(t)) = 𝐂[ u(t) v(t) ]

o también de manera más práctica, si 𝐰1,𝐰2 son las columnas de 𝐂, es decir, 𝐂 = ( 𝐰1 𝐰2 ) , entonces 𝐡(t) = u(t)𝐰1 + v(t)𝐰2.

En la Figura 5 podemos ver el proceso del cambio de variable en el caso de una elipse.

PIC

(a) Cónica con rotación

   

PIC

(b) (x.y) = 𝐂 [ u v ]

   

PIC

(c) (uh)2 a2 + (vk)2 b2 = 1

Figura 5: Cambio de variable .
¿Como escoger el cambio de variable?.

El cambio de variable depende del método de diagonalización que escojamos y depende de los coeficientes A,B y C de la cónica. Estos cambios de variable no siempre son excluyentes.

Cambio de variable con diagonalización por completación de cuadrados.

Supongamos que la ecuación de la cónica propia, con rotación, es

Ax2 + 𝐵𝑥𝑦 + Cy2 + 𝐷𝑥 + 𝐸𝑦 + F = 0 (2)

Podemos completar cuadrados de varias maneras

1.
Si A = C,
Ax2 + 𝐵𝑥𝑦 + Ay2 + Dx + Ey + F = (2A B) (x y 2 )2 + (2A + B) (x + y 2 )2 + Dx + Ey + F
2.
Si A0,
Ax2 + 𝐵𝑥𝑦 + Cy2 + 𝐷𝑥 + 𝐸𝑦 + F = A (x + B 2Ay)2 + (C B2 4A)y2 + 𝐷𝑥 + 𝐸𝑦 + F = 0
3.
Si C0,
Ax2 + 2𝐵𝑥𝑦 + Cy2 + 𝐷𝑥 + 𝐸𝑦 + F = x2A B2 4C + C𝐵𝑥 2C + y2 + 𝐷𝑥 + 𝐸𝑦 + F = 0

Ahora, deducimos un cambio de variable [ x y ] = 𝐂 [ u v ] que diagonaliza la matriz asociada de la cónica: Nos queda una cónica del mismo tipo, pero no necesariamente congruente. Pero es suficiente para parametrizar la cónica original.

1.)
Si A = C, la completación de cuadrados sugiere el cambio de variable u = x y 2 y v = x + y 2 , es decir, x = u + v y = u v.

Si A = C, el cambio de variable es (x,y) = 𝐂 [ u v ] con 𝐂 = ( 1 1 1 1 )

Ax2 + 𝐵𝑥𝑦 + Ay2 + 𝐷𝑥 + 𝐸𝑦 + F = 0 se transforma en (2A + B)u2 + (2A B)v2 + D(u + v) + E(u v) + F = 0

Ahora se puede verificar (en el caso de una cónica propia), que 4(2A + B)(2A B) tiene el mismo signo que B2 4𝐴𝐶 o se anulan simúltaneamente. Es decir, son cónicas del mismo tipo.

Una representación de este caso se puede observar en el ejemplo ??.

Ejemplo 5.

Determine una parametrización trigonométrica para la cónica x2 + 𝑥𝑦 + y2 1 = 0.

Solución: Como A = 1 = C y B = 1, podemos usar (x,y) = [ 1 1 1 1 ] [ u v ]

PIC

Cónica

[ x y ] = [ 1 1 1 1 ] [ u v ]

PIC

u2 1 /3 + v2 = 1

Figura 6: Cambio de variable (x,y) = 𝐂 [ u v ]
(a)
Parametrización en el plano 𝐮𝐯. Haciendo la sustitución nos queda

(2A + B)u2 + (2A B)v2 + D(u + v) + E(u v) + F = 0 3u2 + v2 1 = 0 u2 1 /3 + v2 = 1 𝐫𝑢𝑣(t) = (1 / 3 cos(t),sen(t)),t [0,2π]

(b)
Parametrización en el plano 𝐱𝐲
𝐫𝑥𝑦(t) = [ 1 1 1 1 ] [ 1 / 3 cos(t) sen(t) ] = (sen(t) + cos(t) 3 ,cos(t) 3 sen(t)),t [0,2π]

2.)
Si A0, la completación de cuadrados sugiere la sustitución u = x + B 2Ay y y = v, es decir, x = u B 2Av y y = v.

Si A0, el cambio de variable es (x,y) = 𝐂 [ u v ] con 𝐂 = ( 1 B /2A 0 1 )

Ax2 + 𝐵𝑥𝑦 + Ay2 + 𝐷𝑥 + 𝐸𝑦 + F = 0

se transforma en:

Au2 + v2 (C B2 4A) + 𝐷𝑢 + v (E 𝐵𝐷 2A ) + F = 0.

Ahora se puede verificar (en el caso de una cónica propia), que 4A (C B2 4A) tiene el mismo signo que B2 4𝐴𝐶 o se anulan simúltaneamente. Es decir, son cónicas del mismo tipo.

Una representación de este caso se puede observar en el ejemplo 2..

Ejemplo 6.

Determine una parametrización trigonométrica para la cónica x2 4 2𝑥𝑦 + x y 1 = 0.

solu Como A = 14 y B = 20, podemos usar (x,y) = [ 1 4 0 1 ] [ u v ]

PIC

Cónica 𝐂

[ x y ] = [ 1 4 0 1 ] [ u v ]

PIC

(u h)2 a2 (v k)2 b2 = 1

Figura 7: Cambio de variable (x.y) = 𝐂 [ u v ]
(a)
Parametrización en el plano 𝐮𝐯. Haciendo la sustitución x = u + 4v y y = v nos queda

Au2 + v2 (C B2 4A) + 𝐷𝑢 + v (E 𝐵𝐷 2A ) + F = 0 u2 4 4v2 + u + 3v 1 = 0 (u + 2)2 23 /4 (v 3 /8) 2 23 /64 = 1 { 𝐫𝑢𝑣1(t) = (2 + 23 4 cosh(t),3 8 + 23 / 64senh(t)) 𝐫𝑢𝑣2(t) = (2 23 4 cosh(t),3 8 + 23 / 64senh(t)) ,t

(b)
Parametrización en el plano 𝐱𝐲
𝐫𝑥𝑦(t) = [ 1 4 0 1 ] [ 2 ±23 4 cosh(t) 3 /8 + 23 64senh(t) ] = (2 + 23 2 senh(t) + 3 /2 ±23 /2 cosh(t),3 8 + 23 8 senh(t)),t
3.)
Si C0 la completación de cuadrados sugiere la sustitución u = x y v = y + B 2C x, es decir, x = u y y = v B 2C u

El cambio de variable es (x,y) = 𝐂 [ u v ] con 𝐂 = [ 1 0 B 2C 1 ]

Ax2 + 𝐵𝑥𝑦 + Ay2 + 𝐷𝑥 + 𝐸𝑦 + F = 0 se transforma en (A B2 4C)u2 + Cv2 + u (D 𝐵𝐸 2C ) + 𝐸𝑣 + F = 0

Ahora se puede verificar (en el caso de una cónica propia), que 4 (A B2 4C ) C tiene el mismo signo que B2 4𝐴𝐶 o se anulan simúltaneamente. Es decir, son cónicas del mismo tipo.

Una representación de este caso se puede observar en el ejemplo 2..

Ejemplo 7.

Determine una parametrización trigonométrica para la cónica de ecuación

x2 + y2 + 2𝑥𝑦 + 2x 2y = 0

Solución: Aquí podemos usar el cambio de variable 1.) Pero lo haremos con el cambio de variable 3.). Como C = 1 y B = 2 podemos usar (x,y) = [ 1 0 1 0 ] [ u v ]

PIC

Cónica

[ x y ] = [ 1 0 1 0 ] [ u v ]

PIC

22u + v2 2v = 0

Figura 8: Cambio de variable (x,y) = 𝐂 [ u v ]
1.
Parametrización en el plano 𝐮𝐯. Haciendo la sustitución x = u y y = v u, en el plano 𝐮𝐯 nos queda

(A B2 4C)u2 + Cv2 + u (D 𝐵𝐸 2C ) + 𝐸𝑣 + F = 0 v2 + 22u 2v = 0 despejamos u y tomamos t = v 𝐫𝑢𝑣(t) = (2t t2 22 ,t),t

2.
Parametrización en el plano 𝐱𝐲
𝐫𝑥𝑦(t) = [ 1 0 1 0 ] [ 2t t2 22 t ] = (2t t2 22 ,t 2t t2 22 ),t

En el ejemplo 2. se presenta una forma para determinar una parametrización trigonométrica de la curva de intersección entre un cono y un plano. En este caso, la proyección de la curva de intersección presenta rotación, por lo que, como el propósito es solo parametrizar, podemos usar diagonalización no ortogonal.

Ejemplo 8.

Determinar una parametrización trigonométrica de la curva de intersección entre el cono 𝐒1 : x2 + y2 = z2 y el plano 𝐒2 :4z + y + 2x = 2.

Solución: Como antes, parametrizamos la proyección de la curva con (x(t),y(t),0) y como z = 2 y 2x 4 entonces z(t) = 2 y(t) 2x(t) 4

La proyección de la curva en el plano 𝐱𝐲 es una curva 𝐂𝑥𝑦 con ecuación

3x2 4 𝑥𝑦 4 + x 2 + 15y2 16 + y 4 1 4 = 0

Según el Teorema A (en el apéndice A), esta curva corresponde a un elipse (ver figura 4).

Bien, como A = 3 40 entonces podemos usar 𝐂 = [ 1 1 6 0 1 ]

Aplicando el cambio de variable (x,y) = 𝐂 [ u v ] nos queda

Au2 + v2 (C B2 4A) + 𝐷𝑢 + v (E 𝐵𝐷 2A ) + F = 0 3u2 4 + 11v2 12 + u 2 + v 3 1 4 = 0 33 16 (u + 1 3) 2 + 121 48 (v + 2 11) 2 = 1

1.
Parametrización en el plano 𝐮𝐯. Como h = 1 3,k = 2 11,a = 16 33, y b = 48 121, entonces
𝐫𝐮𝐯(t) = (h + asent,k + bcost)
2.
Parametrización en el plano 𝐱𝐲. La parametrización de la proyección de la curva de intersección es
(x(t),y(t)) = 𝐂 [ h + asent k + bcost ]
3.
La parametrización de la curva de intersección escoger
𝐫(t) = (x(t),y(t),1 4 (2 y(t) 2x(t)))

El código en Mathematica es

ClearAll[x, y, u, v, h, k, a, b, mC, h, k, a, b, rxy, r, cv, n,e1,e2] 
mC = {{1, 1/6}, {0, 1}}; 
h = -1/3; k = -2/11; a = Sqrt[16/33]; b = Sqrt[48/121]; 
u[t_] := h + a Cos[t]; v[t_] := k + b Sin[t]; 
(*Parametrizaciones*) 
cv[t_] := mC.{u[t], v[t]}; 
x[t_] := cv[t][[1]]; y[t_] := cv[t][[2]]; 
rxy[t_] := {x[t], y[t], 0}; 
r[t_] := {x[t], y[t], 1/4 (2 - y[t] - 2 x[t])}; 
(*Una base ortonormal del plano 4z+y+2x=2*) 
n = {2, 1, 4}; 
{e1, e2} = Orthogonalize@{{0, n[[3]], -n[[2]]}, {-n[[3]], 0, n[[1]]}}; 
 
Graphics3D[{ 
  (*Proyecci\’on de la curva*) 
  First@ParametricPlot3D[rxy[t], {t, 0, 2 Pi}, Mesh -> None, 
                 PlotStyle -> {Blue, AbsoluteThickness[2]}], 
  (*Curva de intersecci\’on*) 
  First@ParametricPlot3D[r[t], {t, 0, 2 Pi}, Mesh -> None, 
                 PlotStyle -> {Brown, AbsoluteThickness[2]}], 
  (*Cono*) 
  First@ContourPlot3D[x^2 + y^2 == z^2, {x, -2, 1.5}, {y, -2, 2}, {z, 0, 1.5}, 
              PlotPoints -> 40, Mesh -> None, 
              ContourStyle -> 
              Directive[Orange, Opacity[0.2], Specularity[White, 30]], 
              PerformanceGoal -> "Quality"], 
  (*Plano*) 
  First@ParametricPlot3D[{0,0,1/2}+t*e1+s*e2,{t,-2, 2}, {s,-1,2}, 
                 Mesh -> None, PlotStyle -> {Gray, Opacity[0.2]}] 
  }, Boxed -> False, PlotRange -> All, ImageSize -> 400, 
    Method -> {"ShrinkWrap" -> True}]

3. Parametrización racional de una cónica usando un "haz de rectas".

Una "parametrización racional" de una curva es una parametrización en terminos de cocientes de polinomios. Estas parametrizaciones son muy útiles en muchos contextos, como el de gráficos por computadora, modelado geométrico, ecuaciones diofánticas, cambios de variable en integración, etc.
Las cónicas tienen la propiedad de que cada par de puntos en ella, se pueden obtener por intersección con una línea recta. La idea básica consiste en utilizar un "lápiz o haz de líneas rectas" que pasan por un punto (a,b) de la curva, de manera que al calcular el otro punto de intersección de cada recta genérica del "haz", con la curva, se determina una parametrización de la curva (excepto por uno o dos puntos). Cada línea del "haz" depende de un parámetro t que será el parámetro de la parametrización de la curva (ver figura 9).

PIC

Figura 9: Parametrización de una cónica con un "haz de rectas".

La idea general es como sigue:

1.
Tenemos una cónica 𝐂 (no degenerada) de ecuación Ax2 + 𝐵𝑥𝑦 + Cy2 + 𝐷𝑥 + 𝐸𝑦 + F = 0 y (a,b) 𝐂
2.
El haz de rectas que pasan por (a,b) tienen ecuación general y = t(x a) + b donde t es un parámetro que hace variar las rectas.
3.
Sustituimos en la ecuación de la cónica Ax2 + 𝐵𝑥𝑦 + Cy2 + 𝐷𝑥 + 𝐸𝑦 + F = 0 y queda

Ax2 + 𝐵𝑥 (t(x a) + b) + C (t(x a) + b)2 + 𝐷𝑥 + E (t(x a) + b) + F = 0.

Es decir, recolectando, queda una ecuación cuadrática α(t)x2 + β(t)x + γ(t) = 0.

4.
La ecuación α(t)x2 + β(t)x + γ(t) = 0 tiene a lo sumo dos soluciones. Una de las soluciones es x = a pues el punto (a,b) satisface la ecuación de la cónica. Entonces el resultado de la división de la cuadrática por (x a) tiene que ser un polinomio en x de grado 1, es decir, q(t)x p(t)

α(t)x2 + β(t)x + γ(t) x a q(t)x p(t) Residuo:  0

α(t)x2 + β(t)x + γ(t) = (x a)(q(t)x p(t)) = 0

Por tanto las soluciones son x = a y x = p(t) q(t) donde p(t) y q(t) son polinomios en t de grado a lo sumo 2.

5.
Sustituimos x(t) = p(t) q(t) en y = t(x a) + b y obtenemos
𝐂 : 𝐫(t) = (p(t) q(t),p1(t) q1(t)),t {ceros de q(t) y q1(t)

Observe que si tomamos valores racionales de t, podemos obtener puntos con coordenadas racionales en la cónica.

En el ejemplo 3. se presenta una parametrización racional de una cónica usando un "haz de rectas".

Ejemplo 9.

Determine una parametrización racional de la cónica 𝐂 : 𝑥𝑦 + y + 2 = 0 usando el punto (2,2) 𝐂 (ver figura 10).

PIC

Figura 10: Parametrización racional de 𝐂

Solución:

1.
Para construir el haz de rectas podemos escoger las rectas que pasan por (2,2) y (0,t) por ejemplo. Estas rectas tiene ecuación y(t) = t2 2 (x + 2) + 2
2.
Sustituimos y(t) en 𝑥𝑦 + y + 2 = 0 : x (t2 2 (x + 2) + 2) + t2 2 (x + 2) + 2 + 2 = 0
3.
Polinomio de grado 2 en x: (t 2 1)x2 + (3t 2 1)x + t + 2 = 0
4.
Este polinomio tiene dos raíces, una es x = 2. Para obtener la otra raíz lo que hacemos es dividir este polinomio por x + 2 para obtener el otro factor (el residuo es cero) y luego despejamos x

(t 2 1)x2 + (3t 2 1)x + t + 2 x + 2 (t 2 1)x + t 2 + 1 Residuo:  0

5.
( t 2 1)x + t 2 + 1 = 0x(t) = t 2 t 2
6.
Sustituimos x(t) en y(t) = t2 2 (x + 2) + 2 y obtenemos la parametrización racional,

𝐂 :𝐫(t) = (t + 2 2 t,t 2 2 ),t {2}

Con este método es que se obtienen las parametrizaciones racionales de la Tabla 1.

Hay algoritmos más sofisticados para obtener parametrizaciones racionales de ciertas curvas, pero no todas las curvas tienen una parametrización racional [6], [5].

4. Bibliografía

[1]   David Angell. "Beginning algebraic number theory". https://web.maths.unsw.edu.au/~angell/articles/algnumtheory.pdf

[2]   D. Marsh. (2005) “Applied Geometry for Computer Graphics and CAD”. 2nd ed. Springer.

[3]   Sergei Treil. "Linear Algebra Done Wrong". https://www.math.brown.edu/streil/papers/LADW/LADW-2014-09.pdf. 2022.

[4]   Rick Miranda. (1983). "Integration: why you can and why you can’t". Pi Mu Epsilon Journal, Vol. 7, No. 9 (fall 1983), pp. 557-566.
https://www.math.colostate.edu/~miranda/preprints/IntegrationWhyYouCan-Can%27t.pdf

[5]   Igor R. Shafarevich. (2013). "Basic Algebraic Geometry 1, Varieties in Projective Space". Springer, 2013 edition.

[6]   John Stillwell. (2008). "Naive Lie theory" Springer, Springer; 2008a edition.

Apéndice

A Clasificación de cónicas

Hay casos en los que esta ecuación general (1) no tiene soluciones (no hay lugar geométrico) o el conjunto solución es por ejemplo un punto o dos rectas. Estos casos especiales se llaman "cónicas degeneradas".

Usando la teoría de formas cuadráticas podemos obtener un criterio para clasificar las cónicas a partir de su ecuación general, [2, pp. 112-114]

Teorema 2 ().

Sea C la cónica de ecuación Ax2 + B𝑥𝑦 + Cy2 + Dx + Ey + F = 0.
Si Δ = 4𝐴𝐶𝐹 AE2 B2F + 𝐵𝐷𝐸 CD2, entonces

1.
si B2 4𝐴𝐶 = 0 y Δ0, C es una parábola.
2.
si B2 4𝐴𝐶 < 0 y Δ0, C es una elipse.
3.
si B2 4𝐴𝐶 > 0 y Δ0, C es una hipérbola.

B Acerca del seno hipérbolico y el coseno hipérbolico.

Identidades básicas.

   
senh(t) = et e t 2 cosh(t) + senh(t) = et
   
cosh(t) = et + e t 2 cosh(t) senh(t) = e t
   
tanh(t) = senh(t) cosh(t) cosh2(t) senh2(t) = 1
   
Inversas
   
arsenh(t) = ln (t + t2 + 1),t arcosh(t) = ln (t + t2 1), t 1
   

Crecimiento y signo

PIC

Figura 11: Representación de y = cosh(t).

PIC

Figura 12: Representación de y = sinh(t).
cosh(t) 1 para todo t . Además, cosh(t) si t 0 y cosh(t) si t 0
   
senh(t) es creciente, además senh(t) 0 si t 0y senh(t) 0 si t 0.
   

Ecuaciones. Como se observa en la Figura 11 y la Figura 12, se usan inversas para obtener la o las soluciones de una ecuación.

{ cosh(t) = a { t = acosh(a) t = acosh(a)  si a 1 senh(t) = at = asenh(a)

Simetría. Usando la parametrizaciones 𝐫1 y 𝐫2 para cada rama, 𝐫i(0) es el vértice respectivo y 𝐫i(t) y 𝐫i(t) son simétricos respecto al eje focal (ver Figura 13).

PIC

Figura 13: Simetría.