CM

Published on January 2018 | Categories: Documents | Downloads: 51 | Comments: 0 | Views: 905
of 113
Download PDF   Embed   Report

Comments

Content

Notas de Geometr´ıa Diferencial con aplicaciones Antonio Vald´es 29 de enero de 2014

2

´Indice general 1. Introducci´ on al Sage 1.1. Instalaci´on . . . . . . . . . . . . . . . . . . . 1.2. Primeros pasos con Sage . . . . . . . . . . . 1.2.1. Sage como una calculadora . . . . . . 1.2.2. Sage como int´erprete de Python . . . 1.2.3. C´alculo diferencial y manipulaci´on de 1.2.4. Resolviendo ecuaciones . . . . . . . . 1.2.5. C´alculo diferencial . . . . . . . . . . 1.2.6. Representaci´on gr´afica de funciones . ´ 1.2.7. Algebra lineal . . . . . . . . . . . . .

. . . . . . . . . . . . . . . . . . . . . . . . . . . . expresiones . . . . . . . . . . . . . . . . . . . . . . . . . . . .

. . . . . . . . .

. . . . . . . . .

. . . . . . . . .

5 5 6 6 7 8 10 12 13 14

2. Curvas parametrizadas 19 2.1. Ejemplos de curvas planas . . . . . . . . . . . . . . . . . . . . 20 2.1.1. Circunferencia . . . . . . . . . . . . . . . . . . . . . . . 22 2.1.2. Curvas polin´omicas . . . . . . . . . . . . . . . . . . . . 23 2.2. Curvas de B´ezier . . . . . . . . . . . . . . . . . . . . . . . . . 24 2.2.1. Afinidades . . . . . . . . . . . . . . . . . . . . . . . . . 24 2.2.2. Combinaciones baric´entricas . . . . . . . . . . . . . . . 26 2.2.3. Algoritmo de Jarvis para el c´alculo de la envoltura convexa 29 2.2.4. Polinomios de Bernstein . . . . . . . . . . . . . . . . . 30 2.2.5. Curvas de B´ezier . . . . . . . . . . . . . . . . . . . . . 31 2.2.6. El algoritmo de De Casteljau . . . . . . . . . . . . . . 33 2.3. Las c´onicas como curvas racionales . . . . . . . . . . . . . . . 34 3. Curvas regulares. Curvatura. 39 3.1. Curvas regulares. Longitud de arco. . . . . . . . . . . . . . . . 39 3.2. Definici´on de curvatura . . . . . . . . . . . . . . . . . . . . . . 49 3.2.1. Curvatura de una curva plana . . . . . . . . . . . . . . 50 3.2.2. Curvatura de una curva no parametrizada por la longitud de arco 52 3.2.3. Invariancia de la curvatura . . . . . . . . . . . . . . . . 54 3.2.4. Circunferencia osculadora . . . . . . . . . . . . . . . . 55 3

4

´INDICE GENERAL 3.3. Envolvente de una familia de curvas planas parametrizadas 3.4. Ecuaciones impl´ıcitas de curvas en el plano . . . . . . . . . 3.4.1. C´onicas degeneradas . . . . . . . . . . . . . . . . . 3.4.2. Curvatura de una curva definida impl´ıcitamente . .

4. Ecuaciones de Frenet 4.1. Otra interpretaci´on de la curvatura de una curva plana 4.2. Ecuaciones de Frenet planas . . . . . . . . . . . . . . . 4.3. El teorema fundamental de la teor´ıa de curvas planas . 4.4. Curvas espaciales. Plano osculador . . . . . . . . . . . 4.5. Binormal. F´ormulas de Frenet y consecuencias. . . . . . 4.6. Forma can´onica local de una curva . . . . . . . . . . .

. . . . . .

. . . . . .

. . . .

. . . .

58 61 65 67

. . . . . .

. . . . . .

69 69 71 73 76 78 80

5. Superficies 5.1. Superficies de B´ezier . . . . . . . . . . . . . . . . . . . . . . . 5.2. Superficies parametrizadas regulares . . . . . . . . . . . . . . . 5.2.1. El plano tangente es el que mejor aproxima la superficie 5.3. Ejemplos de superficies . . . . . . . . . . . . . . . . . . . . . . 5.3.1. Grafos . . . . . . . . . . . . . . . . . . . . . . . . . . . 5.3.2. Superficies de rotaci´on . . . . . . . . . . . . . . . . . . 5.4. Rotaciones y cuaterniones . . . . . . . . . . . . . . . . . . . . 5.5. Curvas sobre superficies. Geod´esicas . . . . . . . . . . . . . . . 5.5.1. Geod´esicas . . . . . . . . . . . . . . . . . . . . . . . . .

83 83 86 en el punto 86 88 88 88 89 92 93

6. Las formas fundamentales 97 6.1. Primera forma fundamental . . . . . . . . . . . . . . . . . . . 97 6.1.1. Longitudes de curvas . . . . . . . . . . . . . . . . . . . 97 ´ 6.1.2. Angulos . . . . . . . . . . . . . . . . . . . . . . . . . . 98 ´ 6.1.3. Areas . . . . . . . . . . . . . . . . . . . . . . . . . . . 99 6.2. La segunda forma fundamental . . . . . . . . . . . . . . . . . 99 6.2.1. Curvatura normal y curvatura geod´esica . . . . . . . . 102 6.2.2. Secciones normales . . . . . . . . . . . . . . . . . . . . 103 6.2.3. Diagonalizaci´on de la segunda forma fundamental . . . 104 6.2.4. C´alculos expl´ıcitos . . . . . . . . . . . . . . . . . . . . 105 6.3. L´ıneas de curvatura y l´ıneas asint´oticas . . . . . . . . . . . . . 107 6.4. Una interpretaci´on geom´etrica de la curvatura de Gauss . . . . 109 6.5. El teorema egregio de Gauss . . . . . . . . . . . . . . . . . . . 111 Esta obra se difunde bajo la licencia Creative Commons ReconocimientoCompartirIgual 3.0. Puede leer las condiciones de la licencia en http://creativecommons.org/licenses/by-sa/3.0/deed.es :

Cap´ıtulo 1 Introducci´ on al Sage Esta es una breve introducci´on al Sage en la que s´olo se tratar´an aspectos del mismo que son relevantes para nuestro curso. Sage es realmente inmenso, cosa que se puede comprobar hojeando los miles de p´aginas de los que consta el manual de referencia http://www.sagemath.org/doc/reference/, que es el instrumento esencial para encontrar la parte de Sage relevante para las matem´aticas que te interesen en cada momento. Tambi´en existe mucha informaci´on en la red, as´ı que a menudo podr´as encontrar lo que necesites haciendo una b´ usqueda en Google (en ingl´es, preferiblemente) a la que a˜ nadas el t´ermino “sagemath”.

1.1.

Instalaci´ on

Sage est´a ya instalado en los ordenadores del laboratorio. Cuando el ordenador arranca, hay que elegir “Ubuntu” para que arranque el sistema y, una vez estemos dentro, basta con hacer doble click en el icono de Sage que aparece en el escritorio. La forma m´as sencilla de utilizar Sage fuera del laboratorio es hacerse una cuenta en https://cloud.sagemath.com/. Es un servicio experimental gratuito para el cual es necesario un acceso a internet. Es una opci´on adecuada para hacer c´alculos sencillos, pero probablemente no lo sea para c´alculos de coste computacional alto. Si se quiere instalar Sage en tu propio ordenador, la mejor opci´on es instalar una distribuci´on de GNU/Linux. Si quieres instalar GNU/Linux se recomienda empezar instalando Ubuntu. Si el ordenador es muy antiguo, es posible instalar versiones m´as ligeras como por ejemplo Xubuntu. Las instrucciones para instalar Ubuntu se pueden encontrar en http://www.ubuntu.com/download/desktop/install-desktop-latest. 5

´ AL SAGE CAP´ITULO 1. INTRODUCCION

6

Una vez instalado GNU/Linux, se puede descargar Sage de http://www.sagemath.org/download-linux.html. La instalaci´on es muy sencilla: basta con descomprimir el fichero y ejecutar el fichero sage. Es posible tambi´en hacer un Live USB para arrancar Sage (instrucciones aqu´ı, sin embargo hay que saber que con el Live USB se tiene una versi´on de Sage algo m´as antigua aunque probablemente suficiente para nuestros prop´ositos), instalar Sage en MacOS o usar una m´aquina virtual para hacer que Sage corra en Windows. Instrucciones detalladas se pueden encontrar en http://www.sagemath.org/doc/installation/.

1.2.

Primeros pasos con Sage

Se recomienda encarecidamente hacer el tutorial de Sage, que est´a en cualquier instalaci´on local o en la URL http://www.sagemath.org/doc/tutorial/. Es importante saber que se puede obtener ayuda sobre cualquier funci´on, m´etodo o constante escribiendo su nombre seguido de un signo de interrogaci´on. Por ejemplo: sage: matrix? devolver´a la p´agina de ayuda del constructor de matrices. Si se tiene un objeto (por ejemplo, foo) y un m´etodo met, se puede obtener la ayuda correspondiente haciendo sage: foo.met? Finalmente, es muy u ´til usar el tabulador (quiz´a despu´es de haber escrito algunas letras) para completar las posibles funciones de las que sage dispone. Tambi´en, una vez definido un objeto, escribiendo foo. y pulsando el tabulador, veremos todos los m´etodos que son aplicables a foo. Sage se puede utilizar en la consola o en modo “Hoja de c´alculo”, usando un navegador. Para abrir el navegador basta con dar la instrucci´on notebook() desde Sage o arrancarlo como sage -n. Tambi´en es posible hacer scripts en Python que usen Sage y escribir programas en Sage. Aunque esta introducci´on est´a enfocada principalmente a la consola, en general se recomienda utilizar Sage en el navegador.

1.2.1.

Sage como una calculadora

Con Sage se pueden hacer operaciones como en una calculadora. sage: 1+1

1.2. PRIMEROS PASOS CON SAGE

7

2 sage: 3^20 3486784401 Sage procurar´a, en la medida de lo posible, mantenerse en aritm´etica exacta. Por ejemplo: sage: 5/3 + 2/5 31 15 Si deseamos soluciones aproximadas, podemos obtenerlas por ejemplo usando la funci´on numerical approximation, que se puede abreviar simplemente por n: sage: cos(pi/5) 1 1√ 5+ 4 4 sage: n(cos(pi/5), digits=2) 0,81

1.2.2.

Sage como int´ erprete de Python

Sage utiliza el lenguage Python, de forma que todo lo que se puede hacer en Python se puede hacer en Sage. Por ejemplo, calculemos la suma de los cuadrados de los naturales menores que 1000 que no son m´ ultiplos de 3 ni de 7: sage: sum(x**2 for x in range(1000) if x%3!=0 and x%7!=0) 190046194 Sin embargo, en Sage podemos escribir x^2 en lugar de x**2, si lo deseamos. En Sage se pueden escribir programas en Python, con el cuidado habitual de seguir las reglas de indentaci´on:

´ AL SAGE CAP´ITULO 1. INTRODUCCION

8

sage: def slowFib(n): ... if n == 0: return 0 ... elif n == 1: return 1 ... else: return slowFib(n-1) + slowFib(n-2) sage: slowFib(10)

55

1.2.3.

C´ alculo diferencial y manipulaci´ on de expresiones

En nuestra asignatura nos interesan especialmente las capacidades de c´alculo simb´olico del Sage. Para utilizar variables simb´olicas, es necesario definirlas previamente: sage: var(’x y’) (x, y) sage: A = 1 + 1 / (x + 1/y)^2 + 3*x*y / (1 + y^2); Podemos simplificarla: sage: A.simplify_full() (x2 + 1)y 4 + (3 x3 + 2 x)y 3 + (7 x2 + 2)y 2 + 5 xy + 1 x2 y 4 + 2 xy 3 + (x2 + 1)y 2 + 2 xy + 1 En ocasiones, es preferible factorizar para obtener expresiones m´as simples: sage: A.simplify_full().factor() 3 x3 y 3 + x2 y 4 + 7 x2 y 2 + 2 xy 3 + y 4 + 5 xy + 2 y 2 + 1 (xy + 1)2 (y 2 + 1) En muchas ocasiones, el m´etodo simplify va a ser menos adecuado para nosotros que simplify full, que aplica una sucesi´on de simplificaciones sobre la expresi´on. Para saber exactamente qu´e simplificaciones se hacen, puedes acceder a la ayuda haciendo

9

1.2. PRIMEROS PASOS CON SAGE sage: f.simplify_full?

Para evaluar la expresi´on en x = 2, y = 3 hay varias opciones, por ejemplo: sage: A(x=2, y=3) 731 245 Se puede usar tambi´en el m´etodo subs, que permite hacer sustituciones. Una de las posibilidades es que su argumento sea un diccionario. Para ilustrarlo, hagamos por ejemplo un cambio a coordenadas polares: sage: var(’t’) t sage: A.subs({x:cos(t), y:sin(t)}) 3 cos (t) sin (t) + sin (t)2 + 1

1 1 sin(t)

+ cos (t)

2 + 1

Veamos c´omo obtener el numerador y denominador de una expresi´on: sage: A = x*y/(1+x^2+y^2) sage: A.numerator()

xy sage: A.denominator() x2 + y 2 + 1 Se puede acceder a las distintas partes de una expresi´on mediante el m´etodo operands, mientras que operator nos devuelve la funci´on que une las partes: sage: A.operands() 

1 , x, y 2 x + y2 + 1



´ AL SAGE CAP´ITULO 1. INTRODUCCION

10

sage: A.operator() <function operator.mul> Con esta informaci´on, se puede recuperar la expresi´on original: sage: reduce(A.operator(), A.operands())

x2

xy + y2 + 1

Recordemos que reduce(f, lista es una funci´on de Python que aplica la funci´on f a los dos primeros elementos de la lista, despu´es el resultado al tercero y as´ı sucesivamente hasta que se obtiene un u ´nico valor: sage: reduce(lambda x, y: x*y, range(1,20)) #El resultado es 1*2*...*19 121645100408832000 Hemos usado una λ-funci´on, lambda x: x.simplify full(), que es una forma pr´actica y breve de Python (heredada de Lisp) para definir una funci´on “al vuelo”.

1.2.4.

Resolviendo ecuaciones

En Sage utilizaremos el s´ımbolo == para definir ecuaciones: sage: eq1 = x + y + 1 == 0 #Ejemplo de ecuaci´ on sage: eq2 = 3*x -5 == 4*y Podemos usar la funci´on solve para resolver el sistema: sage: solve([eq1, eq2], [x, y]) 

    8 1 ,y = − x= 7 7

Si queremos utilizar c´omodamente las soluciones, resulta pr´actico utilizar un diccionario como soluci´on: sage: solu = solve([eq1, eq2], [x, y], solution_dict=True) sage: solu

11

1.2. PRIMEROS PASOS CON SAGE

  1 8 x : ,y : − 7 7 sage: solu = solu[0] sage: bool(3*solu[x] - 5 == 4*solu[y]) True Encontremos los cortes de una elipse y una recta: sage: solve([x^2 + 3*y^2 == 1, x - y - 1 == 0], [x, y]) 

     1 1 ,y = − , [x = 1, y = 0] x= 2 2 Cuando no es posible encontrar soluciones exactas, podemos intentar encontrar una soluci´on aproximada: sage: f = x^5 - 4*x^3 + 1 sage: solve(f, x) 

0 = x5 − 4 x3 + 1



Intentar resolverlo exactamente no aporta ninguna soluci´on. Sin embargo un dibujo de la funci´on nos muestra que una de las soluciones est´a entre 0 y 1: sage: plot(f, -2, 2)

6

4

2

-2

-1

1

-2

-4

2

´ AL SAGE CAP´ITULO 1. INTRODUCCION

12

Pidamos a Sage que encuentre esa soluci´on: sage: find_root(f, 0, 1) 0,6541750137 Comprobemos que la soluci´on es correcta: sage: r = find_root(f, 0, 1) sage: f(x=r) 1,52655665886 × 10−16 No hay que obtener la impresi´on incorrecta de que Sage (ni ning´ un otro procedimiento) puede resolverlo todo: los sistemas de ecuaciones no lineales en varias variables pueden ser enormemente complicados de tratar, a´ un num´ericamente.

1.2.5.

C´ alculo diferencial

El m´etodo diff permite calcular derivadas de cualquier orden: sage: f = x*y / (1 + y + x^2) sage: f.diff(x, 2, y, 3) 480 x3 y 108 x 144 xy 288 x3 5 − 6 − 4 + 2 2 2 2 (x + y + 1) (x + y + 1) (x + y + 1) (x + y + 1)5 N´otese que hemos derivado haciendo actuar el m´etodo diff sobre la expresi´on f. Tambi´en podr´ıamos haber usado diff como una funci´on. Sin embargo, no todo m´etodo tiene su equivalente como funci´on ni viceversa. En la mayor parte de las ocasiones, haciendo uso de que en Python RtodoRes un objeto, 1 1 usaremos m´etodos para manipularlos. Calculemos ahora x=0 y=0 f dydx: sage: f.integrate(x, 0, 1).integrate(y, 0, 1)

3 1 log (3) + log (2) + 4 4 En ocasiones no es posible calcular integrales de forma exacta. Se puede hacer entonces una integraci´on simb´olica: −

13

1.2. PRIMEROS PASOS CON SAGE sage: f = sin(x)/x sage: f.integrate(x, 0, 1) Z

1 0

sin (x) dx x

sage: f.nintegrate(x, 0, 1) 0,946083070367, 1,05036320793 × 10−14 , 21, 0



La soluci´on de la integral es 0,946083070367, el error estimado es 1,05036320793× 10−14 , se han hecho 21 evaluaciones de la funci´on a integrar y el u ´ltimo 0 significa que no se han encontrado errores en la evaluaci´on.

1.2.6.

Representaci´ on gr´ afica de funciones

En Sage es sencillo hacer representaciones gr´aficas: plot(sin(1/x), (x, 0, pi)) 1

0.5

0.5

1

1.5

2

2.5

3

-0.5

-1

Veamos ahora c´omo dibujar el grafo de una funci´on f : R2 → R: var(’x y’) f = x*y / (1 + x^2 + y^2) plot3d(f, (x, -1, 1), (y, -1, 1))

´ AL SAGE CAP´ITULO 1. INTRODUCCION

14

A lo largo del curso estaremos interesados en hacer dibujos de curvas parametrizadas y tambi´en de superficies. Veamos, por ejemplo, una representaci´on gr´afica de una h´elice: sage: var(’t’) sage: parametric_plot3d([cos(t), sin(t), t], (t, -2*pi, 2*pi),\ thickness=5)

N´otese que hemos usado “\” para romper una l´ınea demasiado larga y continuar en la siguiente. El superponer varios dibujos es tambi´en f´acil. Aprovechemos para introducir gr´aficas definidas mediante ecuaciones impl´ıcitas: dib_elipse = implicit_plot(x^2 + 3*y^2 == 1, (x, -2, 2), (y, -1, 1)) dib_recta = implicit_plot(x - y - 1 == 0, (x, -2, 2), (y, -1, 1),\ color=’red’) dib_recta + dib_elipse 1

0.5

0

-0.5

-1 -2

1.2.7.

-1

0

1

2

´ Algebra lineal

Veamos c´omo definir vectores y calcular productos escalares y vectoriales: sage: U = vector([-1, 2, 0]) sage: V = vector([3, 1, -5]) sage: U.dot_product(V)

#producto escalar

15

1.2. PRIMEROS PASOS CON SAGE

−1 sage: U.norm()

#norma del vector √ 5

sage: U.cross_product(V) #producto vectorial (−10, −5, −7) En nuestro curso necesitaremos llevar a cabo algunas operaciones sencillas con matrices: multiplicarlas, invertirlas, calcular determinantes, rangos y autovectores. Veamos algunos ejemplos: sage: var(’a b c d’)

(a, b, c, d) sage: M = matrix([[a, b], [c, d]]) sage: M 

a b c d



sage: M.det() −bc + ad sage: M.inverse() 1 a



b

bc a2

(

bc −d a

)

c −d) a( bc a

) − bc1−d

a(

bc −d a a

!

No es posible hacer una simplificaci´on completa de una matriz aplicando directamente el m´etodo .simplify full(). Sin embargo, podemos simplificar aplicando el m´etodo a las coordenadas de la matriz: sage: invM = M.inverse() sage: matrix(map(lambda x: x.simplify_full(), invM))

´ AL SAGE CAP´ITULO 1. INTRODUCCION

16



d − bc−ad c bc−ad

b bc−ad a − bc−ad



Ha sido necesario volver a llamar al constructor matrix porque map devuelve una lista, no una matriz. Veamos un ejemplo de c´alculo de autovectores y autovalores. Primero construyamos una matriz que tenga autovalores sencillos: sage: P = matrix([[1, -1, 1], [1, 1, 1], [0, 1, 1]]) sage: D = matrix([[-1, 0, 0], [0, 2, 0], [0, 0, 0]]) sage: M = P*D*P.inverse() sage: M  1 −2 1  −1 0 1  −1 1 0 

sage: M.characteristic_polynomial(’x’) x3 − x2 − 2x sage: M.eigenvalues() [2, 0, −1] sage: M.eigenvectors_right() [(2, [(1, −1, −1)] , 1) , (0, [(1, 1, 1)] , 1) , (−1, [(1, 1, 0)] , 1)] Naturalmente, la u ´ltima respuesta significa que el vector (1, −1, −1) es un autovector de M por la derecha correspondiente al autovalor 2 y con multiplicidad 1, etc. Todos los ejercicios de este cap´ıtulo deben ser resueltos utilizando Sage: Ejercicio 1. Calc´ ulese el numerador y el denominador de x2

3 x − . +1 x−5

Ejercicio 2. Obt´engase una lista con todos los n´ umeros no primos menores o iguales que 1000.

17

1.2. PRIMEROS PASOS CON SAGE

Ejercicio 3. Dib´ ujense las rectas 2x − y − 1 = 0, x + y = 2 junto con su punto de corte. Ejercicio 4. Resu´elvase el sistema de ecuaciones 2x − y + z = 3y + 2z = 0, proporcionando un vector unitario que genere la recta de soluciones. Ejercicio 5. Factor´ıcese el polinomio p(x) = x4 − 2x3 − 15x2 − 4x + 20. Ejercicio 6. Encu´entrense todas las soluciones de la ecuaci´ on cos(x) − x2 + 1 = 0. Ejercicio 7. Calc´ ulese el gradiente de la funci´on f (x, y) =

2xy . x2 +y 2

Ejercicio 8. Calc´ ulese la integral de la funci´on f (x) = x sin(x2 ). Ejercicio 9. Dib´ ujese la funci´on f (x) = sin(x) junto con su recta tangente en x = π/4. Ejercicio 10. Dada la matrix   1 −1 3 A = 0 −2 0  1 1 −1

calc´ ulense las soluciones de la ecuaci´ on AX = b, siendo b = (2, 1, −1)⊤ . Calc´ ulese tambi´en la inversa de A, su polinomio caracter´ıstico y sus autovalores.

18

´ AL SAGE CAP´ITULO 1. INTRODUCCION

Cap´ıtulo 2 Curvas parametrizadas Comenzaremos adoptando la siguiente idea intuitiva de curva plana: es el conjunto de puntos del plano recorrido por un punto que se mueve en el tiempo. Esto implica que las coordenadas del punto m´ovil x(t) = (x(t), y(t)) deben ser funciones x = x(t), y = y(t) de la variable temporal. Dado que las t´ecnicas que usaremos para estudiar dichas curvas usar´an herramientas del c´alculo diferencial, exigiremos que las coordenadas sean funciones diferenciables con tantas derivadas como sea necesario. = ( dx , dy ). N´otese que conociendo La velocidad de la curva es el vector dx dt dt dt la velocidad v = v(t) de la curva y su posici´on x0 en un instante t = t0 , es posible recuperar la misma: Z t x(t) = x0 + v(t)dt. (2.1) t0

Ejercicio 11. Una curva x = x(t) tiene por velocidad v(t) = (− sin t, cos t). Calc´ ulense sus ecuaciones y dib´ ujese. Ejercicio 12. Calc´ ulense las ecuaciones de las curvas cuya velocidad est´a dada por v(t) = (−a sin t, b cos t) siendo a, b > 0. Ejercicio 13. La curva x(t) = (t2 + 3t + 1, 2t2 − t − 5) es una par´ abola. Encu´entrese su ecuaci´ on impl´ıcita. Ejercicio 14. Demu´estrese que una curva espacial cuyas componentes sean polinomios de grado ≤ 2 est´ a contenida en un plano. (Indicaci´ on: el espacio de polinomios de grado ≤ 2 tiene dimensi´ on tres, y los tres polinomios, junto con el polinomio constante 1 son cuatro polinomios). Ejercicio 15. — Demu´estrese que, si x(t) e y(t) son dos polinomios de grado dos tales que {x(t), y(t), 1} son linealmente independientes, entonces 19

CAP´ITULO 2. CURVAS PARAMETRIZADAS

20

la curva es una par´ abola (Indicaci´ on: hay un cambio af´ın de coordenadas que transforma la curva original en (t, t2 ), que es manifiestamente una par´ abola. Para encontrarlo basta con escribir los polinomios en t´erminos de la base can´onica.) Ejercicio 16. Demu´estrese que una curva espacial cuyas componentes son polinomios independientes de grado ≤ 2 es una par´ abola. ¿Y si son dependientes? Ejercicio 17. De una curva x = x(t) se sabe que x(0,5) = (1, −1) y su 2 velocidad es v(t) = (t, e−t ). Calc´ ulese x(1) con dos d´ıgitos de precisi´ on.

2.1.

Ejemplos de curvas planas

Antes de comenzar el estudio te´orico de las curvas, es conveniente disponer de ejemplos suficientes que muestren su inter´es. Un punto Existe la posibilidad de que el movimiento del punto tenga velocidad nula, es decir, v(t) = 0 = (0, 0). En ese caso la integral de la ecuaci´on (2.1) es nula y por tanto x(t) = x0 . Rectas Cuando un punto se desplaza con velocidad constante, v(t) = v0 entonces Z t v0 dt = v0 (t − t0 ), t0

as´ı que la curva resultante es la recta Z t x(t) = x0 + v(t)dt = x0 + v0 (t − t0 ). t0

En coordenadas, si x0 = (x0 , y0 ) y v0 = (u, v), resulta ( x(t) = x0 + u(t − t0 ) y(t) = y0 + v(t − t0 ) No obstante, una recta puede ser recorrida a velocidad no constante.

2.1. EJEMPLOS DE CURVAS PLANAS

21

Ejercicio 18. Escr´ıbase una curva parametrizada que recorra una recta pero a velocidad no constante. Dib´ ujese, usando Sage. Ejercicio 19. Demu´estrese que los puntos de una curva est´ an contenidos en una recta, es decir, satisfacen una ecuaci´ on lineal del tipo ax + by + c = 0, si y s´ olo si su velocidad v(t) es de la forma v(t) = f (t)v0 para alguna funci´ on diferenciable f (t) y un vector constante v0 . La aceleraci´on de una curva es el vector a(t) =

dv d2 x = 2. dt dt

Tambi´en emplearemos la notaci´on funcional seg´ un la cual dx dt d2 x ′′ x (t) = 2 , dt x′ (t) =

y as´ı sucesivamente. Ejercicio 20. Demu´estrese que una curva x = x(t) con velocidad nunca nula es una recta si y s´ olo si su aceleraci´ on es proporcional a su velocidad (se entiende que los coeficientes de proporcionalidad son funciones de la variable temporal). ¿Es este resultado razonable, intuitivamente? Sugerencia: Para demostrar que la condici´ on es necesaria, puede utilizarse el ejercicio anterior. Para la suficiencia, la ecuaci´ on diferencial x′′ (t) = λ(t)x′ (t) es f´ acil de resolver. Ejercicio 21. Demu´estrese que la funci´on ( 2 e−1/t si t 6= 0 f (t) = 0 si t = 0 es diferenciable (de clase infinito). Indicaci´ on: demu´estrese por inducci´on que las derivadas de f son de la forma ( 2 Rn (t)e−1/t si t 6= 0 (n) f (t) = 0 si t = 0, siendo Rn una funci´on racional y calc´ ulese l´ımt→0 f (n) (t).

CAP´ITULO 2. CURVAS PARAMETRIZADAS

22

Ejercicio 22. Usando la funci´on definida en el ejercicio anterior, escr´ıbase la parametrizaci´on de una curva diferenciable que est´e contenida en el eje OX si t < 0 y en el eje OY si t > 0. ¿Es su aceleraci´ on proporcional a su velocidad? ¿Contradice esto el ejercicio 20? Ejercicio 23. Supongamos que todas las rectas tangentes a una curva cuya velocidad es nunca nula pasan por un punto fijo. Demu´estrese que la curva est´ a contenida en una recta. (Indicaci´ on: Si llamamos al punto x0 , tendremos ′ que x(t) + λ(t)x (t) = x0 para cierta funci´on λ. Der´ıvese esta relaci´on y util´ıcese el ejercicio 20. Ejercicio 24. Dar un ejemplo que muestre que el resultado del ejercicio anterior no es cierto si suponemos que la velocidad de la curva puede anularse.

2.1.1.

Circunferencia

Una circunferencia de centro p0 = (x0 , y0 ) y radio r puede parametrizarse como x(t) = p0 + r(cos t, sin t). Ejercicio 25. Demu´estrese que los puntos de esta curva distan r del centro. Dibujemos la circunferencia con centro (1, 2) y radio 3: t = x0, r = p =

var(’t’) y0 = 1, 2 3 parametric_plot((x0 + r * cos(t), y0 + r*sin(t)), (t, 0 ,2 * pi)) p.set_aspect_ratio(1) 5

4

3

2

1

-2

-1

1

-1

2

3

4

23

2.1. EJEMPLOS DE CURVAS PLANAS

2.1.2.

Curvas polin´ omicas

Una clase especialmente importante de curvas son las curvas polin´omicas, es decir, tales que x(t) e y(t) son polinomios en la variable t, digamos de grado ≤ n. Tienen la ventaja de que su c´alculo es muy r´apido y est´an controladas por n + 1 par´ametros. Sin embargo, la relaci´on entre sus P coeficientes y la curva final es dif´ıcil de controlar. Consideremos una curva ni=0 vi ti , siendo vi puntos de R2 . El siguiente dibujo muestra un ejemplo de la curva resultante con el pol´ıgono formado por los puntos vi . m = 3 v = [vector([-2, 2]), vector([0, 1]), vector([1, 2]),vector([2, 1])] curva = sum(v[i]*t^i for i in range(m + 1)) curva_plot = parametric_plot(curva, (t, 0, 1), color = ’green’) v_plot = line([v[i].list() for i in range(m + 1)], marker = ’o’, color = ’red’) dib = curva_plot + v_plot dib.set_aspect_ratio(1) 6

5

4

3

2

1 -2

-1

1

0

2

Una raz´on que explica la mala relaci´on entre la curva y el pol´ıgono dado por sus coeficientes es que no es geom´etrica: si consideramos una afinidad f : R2 → R2 , es decir, una aplicaci´on de la forma f (x) = Ax + b, siendo A una matriz 2 × 2 y b un punto de R2 , resulta que ! n n X X vi ti 6= f f (vi )ti , i=0

i=0

CAP´ITULO 2. CURVAS PARAMETRIZADAS

24

salvo para valores excepcionales de t. Con palabras, no coincide la transformada de la curva con la curva definida por el pol´ıgono transformado. Ejercicio 26. Demu´estrese, con un ejemplo, que no es cierto que la curva transformada mediante una afinidad coincida con la curva asociada al pol´ıgono transformado. ¿Cu´ales son las curvas que tienen esta propiedad para toda afinidad? Para entender la raz´on por la que las curvas polin´omicas escritas en la base est´andar tienen un mal comportamiento, debemos recordar lo que son las combinaciones afines (tambi´en llamadas baric´entricas) de puntos.

2.2.

Curvas de B´ ezier

Para definir y entender las curvas de B´ezier, es necesario introducir algunos conceptos de geometr´ıa af´ın.

2.2.1.

Afinidades

Las transformaciones afines del plano son de la forma f (x) = Ax + b, siendo A=



a00 a01 a10 a11



una matriz no singular, i.e., det A 6= 0 y   b b= 0 b1 representa una traslaci´on. Algunos casos particularmente importantes de afinidades son los siguientes: Traslaciones Cuando A = Id, la afinidad se reduce a f (x) = x + b, y se llama traslaci´on de vector b.

´ 2.2. CURVAS DE BEZIER

25

Isometr´ıas En el caso de que la matriz A sea ortogonal, es decir, AA⊤ = Id, la afinidad se llama isometr´ıa. Isometr´ıa significa “igual medida”, porque dichas transformaciones preservan las distancias, es decir, para cualesquiera puntos x, y se tiene que d(f (x), f (y)) = d(x, y). Veamos que el que una afinidad preserve las distancias es equivalente a que la matriz A sea ortogonal. Para ello observemos que d(x, y)2 = kx − yk2 = (x − y)⊤ (x − y) = (x − y)⊤ Id(x − y),

(2.2)

mientras que d(f (x), f (y))2 = kAx + b − (Ay + b)k2 = (Ax − Ay)⊤ (Ax − Ay)

(2.3)

= (x − y)⊤ A⊤ A(x − y).

Concluimos por tanto que las expresiones en (2.2) y (2.3) son iguales para todo x, y si y solo si A⊤ A = Id, (2.4) es decir, A es una matriz ortogonal. Ejercicio 27. Demu´estrese esta u ´ltima afirmaci´on. Indicaci´ on: h´agase y = 0. Las derivadas parciales de segundo orden respecto a x de ambas expresiones deben ser iguales y determinan los coeficientes de las matrices. Tomando determinantes en (2.4), vemos que |A|2 = 1, con lo que |A| = ±1. Ejercicio 28. Demu´estrese con un ejemplo que no toda matriz con determinante ±1 tiene por qu´e ser una matriz ortogonal. Las isometr´ıas cuya matriz A tienen determinante +1 se llaman movimientos. Ejercicio 29. Demu´estrese que un movimiento plano tiene una matriz de la forma   cos θ − sin θ . A= sin θ cos θ

Indicaci´ on: anal´ıcese el sistema de ecuaciones que resulta de que A⊤ A = Id.

CAP´ITULO 2. CURVAS PARAMETRIZADAS

26

Ejercicio 30. Demu´estrese que una isometr´ıa plana que no es un movimiento tiene una matriz de la forma   cos θ sin θ . A= sin θ − cos θ

Ejercicio 31. Calc´ ulense todos los movimientos que transforman el segmento √ √ [(1, 1), (1, −1)] en el segmento [(1, 1/ 2), ( 2 + 1, 0)].

2.2.2.

Combinaciones baric´ entricas

Las combinaciones lineales usuales no son compatibles con las afinidades debido esencialmente a la existencia de traslaciones. En efecto, consideremos puntos x, y y escalares λ, µ ∈ R. Sea f la traslaci´on de vector b. Entonces f (λx + µy) = λx + µy + b,

mientras que λf (x) + µf (y) = λx + µy + (λ + µ)b. Si queremos que ambos resultados sean iguales para cualquier traslaci´on, debemos exigir que λ + µ = 1. De forma P m´as general, dados puntos xi y escalares λi ∈ R, i = 0, . . . , m, tales que m i=0 λi = 1, diremos que m X

λi x i

i=0

es la combinaci´on baric´entrica de los puntos xi con pesos λi . Ejercicio 32. Demu´estrese que las afinidades respetan las combinaciones afines de puntos, es decir, ! m m X X λi x i = f λi f (xi ) i=0

si f es una afinidad y

i=0

m X

λi = 1.

i=0

Ejercicio 33. Demu´estrese que las afinidades transforman rectas en rectas. Ejercicio 34. Dada la afinidad f (x) =



   1 1 −1 , x+ 2 0 1

encu´entrese la imagen de la recta y = x + 1.

´ 2.2. CURVAS DE BEZIER

27

Convexidad El segmento determinado por los puntos x y y est´a dado por → : µ ∈ [0, 1]}. [x, y] := {x + µ− xy

→ = y − x, es Puesto que − xy

→ = (1 − µ)x + µy. x + µ− xy

Ejercicio 35. Demu´estrese que el segmento [x, y] es el conjunto de combinaciones baric´entricas λx + µy con λ + µ = 1 y λ, µ ≥ 0. Los pesos λ, µ fijan en qu´e proporci´on divide el punto p = λx + µy al segmento [x, y]. En efecto, si escribimos → p = x + µ− xy, o, equivalentemente

− → = µ− → = µ(− →+− → xp xy xp py),

de donde se obtiene que − →= xp

µ − → = µ− → py py. 1−µ λ

Llamaremos a la cantidad µ/λ raz´on simple de los puntos x, p y y y la denotaremos como (x, p, y). Esto a veces se denota de la siguiente forma: − →:− → = µ : λ. xp py Ejercicio 36. Demu´estrese que las afinidades preservan la raz´on simple, es decir, si f es una afinidad y x, p, y son puntos alineados, entonces (x, p, y) = (f (x), f (y), f (p)). Ejercicio 37. Dados los puntos x = (1, 1)⊤ y y = (2, −1)⊤ , encu´entrese un punto p cuya raz´on simple (x, p, y) = t, siendo t ∈ R. Ejercicio 38. Demu´estrese con un ejemplo que las afinidades no preservan ni distancias, ni ´ angulos, ni ´ areas. Ejercicio 39. Diremos que los puntos p0 , . . . , pn forman una referencia af´ın de Rn si los vectores p0 p1 , . . . , p0 pn forman una base de Rn . Demu´estrese que, en este caso, cualquier punto ´nica como una Pn p se escribe de forma u combinaci´on baric´entrica p = i=0 λi pi .

CAP´ITULO 2. CURVAS PARAMETRIZADAS

28

Ejercicio 40. Demu´estrese que los puntos (0, 1, 1), (1, −1, 1), (1, 1, 0), (1, 1, 1) son una referencia af´ın y calc´ ulese las coordenadas baric´entricas del punto (x, y, z). Ejercicio 41. Progr´amese en Sage una funci´on que dibuje el segmento [x, y] y, dada una raz´on simple r dibuje el punto correspondiente p y los segmentos [x, p] y [p, y]. Sugerencia: u ´sese la funci´on @interact para poder modificar interactivamente el valor de r. Ejercicio 42. Progr´amese en Sage una funci´on que dibuje el segmento [x, y] y, dado un escalar λ dibuje el punto correspondiente p = λx + µy, λ + µ = 1. Sugerencia: u ´sese la funci´on @interact para poder modificar interactivamente el valor de λ. Veamos ahora c´omo interpretar una combinaci´on af´ın de tres puntos, λx + µy + νz, λ + µ + ν = 1. Escribiendo la combinaci´on de la siguiente forma:   µ ν λx + µy + νz = λx + (µ + ν) y+ z , (2.5) µ+ν µ+ν vemos que el punto resultante se puede escribir como combinaci´on af´ın de dos puntos, por una parte el punto x y un punto de la recta yz. Veamos a continuaci´on que el conjunto de combinaciones baric´entricas de tres puntos x, y, z con coeficientes no negativos coinciden con el tri´angulo △

cerrado xyz. Sea una combinaci´on baric´entrica, λx + µy + νz, λ + µ + ν = 1 y λ, µ, ν ≥ 0. En primer lugar, notemos que si es µ + ν = 0, al ser λ, µ ≥ 0 entonces λ = µ = 0 y la combinaci´on simplemente es el punto x. Supongamos µ ν pues que que µ+ν 6= 0. Entonces los coeficientes de µ+ν y + µ+ν z son tambi´en positivos y sabemos que esto significa que dicho punto est´a en el segmento [y, z]. Como adem´as el coeficiente λ ≥ 0 entonces el punto λx + µy + νz △

est´a en el tri´angulo cerrado xyz.



Para demostrar el otro contenido, basta con observar que xyz se puede escribir como uni´on de los segmentos [x, p] con p ∈ [y, z]. △

Ejercicio 43. Usando Sage, dib´ ujese el tri´ angulo xyz y el punto λx+µy+νz, λ+µ+ν = 1. Usando @interact, perm´ıtase mover los coeficientes y observar el resultado. Expl´ıquese, seg´ un el signo de los coeficientes, d´onde debe hallarse el punto respecto al tri´ angulo. Un conjunto S se dira convexo si para cualesquiera x, y ∈ S, se tiene que [x, y] ⊂ S. Dado un conjunto S ⊂ Rn , definimos la envoltura convexa de S, hull(S), como el menor de los conjuntos convexos que contienen a S.

´ 2.2. CURVAS DE BEZIER

29

Ejercicio 44. Dado S ⊂ Rn , se tiene que ) ( m m X X λi xi : λi ≥ 0, xi ∈ S, λi = 1, m ∈ N . hull(S) = i=0

i=0

Indicaci´ on: Demu´estrese que el conjunto dado por el lado derecho de la igualdad contiene a S, es convexo y que est´ a contenido en cualquier conjunto convexo que contenga a S. Para esto u ´ltimo, se puede hacer un argumento inductivo basado en (2.5). Ejercicio 45. Dib´ ujese una aproximaci´ on a la envoltura convexa de un conjunto finito dibujando un n´ umero suficientemente grande de puntos de la misma. ¿Se obtienen as´ı buenos resultados? ¿Qu´e problemas puede tener este m´etodo?

2.2.3.

Algoritmo de Jarvis para el c´ alculo de la envoltura convexa

Terminemos viendo un algoritmo sencillo, aunque no el m´as eficiente, para calcular la envoltura convexa de un conjunto finito S de puntos en el plano. Por simplicidad, supondremos que los puntos est´an situados en posici´on general, en el sentido de que no hay tres puntos de S alineados. La idea consiste en comenzar el c´alculo de la envoltura convexa con un punto p0 ∈ S con la menor abscisa posible. A continuaci´on consideramos todos los segmentos de la forma p0 p con p ∈ S, hasta que encontramos uno, p1 , con la propiedad de que p0 p1 deja a todos los otros puntos de S a su derecha. Incorporamos p1 a la frontera de la envoltura convexa y repetimos el procedimiento con ´el y as´ı sucesivamente hasta alcanzar un punto pm = p0 , en cuyo caso hemos terminado. N´otese que para decidir si un punto z est´a a la izquierda del segmento → − → forman una base positivamente [x, y] basta con ver si los vectores {− xy, xz} orientada, es decir, si su determinante es positivo. Ejercicio 46. Implem´entese, usando Sage, el algoritmo de Jarvis. Ejercicio 47. ¿Sabr´ıas implementar el algoritmo si los puntos no est´ an necesariamente en posici´on general? Dado que por cada punto de la frontera de la envoltura convexa hemos de hacer un n´ umero constante de operaciones por cada punto de S que no es ´el mismo, el coste del algoritmo de Jarvis es O(kn), siendo n el n´ umero de puntos de S y k el n´ umero de puntos de ∂ hull S.

30

CAP´ITULO 2. CURVAS PARAMETRIZADAS

Algorithm 1 Algoritmo de Jarvis para el c´alculo de la envoltura convexa Input: Un conjunto finito S de puntos del plano en posici´on general. Output: El pol´ıgono frontera de hull S, ∂ hull S, dada como una lista ordenada de puntos. p0 ← punto de S con menor abscisa. pinic ← p0 Inicializa ∂ hull S con p0 . repeat Para cada p ∈ S if p 6= pinic and [pinic , p] deja a los dem´as puntos de S a su derecha then A˜ nade p a la lista de puntos de ∂ hull S pinic ← p end if until p = p0

2.2.4.

Polinomios de Bernstein

Una forma de mejorar la relaci´on entre curva y pol´ıgono ser´ıa sustituir la base de polinomios {1, t, . . . , tn } por otra base {P0n (t), . . . , Pnn (t)} que satisfaga la condici´on X Pin (t) = 1, i

de forma que la curva est´e dada por una combinaci´on baric´entrica de puntos del pol´ıgono de control. La forma m´as sencilla de hacer esto es escribir 1 = (t + (1 − t))n y desarrollar esta identidad mediante el binomio de Newton. Resulta as´ı         n n n n−1 n n−2 n 2 1= t + t (1 − t) + t (1 − t) + · · · (1 − t)n . 0 1 2 n Este desarrollo sugiere definir los polinomios   n i n t (1 − t)n−i , Bi (t) = i que se llaman polinomios de Bernstein. Obviamente se trata de n + 1 polinomios de grado n, de forma que ser´a suficiente comprobar que son linealmente independientes para ver que forman una base.

´ 2.2. CURVAS DE BEZIER

31

Ejercicio 48. Dib´ ujense los polinomios de Bernstein con t ∈ [0, 1] para n = 1, . . . , 10. Ejercicio 49. Demu´estrese que los polinomios de Bernstein se pueden determinar recursivamente mediante las f´ormulas n−1 Bin (t) = (1 − t)Bin−1 (t) + tBi−1 (t),

B00 (t) = 1.

Indicaci´ on: Recu´erdese que los coeficientes binomiales satisfacen la relaci´on recursiva       n−1 n−1 n . + = i−1 i i

Ejercicio 50. Calc´ ulese el m´aximo de Bin (t) con t ∈ [0, 1].

Ejercicio 51. Demu´estrese que los polinomios de Bernstein de grado n son linealmente independientes y, por tanto, forman una base. Indicaci´ on: Obs´ervese que la matriz de sus coeficientes respecto a la base can´onica tiene forma triangular.

2.2.5.

Curvas de B´ ezier

Estamos ya en condiciones de introducir las curvas de B´ezier. Dado que los polinomios de B´ernstein son una base, cualquier curva polin´omica x = x(t) se podr´a escribir como n X x(t) = bi Bin (t). i=0

Una curva polin´omica escrita de esta forma y parametrizada con t ∈ [0, 1] se llama curva de B´ezier y los coeficientes (ordenados) bi se llaman pol´ıgono de control de la curva. Veamos un ejemplo de curva de B´ezier junto con el pol´ıgono de control correspondiente: B(n, i, t) = binomial(n, i) * t^i * (1 - t)^(n - i) b = [vector([-2, 2]), vector([0, 1]), vector([1, 2]),vector([2, 1])] m = len(b) - 1 curva = sum(b[i]*B(m,i,t) for i in range(m + 1)) curva_plot = parametric_plot(curva, (t, 0, 1), color = ’green’) b_plot = line([b[i].list() for i in range(m + 1)], marker = ’o’, color = ’red’) dib = curva_plot + b_plot dib.set_aspect_ratio(1)

CAP´ITULO 2. CURVAS PARAMETRIZADAS

32 2 1.8 1.6 1.4

1.2 1 -2

-1

1

0

2

Ejercicio 52. Demu´estrese que una curva de B´ezier satisface x(0) = b0 , x(1) = bn . Ejercicio 53. Demu´estrese que las curvas de B´ezier son invariantes bajo transformaciones afines, es decir, ! n n X X bi Bin (t) = f (bi )Bin (t) f i=0

i=0

si f es una transformaci´ on af´ın.

P Ejercicio 54. Demu´estrese que una curva de B´ezier x(t) = ni=0 bi Bin (t) permanece dentro de la envoltura convexa del pol´ıgono de control si t ∈ [0, 1]. Ejercicio 55. Simetr´ıa de las curvas de B´ezier: pru´ebese que la imagen de una curva de B´ezier con pol´ıgono de control {b0 , . . . , bn } coincide con la de la curva que tiene el pol´ıgono de control invertido, es decir, {bn , . . . , b0 }. Ejercicio 56. Invariancia por combinaciones afines: dadas curvas de B´ezier con pol´ıgonos de control {bi } y {ci }, que llamaremos b(t) y c(t), respectivamente, entonces la curva que tiene pol´ıgono de control {αbi + βci }, siendo α + β = 1, est´ a dada por αb(t) + βc(t). Ejercicio 57. Los polinomios de Bernstein satisfacen la identidad m X j m B (t) = t. m j j=0

Es decir, las coordenadas de t en la base de polinomios de Bernstein son {j/m}. Demu´estrese esto. Ejercicio 58. Demu´estrese que si los v´ertices del pol´ıgono de control de una curva de B´ezier se situan uniformemente distribuidos en el segmento [p, q], es decir,   j j bj = 1 − p + q, m m entonces la curva de B´ezier es la recta que une p con q parametrizada de forma natural.

´ 2.2. CURVAS DE BEZIER

33

Ejercicio 59. Demu´estrese que las derivadas de los polinomios de Bernstein satisfacen la relaci´on  d n n−1 Bi (t) = n Bi−1 (t) − Bin−1 (t) . dt Pn n Ejercicio 60. La derivada de una curva de B´ezier b(t) = i=0 bi Bi (t) viene dada mediante la f´ormula ′

b (t) =

n−1 X

n∆bi Bin−1 (t),

i=0

donde ∆bi = bi+1 − bi .

Ejercicio 61. Demu´estrese que una curva de B´ezier es tangente al pol´ıgono de control en los extremos del mismo. Ejercicio 62. Calc´ ulese la curva de B´ezier cuyo pol´ıgono de control es (1, 2, 0), (0, −1, 1), (3, 0, 1), (1, −1, 2), (0, 1, 5). Ejercicio 63. Escr´ıbase la curva x(t) = (t3 − 2t + 1, 4t − 7) con t ∈ [0, 1] como una curva de B´ezier.

Ejercicio 64. Escr´ıbase la curva x(t) = (t2 + 5t − 1, 4t2 − 1) con t ∈ [−3, 2] como una curva de B´ezier. (Indicaci´ on: comp´ongase la curva con una transformaci´ on af´ın que lleve el intervalo [−3, 2] en el [0, 1]). Ejercicio 65. Escr´ıbase la curva de B´ezier del ejercicio anterior como una curva de B´ezier de grado tres, es decir, con un pol´ıgono de control de cuatro v´ertices.

2.2.6.

El algoritmo de De Casteljau

Veamos en esta secci´on un m´etodo alternativo para generar curvas de B´ezier. Aunque computacionalmente es menos eficaz que la forma de Bernstein, es un procedimiento geom´etrico m´as sugerente. Se trata, b´asicamente, de hacer una interpolaci´on lineal iterada de los v´ertices del pol´ıgono de control. Ejercicio 66. Demu´estrese que la curva resultante del algoritmo de De Casteljau coincide con la curva de B´ezier correspondiente. Indicaci´ on: u ´sese inducci´on y el ejercicio (49). Ejercicio 67. Usando la funci´on animate de Sage, implem´entese el algoritmo de De Casteljau. Ejercicio 68. Calc´ ulense las curvas de B´ezier de grado 3 que en sus extremos (1, 2) y (1, 5) tienen velocidades (6, −9) y (−3, 9), respectivamente.

34

CAP´ITULO 2. CURVAS PARAMETRIZADAS

Algorithm 2 Algoritmo de De Casteljau Input: El pol´ıgono de control {b0 , . . . , bn } de la curva y t ∈ R. Output: La curva de B´ezier asociada b = b(t). b0i (t) = bi . for r = 1 → n do for i = 1 → n − r do r−1 bri (t) = (1 − t)bir−1 (t) + tbi+1 (t) end for end for return b(t) = bn0 (t).

2.3.

Las c´ onicas como curvas racionales

Las curvas polin´omicas no son suficientes para obtener algunos tipos de curvas importantes, como por ejemplo las c´onicas. Ejercicio 69. Demu´estrese que no existe ninguna curva polin´ omica no constante cuya imagen se encuentre contenida en la circunferencia x2 + y 2 = 1. (Indicaci´ on: si x(t), y(t) son polinomios, ¿qu´e sucede con el t´ermino de grado m´as alto de x(t)2 + y(t)2 − 1?) Una curva racional es la afinizaci´on de una curva polin´omica proyectiva. Es decir, si x(t) = (p(t) : q(t) : r(t)) es una curva polin´omica proyectiva, entonces   p(t) q(t) x(t) = , r(t) r(t) es la correspondiente curva racional af´ın (correspondiente a la afinizaci´on que supone que z = 0 es la recta del infinito). Las curvas racionales no aumentan excesivamente el coste computacional y permiten dibujar c´onicas. Para parametrizar una c´onica, que suponemos no degenerada, tomemos un punto x0 = (x0 , y0 ) de la misma y consideremos el haz de rectas de base x0 . Para ello tomemos dos rectas r0 (x, y) = 0 y r1 (x, y) = 0 no proporcionales cualesquiera que pasen por un punto x0 , (por ejemplo, si x0 = (x0 , y0 ) podemos tomar las rectas r0 : x − x0 = 0 y r1 : y − y0 = 0). El haz de rectas de base x0 es, por definici´on, el conjunto de todas las rectas que pasan por x0 y coincide con las combinaciones lineales (no triviales) de r0 y r1 , es decir, el haz de rectas es el conjunto de rectas de la forma ar0 (x, y) + br1 (x, y) = 0

´ 2.3. LAS CONICAS COMO CURVAS RACIONALES

35

donde los coeficientes a, b no son simult´aneamente nulas. Ejercicio 70. Demu´estrese esta u ´ltima afirmaci´on. Como cualquier m´ ultiplo no nulo de la ecuaci´on de una recta produce la misma recta, podemos usar el par´ametro t = b/a en tanto que a 6= 0 y describir el haz como las rectas como r0 (x, y) + t r1 (x, y) = 0. Si hacemos esto, tenemos el inconveniente de que perdemos la recta r1 (x, y) = 0, a la que nos podemos acercar tanto como queramos haciendo que t → ∞. Cada recta de dicho haz corta a la c´onica en el punto x0 y en otro punto adicional, x = x(t). Veamos que as´ı se obtiene una parametrizaci´on racional de la curva. Para ello, conviene simplificar los c´alculos haciendo el cambio af´ın de coordenadas ( X = r0 (x, y) Y = r1 (x, y)

4

r0 :2x + y−7 =0 3.5

r1 :−x +2y−4 =0 3

2.5

r = r0 −5r1

El que las rectas se cortan en un punto finito x0 nos asegura que las ecuaciones del cambio son invertibles Figura 2.1: Haz engendrado por r0 y (¿por qu´e?). En las nuevas coordena- r1 das, el punto de corte es (0, 0) y el haz de rectas es X + tY = 0. La c´onica, dado que tiene que pasar por el origen de coordenadas, no tendr´a t´ermino independiente, es decir, su ecuaci´on ser´a de la forma 2

1

1.5

2

2.5

3

f (X, Y ) = aX 2 + bY 2 + 2cXY + 2dX + 2eY = 0. Haciendo X = −tY , obtenemos la ecuaci´on para calcular la ordenada de los puntos de intersecci´on:  f (−tY, Y ) = (at2 − 2ct + b)Y + 2(−dt + e) Y = 0,

que tiene las soluciones Y = 0 e

Y =

2(dt − e) . at2 − 2ct + b

CAP´ITULO 2. CURVAS PARAMETRIZADAS

36

Como la ordenada es X = −tY , la parametrizaci´on que se obtiene es pues   2(dt − e) −2(dt − e)t , c(t) = at2 − 2ct + b at2 − 2ct + b Deshaciendo el cambio af´ın de coordenadas, obtenemos la parametrizaci´on racional deseada (¿por qu´e?). Veamos un ejemplo. Consideremos la elipse f (x, y) = 0 siendo f(x, y) = 4*x^2 + 9*y^2 - 8*x + 18*y + 12 elipse_dib = implicit_plot(f(x, y), (x, 0, 2), (y, -2, 0)) elipse_dib.set_aspect_ratio(1) 0

-0.5

-1

-1.5

-2 0

0.5

1

1.5

2

Encontremos primero un punto de la misma. Para ello, cortemos la c´onica con una recta de la forma x = c. Como los puntos de corte pueden ser imaginarios, as´ı que elegiremos c de forma que sean puntos reales. Las soluciones de f (c, y) = 0 son 

 1√ 1√ 2 2 y=− −4 c + 8 c − 3 − 1, y = −4 c + 8 c − 3 − 1 . 3 3

Para que el polinomio −4 c2 + 8 c − 3 no sea negativo, debe ser 1/2 ≤ c ≤ 3/2. Si tomamos, por ejemplo, c = 1/2 obtenemos el punto x0 = (1/2, −1). Ahora que tenemos un punto de la c´onica, parametricemos el haz de rectas que pasan por ´el usando por ejemplo las rectas x − 1/2 = 0 e y + 1 = 0, de forma que una recta gen´erica es de la forma (x − 1/2) + t(y + 1) = 0. Cort´emosla con la elipse: sage: var(’t’)

´ 2.3. LAS CONICAS COMO CURVAS RACIONALES

37

t sage: sols = solve([x-1 +t*(y+4/3)==0, f(x,y)==0], x, y) Se observa que la soluci´on no constante viene dada por la curva racional x=

2 (8 t2 + 9) 4 t2 − 6 t + 9 , y = − 4 t2 + 9 3 (4 t2 + 9)

como se quer´ıa. curva = map(lambda X: X.rhs(),sols[1]) elip_param = parametric_plot(curva, (t, -10, 10)) elip_param.set_aspect_ratio(1) 0.6

0.8

1

1.2

1.4

-0.7

-0.8

-0.9

-1

-1.1

-1.2

-1.3

Ejercicio 71. Como se observa, al dibujo de la elipse como curva parametrizada le falta un trozo. ¿Por qu´e? ¿C´omo podr´ıas hacer para pintar la curva entera? Ejercicio 72. Apl´ıquese un procedimiento parecido al anterior para obtener parametrizaciones de la hip´erbola de ecuaci´ on (x + y − 1)(x − y + 1) = 1. Ejercicio 73. H´agase el ejercicio anterior haciendo el cambio af´ın de coordenadas que sugiere la teor´ıa y deshaci´endolo al final. Ejercicio 74. Usando Sage, encu´entrese la matriz de la c´onica recorrida por la curva parametrizada como   2t − 1 t+1 , t2 + t + 1 t2 + t + 1 Ejercicio 75. Progr´amese, usando Sage, un procedimiento que dada una c´onica devuelva un punto real de la misma, si tal punto existe.

38

CAP´ITULO 2. CURVAS PARAMETRIZADAS

Cap´ıtulo 3 Curvas regulares. Curvatura. 3.1.

Curvas regulares. Longitud de arco.

Sea I ⊂ R in intervalo abierto y sea x : I → Rn una curva diferenciable. Un punto t ∈ I se dice singular cuando x′ (t) = 0. Una curva es regular cuando no tiene puntos singulares. Ejemplo 1. La curva x(t) = (t2 , t3 ) tiene en t = 0 un punto singular. 1

0.5

0.2

0.4

0.6

0.8

1

-0.5

-1

Llamaremos recta tangente a la curva en un punto regular t a la recta definida por x(t) y el vector x′ (t). En los puntos singulares no existe recta tangente, seg´ un esta definici´on. 39

CAP´ITULO 3. CURVAS REGULARES. CURVATURA.

40

Ejercicio 76. Demu´estrese que la curva x(t) = (t3 − 4t, t2 − 4) es regular. Sin embargo, en el (0, 0) tiene dos posibles rectas tangentes. ¿Contradice esto 2

1

-4

-2

2

4

-1

-2

-3

la regularidad?

-4

Ejercicio 77. Encu´entrese una curva x : R → R2 tal que x(R) sea el eje OX, pero que no sea una curva regular. Ejercicio 78 (*). Sea p un punto de una circunferencia que rueda sobre una recta sin resbalar. Parametr´ıcese la curva que p describe y encu´entrese sus puntos singulares. Dib´ ujese. Una tal curva se llama cicloide. Ejercicio 79. La hipocicloide es la curva generada por la traza de un punto p sobre una circunferencia de radio r que rueda sin deslizarse por el interior de un c´ırculo de radio R > r. Encu´entrense sus ecuaciones. Ejercicio 80. En el ejercicio anterior, dib´ ujese la curva y encu´entrense sus puntos singulares en el caso de que r = 1, R = 3. Ejercicio 81. Est´ udiese si la curva x(t) = (t2 + 1, t3 − t, t − 1) es regular. Calc´ ulese su pol´ıgono de control. Dib´ ujese. La longitud de una curva restringida a un intervalo cerrado [a, b] (¿por qu´e esta restricci´on?) se puede aproximar del siguiente modo: dada una partici´on P = {t0 = a < t1 < · · · < tk = b}, la longitud asociada a la partici´on se define como k−1 X L(P) = kx(ti+1 ) − x(ti )k i=0

Ejercicio 82 (*). Escr´ıbase un programa en Sage que calcule L(P) siendo la partici´on la dada al dividir el intervalo en n trozos iguales. Se llama longitud de la curva al n´ umero L(x) = sup{L(P) : P es partici´on del intervalo [a, b]}. N´otese que para que este n´ umero no sea infinito, el conjunto debe estar acotado, lo que requiere demostraci´on.

3.1. CURVAS REGULARES. LONGITUD DE ARCO.

41

Teorema 1. Dada una curva diferenciable x, se tiene que

L(x) =

Z

b a

kx′ (t)kdt.

Demostraci´ on. Queremos ver que dos supremos son iguales, por una parte L(x) = sup{L(P) : P es partici´on del intervalo [a, b]}. y por otra parte Z

b a

kx′ (t)kdt = sup {I(P) : P es partici´on de [a, b]} ,

donde I(P) denota la suma inferior de Riemann

I(P) =

k−1 X i=0

mi (ti+1 − ti )

y, en la u ´ltima ecuaci´on, P = {ti }ki=0 y mi = m´ın{kx′ (t)k : t ∈ [ti , ti+1 ]}. El que la integral se puede escribir de esta forma es simplemente parte de la definici´on de integral de Riemann. Sabemos que, por ser kx′ (t)k continua en [a, b], su integral es finita. Sin embargo esto no est´a tan claro para L(x), as´ı que asegur´emonos de que L(x) < ∞. Para ello, observemos que, dada una partici´on P, por el teorema de Lagrange del valor medio, podemos escribir xh (ti+1 ) − xh (ti ) = x′h (αih )(ti+1 − ti ) para ciertos αih ∈ (ti , ti+1 ) y cada h = 1, . . . , n. Pero nuevamente por tratarse x′h de un n´ umero finito de funciones continuas en un compacto, est´an uniformemente acotadas, es decir, existe M ≥ 0 tal que |x′h (t)| ≤ M para toda

CAP´ITULO 3. CURVAS REGULARES. CURVATURA.

42

componente h y t ∈ [a, b]. Por tanto L(P) =

k−1 X i=0

v k−1 uX X u n t (xh (ti+1 ) − xh (ti ))2 kx(ti+1 ) − x(ti )k = i=0

h=1

v k−1 uX X u n 2 t x′h (αih )(ti+1 − ti ) = i=0

h=1

v u n k−1 X uX 2 t = x′h (αih ) (ti+1 − ti ) i=0

(3.1)

h=1

v u n k−1 X uX (ti+1 − ti ) ≤t M2

=



h=1

i=0

nM (b − a).

Como esta acotaci´on no depende de la partici´on empleada, concluimos que L(x) ≤



nM (b − a) < ∞.

Por definici´on de supremo y por ser la longitud asociada a una partici´on ≥ 0, resulta que para todo ǫ > 0 existe una partici´on P tal que 0 ≤ L(x) − L(P) < ǫ. An´alogamente, por ser kx′ (t)k ≥ 0 y tambi´en por la definici´on de supremo, existe una partici´on P ′ tal que 0≤

Z

b a

kx′ (t)kdt − I(P ′ ) < ǫ.

Veamos que podemos conseguir un refinamiento com´ un de P y P ′ que satisface las mismas desigualdades. Para ello, notemos que refinar P (o P ′ ) incrementa L (o I). Para ver esto, basta con a˜ nadir un punto a P, obteniendo ¯ P. Por simplicidad, suponemos que a˜ nadimos un punto al final, de forma ¯ que P = {t0 , t1 , . . . , tk , tk+1 }. Se argumentar´ıa de forma similar si se a˜ nade

3.1. CURVAS REGULARES. LONGITUD DE ARCO.

43

en cualquier otra posici´on. Tenemos que L(P) = = = ≤

k−1 X

i=0 k−2 X

i=0 k−2 X

i=0 k−2 X i=0

kx(ti+1 ) − x(ti )k kx(ti+1 ) − x(ti )k + kx(tk ) − x(tk−1 )k kx(ti+1 ) − x(ti )k + kx(tk ) − x(tk+1 ) + x(tk+1 ) − x(tk−1 )k kx(ti+1 ) − x(ti )k + kx(tk ) − x(tk+1 )k + kx(tk+1 ) − x(tk−1 )k

¯ = L(P).

Un argumento inductivo muestra que, en general, si P¯ es m´as fina que P (y ¯ entonces escribiremos P ≤ P) ¯ L(P) ≤ L(P).

Se demuestra de forma parecida que si P ′ ≤ P¯ ′

I(P ′ ) ≤ I(P¯ ′ )

Ejercicio 83 (*). Pru´ebese esta u ´ltima afirmaci´on. De forma que podemos refinar P y P ′ obteniendo una partici´on P ′′ m´as fina que ambas y tan fina como deseemos. Para esta partici´on ocurrir´a que 0 ≤ L(x) − L(P ′′ ) ≤ L(x) − L(P) < ǫ 0 ≤ I(x) − I(P ′′ ) ≤ I(x) − I(P ′ ) < ǫ. En adelante escribiremos P ′′ = P, por simplicidad. Veamos cu´al es la distancia entre L(P) y I(P): k−1 X |L(P) − I(P)| = kx(ti+1 ) − x(ti )k − I(P) i=0 v X n k−1 k−1 uX X u  2 j t x′j (αi ) (ti+1 − ti ) − mi (ti+1 − ti ) = i=0 j=1 i=0  v X u n k−1 uX 2 j ′ t xj (αi ) − mi  (ti+1 − ti ) . = i=0 j=1

CAP´ITULO 3. CURVAS REGULARES. CURVATURA.

44

Tenemos que analizar si podemos hacer tan peque˜ na como queramos la diferencia v uX u n 2 t x′j (αij ) − mi . j=1

Para ello, definamos la funci´on

F : [a, b]n → R como F (t) = k(x′1 (t1 ), . . . , x′n (tn ))k.

donde hemos llamado, para abreviar, t = (t1 , . . . , tn ). Esta es una funci´on cont´ınua y [a, b]n es compacto, as´ı que es uniformemente continua. Esto significa que para todo ǫ > 0 existe δ > 0 tal que si kt − t′ k < δ entonces |F (t) − F (t′ )| < ǫ/(b − a). Por otra parte, como kx′ k es continua y [ti , ti+1 ] es compacto, el m´ınimo mi se alcanza en alg´ un punto del intervalo [ti , ti+1 ], digamos βi . Ahora refinemos √ la partici´on P de forma que |P| = m´ın{ti+1 − ti , i = 0, . . . , k − 1} < δ/ n. Entonces v uX u n  2 j ′ = F (αi1 , . . . , αin ) − F (βi , . . . , βi ) < ǫ/(b − a) t x (α ) − m i j i j=1

ya que

v v uX uX n u u n √ j 1 n 2 t k(αi , . . . , αi ) − (βi , . . . , βi )k = (αi − βi ) ≤ t (δ/ n)2 = δ. j=1

j=1

As´ı que ya tenemos que

 v X uX n k−1 u 2 j ′ t |L(P) − I(P)| = xj (αi ) − mi  (ti+1 − ti ) i=0 j=1 v k−1 uX X u n 2 j ′ t ≤ xj (αi ) − mi (ti+1 − ti ) i=0 j=1 n

ǫ X ǫ ≤ (ti+1 − ti ) = (b − a) = ǫ. b − a j=1 b−a

45

3.1. CURVAS REGULARES. LONGITUD DE ARCO. Finalmente, vemos que Z b ′ L(x) − kx (t)kdt a Z b ′ kx (t)kdt = L(x) − L(P) + L(P) − I(P) + I(P) − a Z b ′ kx (t)kdt ≤ |L(x) − L(P)| + |L(P) − I(P)| + I(P) − a

≤ ǫ + ǫ + ǫ = 3ǫ.

Como esto ocurre para todo ǫ > 0, concluimos que Z b kx′ (t)kdt, L(x) = a

como quer´ıamos demostrar. Observaci´ on 1. Desde un punto de vista intuitivo, el teorema anterior puede interpretarse as´ı: en un tiempo infinitesimal dt la curva recorre un arco de longitud ds. 1

ds2

0.9

dx

= dx2 + dy2

dy

0.8

0.7

0.6

0

0.1

0.2

0.3

0.4

0.5

0.6

0.7

0.8

La infinitesimalidad del arco nos permite suponer que ds en realidad es la hipotenusa de un tri´ angulo rect´ angulo de lados dx y dy. El teorema de Pit´ agoras aplicado a los elementos infinitesimales asegura que p ds = dx2 + dy 2 con lo que, integrando,

Z Z p dx2 + dy 2 = s=

s

dx dt

2

+



dy dt

2

dt.

46

CAP´ITULO 3. CURVAS REGULARES. CURVATURA.

Ejercicio 84. Este argumento es intuitivamente interesante, pero no es riguroso. ¿Qu´e deficiencias le ves? Ejercicio 85. Calc´ ulese la longitud del segmento de recta determinado por los puntos p y q. Ejercicio 86. Demu´estrese que la longitud de un arco de circunferencia de a´ngulo θ y radio r es θr. Ejercicio 87. Calc´ ulese una aproximaci´ on num´erica a la longitud de una elipse de semiejes a = 2 y b = 3 con un error menor que 10−3 . Ejercicio 88 (*). Una curva se desplaza en un cilindro circular de radio a, girando alrededor del cilindro con velocidad angular ω y subiendo a lo largo del mismo con velocidad b. Calc´ ulese su longitud que recorre cuando da una vuelta completa alrededor del cilindro. Ejercicio 89 (*) (La distancia m´as corta entre dos puntos es la l´ınea recta). Demu´estrese que si la longitud de una curva x : [a, b] → Rn satisface la desigualdad kx(b) − x(a)k ≥ L(x),

entonces la curva es una recta (en el sentido de que est´ a contenida en el segmento [x(a), x(b)] y su longitud es kx(b) − x(a)k. (Indicaci´ on: si en alg´ un instante t1 la curva no est´ a contenida en el segmento, consid´erese la partici´on {t0 = a, t1 , t2 = b}. La longitud de arco de una curva es la funci´on s : I → R definida por Z t s(t) = kx′ (u)kdu, t0

siendo t0 ∈ I. Intuitivamente, la longitud de arco mide la distancia sobre la curva entre el punto parametrizado por t0 y el parametrizado por t, teniendo en cuenta que el signo viene dado por el sentido de recorrido de la curva. Cuando se cambia t0 se obtiene una nueva longitud de arco que difiere de la anterior en una constante. Ejercicio 90. Demu´estrese esta u ´ltima afirmaci´on. Ejercicio 91. Encu´entrese la longitud de arco de la cicloide, definida en el ejercicio 78. Ejercicio 92. Consid´erese la curva x(t) = aebt (cos t, sin t), siendo a > 0 y b < 0. Dib´ ujese para algunos valores de a y b. Demu´estrese que la curva tiende al origen de coordenadas cuando t → ∞, dando infinitas vueltas alrededor de ´el, pero que, sin embargo, el l´ımite de su longitud de arco l´ımt→∞ s(t) < ∞.

3.1. CURVAS REGULARES. LONGITUD DE ARCO.

47

Una curva se dice que est´a parametrizada por la longitud de arco cuando est´a recorrida con velocidad uno. Equivalentemente, s(t) = t − t0 . Ejercicio 93 (*). Demu´estrese la equivalencia. ¿Qu´e significado geom´etrico tiene el que s(t) = t − t0 ? Las curvas parametrizadas por longitud de arco tienen la ventaja de que est´an libres de los elementos no geom´etricos que suponen los cambios de velocidad. Todas las curvas regulares pueden reparametrizarse por la longitud de arco. Definamos formalmente lo que es una reparametrizaci´on: dos curvas x : I → Rn , y : J → Rn se dicen que son reparametrizaciones una de la otra cuando existe un difeomorfismo f : I → J tal que x = y ◦ f . Como f es difeomorfismo, es sobreyectiva, lo que asegura que la imagen de la curva reparametrizada es la misma que la de la original, es decir y(I) = (x ◦ f )(I). Por otra parte, por tener f una inversa diferenciable, f ′ nunca se anula (¿por qu´e?). Por tanto, y′ (t) = x′ (f (t))f ′ (t) y, si la curva x es regular, entonces podemos asegurar que su reparametrizaci´on y tambi´en es regular. N´otese que esta f´ormula nos dice que la velocidad a la que es recorrida la curva reparametrizada es la de la original escalada por el factor f ′ (t). Ejercicio 94 (*). Dar un ejemplo que muestre que el que dos curvas tengan la misma imagen no implica que sean equivalentes por reparametrizaciones. Ejercicio 95. Diremos que dos curvas son equivalentes cuando una es una reparametrizaci´on de la otra. Demu´estrese as´ı tenemos una relaci´on de equivalencia en el conjunto de curvas diferenciables. Ejercicio 96 (*). Se consideran las parametrizaciones de la circunferencia dadas por x(θ) = ¯ (t) = x



(cos θ, sin θ),  1 − t2 2t . , 1 + t2 1 + t2

Encu´entrese un cambio de par´ ametros que las relacione. Ejercicio 97. Dada la concoide de Nicomedes en coordenadas polares: r = a + c, a, c 6= 0, dib´ ujese y encu´entrese sus ecuaciones en coordenadas cos θ rectangulares.

48

CAP´ITULO 3. CURVAS REGULARES. CURVATURA.

Teorema 2. Dada una curva regular, existe siempre una reparametrizaci´on de la misma por la longitud de arco. Demostraci´ on. La funci´on longitud de arco s : I → R es una funci´on diferenciable cuya derivada es s′ (t) = kx′ (t)k > 0. Por tanto define un difeomorfismo sobre su imagen J = s(I). La curva y : J → Rn dada por y = x ◦ s−1 es una reparametrizaci´on de x y, debido a la regla de la cadena, y′ (ξ) =x′ (s−1 (ξ))(s−1 )′ (ξ) 1 =x′ (s−1 (ξ)) ′ −1 s (s (ξ)) ′ −1 x (s (ξ)) , = ′ −1 kx (s (ξ))k

(3.2)

con lo que claramente y′ tiene norma uno, es decir, es una curva parametrizada por la longitud de arco. Intentaremos descargar un poco la notaci´on: la inversa de la longitud de arco s = s(t) se llamar´a simplemente t = t(s). Esto conlleva la confusi´on de denotar por las letras s y t simult´aneamente a funciones y variables, pero esperamos que est´e claro por el contexto de qu´e se trata en cada caso. M´as a´ un, denotaremos por el mismo s´ımbolo x tanto a la curva como a su reparametrizaci´on. Con esta notaci´on, las ecuaciones (3.2) quedar´an como dx dt dx 1 dx 1 dx = = = . ds dt ds dt ds/dt dt kdx/dtk

(3.3)

En adelante, cuando llamemos s al par´ametro de una curva, supondremos que est´a parametrizada por la longitud de arco. En general, no es posible calcular expl´ıcitamente una reparametrizaci´on al arco: por una parte no suele poder calcularse expl´ıcitamente s(t) y a´ un menos su inversa. No obstante, su existencia te´orica es de gran ayuda para clarificar definiciones y demostraciones, como enseguida veremos. Ejercicio 98. Reparametr´ıcese la curva x(t) = (a cos t, a sin t, bt) por la longitud de arco. Ejercicio 99. Reparametr´ıcese la espiral logar´ıtmica x(t) = aebt (cos t, sin t) por la longitud de arco. Ejercicio 100. Reparametr´ıcese la cicloide por la longitud de arco en la zona comprendida entre dos puntos singulares, excluyendo los mismos.

´ DE CURVATURA 3.2. DEFINICION Ejercicio 101. Dada la curva x(t) = ( 21 t, 2t1 , por la longitud de arco.

49 √

2 2

log t), t > 0, parametr´ıcese

Ejercicio 102. Demu´estrese que la longitud de una curva dada por una funci´on y = f (x) con x ∈ [a, b] est´ a dada por Z bp L= 1 + f ′ (x)2 dx. a

Ejercicio 103. Demu´estrese que la longitud de una curva dada en coordenadas polares por ρ = ρ(θ) est´ a dada por Z bp L= ρ2 + (ρ′ )2 dθ. a

Ejercicio 104 (*). Parametr´ıcese la elipse de ecuaci´ on 4x2 + 9y 2 = 1 usando funciones trigonom´etricas y calc´ ulese num´ericamente, usando Sage, una reparametrizaci´on de la misma por la longitud de arco. Dib´ ujese los puntos de la misma equiespaciados por la longitud de arco, como comprobaci´ on de la correcci´on de los c´alculos.

Ejercicio 105 (*). Dib´ ujense los puntos de la hip´erbola xy = 1 que cortan a las rectas y = 10x e y = x/10 en el primer cuadrante, junto con la secuencia de puntos sobre la hip´erbola que los une y que se hayan equiespaciados sobre la hip´erbola a distancia igual a 0,03 unidades.

3.2.

Definici´ on de curvatura

El vector tangente o velocidad x′ (s) de una curva parametrizada por el arco es un vector unitario y la velocidad con la que ´este var´ıa indica cu´anto se est´a curvando una curva. Con m´as precisi´on, definimos pues el vector curvatura k(s) = x′′ (s) y la funci´on curvatura est´a dada por k(s) = kk(s)k. Ejemplo 2. Con lo que sabemos hasta ahora, s´ olo podemos calcular la curvatura de curvas parametrizadas al arco. Por ejemplo, sea v un vector unitario y sea x(s) = p + sv la recta que pasa por p y tiene vector director v. Como x′ (s) = v, es kx′ (s)k = 1, con lo que la recta est´ a efectivamente parametrizada por la longitud de arco. Su vector de curvatura es k′′ (s) = 0 y su funci´ on curvatura es pues k(s) = 0, como cab´ıa esperar. Observemos que, derivando la ecuaci´on   dx dx , =1 ds ds

CAP´ITULO 3. CURVAS REGULARES. CURVATURA.

50 obtenemos que

2



d2 x dx , ds2 ds



= 0,

es decir, la aceleraci´on de una curva parametrizada por el arco es perpendicular a la velocidad.

Ejercicio 106. Esta afirmaci´on es, intuitivamente, evidente si se observa que los vectores tangentes recorren una curva sobre la circunferencia o esfera unidad. Expl´ıquese por qu´e.

3.2.1.

Curvatura de una curva plana

Llamemos t(s) = x′ (s). Este vector, si la curva es plana, puede ser completado con un vector unitario n(s) de forma que {t(s), n(s)} formen una base ortonormal positivamente orientada. Como t(s) = (x′ (s), y ′ (s)), vemos que el vector normal es n(s) = (−y ′ (s), x′ (s)). La interpretaci´on geom´etrica del producto escalar hn(s), x′′ (s)i es la proyecci´on del vector de curvatura sobre el vector normal. Como ambos vectores son paralelos y n(s) es unitario, es

| hn(s), x′′ (s)i | = kx′′ (s)k. El signo de hn(s), x′′ (s)i es positivo si x′′ (s) tiene la misma direcci´on que n(s) y negativo en caso contrario. Como x′′ (s) es la velocidad con la que se mueve el vector tangente t(s) y el vector normal completa el tangente formando una base ortonormal positivamente orientada, es decir, siempre se encuentra a la izquierda del vector tangente, deducimos que hn(s), x′′ (s)i es positivo cuando la curva est´a girando “a izquierdas” respecto al vector de velocidad y negativo cuando gira “a derechas” (v´ease el siguiente dibujo).

´ DE CURVATURA 3.2. DEFINICION

51 dt/ds t(s)

k(s) >0

n(s)

k(s) <0 t(s)

dt/ds

n(s)

En adelante es a esta funci´on a la que llamaremos curvatura, cuando la curva es plana: ′ x (s) x′′ (s) ′′ ′ ′′ ′ ′′ . k(s) = hn(s), x (s)i = −y (s)x (s) + x (s)y (s) = ′ y (s) y ′′ (s) Por tanto, la curvatura de una curva plana tiene m´as informaci´on que la de una curva en dimensi´on mayor, puesto que su signo nos informa sobre el sentido de giro de la curva. El vector de curvatura se puede escribir como

dt = k(s)n(s), ds siendo ´esta la primera de las f´ormulas de Frenet, que ser´an el objeto del siguiente cap´ıtulo. En el caso de que la curva sea espacial, seguir´a siendo k(s) = kx′′ (s)k. k(s) =

Ejercicio 107. Parametr´ıcese la circunferencia de radio r por la longitud de arco y calc´ ulese su vector de curvatura y su curvatura.

Ejercicio 108. Si una curva x = x(s) se reparametriza haciendo s¯ = s¯(s), demu´estrese que la curva x(¯ s)) estar´a parametrizada por la longitud de arco si y s´ olo si s¯ = ±s + s0 .

CAP´ITULO 3. CURVAS REGULARES. CURVATURA.

52

Ejercicio 109 (*). ¿Cu´al es el efecto sobre el recorrido de la curva de la reparametrizaci´on s¯ = −s + s0 ? ¿Y sobre el vector de la curvatura? ¿Y sobre la curvatura? (dist´ınganse los casos plano y espacial)

3.2.2.

Curvatura de una curva no parametrizada por la longitud de arco

Nuestro siguiente objetivo es calcular la curvatura de una curva no parametrizada por la longitud de arco sin necesidad de reparametrizarla, es decir, dada la curva x = x(t) calcular k(s(t)) (que llamaremos, abusando de la notaci´on, k(t). Esto se puede hacer sin necesidad de conocer expl´ıcitamente la reparametrizaci´on por el arco gracias a la regla de la cadena y el teorema de la funci´on inversa en una variable. F´ ormula de la curvatura de una curva plana no parametrizada por el arco. Derivando la ecuacion (3.3) obtenemos que   d2 x d dx dt = ds2 ds dt ds   d dx dt dx d2 t = + ds dt ds dt ds2   2 d2 x dt dx d2 t = 2 + dt ds dt ds2

(3.4)

En el lado derecho de esta ecuaci´on conocemos todos los t´erminos salvo dt2 /ds2 , que podr´ıamos tambi´en calcular. No obstante, si s´olo queremos calcular k(t) realmente no es necesario el c´alculo expl´ıcito de dt2 /ds2 . Basta con calcular   dx d2 x , k(s(t)) = det ds ds2 !  2 dx d✚2✚ dx dt d2 x dt t = det + ✚ 2 , dt ds dt2 ds dt ds ✚     3 dt dx d2 x (3.5) , = det dt dt2 ds   dx d2 x det , dt dt2 , =

3

dx



dt

´ DE CURVATURA 3.2. DEFINICION

53

donde hemos usado que 1 1 dt = = ′ . ds ds/dt kx (t(s))k F´ ormula de la curvatura de una curva espacial no parametrizada por el arco. En el caso de que la curva est´e parametrizada por el arco y sea espacial, sabemos que su curvatura est´a definida por

2

d x

k(s) =

ds2 .

Podr´ıamos usar directamente la expresi´on (3.4), pero es m´as eficiente observar que, como x′ (s) es un vector unitario y ya hemos observado que es perpendicular a x′′ (s), entonces

2



2



d x dx d2 x dx



sin(π/2) = d x ,

= ×

ds2

ds2 ds ds2 ds lo que nos va a permitir hacer una simplificaci´on an´aloga a la que hemos hecho en el caso plano, ahorr´andonos el c´alculo de dt2 /ds2 .

Ejercicio 110 (*). Demu´estrese que la curvatura de una curva espacial viene dada por la f´ ormula

dx d2 x

dt × dt2 k(t) =

3

dx



dt

Ejemplo 3. Calculemos la curvatura de la curva x(t) = (cos t, sin 2t): var(’t’) x(t) = (cos(t), sin(2*t)) dx = x.diff(t); d2x = x.diff(t,t) result = matrix([dx,d2x]).det() / dx.norm()^3 curv = result.full_simplify().trig_expand() que da como resultado una funci´on curvatura  2 2 cos (t)3 − 3 cos (t) k(t) = −  3 2 2 2 2 4 2 cos (t) − 1 + |sin (t)|

54

CAP´ITULO 3. CURVAS REGULARES. CURVATURA.

Ejercicio 111 (*). Dib´ ujese la funci´on de curvatura y el vector curvatura a lo largo de la curva del ejercicio anterior. Interpr´etese geom´etricamente el signo de la curvatura, sus m´aximos, m´ınimos y ceros. Ejercicio 112. La nefroide es la curva de ecuaciones x(t) = a(3 cos t − cos 3t), y(t) = a(3 sin t − sin 3t). Dib´ ujese y calc´ ulese su curvatura en los puntos regulares. Ejercicio 113 (*). Encu´entrese la ecuaci´ on param´etrica de la curva x3 + y 3 − 3axy = 0 cort´andola con el haz de rectas que pasa por el origen. ¿Por qu´e funciona este procedimiento? Calc´ ulese la curvatura de la curva. Dib´ ujese. Esta curva se llama folium de Descartes. Ejercicio 114. Parametr´ıcese la elipse de ecuaci´ on x2 /a2 +y 2 /b2 = 1 usando funciones trigonom´etricas y calc´ ulese su curvatura en los v´ertices. Ejercicio 115. Encu´entrese el vector de curvatura y la curvatura de la curva x(t) = (t2 + 3t − 1, t2 + t + 1, −t2 + 5). Ejercicio 116. Calc´ ulese la curvatura de una curva dada por un grafo y = f (x). Ejercicio 117 (*). Demu´estrese que el vector normal a una curva plana satisface la ecuaci´ on dn = −kt. ds Ejercicio 118. Demu´estrese que una curva regular tiene curvatura nula si y s´ olo si es una recta.

3.2.3.

Invariancia de la curvatura

Un movimiento es una aplicaci´on M : Rn → Rn de la forma M (p) = Rp + b siendo R una matriz ortogonal, es decir, tal que RR⊤ = Id y det R > 0 y b ∈ Rn . Si transformamos una curva mediante un movimiento entonces la curvatura de la curva transformada coincide con el de la original. Ejercicio 119 (**). Demu´estrese este resultado. N´otese que el resultado es cierto tanto para la curvatura de curvas espaciales (sin signo) como para la curvatura con signo de curvas planas.

´ DE CURVATURA 3.2. DEFINICION

55

Ejercicio 120 (*). Si en lugar de movimientos consideramos isometr´ıas (es decir, permitimos que el det R sea negativo), ¿sigue siendo cierto el resultado del ejercicio anterior? (Indicaci´ on: consid´erese por separado los casos n = 2 y n = 3). Adem´as de la invariancia por movimientos, la expresi´on de la curvatura dada por la expresi´on (3.5) es invariante por reparametrizaciones que no invierten la orientaci´on. Con m´as precisi´on, si t = f (t¯) es un cambio de ¯ est´a dada por x ¯ (t¯) = x(f (t¯)) entonces kx¯ (t¯) = kx (f (t¯)). par´ametros y x Ejercicio 121 (*). Demu´estrese este resultado. N´otese que se trata tan s´ olo de una aplicaci´ on de la regla de la cadena y que, observando que cuando un determinante tiene dos columnas proporcionales entonces se anula, los c´alculos necesarios se pueden minimizar. Ejercicio 122. ¿Y qu´e sucede en el caso de que se admitan cambios de par´ ametro que inviertan la orientaci´ on?

3.2.4.

Circunferencia osculadora

A continuaci´on describiremos otros m´etodos para estudiar la curvatura de una curva. En primer lugar, supongamos que la curva es plana. Si queremos estudiar la curvatura en un punto x(s0 ), podemos tomar dos puntos cercanos, digamos x(s1 ) y x(s2 ) y considerar la circunferencia determinada por los tres puntos. Veamos un ejemplo, usando Sage, en el que calculamos la circunferencia que pasa por tres puntos concretos de la curva x(t) = (1 − t2 , t3 − t): var(’t’) X = vector([1-t^2, t^3-t]).function(t) var(’x x0 x1 x2 y y0 y1 y2 A B C ’) circ = x^2 + y^2 - 2*A*x - 2*B*y + C solus = solve([circ(x=x0, y=y0),circ(x=x1, y=y1),circ(x=x2, y=y2)],\ [A,B,C], solution_dict= True) centro = vector([solus[0][A],solus[0][B]]) r_square = centro[0]^2 + centro[1]^2 - solus[0][C] t0 = 0.3; t1 = 0.6; t2 = -0.2 c =

centro.subs(x0 = X(t0)[0], y0 = X(t0)[1], \ x1 = X(t1)[0], y1 = X(t1)[1], \

56

CAP´ITULO 3. CURVAS REGULARES. CURVATURA.

x2 = X(t2)[0], y2 = X(t2)[1]) r = r_square.subs(x0 = X(t0)[0], y0 = X(t0)[1], x1 = X(t1)[0], y1 = X(t1)[1], x2 = X(t2)[0], y2 = X(t2)[1]).sqrt( curva = parametric_plot(X, (t, -1.5, 1.5)) circulo = circle(c, r) puntos = point([X(t0), X(t1), X(t2)], color = ’red’) dib = curva + circulo + puntos dib.set_aspect_ratio(1)

1.5

1

0.5

-1

-0.5

0.5

1

-0.5

-1

-1.5

Dado que que las circunferencias son curvas de curvatura constante y determinada por el radio, el l´ımite l´ıms1 ,s2 →s0 r(s0 , s1 , s2 ) = r(s0 ), de existir, ser´ıa otra forma posible de medir la curvatura de la curva en un punto. A la circunferencia correspondiente se la denomina circunferencia osculatriz. Para ello consideremos la ecuaci´on del c´ırculo que pasa por los tres puntos x(si ), que ser´a de la forma hx − c(s0 , s1 , s2 ), x − c(s0 , s1 , s2 )i = r(s0 , s1 , s2 )2 , donde x denota un punto gen´erico del plano. Si llamamos f (s) = hx(s) − c(s0 , s1 , s2 ), x(s) − c(s0 , s1 , s2 )i − r(s0 , s1 , s2 )2 , resulta que f es una funci´on diferenciable que se anula en s0 , s1 y s2 . Por el teorema de Rolle, existen puntos intermedios ξ0 , ξ1 tales que f ′ (ξ0 ) = f ′ (ξ1 ) =

´ DE CURVATURA 3.2. DEFINICION

57

0. Una nueva aplicaci´on del teorema de Rolle asegura que existe un nuevo punto intermedio η0 tal que f ′′ (η0 ) = 0. Dado que es 1 ′ f (s) = hx′ (s), x(s) − c(s0 , s1 , s2 )i , 2 1 ′′ f (s) = hx′′ (s), x(s) − c(s0 , s1 , s2 )i + hx′ (s), x′ (s)i 2 = hx′′ (s), x(s) − c(s0 , s1 , s2 )i + 1,

(3.6)

hx′ (s), x(s) − c(s)i = 0, hx′′ (s), x(s) − c(s)i + 1 = 0.

(3.7)

c(s) = x(s) + λn(s).

(3.8)

tomando el l´ımite cuando s1 , s2 → s0 y, por simplificar la notaci´on, denotando s0 por la letra s se tiene que ξ0 , ξ1 → s y η0 → s, con lo que

De la primera de estas dos u ´ltimas ecuaciones se deduce, dado que suponemos que la curva es plana, que

La sustituci´on de (3.8) en la segunda ecuaci´on de (3.7) nos lleva a hx′′ (s), −λn(s)i + 1 = 0, con lo que λ= Por tanto tenemos que

1 1 = . hx′′ (s), n(s)i k(s)

c(s) = x(s) +

1 n(s), k(s)

(3.9)

La curva recorrida por los centros de las circunferencias osculatrices se llama evoluta de la curva x(s). El radio de la curva es 1 1 n(s) = . r(s) = kx(s) − c(s)k = k(s) |k(s)|

Por tanto el radio de la circunferencia osculatriz no aporta una nueva medida de la curvatura de la curva. N´otese que si en un punto la curvatura es nula entonces la circunferencia osculatriz tiene radio infinito. En el caso de que la curva no est´e parametrizada por la longitud de arco, es posible calcular la evoluta. Para ello definimos el vector tangente a la curva como x′ (t) t(t) = ′ . kx (t)k El vector normal n(t) est´a u ´nicamente determinado por ser {t, n} una base ortonormal positivamente orientada y haciendo uso de (3.5), se calcula c(t).

58

CAP´ITULO 3. CURVAS REGULARES. CURVATURA.

Ejercicio 123. Calc´ ulese la evoluta de la cicloide. Dib´ ujese. Ejercicio 124. Escr´ıbase el c´odigo necesario para obtener el dibujo de la evoluta de la elipse de ecuaci´ on x2 /4 + y 2 /9 = 1, que incluimos a continuaci´on: 3

2

1

-2

-1

1

2

-1

-2

-3

¿Sabes interpretar el dibujo? ¿Eres capaz de visualizar las circunferencias osculatrices? Ejercicio 125 (*). Demu´estrese que las rectas tangentes a la evoluta de una curva son perpendiculares a la curva inicial. Esto se suele expresar diciendo que la evoluta es envolvente de las rectas normales a la curva.

3.3.

Envolvente de una familia de curvas planas parametrizadas

Supongamos que tenemos una familia uniparam´etrica de curvas, x = x(p, t), entendiendo que t es el par´ametro usual de la curva y p describe las distintas curvas de la familia. Ejemplo 4. Consideremos la familia de circunferencias x(p, t) = (p+p cos t, p+ p sin t), t ∈ [0, 2π], p ∈ R. Se trata de circunferencias con centro en la bisectriz y = x y tangentes a los ejes coordenados: var(’t p’) X = vector([p + p*cos(t), p + p*sin(t)]) a, b = -10, 10 dibs = [parametric_plot(X(p=var), (t, 0, 2*pi))\ for var in srange(a, b, 0.5)] circs = sum(dibs)

3.3. ENVOLVENTE DE UNA FAMILIA DE CURVAS PLANAS PARAMETRIZADAS59

15

10

5

-20

-15

-10

-5

5

10

15

-5

-10

-15

-20

Pediremos que las curvas sean todas regulares y que adem´as la variaci´on con respecto a p tambi´en sea regular, es decir: xt 6= 0 6= xp . Una envolvente de esta familia de curvas es una curva con la propiedad de ser tangente en cada punto a alguna de las curvas de la familia. Para facilitar las cosas, exigiremos que exista una aplicaci´on diferenciable u 7→ (p(u), t(u)) de forma que la envolvente sea e(u) = x(p(u), t(u)). Debe ocurrir que e′ (u) = xp (p(u), t(u))p′ (u) + xt (p(u), t(u))t′ (u) debe, por ser tangente a la curva, ser proporcional a xt . De aqu´ı concluimos que xp y xt deben ser linealmente dependientes. Rec´ıprocamente, si xp y xt son dependientes en los puntos de la forma (p(u), t(u)), como suponemos que xt 6= 0 podemos escribir xp = αxt y entonces e′ (u) = xp p′ + xt t′ = (t′ + αp′ )xt y deducimos que e es una envolvente. Hemos demostrado el siguiente resultado: Teorema 3. Con las hip´ otesis de regularidad enunciadas m´as arriba, la condici´ on necesaria y suficiente para que e(u) = x(p(u), t(u)) sea una envolvente es que los vectores xp y xt sean linealmente dependientes en los puntos de la forma (p(u), t(u)) o, equivalentemente, que  det xt xp = 0.

Ejemplo 5. El teorema que acabamos de ver proporciona un procedimiento para calcular envolventes. Por ejemplo, consideremos la familia de elipses x(p, t) = (p cos t, (1 − p) sin t).

CAP´ITULO 3. CURVAS REGULARES. CURVATURA.

60

1

0.5

-1

-0.5

0.5

1

-0.5

-1

Figura 3.1: La familia de elipses (en azul) junto con su envolvente. Entonces xp = (cos t, − sin t) y xt = (−p sin t, (1 − p) cos t), y el determinante  det xt xp = cos2 t − p,

con lo que la envolvente es e(u) = x(cos2 t, t) = (cos3 t, sin3 t). Esta curva se llama astroide.

Ejercicio 126. Compru´ebese que la astroide mencionada m´as arriba coincide con los puntos del plano que satisfacen la ecuaci´ on x2/3 + y 2/3 = 1. Ejercicio 127. Calc´ ulense las envolventes de la familia de c´ırculos de radio 1 con centro en el eje OX. Ejercicio 128. Demu´estrese que las envolventes de la familia de√c´ırculos x(p, t) = (p + a cos t, −p + a sin t) con a > 0 son las rectas x + y ± 2a = 0. Ejercicio 129. Calc´ ulese la envolvente de los c´ırculos x(p, t) = (p+p cos t, p+ p sin t). Ejercicio 130. En cada caso, calcul´ese la envolvente de la familia de curvas: 1. x(p, t) = (t, p − p2 t) 2. x(p, t) = (t, p2 − pt)

3.4. ECUACIONES IMPL´ICITAS DE CURVAS EN EL PLANO

61

3. x(p, t) = (t, p2 − p2 t) 4. x(p, t) = (t, 27p − p3 t2 ) 5. x(p, t) = (t, p3 + pt) Ejercicio 131. En cada caso, calcul´ese la envolvente de la familia de curvas: 1. x(p, t) = (t, p − p2 t) 2. x(p, t) = (p + t, pt2 ) 3. x(p, t) = (p + t2 , pt) 4. x(p, t) = (t2 + p, pt) Ejercicio 132. Calcular la envolvente de la familia de curvas x(p, t) = (t2 − p2 + 2t, 2tp − 2p). La curva que resulta se llama deltoide. Ejercicio 133. Calc´ ulese la envolvente de los c´ırculos que tienen sus centros sobre el c´ırculo unidad y son tangentes al eje OX. La curva resultante se llama nefroide. Ejercicio 134. Calc´ ulese la envolvente de la familia de rectas de ecuaciones x(p, t) = ((1−t) cos p+t cos np, (1−t) sen p+t sen np) siendo n 6= 1 un n´ umero real fijo. Las curvas resultantes son epicicloides, si n > 1 y hipocicloides, si n < 1. En el caso particular n = 2 la curva se llama cardioide. Ejercicio 135. Demostrar que la envolvente de la familia de c´ırculos x(p, t) = (cos p + sen p cos t, sen p sen t) es una elipse. Ejercicio 136. Un c´ırculo de radio 1 rueda a lo largo de la recta y = −1. Consid´erese el di´ametro del c´ırculo, que fijamos como el eje OX cuando el centro est´ a en el origen. Al rodar el c´ırculo el di´ametro describe una familia de rectas. Demu´estrese que la envolvente de dicha familia es una cicloide.

3.4.

Ecuaciones impl´ıcitas de curvas en el plano

Comencemos con algunas ideas intuitivas: una curva x = x(t) es un objeto geom´etrico descrito por un grado de libertad, el par´ametro temporal t. Una idea alternativa para generar un objeto con un grado de libertad es partir del plano, cuyos puntos x = (x, y), est´an descritos por dos coordenadas y eliminar una de las libertades forzando el que se satisfaga una ecuaci´on f (x, y) = 0.

(3.10)

62

CAP´ITULO 3. CURVAS REGULARES. CURVATURA.

La curva ser´a el conjunto de puntos que satisfacen una tal la ecuaci´on. Pero no cualquier funci´on es adecuada. Desde luego, pediremos que la funci´on sea diferenciable, pero esto no es suficiente: por ejemplo, podr´ıa ser f ≡ 0, con lo que la curva ser´ıa todo el plano, cosa que no deseamos o f ≡ 1, entonces ning´ un punto satisface la ecuaci´on y la curva ser´ıa vac´ıa, cosa que tampoco es deseable, as´ı que vemos que debemos exigir m´as condiciones a la ecuaci´on f. El que una variable determine la otra se puede expresar matem´aticamente pidiendo que exista una funci´on y = y(x) (o bien x = x(y)) que cumpla que f (x, y) = 0 si y s´olo si y = y(x). Esta funci´on es la que ’despeja’ y en funci´on de x. Si existe, deberemos tener que f (x, y(x)) ≡ 0. Si derivamos esta identidad respecto a x obtenemos que ∂f ∂f dy + = 0. ∂x ∂y dx De aqu´ı se deduce que ∂f dy ∂x = − ∂f . dx ∂y

Como queremos que dy/dx sea diferenciable, es claro que una condici´on necesaria para que exista una tal y = y(x) es que ∂f 6= 0. ∂y Esta condici´on tambi´en se puede interpretar geom´etricamente. Para ello debemos entender geom´etricamente el vector gradiente ∇f =



∂f ∂f , ∂x ∂y

⊤

.

Si asumimos anticipadamente que las soluciones de la ecuaci´on (3.10) son (localmente) una curva parametrizada regular de la forma x = (x(t), y(t)), tendremos que f (x(t), y(t)) ≡ 0. Derivando esta ecuaci´on obtenemos ∂f ′ ∂f ′ x (t) + y (t) = 0, ∂x ∂y

3.4. ECUACIONES IMPL´ICITAS DE CURVAS EN EL PLANO

63

o, equivalentemente, h∇f (x(t)), x′ (t)i = 0. Es decir, el gradiente de f en los puntos de la curva es perpendicular a la misma. Por tanto, la condici´on ∂f 6= 0 se puede interpretar como la ∂y condici´on de que la curva no tenga una normal paralela al eje de abscisas, es decir, que la curva, vista como una funci´on desde el eje OX, no tenga pendiente infinita. Si tuviese pendiente infinita, la funci´on y = y(x) no ser´ıa una funci´on diferenciable. El teorema de la funci´on impl´ıcita nos asegura que esta condici´on necesaria para que f (x, y) = 0 defina una curva es tambi´en suficiente. Teorema 4 (Teorema de la funcion implicita en el plano). Sea f una funci´ on ∂f diferenciable definida en un abierto del plano y supongamos que ∂y (x0 , y0 ) 6= 0 y f (x0 , y0 ) = 0. Entonces existe un entorno abierto U ⊂ R2 del punto (x0 , y0 ), un entorno abierto V de x0 y una funci´on diferenciable y = y(x) definida en V tales que {(x, y) ∈ U : f (x, y) = 0} = {(x, y(x)) ∈ R2 : x ∈ V }. Puede hacerse un enunciado an´alogo en el caso de que sea ∂f (x0 , y0 ) 6= 0. ∂x En adelante diremos que un subconjunto no vac´ıo C ⊂ R2 es una curva plana definida impl´ıcitamente cuando existe una funci´on diferenciable f : U → R tal que C = {(x, y) ∈ U : f (x, y) = 0} y adem´as ∇f (x, y) 6= 0 para todo (x, y) ∈ C. Si C es una curva definida impl´ıcitamente, el teorema de la funci´on impl´ıcita nos asegura que cada punto de la misma tiene un entorno que se puede parametrizar usando un grafo y, en particular, se puede ver localmente como una curva parametrizada regular. La noci´on de curva parametrizada que hemos empleado en las secciones anteriores est´a relacionada, pero no coincide con la de curva definida impl´ıcitamente. V´eamoslo en los siguiente ejercicios: Ejercicio 137 (*). Sea x : I ⊂ R → R2 una curva parametrizada regular. Demu´estrese que cada t0 tiene un entorno abierto V ⊂ R de forma que {x(t) : t ∈ V } es un grafo. (Indicaci´ on: una de las dos funciones coordenadas de la curva es invertible en t0 .) Ejercicio 138 (*). Demu´estrese que un grafo es una curva definida impl´ıcitamente. Ejercicio 139. Sea C = {(x, y) : f (x, y) = (x − y)2 = 0}. ¿Es C una curva definida impl´ıcitamente?

64

CAP´ITULO 3. CURVAS REGULARES. CURVATURA.

Ejercicio 140 (*). Demu´estrese que la curva parametrizada dada por x(t) = (cos t, sen 2t) con t ∈ R no es una curva impl´ıcita. (Indicaci´ on: si existe una funci´on f (x, y) = 0 que la defina, su gradiente en el origen deber´ıa ser perpendicular a los dos vectores independientes tangentes a las dos ramas de la curva que pasan por el origen). Por otra parte, no toda curva definida impl´ıcitamente se puede parametrizar, como muestra el siguiente ejercicio: Ejercicio 141. Encu´entrese una curva definida impl´ıcitamente cuya imagen no sea una curva parametrizada. (Indicaci´ on: b´ usquese una ecuaci´ on cuya soluci´on tenga m´as de una componente conexa). Sin embargo, si se requiere que la curva definida impl´ıcitamente sea conexa, se puede demostrar que siempre se puede parametrizar. Ejercicio 142. En los siguientes ejemplos, est´ udiese si la ecuaci´ on f (x, y) = 0 define impl´ıcitamente una curva: 1. ax + by + c, 2. x2 − y 2 , 3. x2 − y 2 − 9 4. xy − 1. 5. 5x2 − 4xy + 8y 2 + 18x − 36y + 9 = 0. Dib´ ujese el lugar de los ceros y anal´ıcese en qu´e zonas de los lugares resultantes la curva puede escribirse localmente como un grafo, construy´endolo expl´ıcitamente. Ejemplo: c´ onicas Sea f un punto fijo, que llamaremos foco, y sea D una recta fija, que llamaremos directriz. Fijemos tambi´en una constante positiva e que llamaremos excentricidad. La c´onica dada por estos elementos es el lugar geom´etrico de los puntos del plano x que satisfacen que la raz´on de la distancia d(x, f ) a la distancia de x a la recta D es e. Si e = 1 la c´onica se dice una par´abola, si e > 1 se trata de una hip´erbola y si e < 1 diremos que la c´onica es una elipse. Ejercicio 143. Encu´entrese la ecuaci´ on de la c´onica de excentricidad 2 cuya directriz es la recta 4x − 3y + 2 = 0 y cuyo foco es (1, −1). Ejercicio 144. Demu´estrese que la ecuaci´ on de una c´onica, seg´ un se ha definido, es de la forma f (x, y) = ax2 + 2hxy + by 2 + 2gx + 2f y + c = 0.

(3.11)

Ejercicio 145. Escr´ıbase la ecuaci´ on de una c´onica de forma matricial:   x  f (x, y) = x y 1 C y  , 1

3.4. ECUACIONES IMPL´ICITAS DE CURVAS EN EL PLANO

65

siendo C una matriz sim´etrica. Demu´estrese que, si llamamos a las columnas de C = (a0 , a1 , a2 ) y es X = (x, y, 1)⊤ es ⊤ ∇f = 2(a⊤ 0 X, a1 X).

Ejercicio 146 (*). Usando el ejercicio anterior, concl´ uyase que si C es re2 gular y ab − h 6= 0, entonces la c´onica es una curva definida impl´ıcitamente.

3.4.1.

C´ onicas degeneradas

En adelante, supondremos que la matriz de las c´onicas C 6= 0. Una c´onica se dice degenerada cuando es un par de rectas o una recta doble. Para ver cu´ando una c´onica es degenerada, partimos de su ecuaci´on (3.11) y, suponemos que a 6= 0. Si a = 0, se argumenta de forma similar con b, y si a = b = 0, la c´onica es un producto de rectas si y s´olo si se trata de la la multiplicamos por a y completamos el cuadrado como sigue: a2 x2 + 2ax(hy + g) + (hy + g)2 = −aby 2 − 2af y − ac + (hy + g)2 = (h2 − ab)y 2 + 2(hg − af )y + g 2 − ac de forma que (ax + hy + g) = ±

p

(h2 − ab)y 2 + 2(hg − af )y + g 2 − ac.

Por tanto, para que la c´onica se descomponga en dos rectas, el polinomio (h2 − ab)y 2 + 2(hg − af )y + g 2 − ac ha de ser un cuadrado perfecto. Ejercicio 147. Compru´ebese que la condici´ on para que el polinomio αy 2 + 2βy + γ sea un cuadrado perfecto es que su discriminante se anule, es decir, β 2 − αγ = 0. Teniendo en cuenta el ejercicio anterior, la condici´on para que la c´onica se descomponga en dos rectas es que 0 = (hg − af )2 − (h2 − ab)(g 2 − ac) = a(af 2 − 2hgf + h2 c + bg 2 − bca) = −a det C. Por tanto, una c´onica es degenerada si y s´olo si det C = 0. Ejercicio 148. Clasifica las c´onicas con las ecuaciones siguientes y determina, cuando existan, su centro y as´ıntotas: 1. X 2 + Y 2 − 4X + 4Y + 6 = 0

66

CAP´ITULO 3. CURVAS REGULARES. CURVATURA. 2. 3X 2 − 2XY − Y 2 − 2Y = 0 3. X − XY + Y = 0 4. 5 + 5X 2 + 5Y 2 + 4X + 4Y + 4XY = 0 5. 3 − 5X 2 − 3Y 2 − 4X + 4Y − 4XY = 0 6. 2 + X 2 + Y 2 − 4X − 2Y − 2XY = 0 7. 2X 2 − 13XY + 15Y 2 + 5X − 11Y + 7 = 0

Ejercicio 149. Clasifica af´ınmente, seg´ un los valores de t, la c´onica de ecuaci´on 1 − 2X + 2Y + tX 2 + 2XY + tY 2 = 0. Ejercicio 150. Encu´entrese el eje, el v´ertice y la tangente en el v´ertice de la par´ abola de ecuaci´ on 4x2 − 12xy + 9y 2 − 3x − 2y + 4 = 0. Mostrar un sistema de coordenadas en el que tenga ecuaci´ on can´onica y encontrar la misma. Ejercicio 151. Encu´entrese los ejes, centro (si existe), v´ertices, y focos de las siguientes c´onicas: 8x2 −4xy+5y 2 −36x+18y+9 = 0, 2x2 −4xy−y 2 +7x− 2y + 3 = 0., 5x2 − 4xy + 8y 2 + 18x − 36y + 9 = 0, x2 − 4xy + 4y 2 + 5y − 9 = 0, 2x2 + 3xy − 2y 2 − 11x − 2y + 12 = 0, x2 + 2xy + y 2 − 2x − 2y − 3 = 0, x2 −4xy−2y 2 −2x+7y−3 = 0. Compru´ebese que los ejes son perpendiculares. Mostrar un sistema de coordenadas en el que tengan ecuaci´ on can´onica y encontrar las mismas. Ejercicio 152. Demu´estrese que los ejes de una hip´erbola bisecan las as´ıntotas. Encu´entrese las as´ıntotas de la hip´erbola 2x2 − 3xy − 2y 2 + 3x − y + 8 = 0 y, a partir de ellas, sus ejes. Sabiendo las ecuaciones de los ejes, dar las ecuaciones de un cambio de coordenadas que la reduzca a su forma can´onica. Ejercicio 153. Escr´ıbase la ecuaci´ on de todas las hip´erbolas que tienen como as´ıntotas los ejes coordenados. Ejercicio 154. Encu´entrese la ecuaci´ on de la hip´erbola que tiene por as´ıntotas 2x − 3y + 1 = 0 y x + y − 3 = 0 y pasa por el punto (1, −2). Ejercicio 155. Demu´estrese que los centros de las hip´erbolas cuyas as´ıntotas son paralelas a los ejes coordenados y que pasan por los puntos (2, 5) y (3, 2) est´ an sobre la recta 3x − y − 4 = 0. Encu´entrese la ecuaci´ on de la hip´erbola de dicho conjunto que pasa por el punto (−2, 3).

3.4. ECUACIONES IMPL´ICITAS DE CURVAS EN EL PLANO

67

Ejercicio 156. Demu´estrese que los centros de las c´onicas de ecuaciones ax2 + 2dxy + by 2 + 2etx + 2f ty + c = 0 se encuentran sobre una recta que pasa por el origen de coordenadas. Ejercicio 157. ¿Para qu´e valores de c es x2 − 3xy + 2y 2 + x − y + c = 0 la ecuaci´ on de un par de rectas que se cortan? Encu´entrese dichas rectas. Ejercicio 158. Demu´estrese que toda hip´erbola se puede escribir de la forma (ax+by+c)(a′ x+b′ y+c′ ) = k, y que entonces ax+by+c = 0 y a′ x+b′ y+c′ = 0 son sus as´ıntotas.

3.4.2.

Curvatura de una curva definida impl´ıcitamente

Puesto que una curva definida impl´ıcitamente est´a dada localmente por parametrizaciones (grafos), tiene sentido definir su curvatura en un punto. El siguiente teorema muestra c´omo calcular la curvatura a partir de la funci´on que define impl´ıcitamente la curva. Antes introduzcamos la siguiente notaci´on: llamaremos J : R2 → R2 a la aplicaci´on lineal que gira un vector 90 grados en sentido antihorario, es decir, J(x, y)⊤ = (−y, x)⊤ . Denotaremos por H(f ) a la matriz Hessiana de f definida por   fxx fxy . H(f ) = fyx fyy Teorema 5. La curvatura de la curva definida impl´ıcitamente por la ecuaci´on f (x, y) = 0 y recorrida en el sentido que hace que ∇f tenga el mismo sentido que el vector normal n est´ a dada por k(x, y) = −

J(∇f )⊤ H(f )J(∇f ) . k∇f k3

Demostraci´ on. Sea x0 = (x0 , y0 ) un punto de la curva definida impl´ıcitamente y sea x = x(s) una parametrizaci´on por la longitud de arco definida en un entorno de x0 con x(s0 ) = x0 . Puesto que f (x(s), y(s)) ≡ 0, derivando esta ecuaci´on tenemos que fx x′ + fy y ′ = 0, con lo que obtenemos que (x′ , y ′ ) = ±

(−fy , fx ) . (fx2 + fy2 )1/2

(3.12)

CAP´ITULO 3. CURVAS REGULARES. CURVATURA.

68

Como suponemos que n = (−y ′ , x′ ) tiene el mismo sentido que (fx , fy ), el signo en la u ´ltima ecuaci´on debe ser negativo. Derivando nuevamente (3.12) se obtiene fxx (x′ )2 + 2fxy x′ y ′ + fyy (y ′ )2 + fx x′′ + fy y ′′ = 0. Por tanto k = x′ y ′′ − y ′ x′′ =

fx x′′ + fy y ′′ (fx2 + fy2 )1/2

fxx (x′ )2 + 2fxy x′ y ′ + fyy (y ′ )2 (fx2 + fy2 )1/2    ′  fxx fxy x ′ ′ x y y′ fyx fyy =− (fx2 + fy2 )1/2     fxx fxy −fy −fy fx fx fyx fyy J(∇f )⊤ H(f )J(∇f ) =− =− (fx2 + fy2 )3/2 k∇f k3 =−

(3.13)

Ejercicio 159. Calc´ ulese la curvatura de la elipse y de la hip´erbola de semiejes a y b en mediante los dos m´etodos que conocemos: parametrizando las curvas y mediante su ecuaci´ on impl´ıcita. Compru´ebese que se obtiene el mismo resultado. Ejercicio 160. Se denomina v´ertice de una curva a los puntos cr´ıticos de la curvatura. Encu´entrense los v´ertices de la elipse x2 /a2 + y 2 /b2 = 1. Ejercicio 161. Calc´ ulese la curvatura de la curva de ecuaci´ on y 2 (a2 − x2 ) = (x2 + 2ay − a2 )2 . Dib´ ujese para algunos valores de a. Ejercicio 162 (*). Demu´estrese que la curvatura de una curva definida impl´ıcitamente por la ecuaci´ on f (x, y) = 0 viene dada por   ∇f k = div , k∇f k donde el operador divergencia se define de forma general como n X ∂F , div(F ) = ∂xi i=1

siendo F : Rn → Rn una aplicaci´ on diferenciable.

Cap´ıtulo 4 Ecuaciones de Frenet 4.1.

Otra interpretaci´ on de la curvatura de una curva plana

Una curva plana se puede recuperar a partir de su curvatura k = k(s). En efecto, si t(s) es el vector tangente a la curva, dado que es un vector unitario, podremos escribirlo como t(s) = (cos(θ(s)), sin(θ(s)))

(4.1)

para una cierta funci´on diferenciable θ = θ(s). El que esta funci´on est´a definida en todo el dominio de la curvatura no es completamente evidente y es el contenido del siguiente teorema: Teorema 6. Sea I ⊂ R un intervalo conexo y sea f : I → S 1 una funci´ on diferenciable. Entonces existe una funci´on diferenciable θ : I → R tal que f (s) = (cos θ(s), sin θ(s)). Adem´ as, cualquier otra funci´on θ¯ : I → R con la misma propiedad difiere de θ en un m´ ultiplo entero de 2π. Demostraci´ on. Si f (s) = (u(s), v(s)) = (cos(θ(s)), sin(θ(s)), debe ser f ′ (s) = θ′ (s)(−v(s), u(s)), por lo que hf ′ (s), (−v(s), u(s))i = θ′ (s).

Esto demuestra que θ est´a u ´nicamente determinada salvo constantes y podemos recuperarla mediante la f´ormula Z s θ(s) = (v ′ (ξ)u(ξ) − v(ξ)u′ (ξ))dξ + θ0 s0

69

70

CAP´ITULO 4. ECUACIONES DE FRENET

estando dado θ0 por la condici´on f (s0 ) = (cos θ0 , sin θ0 ), lo que le determina u ´nicamente salvo m´ ultiplos enteros de 2π. Hemos de comprobar ahora que la funci´on θ que acabamos de definir cumple que (u(s), v(s)) = (cos(θ(s)), sin(θ(s)). Veamos para ello que la funci´on F (s) = (u(s) − cos θ(s))2 + (v(s) − sin θ(s))2 = 2 − 2(u cos θ + v sin θ) es id´enticamente cero, lo que es equivalente a demostrar que u cos θ + v sin θ = 1, (de hecho, pod´ıamos haber empezado por aqu´ı: dos vectores de la circunferencia unidad coinciden si y s´olo si su producto escalar es 1). Derivando (u cos θ + v sin θ)′ = u′ cos θ + v ′ sin θ − uθ′ sin θ + vθ′ cos θ = u′ cos θ + v ′ sin θ − u(v ′ u − vu′ ) sin θ + v(v ′ u − vu′ ) cos θ. Usemos ahora que, como u2 + v 2 = 1 es uu′ + vv ′ = 0, de forma que (u cos θ + v sin θ)′ = u′ cos θ + v ′ sin θ − (v ′ u2 − vuu′ ) sin θ + (vv ′ u − v 2 u′ ) cos θ = u′ cos θ + v ′ sin θ − v ′ (u2 + v 2 ) sin θ − u′ (v 2 + u2 ) cos θ = u′ cos θ + v ′ sin θ − v ′ sin θ − u′ cos θ = 0. Por tanto la funci´on F es constante, pero como F (s0 ) = 0, concluimos que F ≡ 0. El resultado que acabamos de demostrar es tambi´en cierto en el caso de que f sea u ´nicamente continua, pero la demostraci´on es m´as complicada. Ejercicio 163. Util´ıcese el teorema que acabamos de ver para encontrar θ(s) en el caso de que   2s 1 − s2 . , f (s) = 1 + s2 1 + s2 Usemos el teorema anterior en el caso de que f = t se tendr´a entonces que dt = θ′ (s)(− sin(θ(s)), cos(θ(s)). (4.2) ds Como (− sin(θ(s)), cos(θ(s)) = n(s), resulta que k(s) = θ′ (s).

4.2. ECUACIONES DE FRENET PLANAS

71

Ejercicio 164. Encu´entrese la curvatura de la curva x(t) = et (cos t, sin t) calculando θ′ (s). Por tanto si tenemos k(s) podemos recuperar θ(s) integrando: Z s θ(s) = k(ξ)dξ + θ0 s0

y la curva x = x(s) mediante una segunda integraci´on: Z s  Z s x(s) = cos θ(ξ)dξ + x0 , sin θ(ξ)dξ + y0 . s0

s0

Ejercicio 165. Esto demuestra que a partir de la curvatura recuperamos la curva salvo una traslaci´on y un giro. ¿Por qu´e? Ejercicio 166. Usando la expresi´ on que acabamos de calcular, demu´estrese de nuevo que las rectas son las curvas de curvatura nula y las circunferencias las de curvatura constante. Ejercicio 167. Calc´ ulese una curva cuya curvatura est´ a dada por k(s) = √ −1 , s ∈ (−1, 1). 1−s2

4.2.

Ecuaciones de Frenet planas

Sabemos ya que k(s) = x′′ (s) =

dt = k(s)n(s) ds

Esta es conocida como la primera de las ecuaciones de Frenet. Calculemos . En primer lugar, derivando hn, ni = 1 obtenemos que ahora dn ds   dn 0=2 ,n , ds con lo que

dn ds

no tiene componente n y es pues de la forma dn = µ(s)t(s). ds

Para calcular el coeficiente, observamos que µ(s) = dn , t y derivando ds ht, ni = 0 vemos que         dn dn dn dt , n = t, + + hk(s), ni = t, + k(s). 0 = t, ds ds ds ds

CAP´ITULO 4. ECUACIONES DE FRENET

72

Por tanto dn = −k(s)t(s), que es la segunda de las f´ormulas de Frenet. ds Tenemos pues  dt   = k(s)n(s) ds (4.3)   dn = −k(s)t(s) ds De forma matricial, estas ecuaciones se pueden escribir como 

dt ds

dn ds



 = t n



 0 −k(s) . k(s) 0

La idea tras las f´ormulas de Frenet es la siguiente: hemos calculado la velocidad con la que var´ıa una referencia asociada geom´etricamente a la curva ({t, n}) respecto a un par´ametro tambi´en geom´etrico (la longitud de arco) y la hemos expresado en t´erminos de la propia referencia. Obtenemos as´ı una matriz de coordenadas que nos proporciona funciones (en este caso, esencialmente una u ´nica funci´on, la curvatura) con informaci´on geom´etrica sobre la curva. Estas funciones se llaman invariantes diferenciales esta t´ecnica, que es de utilidad en muchos otros contextos, se llama m´etodo de la referencia m´ovil de Cartan.

dt/ds

t(s)

n(s)

dn/ds

Diedro de Frenet

4.3. EL TEOREMA FUNDAMENTAL DE LA TEOR´IA DE CURVAS PLANAS73

4.3.

El teorema fundamental de la teor´ıa de curvas planas

Las ecuaciones de Frenet (4.3) pueden verse como un sistema de cuatro ecuaciones diferenciales ordinarias (ya que cada ecuaci´on tiene dos componentes), las inc´ognitas son t, n y la curvatura k(s) es un dato. Dado que es un sistema lineal de ecuaciones diferenciales ordinarias, tenemos garantizada la existencia de una soluci´on para cada condici´on inicial que demos, t(s0 ) = t0 , n(s0 ) = n0 . La soluci´on, precisamente por la linealidad, se extiende a todo el intervalo en el que est´e definida k. Esto nos deja a un paso de recuperar la curva a partir de la curvatura, ya que ser´a

x(s) =

Z

s

t(ξ)dξ + x0 , s0

como sabemos. Como ejemplo, calculemos la curva cuya curvatura est´a dada por k(s) = (s + 4)s(s − 4)(s − 6) en el intervalo [−5, 7]:

var(’s x0 x1 t0 t1 n0 n1’) k = (s+4)*s*(s-4)*(s-6) ecuaciones = [t0, t1, k*n0, k*n1, -k*t0, -k*t1] variables = [x0, x1, t0, t1, n0, n1] condiciones_iniciales = [0, 0, 0, 1, 0, 0, 1] intervalo = [-5, 7] P = desolve_system_rk4(ecuaciones, variables, ics = condiciones_iniciales, ivar = s, end_points = intervalo, step = 0.001 ) curva = [[x0, x1] for [s, x0, x1, t0, t1, n0, n1] in P] dib_curva = line(curva) dib_curva.set_aspect_ratio(1)

CAP´ITULO 4. ECUACIONES DE FRENET

74

0.4

0.3

0.2

0.1

-0.1

0.1

0.2

0.3

-0.1

Ejercicio 168. ¿C´omo est´ a recorrida la curva? ¿Puedes se˜ nalar cu´ales son los puntos en los que se anula la curvatura? ¿A qu´e se deben los “remolinos”? ¿Puedes explicar el sentido de giro alrededor de los mismos? Este argumento no es una demostraci´on: resolviendo el sistema anterior obtenemos una curva x(s) pero, ¿su curvatura ser´a k(s)? De hecho, no est´a claro que la curva que se obtiene por este procedimiento est´e siquiera parametrizada al arco, porque podr´ıa suceder que t, n no constituyesen una base ortonormal para cada s. Veamos que esto es as´ı. En primer lugar, demostremos que los vectores tangente y normal que se obtienen forman una base

4.3. EL TEOREMA FUNDAMENTAL DE LA TEOR´IA DE CURVAS PLANAS75 ortonormal. Para ello calculamos ht(s), t(s)i′ = 2 ht′ (s), t(s)i = 2k(s) hn(s), t(s)i

hn(s), n(s)i′ = 2 hn′ (s), n(s)i = −2k(s) ht(s), n(s)i

ht(s), n(s)i′ = ht′ (s), n(s)i + ht(s), n′ (s)i = k(s)(hn(s), n(s)i − ht(s), t(s)i)

La manipulaci´on directa de estas ecuaciones no permite concluir que {t, n} es una base ortonormal. Sin embargo, podemos observar que las tres ecuaciones anteriores constituyen un sistema de tres ecuaciones diferenciales ordinarias con tres inc´ognitas y que, por tanto, tiene soluci´on u ´nica para cada condici´on inicial que demos. Dado que ht, ti = hn, ni = 1, ht, ni = 0 es soluci´on de las mismas y compatible con la condici´on inicial de ortonormalidad, concluimos la ortonormalidad de la base. Ofrezcamos ahora una demostraci´on matricial de este hecho, que har´a transparente que lo que realmente se est´a usando es la antisimetr´ıa de las ecuaciones de Frenet. N´otese que la matriz  F= t n satisface que



F F=



t⊤ n⊤



 t n =



 ht, ti ht, ni . hn, ti hn, ni

La matriz G = F⊤ F es una matriz sim´etrica. Tenemos adem´as, debido a las f´ormulas de Frenet, que     0 −k(s) ′ ′ ′ . F = t n = t n k(s) 0   0 −k(s) , de forma que las f´ormulas de Frenet se escriLlamemos A = k(s) 0 ben simplemente como F′ = FA. N´otese que A es una matriz antisim´etrica, es decir, A⊤ = −A. Veamos que G es constante. Para ello derivamos G′ = (F⊤ F)′ = F′⊤ F + F⊤ F′ = (FA)⊤ F + F⊤ (FA) = −AG + GA. Ahora observemos que G′ = −AG + GA puede interpretarse como un sistema lineal de ecuaciones diferenciales (en las coordenadas de G con tantas ecuaciones como inc´ognitas: tres. Por tanto para cada condici´on inicial hay una soluci´on u ´nica. Nosotros tenemos como

76

CAP´ITULO 4. ECUACIONES DE FRENET

condici´on inicial la dada por una elecci´on de base ortonormal para alg´ un s = s0 , de forma que G(s0 ) = Id. Como la matriz constante Id es una soluci´on del sistema de ecuaciones compatible con la condici´on inicial, concluimos que necesariamente G(s) = Id para todo s, es decir, {t(s), n(s)} son una base ortonormal. Siendo t unitario, la curva que se obtiene integr´andolo est´a parametrizada al arco y el hecho de que se satisfacen las f´ormulas de Frenet F′ = FA fuerza a k(s) a ser su curvatura. Adem´as, la unicidad de la soluci´on de las ecuaciones de Frenet asegura que tenemos tantas curvas como posibles condiciones iniciales, pero todas ellas difieren s´olo movimientos que llevan unas condiciones iniciales en otras. Ejercicio 169. Det´ allese la expresi´ on matricial de dichos movimientos.

4.4.

Curvas espaciales. Plano osculador

En el caso de que la curva sea espacial, no hay un vector normal a t u ´nicamente determinado. Sin embargo, en tanto que x′′ (s) 6= 0 podemos definir un vector privilegiado n(s) =

x′′ (s) kx′′ (s)k

que llamaremos vector normal a la curva. La curvatura es ahora k(s) = kx′′ (s)k y satisface, por definici´on, la primera ecuaci´on de Frenet: dt = k(s)n(s). ds

(4.4)

En adelante, supondremos que la curva cumple la condici´on adicional de regularidad x′′ (s) 6= 0 para todo s. Sin embargo, el s´olo conocimiento de la curvatura no determina la curva espacial salvo movimientos, como se demuestra en el siguiente ejercicio: Ejercicio 170. Calc´ ulese la curvatura de la h´elice x(t) = (a cos t, a sin t, bt) Concl´ uyase que la curvatura no es suficiente para determinar una curva espacial salvo movimientos.

4.4. CURVAS ESPACIALES. PLANO OSCULADOR

77

La curvatura apareci´o como la velocidad del movimiento de la recta tangente a una curva parametrizada por el arco. La recta tangente en un punto es la recta que mejor se ajusta a la curva en dicho punto. De forma an´aloga, se puede considerar el plano que mejor se ajusta a una curva espacial en un punto. Para calcularlo, procederemos de forma semejante a c´omo calculamos la circunferencia osculadora: tomemos tres puntos de la curva, calculemos el plano que pasa por los mismos (suponiendo que no est´an alineados) y veamos cu´al es el l´ımite cuando dos de ellos se aproximan al tercero. Si los puntos son x(si ), i = 0, 1, 2 y llamamos v a un vector director del plano la ecuaci´on del plano es hx − x(s0 ), vi = 0. Como la funci´on f (s) = hx(s) − x(s0 ), vi se anula en s0 , s1 y s2 , habr´a dos puntos intermedios en los que se anule la primera derivada y uno en el que se anula la segunda derivada. Tomando l´ımites cuando s1 , s2 → s0 tenemos que f (s0 ) = hx(s0 ) − x(s0 ), vi = 0, lo que es trivial, y f ′ (s0 ) = hx′ (s0 ), vi = 0 f ′′ (s0 ) = hx′′ (s0 ), vi = 0

(4.5)

De la segunda y terceras ecuaciones concluimos que el del plano l´ımite, que llamaremos plano osculador, est´a generado por x′ (s0 ) y x′′ (s0 ). En un ejercicio posterior se demostrar´a que la condici´on x′′ (s) 6= 0 para todo s de hecho garantiza la existencia de dicho plano. Ejercicio 171. Demu´estrese que     dx d2 x dx d2 x L =L , , dt dt2 ds ds2 de forma que, en el caso de que la curva no est´e parametrizada por el arco, podemos escribir el plano osculador como   dx d2 x π = x(t) + L . , dt dt2 Ejercicio 172. Demu´estrese que la ecuaci´ on del plano osculador es   dx d2 x det x − x(t), , 2 . dt dt donde x = (x, y, z)⊤ .

78

CAP´ITULO 4. ECUACIONES DE FRENET

N´otese que no es lo mismo el plano osculador, que pasa por un punto de la curva, que el plano vectorial asociado, que es paralelo al mismo y contiene al origen de coordenadas, es decir, es el plano vectorial engendrado por los vectores t(s) y n(s). No obstante, para no sobrecargar demasiado la terminolog´ıa, llamaremos a ambos planos osculadores y se sobreentender´a por el contexto a cu´al nos estamos refiriendo.

4.5.

Binormal. F´ ormulas de Frenet y consecuencias.

Para analizar la variaci´on del plano osculador vectorial, observemos que existe un vector director del mismo can´onico: el producto vectorial del tangente por el normal. Dicho vector se denomina vector binormal: b(s) = t(s) × n(s). Estudiar la variaci´on del vector binormal es equivalente a estudiar la variaci´on del plano osculador. Es m´as pr´actico, de hecho, analizar simult´aneamente la variaci´on del triedro de Frenet F = (t, n, b). Sabemos, dado que el triedro es una base, que F′ = FA para cierta matriz A. Veamos que esta matriz ha de ser antisim´etrica: como F⊤ F = Id, derivando esta identidad tenemos que 0 = (F⊤ F)′ = F′⊤ F + F⊤ F′ = A⊤ F⊤ F + F⊤ FA = A⊤ + A, con lo que A es efectivamente antisim´etrica, como quer´ıamos demostrar. La ecuaci´on (4.4) nos dice c´omo es la primera columna de A. Teniendo en cuenta la antisimetr´ıa de la misma, concluimos que su forma es   0 −k(s) 0 0  , A = k(s) 0 − 0

donde  es simplemente una coordenada de la matriz que desconocemos. Pero, por la definici´on de la matriz A, es

db = n. ds Por tanto,  es justamente la variaci´on del vector binormal que est´abamos buscando. A esta funci´on le llamaremos torsi´on de la curva y se denotar´a por τ (s) =  = h

db , n(s)i. ds

´ 4.5. BINORMAL. FORMULAS DE FRENET Y CONSECUENCIAS.

79

Seg´ un esta definici´on, la matriz de variaci´on A es   0 −k(s) 0 0 τ (s) A = k(s) 0 −τ (s) 0

y las f´ormulas de Frenet correspondientes, escritas como un sistema de ecuaciones, son  dt   = k(s)n(s)   ds   dn (4.6) = −k(s)t(s) − τ (s)b(s)  ds      db = τ (s)n(s) ds De forma an´aloga a como se trat´o el caso plano, podemos demostrar ahora que dadas dos funciones diferenciables k(s) > 0 y τ (s) arbitrarias, existe, salvo movimientos del espacio, una u ´nica curva diferenciable que las produce como curvatura y torsi´on. Ejercicio 173. Aparentemente, en la resoluci´on del sistema de ecuaciones diferenciales ordinarias no se emplea en absoluto el que la curvatura sea positiva. ¿Qu´e ocurre si se proporciona una funci´on de curvatura que cambia de signo? Ejercicio 174. Demu´estrese cuidadosamente que si se transforma la curva mediante una isometr´ıa del espacio que preserva la orientaci´ on, entonces la nueva curva tiene la misma curvatura y torsi´on. ¿Qu´e ocurre si la isometr´ıa invierte la orientaci´ on? Ejercicio 175. ¿Qu´e ocurre con la curvatura y la torsi´on de una curva cuando se invierte su sentido de recorrido? Ejercicio 176. De forma an´aloga a como se ha hecho en el caso plano, demu´estrese que la torsi´on de una curva no parametrizada por el arco viene dada por la f´ ormula τ (t) = −

det (x′ (t), x′′ (t), x′′′ (t)) kx′ (t) × x′′ (t)k2

Ejercicio 177. Demu´estrese que la torsi´on de una curva es nula si y s´ olo si est´ a contenida en un plano. Ejercicio 178. Demu´estrese que no existe ning´ un movimiento que transforma la h´elice lev´ ogira x(t) = (a cos t, a sin t, bt) en la h´elice dextr´ogira y(t) = (a sin t, a cos t, bt).

80

CAP´ITULO 4. ECUACIONES DE FRENET

Ejercicio 179. Encu´entrese la expresi´ on anal´ıtica de las curvas con curvatura contante k(s) = k0 > 0 y torsi´on tambi´en constante τ (s) = τ0 . Ejercicio 180. De forma an´aloga a como se hizo en el caso de curvas planas, progr´amese el c´odigo necesario para calcular una curva espacial a partir de su torsi´on y su curvatura. Ejercicio 181. Progr´amese el c´odigo necesario para calcular la curvatura y torsi´on de una curva no parametrizada por el arco. Ejercicio 182. Progr´amese el c´odigo necesario para calcular el triedro de Frenet de una curva no parametrizada por el arco. Compru´ebese que las f´ ormulas de Frenet se cumplen para la curva x(t) = (et , 1 + t2 , cos t). Ejercicio 183. Calc´ ulese la recta tangente y el plano osculador a la curva 2 2 2 dada por x + y + z − 1 = x2 + y 2 − x = 0 en el punto (0, 0, 1).

4.6.

Forma can´ onica local de una curva

Recordemos que el desarrollo de Taylor de una curva x = x(s) comenzar´ıa como 1 1 ′′ x (s0 )(s − s0 )2 + x′′′ (s0 )(s − s0 )3 + R(s), 2! 3! (4.7) 3 donde O(s − s0 ) es una funci´on diferenciable que satisface que

x(s) − x(s0 ) = x′ (s0 )(s − s0 ) +

l´ım

s→s0

R(s) = 0. (s − s0 )3

Las derivadas x(i) (s0 ), dado que son vectores, se pueden escribir en t´erminos del triedro de Frenet {t(s0 ), n(s0 ), b(s0 )} y la clave para hacer esto es observar que x′ (s0 ) = t(s0 ) y que las derivadas de orden superior se pueden calcular usando las f´ormulas de Frenet. Concretamente, x′′ (s0 ) = t′ (s0 ) = k(s0 )n(s0 ) x′′′ (s0 ) = (kn)′ (s0 ) = k ′ (s0 )n(s0 ) + k(s0 )n′ (s0 ) = k ′ (s0 )n(s0 ) + k(s0 )(−k(s0 )t(s0 ) − τ (s0 )b(s0 )) = −k(s0 )2 t(s0 ) + k ′ (s0 )n(s0 ) − k(s0 )τ (s0 )b(s0 ).

(4.8)

Si usamos como sistema de coordenadas aquel con origen x(s0 ) y base el triedro de Frenet en s0 , sustituyendo las ecuaciones (4.8) en (4.7) obtenemos

´ 4.6. FORMA CANONICA LOCAL DE UNA CURVA

81

que las componentes de la curva son 1 x(s) = s − k(0)2 s3 + Rx (s) 6 1 1 y(s) = k(0)s2 + k ′ (0)s3 + Ry (s) 2 6 1 z(s) = − k(0)τ (0)s3 + Rz (s) 6 donde, por simplificar el resultado, hemos supuesto s0 = 0. Estas ecuaciones de la curva reciben el nombre de forma can´onica local. Se podr´ıa, si es preciso, a˜ nadir m´as t´erminos al desarrollo que depender´an tambi´en de la curvatura, la torsi´on y sus derivadas. N´otese que, en el caso de que la curva sea anal´ıtica, esta t´ecnica demuestra alternativamente que la curvatura y la torsi´on determina la curva salvo movimientos del espacio. El plano osculador, generado por el tangente y el normal, es en este sistema de referencia el plano xy. Al plano generado por el vector normal y el binormal se le llama plano normal a la curva y es el plano yz. Finalmente, al plano engendrado por el vector tangente y binormal a la curva se le llama plano rectificante, y resulta ser el plano xz. Ejercicio 184. Demu´estrese que si las rectas normales a una curva pasan todas por un punto fijo, entonces la curva est´ a contenida en un c´ırculo. Ejercicio 185. Supongamos que k(s) 6= 0 6= τ (s) para todo s. Demu´estrese que la curva yace en una esfera si y s´ olo si R2 + (R′ )2 T 2 = cte. siendo R = 1/k y T = 1/τ . (Indicaci´ on: sup´ongase la curva parametrizada por el arco. Si el centro de la esfera es c, der´ıvese sucesivamente la ecuaci´ on hx − c, x − ci = r2 para obtener las coordenadas de x−c respecto al triedro de Frenet. Una vez obtenidas, der´ıvese la ecuaci´ on que expresa x−c en el triedro de Frenet para concluir la relaci´on pedida. Para el rec´ıproco, compru´ebese que la derivada de x + Rn − R′ T b es nula y concl´ uyase que existe c de forma que hx − c, x − ci es constante). Ejercicio 186. Demu´estrese que si la curvatura de una curva es no nula en un valor del par´ ametro s0 entonces cualesquiera tres valores de la curva en par´ ametros distintos y suficientemente cercanos a s0 determinan un plano. Ejercicio 187. Demu´estrese que el vector binormal determina la curvatura de la curva, as´ı como el valor absoluto de la torsi´on.

82

CAP´ITULO 4. ECUACIONES DE FRENET

Ejercicio 188. Demu´estrese que el vector normal a una curva con torsi´ on no nula en todos los puntos determina las funciones curvatura y torsi´ on. Ejercicio 189. Sea x = x(s) una curva plana. Demu´estrese que si s0 es tal que la distancia de x(s0 ) al origen es m´axima, entonces la curvatura satisface que |k(s0 )| ≥ 1/kx(s0 )k.

Cap´ıtulo 5 Superficies De la misma forma que una curva estaba descrita por un punto que se mov´ıa a lo largo del tiempo, se puede obtener una superficie mediante una curva que se mueve a lo largo del tiempo. Modelizaremos esta idea definiendo, provisionalmente, una superficie como una aplicaci´on diferenciable x : U ⊂ R2 → R3 ,

siendo U un abierto del plano R2 . Si denotamos por (u, v) ∈ U a un punto gen´erico de U , fijado u la aplicaci´on v 7→ x(u, v) es una curva parametrizada por u, de forma que la superficie est´a siendo barrida por una curva m´ovil. Ejemplo 6. Consideremos dos curvas x1 , x2 : I ⊂ R → R3 . Si interpolamos linealmente cada par de puntos correspondientes x1 (u), x2 (u), obtenemos una aplicaci´ on diferenciable x : I × R → R3 definida por x(u, v) = (1 − v)x1 (u) + vx2 (u). Por ejemplo, sea x1 (u) = (u, 0, 0) y x2 (u) = (u, 1, u). Entonces x(u, v) = (1 − v)(u, 0, 0) + v(u, 1, u) = (u, v, uv), y obtenemos as´ı el paraboloide hiperb´olico z = xy. Ejercicio 190. Demu´estrese que la superficie anterior coincide efectivamente con el hiperboloide z = xy.

5.1.

Superficies de B´ ezier

Supongamos que x = x(u) es una curva de B´ezier de grado n, de forma que n X x(u) = bi Bin (u). i=0

83

84

CAP´ITULO 5. SUPERFICIES

Figura 5.1: Una recta m´ovil barre una superficie, obteni´endose un paraboloide hiperb´olico.

´ 5.1. SUPERFICIES DE BEZIER

85

Una forma razonable de desplazar la curva es hacer que los puntos del pol´ıgono de control bi vayan recorriendo a su vez una curva de B´ezier, de forma que m X bi = bi (v) = bij Bjm (v). j=0

De esta forma la superficie que estas curvas barren ser´a x(u, v) = = =

n X

i=0 n X

bi (v)Bin (u) m X

j=0 i=0 n m XX

!

bij Bjm (v) Bin (u)

(5.1)

bij Bin (u)Bjm (v)

i=0 j=0

A los puntos {bij } se les llama malla de control de la superficie, que ser´a una superficie polin´omica de bigrado (n, m).

Figura 5.2: Dibujo de una superficie de B´ezier bic´ ubica junto con su malla de control. Ejercicio 191. En general, ¿por cu´antos puntos de la malla de control pasa la superficie de B´ezier?

86

CAP´ITULO 5. SUPERFICIES

Ejercicio 192. Escr´ıbase un programa en Sage que, dada una malla de control, calcule y dibuje la superficie correspondiente.

5.2.

Superficies parametrizadas regulares

Una aplicaci´on diferenciable x : U ⊂ R2 → R3 se dice regular si en cada punto (u, v) ∈ U la diferencial dx(u,v) : R2 → R3 tiene rango dos, es decir, los vectores dx(e1 ) = ∂x/∂u = xu y dx(e2 ) = ∂x/∂v = xv son linealmente independientes en cada (u, v). Normalmente simplificaremos la notaci´on no especificando el punto (u, v) a menos que sea necesario. Cuando la aplicaci´on x es regular en todo punto, diremos que es una superficie parametrizada. La condici´on de regularidad nos permite definir el espacio vectorial tangente a la superficie en un punto p = x(u, v): es simplemente la imagen de la diferencial dx(u,v) , es decir, el espacio vectorial generado por xu , xv . El espacio af´ın tangente ser´a simplemente el trasladado de este espacio vectorial al punto p, es decir, x(u, v) + L(xu (u, v), xv (u, v)). Normalmente llamaremos a ambos espacios simplemente espacio tangente, estando claro por el contexto a cu´al de ellos nos estamos refiriendo. Ejercicio 193. Calc´ ulense las ecuaciones param´etricas e impl´ıcitas del plano tangente a la superficie x(u, v) = (u, v, uv) en el punto (1, −1, −1). Ejercicio 194. Encu´entrese los puntos de la superficie x(u, v) = (u − v, u + v, u2 + v 2 ) en los que el plano tangente es paralelo al plano x − y + z = 0. Ejemplo 7. Un plano de R3 es una superficie. En efecto, si x0 es un punto del plano y u, v son dos vectores del mismo linealmente independientes, la aplicaci´ on x(u, v) = x0 + uu + vv, con u, v ∈ R

es una superficie dado que es una aplicaci´ on diferenciable y xu = u, xv = v.

5.2.1.

El plano tangente es el que mejor aproxima la superficie en el punto

El plano tangente est´a caracterizado por ser el que mejor aproxima la superficie en el punto de tangencia. Esto se puede formalizar del siguiente

5.2. SUPERFICIES PARAMETRIZADAS REGULARES

87

Figura 5.3: Plano tangente a una superficie. Se muestran tambi´en los vectores xu , xv que lo generan.

CAP´ITULO 5. SUPERFICIES

88

modo: la distancia (con signo) de un punto x(u, v) al plano que pasa por x(u0 , v0 ) y es normal al vector unitario u est´a dada por d(u, v) = hx(u, v) − x(u0 , v0 ), ui.

(5.2)

Calculemos el desarrollo de Taylor de d(u, v) en (u0 , v0 ). Para ello recordemos que el desarrollo de Taylor de x(u, v) comienza por x(u, v) − x(u0 , v0 ) = xu (u0 , v0 )(u − u0 ) + xv (u0 , v0 )(v − v0 ) + R(u, v). Sustituyendo en 5.2 obtenemos d(u, v) = hxu (u0 , v0 ), ui(u − u0 ) + hxv (u0 , v0 ), ui(v − v0 ) + R(u, v), siendo R(u, v) = hR(u, v), ui el resto, el cual satisface R(u, v) = 0. (u,v)→(u0 ,v0 ) k(u − u0 , v − v0 )k l´ım

Vemos pues que el que u ⊥ xu , xv , es decir, que el plano sea el que hemos definido como plano tangente, asegura que el plano sea el m´as pr´oximo posible a la superficie en el punto x(u0 , v0 ).

5.3.

Ejemplos de superficies

5.3.1.

Grafos

Siempre que tengamos una funci´on diferenciable f : U → R, su grafo produce una superficie parametrizada dada por x(u, v) = (u, v, f (u, v)). Ejercicio 195. Demu´estrese que un grafo es siempre una superficie parametrizada.

5.3.2.

Superficies de rotaci´ on

Dada una curva en el plano XZ con ecuaciones x = x(t), y = 0, z = z(t), podemos obtener una superficie gir´andola alrededor del eje OZ. Como la matriz de un giro de ´angulo θ alrededor del eje OZ es   cos θ − sin θ 0 Rθ =  sin θ cos θ 0 , 0 0 1

5.4. ROTACIONES Y CUATERNIONES

89

la parametrizaci´on debe ser   x(t) x(t, θ) = Rθ  0  . z(t)

Ejercicio 196. Compru´ebese que si la curva (x(t), z(t)) es regular y x(t) > 0 entonces la parametrizaci´on de la superficie de rotaci´ on es regular. Este ejemplo suscita la cuesti´on de c´omo girar un objeto alrededor de una recta arbitraria. Una alternativa ser´ıa hacer un cambio de coordenadas, de forma que la recta se convierta en el eje OZ, pero esto no es lo m´as conveniente. Uno de los m´etodos m´as eficaces consiste en usar cuaterniones.

5.4.

Rotaciones y cuaterniones

Los n´ umeros cuaterni´onicos son una generalizaci´on de los n´ umeros complejos. Por definici´on, los n´ umeros cuaterniones H son el espacio vectorial usual de cuaternas de n´ umeros reales R4 , en el cual introducimos una operaci´on adicional de producto. Para definirlo, usaremos los s´ımbolos formales i, j y k para escribir un cuaterni´on q = (a, b, c, d) ∈ H de la forma q = a + bi + cj + dk. La multiplicaci´on de cuaterniones se lleva a cabo teniendo en cuenta las reglas i2 = j2 = k2 = ijk = −1, de las que puede deducirse el resultado del producto de cualesquiera otras dos unidades imaginarias teniendo en cuenta que el producto es asociativo. Ejercicio 197. Calc´ ulense los productos de cualesquiera dos unidades imaginarias en cualquier orden. Los n´ umeros cuaterni´onicos, dotados de estas dos operaciones, resultan ser un cuerpo no conmutativo. Para calcular el inverso de un cuaterni´on, definamos primero la noci´on de conjugado. Dado el cuaterni´on q = a + bi + cj + dk se define ¯ = a − bi − cj − dk. q Ejercicio 198. La operaci´ on de conjugaci´ on se comporta con respecto al producto de cuaterniones del siguiente modo: ¯. ¯′q qq′ = q

CAP´ITULO 5. SUPERFICIES

90 Ejercicio 199. Compru´ebese que

¯ q = kqk2 = a2 + b2 + c2 + d2 . q¯ q=q El inverso de q 6= 0 est´a dado por q−1 =

¯ q . kqk2

Ejercicio 200. Compru´ebese que, efectivamente, qq−1 = q−1 q = 1. Consideraremos la inmersi´on del espacio R3 dentro de los cuaterniones que identifica un vector v = (u, v, w) con el cuaterni´on imaginario, que denotaremos tambi´en por v, v = ui + vj + wk. Dado un cuaterni´on q 6= 0, introducimos la operaci´on de conjugaci´on que asigna al vector v el cuaterni´on v′ = q−1 vq. Resulta que v′ es tambi´en un cuaterni´on imaginario. En efecto, los cuater¯ = −v. Se tiene que niones imaginarios se caracterizan porque v ¯ ′ = q−1 vq = q ¯v ¯ v

¯ q = −q−1 vq = −v′ , kqk2

con lo que v′ es efectivamente imaginario puro. Afirmamos que, adem´as, kvk = kv′ k. Esto es consecuencia de la identidad, m´as general, kqq′ k = kqkkq′ k. Esta u ´ltima igualdad se demuestra usando la notable identidad (a2 + b2 + c2 + d2 )(x2 + y 2 + z 2 + t2 ) =(−d2 + ax − by − cz)2 + (cd + bx + ay − dz)2 + (−bd + cx + dy + az)2 + (ad + dx − cy + bz)2 , tambi´en llamada identidad de los cuatro cuadrados de Euler. Ejercicio 201. Demu´estrese la identidad de los cuatro cuadrados de Euler y kqq′ k = kqkkq′ k.

91

5.4. ROTACIONES Y CUATERNIONES

Ejercicio 202. Demu´estrese que la aplicaci´ on v 7→ v′ = qvq−1 define una isometr´ıa lineal Rq de R3 . Ejercicio 203. Demu´estrese que la matriz de Rq est´ a dada por   2 2 2 2  Rq = 

a +b −c −d a2 +b2 +c2 +d2 2bc+2ad a2 +b2 +c2 +d2 −2ac+2bd a2 +b2 +c2 +d2

2bc−2ad a2 +b2 +c2 +d2 2 a −b2 +c2 −d2 a2 +b2 +c2 +d2 2ab+2cd a2 +b2 +c2 +d2

2ac+2bd a2 +b2 +c2 +d2 −2ab+2cd a2 +b2 +c2 +d2 2 a −b2 −c2 +d2 a2 +b2 +c2 +d2

 .

Ejercicio 204. Demu´estrese que el determinante de Rq es positivo construyendo una curva de cuaterniones qt = a + t(bi + cj + dk) y observando que det Rq = det Rq0 = 1. Ejercicio 205. Los n´ umeros cuaterni´onicos se pueden descomponer como pares q = (s, v), donde s es la parte real de q y v su parte imaginaria. Con esta notaci´on, demu´estrese que el producto de cuaterniones est´ a dado por la f´ormula (s, v) · (s′ , v′ ) = (ss′ − hv, v′ i, sv′ + s′ v + v × v′ ). Ejercicio 206. Demu´estrese que Rq , siendo q = (cos(θ/2), sin(θ/2)n), representa un giro de ´ angulo θ alrededor del vector unitario n. Ejercicio 207. Escr´ıbase un programa en Sage que dada una curva espacial y un vector, calcule y dibuje la superficie que se obtiene al girarla. Ejercicio 208. Dada una recta en el espacio, determinada por un punto y un vector director, modif´ıquese el programa anterior para que calcule la superficie que se obtiene al girar la curva alrededor de la recta. Ejercicio 209. Usar el programa del ejercicio anterior para girar una recta alrededor de otra. Dib´ ujense los distintos casos, seg´ un la posici´on relativa de ambas. Ejercicio 210. Idem. al ejercicio anterior, pero para una circunferencia y una recta. Ejercicio 211. Idem. al ejercicio anterior, pero para una elipse y una recta. Ejercicio 212. Idem. al ejercicio anterior, pero para una hip´erbola y una recta. Ejercicio 213. Encu´entrese las ecuaciones param´etricas e impl´ıcitas de la superficie que se obtiene al girar la a) circunferencia x2 + y 2 = 1 b) elipse x2 /a2 + y 2 /b2 = 1 c) hip´erbola x2 /a2 − y 2 /b2 = 1 d) par de rectas x2 = 1 e) par de rectas x2 − y 2 = 0 alrededor del eje OX. Dib´ ujense. ¿Qu´e nombres reciben estas superficies?

CAP´ITULO 5. SUPERFICIES

92

Ejercicio 214. Compru´ebese que las superficies parametrizadas del problema anterior son regulares (si la parametrizaci´on es la obvia sugerida por su construcci´on). Ejercicio 215. Encu´entrese las ecuaciones param´etricas e impl´ıcitas de la superficie que se obtiene al girar la circunferencia de radio 1 centrada en el punto (0, 2, 0) y contenida en el plano XY alrededor del eje OX. Esta superficie es un toro de revoluci´ on. Ejercicio 216. Compru´ebese que la superficies parametrizada del problema anterior son regulares (si la parametrizaci´on es la obvia sugerida por su construcci´on). Ejercicio 217. Se considera el segmento L que une los puntos (0, 2, 1) y (0, 2, −1). El segmento comienza a girar alrededor del eje OZ a velocidad angular 1 al mismo tiempo que gira alrededor de su centro con velocidad angular 1/2 en el plano determinado por el eje OZ y el centro del segmento, barriendo de esta forma una superficie. Calc´ ulense unas ecuaciones param´etricas de dicha superficie y demu´estrese que es una parametrizaci´on regular. Esta superficie es un ejemplo de banda de M¨obius.

5.5.

Curvas sobre superficies. Geod´ esicas

Una curva sobre una superficie parametrizada x : U → R3 es, por definici´on, la imagen de una curva diferenciable u = u(t), v = v(t) con t en un intervalo abierto I ⊂ R, es decir, es la aplicaci´on I ∋ t 7→ x(u(t), v(t)) ∈ R3 . Escribiremos a veces simplemente que x = x(t) es una curva de la superficie x = x(u, v), entendi´endose que x(t) = x(u(t), v(t)). El vector tangente a la curva x(t) vendr´a dado, seg´ un la regla de la cadena, por x′ (t) = u′ (t)xu (u(t), v(t)) + v ′ (t)xv (u(t), v(t)), aunque simplemente escribiremos x ′ = u′ x u + v ′ x v . Equivalentemente,  ′ u x (t) = dx ′ . v ′

(5.3)

´ 5.5. CURVAS SOBRE SUPERFICIES. GEODESICAS

93

Figura 5.4: Curva u(t) = cos(t), v(t) = sin(t) sobre la superficie x(u, v) = (u, v, cos(uv)).

5.5.1.

Geod´ esicas

Una curva x = x(t) se dice geod´esica cuando su aceleraci´on x′′ (t) es siempre normal a la superficie. Desde un punto de vista intuitivo, esto significa que no sufre aceleraciones detectables desde la superficie. La condici´on de ser geod´esica se puede escribir como un sistema de ecuaciones diferenciales ordinarias de segundo orden. Veamos c´omo: si derivamos la ecuaci´on (5.3) obtenemos que x′′ = u′′ xu + v ′′ xv + (u′ )2 xuu + 2u′ v ′ xuv + (v ′ )2 xvv y por tanto las condiciones hx′′ , xu i = hx′′ , xv i = 0 se pueden escribir como hx′′ , xu i = u′′ hxu , xu i + v ′′ hxv , xu i + (u′ )2 hxuu , xu i + 2u′ v ′ hxuv , xu i + (v ′ )2 hxvv , xu i = 0 hx′′ , xv i = u′′ hxu , xv i + v ′′ hxv , xv i + (u′ )2 hxuu , xv i + 2u′ v ′ hxuv , xv i + (v ′ )2 hxvv , xv i = 0 (5.4)

CAP´ITULO 5. SUPERFICIES

94

Es quiz´a conveniente introducir la terminolog´ıa cl´asica siguiente: las funciones E = hxu , xu i F = hxu , xv i G = hxv , xv i se llaman coeficientes de la primera forma fundamental de la superficie. Son funciones E, F, G : U → R. N´otese que sin embargo en (5.4) las funciones E, F y G est´an evaluadas en (u(t), v(t)). Con esta notaci´on, podemos escribir las ecuaciones de las geod´esicas de forma matricial:   ′′     A u E F =0 + B v ′′ F G donde A = (u′ )2 hxuu , xu i + 2u′ v ′ hxuv , xu i + (v ′ )2 hxvv , xu i, B = (u′ )2 hxuu , xv i + 2u′ v ′ hxuv , xv i + (v ′ )2 hxvv , xv i. Como estamos suponiendo que la superficie parametrizada es regular, los vectores xu y xv son linealmente independientes. Esto es equivalente que la matriz   E F I= F G sea invertible.

Ejercicio 218. Demu´estrese que efectivamente la matriz I es invertible viendo que kxu × xv k2 = EG − F 2 y teniendo en cuenta que, por hip´ otesis, xu y xv son linealmente independientes. Por tanto las ecuaciones de las geod´esicas producen el sistema de ecuaciones diferenciales ordinarias de segundo orden    ′′  u −1 A , = −I B v ′′ en donde, al estar dado el lado izquierdo por funciones diferenciables en u, v, u′ , v ′ y t concluimos que tiene soluci´on u ´nica en un entorno de t0 para

´ 5.5. CURVAS SOBRE SUPERFICIES. GEODESICAS

95

cada cualquier condici´on inicial que demos u(t0 ) = u0 , v(t0 ) = v0 , y u′ (t0 ) = u′0 , v ′ (t0 ) = v0′ . tales que (u0 , v0 ) ∈ U . Desde el punto de vista geom´etrico, esto significa que existe una u ´nica geod´esica que a tiempo t0 pase por un punto dado de la superficie con velocidad dada. La vida de dicha geod´esica, sin embargo, no tiene porqu´e ser todo R, pudiendo estar reducida s´olo a un intervalo abierto en el cual se encuentra t0 . La velocidad de una geod´esica es constante. En efecto, si llamamos f (t) = hx′ (t), x′ (t)i entonces f ′ (t) = 2hx′′ (t), x′ (t)i = 0, ya que x′′ es perpendicular a la superficie y x′ es un vector tangente a la misma. Por tanto f debe ser constante, es decir, la velocidad de la geod´esica es constante. Ejercicio 219. El rec´ıproco no es cierto: dar un ejemplo de curva sobre una superficie que tenga velocidad constante pero que no sea geod´esica. Ejercicio 220. Demu´estrese que las geod´esicas del plano son las l´ıneas rectas recorridas a velocidad constante. Ejercicio 221. Demu´estrese que las geod´esicas de una esfera son sus c´ırculos m´aximos. Ejercicio 222. Encu´entrense las ecuaciones de las geod´esicas de una superficie parametrizada como un grafo. Ejercicio 223. Dada una superficie de revoluci´ on, est´ udiese si los paralelos y los meridianos son en alg´ un caso geod´esicas. Ejercicio 224. Supongamos que dos superficies son tangentes a lo largo de una curva. Demu´estrese que son geod´esicas para una de ellas si y s´ olo si lo son para la otra. Ejercicio 225. Demu´estrese que si una recta est´ a contenida en una superficie, entonces es una geod´esica de la misma. Ejercicio 226. Consideremos todas las geod´esicas que parten de un punto de una superficie con velocidad de m´odulo r. El conjunto de los puntos que dichas geod´esicas alcanzan a tiempo 1 se llama c´ırculo geod´esico de radio r. Escr´ıbase un programa en Sage que, dada una superficie y un punto de la misma, dibuje el mencionado c´ırculo geod´esico.

96

CAP´ITULO 5. SUPERFICIES

Cap´ıtulo 6 Las formas fundamentales 6.1.

Primera forma fundamental

Hemos definido previamente los coeficientes de la primera forma fundamental como E = hxu , xu i, F = hxu , xv i, G = hxv , xv i. Ejercicio 227. Dada una curva parametrizada por el arco, x = x(s), se considera la superficie tangencial x(s) + λt(s). Calc´ ulese su primera forma fundamental. Con ellos es posible realizar medidas intr´ınsecas en la superficie: calcular longitudes de curvas, ´angulos y ´areas. Ve´amoslo.

6.1.1.

Longitudes de curvas

La longitud de una curva x = x(t), t ∈ [a, b], sobre una superficie parametrizada x = x(u, v) viene dada por Z

b a

hx′ (t), x′ (t)i1/2 dt.

Dado que x′ = u′ xu + v ′ xv , se tiene que hx′ (t), x′ (t)i =hxu , xu i(u′ )2 + 2hxu , xv iu′ v ′ + hxv , xv i(v ′ )2 =E(u′ )2 + 2F u′ v ′ + G(v ′ )2 . 97

CAP´ITULO 6. LAS FORMAS FUNDAMENTALES

98

Desde un punto de vista intuitivo, merece la pena recordar la interpretaci´on cl´asica del elemento de longitud escrito en t´erminos de la primera forma fundamental. Llamando ds2 = hx′ (t), x′ (t)idt2 y escribiendo u′ = du/dt, v ′ = dv/dt, tenemos que ds2 = Edu2 + 2F dudv + Gdv 2 , que se puede ver como una versi´on infinitesimal de un teorema de Pit´agoras distorsionado. La geometr´ıa de la superficie ser´a Eucl´ıdea cuando el elemento de longitud satisfaga ds2 = du2 + dv 2 . Ejercicio 228. Calc´ ulese la primera forma fundamental de las siguientes superficies parametrizadas: 1. x (u, v) = (a sin u cos v, b sin u sin v, c cos u) , (elipsoide). 2. x (u, v) = (au cos v, bu sin v, u2 ) (paraboloide el´ıptico). 3. x (u, v) = (au cosh v, bu sinh v, u2 ) (paraboloide hiperb´olico). 4. x (u, v) = (a sinh u cos v, b sinh u sin v, c cosh u) (hiperboloide de dos hojas). Dib´ ujense las superficies y sobre las mismas las curvas coordenadas u = 0 y v = 0. Calc´ ulese num´ericamente la longitud de la curva u = t, v = t2 , t ∈ [−1, 1] usando los coeficientes de la primera forma fundamental y dib´ ujese tambi´en, tanto la curva sobre la superficie como en el espacio de par´ ametros.

6.1.2.

´ Angulos

Dos vectores u1 , u2 tangentes a la superficie en un punto (u0 , v0 ), podemos escribirlos en funci´on de la base xu (u0 , v0 ), xv (u0 , v0 ) del espacio tangente en dicho punto como ui = u′i xu + vi′ xv 1 . El coseno del ´angulo α formado por los vectores u1 y u2 vendr´a dado por cos α =

hu1 , u2 i , ku1 kku2 k

pero ya sabemos c´omo escribir kui k en t´erminos de sus coordenadas y la primera forma fundamental y respecto al numerador, se tiene que hu1 , u2 i = u′1 u′2 hxu , xu i + (u′1 v2′ + u′2 v1′ )hxu , xu i + v1′ v2′ hxv , xv i = Eu′1 u′2 + F (u′1 v2′ + u′2 v1′ ) + Gv1′ v2′ . 1

Estamos llamando u′i , vi′ a los coeficientes de los vectores aunque no tienen por qu´e provenir de curvas. No obstante, siempre es posible definir curvas de las que provengan, de hecho, una infinidad de ellas. Por ejemplo, si se desea se puede considerar las dadas por ui (t) = u0 + tu′i y vi (t) = v0 + tvi′ .

6.2. LA SEGUNDA FORMA FUNDAMENTAL

6.1.3.

99

´ Areas

El paralelogramo curvil´ıneo delimitado por las l´ıneas coordenadas y de v´ertices los puntos x(u, v), x(u + du, v), x(u, v + dv), x(u + du, v + dv) se puede aproximar, si du y dv son suficientemente peque˜ nos, por el paralelogramo en el espacio dado por los tres primeros v´ertices. Dicho paralelogramo tendr´a como ´area dA = k(x(u + du, v) − x(u, v)) × (x(u, v + dv) − x(u, v))k. Aproximando x(u + du, v) − x(u, v) ∼ xu du y x(u, v + dv) − x(u, v) ∼ xv dv resultar´a que dA = kxu × xv kdudv.

Esta es la motivaci´on intuitiva de la siguiente definici´on: dada una superficie parametrizada x y una regi´on medible R en el dominio de la parametrizaci´on en la que x es inyectiva, definimos el ´area de la superficie x(R) como Z A= kxu × xv kdudv. R

N´otese que como x no tiene por qu´e ser inyectiva globalmente, puede ocurrir que este ´area no se corresponda con el ´area de x(R). En el caso de que haya una biyecci´on entre R y x(R) (salvo, quiz´a, en un conjunto de medida nula), diremos que A es el ´area de x(R). Puesto que kxu × xv k2 = kxu k2 kxv k2 sin2 θ = kxu k2 kxv k2 (1 − cos2 θ) = kxu k2 kxv k2 − hxu , xv i2 = EG − F 2 , vemos que el ´area se puede calcular tambi´en conociendo la primera forma fundamental de la superficie. Ejercicio 229. Calc´ ulese el ´ area del toro que se obtiene al hacer girar la circunferencia de radio a y centro el punto (0, b, 0) alrededor del eje OZ.

6.2.

La segunda forma fundamental

Dado que xu y xv son, por hip´otesis, siempre linealmente independientes, disponemos de un vector normal unitario a la superficie N=

xu × xv . kxu × xv k

100

CAP´ITULO 6. LAS FORMAS FUNDAMENTALES

Para estudiar la forma de la superficie en un entorno del punto x(u0 , v0 ) es natural considerar la altura (con signo) a la que est´a cada punto respecto al plano tangente a la superficie en x(u0 , v0 ), es decir, la funci´on h(u, v) = hx(u, v) − x(u0 , v0 ), N(u0 , v0 )i. Usando el desarrollo de Taylor alrededor de (u0 , v0 ) tenemos que h(u, v) =hxu (u0 , v0 )∆u + xv (u0 , v0 )∆v  1 + xuu (u0 , v0 )∆u2 + 2xuv (u0 , v0 )∆u∆v + xvv (u0 , v0 )∆v 2 , N(u0 , v0 )i 2 + O(∆u2 + ∆v 2 ) 1 = hxuu (u0 , v0 ), N(u0 , v0 )i∆u2 2 + hxuv (u0 , v0 ), N(u0 , v0 )i∆u∆v  + hxvv (u0 , v0 ), N(u0 , v0 )i∆v 2 , + O(∆u2 + ∆v 2 ) donde ∆u = u − u0 y ∆v = v − v0 . Definimos los coeficientes de la segunda forma fundamental como e = hxuu , Ni, f = hxuv , Ni y g = hxvv , Ni. y la segunda forma fundamental en un punto (u, v), II(u,v) , es la aplicaci´on (∆u, ∆v) 7→ II(u,v) (∆u, ∆v) = e(u, v)∆u2 + 2f (u, v)∆u∆v + g(u, v)∆v 2 . Por tanto

1 h(u, v) = II(u0 ,v0 ) (∆u, ∆v) + O(∆u2 + ∆v 2 ), 2 es decir, la segunda forma fundamental es, salvo el factor 1/2, la primera aproximaci´on no nula a la funci´on altura sobre la normal en el punto de par´ametros (u0 , v0 ). Viene dada por una forma cuadr´atica cuyos coeficientes son e, f y g. Si la segunda forma fundamental es definida, lo que sucede cuando su determinante eg − f 2 > 0, como es (salvo el factor 1/2) la aproximaci´on de orden 2 del desarrollo de Taylor de h, tenemos que en un entorno reducido de (u0 , v0 ) la funci´on altura es estrictamente positiva o bien estrictamente negativa. En cualquiera de los dos casos, esto significa que la superficie se encuentra, en un entorno del punto, a un solo lado del plano tangente y s´olo toca al plano tangente en el punto de tangencia. Los puntos de la superficie con segunda forma fundamental definida se llaman puntos el´ıpticos.

6.2. LA SEGUNDA FORMA FUNDAMENTAL

101

Si, por el contrario, es eg −f 2 < 0 entonces la segunda forma fundamental es no definida y no degenerada. Esto significa que en todo entorno reducido del punto en cuesti´on hay puntos en los que la superficie est´a a un lado del plano tangente y puntos en los que se encuentra al lado contrario. Los puntos en los que la segunda forma fundamental es no definida y no degenerada se llaman hiperb´olicos. Finalmente, si la segunda forma fundamental es degenerada, es decir, eg − f 2 = 0, pero no es nula, entonces se dice que el punto es parab´olico. Dado que habr´a una recta a lo largo de la cual la segunda forma fundamental se anula, existe la posibilidad de que t´erminos de orden superior dominen, por lo que no se puede asegurar nada sobre la posici´on de la superficie respecto al plano tangente. El u ´ltimo caso a considerar es el de la anulaci´on completa de la segunda forma fundamental en el punto. En este caso el punto se llama plano (lo cual no significa que la superficie sea un plano). Ejercicio 230. Dar un ejemplo de superficie con un punto parab´ olico en el que la misma se encuentra a un solo lado del plano tangente y corta al mismo en una curva. Idem., estando a un s´ olo lado y cortando al mismo en un s´ olo punto. Idem., estando a los dos lados del plano tangente. Ejercicio 231. Dar un ejemplo de superficie en el que un punto es plano, pero que no sea plana en ning´ un entorno del mismo. Podemos tambi´en considerar una curva x = x(u(t), v(t)) y calcular la altura h = h(t) respecto al punto de par´ametro t0 . Realizando unos c´alculos similares a los anteriores, vemos que h(t0 ) = h′ (t0 ) = 0 y 1 h(t) = II(u′ (t0 ), v ′ (t0 ))(t − t0 )2 + O((t − t0 )2 ). 2 Ejercicio 232. Dar los detalles de la demostraci´ on de esta u ´ltima afirmaci´on. N´otese que II(u′ , v ′ ) resulta ser una forma cuadr´atica en u′ , v ′ . Esto nos ser´a u ´til para responder a preguntas tales como ¿cu´ales son las direcciones en las que la funci´on altura var´ıa con mayor o menor rapidez?, ya que tendremos a nuestra disposici´on la teor´ıa de formas cuadr´aticas. Pero antes de proseguir, observemos que derivando las ecuaciones hxu , Ni = hxv , Ni = 0 obtenemos las relaciones e = hxuu , Ni = −hxu , Nu i f = hxuv , Ni = −hxu , Nv i = −hxv , Nu i g = hxvv , Ni = −hxv , Nv i que, en ocasiones, simplificar´an algunos c´alculos te´oricos.

CAP´ITULO 6. LAS FORMAS FUNDAMENTALES

102

Ejercicio 233. Calc´ ulese la segunda forma fundamental de las siguientes superficies parametrizadas: 1. x (u, v) = (a sin u cos v, b sin u sin v, c cos u) , (elipsoide). 2. x (u, v) = (au cos v, bu sin v, u2 ) (paraboloide el´ıptico). 3. x (u, v) = (au cosh v, bu sinh v, u2 ) (paraboloide hiperb´olico). 4. x (u, v) = (a sinh u cos v, b sinh u sin v, c cosh u) (hiperboloide de dos hojas).

6.2.1.

Curvatura normal y curvatura geod´ esica

Sea x = x(s) = x(u(s), v(s))2 una curva sobre la superficie, que suponemos, para simplificar los c´alculos, que est´a parametrizada por el arco. Su aceleraci´on coincide entonces con su vector curvatura, k(s) = x′′ (s), y en cada punto de la curva podemos descomponer el vector de curvatura en suma de dos vectores, k(s) = kg (s) + kn (s) donde kg (s) es tangente a la superficie y kn (s) es perpendicular a la misma. Llamaremos a kg vector de curvatura geod´esica y a kn , vector de curvatura normal. No es dif´ıcil calcularlos si conocemos el vector normal N(s) = N(u(s), v(s)): tendremos que kn (s) = κn (s)N(s) y, por tanto, hk, Ni = hkn , Ni = κn . En consecuencia, kn = κn N = hx′′ (s), NiN kg = x′′ (s) − hx′′ (s), NiN N´otese que la curva es una geod´esica si y s´olo si su vector de curvatura geod´esica es nulo. Como hx′ , Ni = 0, derivando esta identidad obtenemos que hx′′ , Ni = −hx′ , N′ i, y como N ′ = N u u ′ + N v v ′ y x ′ = x u u′ + x v v ′ , 2

N´ otese que la curva x(s) est´ a parametrizada por el arco, pero no necesariamente la curva en el espacio de par´ ametros (u(s), v(s))

6.2. LA SEGUNDA FORMA FUNDAMENTAL

103

tenemos que κn = −hx′ , N′ i = −hxu , Nu i(u′ )2 + (−hxu , Nv i − hxv , Nu i)u′ v ′ − hxv , Nv i(v ′ )2 = e(u′ )2 + 2f u′ v ′ + g(v ′ )2 . Es decir, la segunda forma fundamental, actuando sobre un vector unitario, puede interpretarse como la curvatura normal de una curva de la superficie con dicho vector tangente. En el caso de que la curva no est´e parametrizada por el arco, esto no es un gran problema: podemos sustituir la coordenadas (u′ , v ′ ) por (u′ /kx′ (t)k, v ′ /kx′ (t)k), que s´ı se corresponden con las del vector unitario de la curva reparametrizada por la longitud de arco. Quedar´ıa entonces que κn =

e(u′ )2 + 2f u′ v ′ + g(v ′ )2 II(u′ , v ′ ) = . kx′ (t)k2 I(u′ , v ′ )

Es decir, la curvatura normal viene dada por el cociente de la segunda forma fundamental entre la primera. Otro resultado interesante que acabamos de demostrar es que la curvatura normal no depende m´as que de la direcci´on tangente a la curva en el punto del que se trate o, dicho de otra forma, todas las curvas en una superficie que tengan la misma recta tangente en un punto tendr´an en dicho punto la misma curvatura normal. Este resultado se conoce como teorema de Meusnier.

6.2.2.

Secciones normales

Como acabamos de ver que la curvatura normal s´olo depende de la recta tangente, dado un punto de una superficie podemos considerar el haz de planos normales a la misma que pasan por dicho punto. Los cortes de dichos planos con la superficie nos proporcionar´an curvas3 con todas las direcciones tangentes posibles. A los cortes de los planos normales a la superficie en un punto con la superficie misma les llamaremos secciones normales. Este haz de planos genera una familia de curvas planas y calcular la curvatura de las mismas en el punto base no es dif´ıcil. Dada una de estas secciones normales a la superficie, el vector de curvatura k debe estar contenido en el plano normal, simplemente por ser la curva plana. Por ser vector de curvatura, debe ser tambi´en normal al vector tangente, as´ı que la u ´nica posibilidad 3

Estas curvas estar´an bien definidas al menos en un entorno del punto, lo que es suficiente para nuestros prop´ositos.

104

CAP´ITULO 6. LAS FORMAS FUNDAMENTALES

es que sea proporcional al vector normal a la superficie, por lo tanto (en el punto base): k = kn = κn N. En consecuencia la curvatura de la secci´on normal kkk = k coincide con el valor absoluto de κn . Ejercicio 234. Como la secciones normales son curvas planas, es posible orientar el plano y dar un signo a la curvatura, que llamaremos tambi´en k. ¿Es posible orientar el plano de forma que k = κn .

6.2.3.

Diagonalizaci´ on de la segunda forma fundamental

Consideremos el c´ırculo unidad c en el plano tangente a la superficie en un punto de par´ametros (u, v). Podemos considerar la segunda forma fundamental restringida a dicho c´ırculo, II(u,v) |c : c → R. que, por ser una aplicaci´on continua definida sobre un conjunto compacto, debe alcanzar su m´aximo en alg´ un punto e1 ∈ c. El vector unitario e1 puede ser completado con otro vector e2 ∈ c de forma que {e1 , e2 } formen una base ortonormal. Veamos cu´al es la expresi´on de la segunda forma fundamental en esta base: si x′ es un vector tangente en el punto de par´ametros (u, v), tendremos que x′ = xe1 + ye2 y ser´a II(x′ ) = αx2 + 2βxy + γy 2 para ciertos coeficientes α, β y γ. Dado que x′ ∈ c ser´a x′ = cos θe1 + sin θe2 y como e1 es un m´aximo, la funci´on f (θ) = α cos2 θ + 2β cos θ sin θ + γ sin2 θ tiene un m´aximo en θ = 0. Por tanto satisface que f ′ (0) = 2β = 0. En consecuencia, II(x′ ) = αx2 + γy 2 . Nuevamente utilizando que e1 es m´aximo vemos que α = II(e1 ) ≥ II(e2 ) = γ. Tenemos pues las desigualdades α(x2 + y 2 ) ≥ αx2 + γy 2 ≥ γ(x2 + y 2 ). Por tanto si x′ ∈ c, i.e., x2 + y 2 = 1, concluimos que α ≥ αx2 + γy 2 ≥ γ,

6.2. LA SEGUNDA FORMA FUNDAMENTAL

105

es decir, la segunda forma fundamental restringida al c´ırculo unidad alcanza su m´ınimo γ en la direcci´on ortogonal a la del m´aximo, es decir, la de e2 . Finalmente notemos que ambos valores coinciden s´olo cuando α = γ, en cuyo caso la segunda forma fundamental resulta ser constante sobre el c´ırculo unidad. Recordemos que si x′ ∈ c entonces II(x′ ) es el valor de la curvatura normal κn en la direcci´on de x′ . Hemos demostrado pues que la curvatura normal en un punto de la superficie o bien es constante o bien tiene dos valores extremos distintos. En la notaci´on anterior, el m´aximo era α y el m´ınimo γ, pero es costumbre denotarlos como κ1 y κ2 , respectivamente. Estos valores se alcanzan en dos direcciones ortogonales y se llaman curvaturas principales. Las direcciones correspondientes se denominan direcciones principales. Si la segunda forma fundamental es constante sobre el c´ırculo unidad, entonces diremos que κ1 = κ2 es este valor com´ un y las direcciones principales no est´an definidas. Los puntos en los que ocurre esto se llaman puntos umbilicales. Tenemos ya los ingredientes necesarios para definir las dos principales medidas de la curvatura de una superficie en un punto: la curvatura de Gauss, K = κ1 κ2 , y la curvatura media

κ1 + κ2 . 2 La primera interpretaci´on de estas dos importantes funciones es algebraica: son, respectivamente, el determinante y la traza de la segunda forma fundamental respecto a la base ortonormal {e1 , e2 }. H=

6.2.4.

C´ alculos expl´ıcitos

Veamos ahora c´omo podemos calcular, en t´erminos de una parametrizaci´on, las direcciones principales y las curvaturas principales. Buscamos una direcci´on tangente x′ = u′ xu + v ′ xv que haga m´axima la curvatura normal κn =

II(u′ , v ′ ) e(u′ )2 + 2f u′ v ′ + g(v ′ )2 = . I(u′ , v ′ ) E(u′ )2 + 2F u′ v ′ + G(v ′ )2

Supongamos que no es direcci´on principal ni la de xu ni la de xv , de forma que si dividimos numerador y denominador por (u′ )2 y llamamos m = v ′ /u′ tenemos que la funci´on κn (m) =

e + 2f m + gm2 E + 2F m + Gm2

CAP´ITULO 6. LAS FORMAS FUNDAMENTALES

106

debe satisfacer, en la direcci´on principal, que dκn /dm = 0. De esto se obtiene que dκn 2(f + gm)(E + 2F m + Gm2 ) − 2(F + Gm)(e + 2f m + gm2 ) = = 0, dm (E + 2F m + Gm2 )2 de donde se obtiene que (f + gm) − (F + Gm)κn = 0, es decir, f u′ + gv ′ − (F u′ + Gv ′ )κn = 0.

(6.1)

eu′ + f v ′ − (Eu′ + F v ′ )κn = 0.

(6.2)

Si hacemos unos c´alculos an´alogos para m = u′ /v ′ obtenemos que

Es decir, poniendo en forma matricial estas dos u ´ltimas ecuaciones:   ′     u E F e f = 0, − κn v′ F G f g o, equivalentemente, si κn es curvatura principal y (u′ , v ′ ) es una direcci´on principal, entonces est´a en el n´ ucleo de II − κn I. Equivalentemente, podemos decir que (u′ , v ′ ) es un autovector de I−1 II con autovalor κn . De aqu´ı resulta que las dos curvaturas principales se obtienen resolviendo la ecuaci´on de segundo grado en κn dada por  det(II − κn I) = EG − F 2 k 2 − (Eg − 2 F f + Ge)k + eg − f 2 = 0,

y los dos autovectores asociados a las soluciones son las coordenadas de las direcciones principales. Como EG − F 2 6= 0, podemos dividir la ecuaci´on por este determinante, obteniendo que las curvaturas principales son soluci´on de k2 −

(Eg − 2 F f + Ge) eg − f 2 k + = 0, (EG − F 2 ) (EG − F 2 )

de donde deducimos que Eg − 2 F f + Ge 2 (EG − F 2 ) eg − f 2 . K= (EG − F 2 ) H=

´ 6.3. L´INEAS DE CURVATURA Y L´INEAS ASINTOTICAS

107

Ejercicio 235. Calc´ ulese la primera y segunda forma fundamental, as´ı como la curvatura de Gauss, la curvatura media y las curvaturas principales de: una esfera, un elipsoide de rotaci´ on, un hiperboloide de rotaci´ on de una hoja, un hiperboloide de rotaci´ on de dos hojas, un paraboloide de rotaci´ on, un cilindro de secci´on circular, un cono circular sin v´ertice, un toro, una catenoide de rotaci´ on, una seudoesfera y un helicoide recto. Ejercicio 236. Demu´estrese que si una superficie est´ a contenida en un plano entonces sus curvaturas principales son nulas. Ejercicio 237. Calc´ ulese las curvaturas principales en los v´ertices del hiperboloide de dos hojas x2 /a2 − y 2 /b2 − z 2 /c2 = 1. Ejercicio 238. Calcular las curvaturas principales de la superficie z = xy en el punto (1, 1, 1). Ejercicio 239. Calcular las curvaturas principales de la superficie x2 /p + y 2 /q = 2z en el punto (0, 0, 0). Ejercicio 240. Demu´estrese que en todo punto de la superficie x = u cos v, y = u sin v, z = λu una de las secciones normales principales es una recta. Ejercicio 241. En la superficie x = u2 +v 2 , y = u2 −v 2 , z = uv se considera el punto P de par´ ametros u = v = 1. Calc´ ulense las curvaturas principales de la superficie en P , las ecuaciones de las secciones normales principales y las rectas tangentes correspondientes en P y la curvatura de la secci´on normal en P que es tangente a la curva v = u2 . Ejercicio 242. Demu´estrese que un punto es umbilical si y s´ olo si e/E = f /F = g/G. Ejercicio 243. Calc´ ulense los puntos umbilicales del elipsoide de semiejes a, b, c.

6.3.

L´ıneas de curvatura y l´ıneas asint´ oticas

En un punto no umbilical de la superficie tenemos pues bien definidas dos direcciones principales, en las que las secciones normales alcanzan sus valores m´aximo y m´ınimo. Las curvas con la propiedad de que su direcci´on tangente es siempre una direcci´on principal se llaman l´ıneas de curvatura. Para calcularlas, usemos las ecuaciones (6.1) y (6.2) para obtener f u′ + gv ′ eu′ + f v ′ = , F u′ + Gv ′ Eu′ + F v ′

CAP´ITULO 6. LAS FORMAS FUNDAMENTALES

108 es decir,

(f E − eF )(u′ )2 + (gE − eG)u′ v ′ + (gF − f G)(v ′ )2 = 0, que son las ecuaciones diferenciales de las l´ıneas de curvatura. Por otra parte, en el caso de que la curvatura de Gauss sea negativa entonces existen dos direcciones independientes en las cuales la curvatura normal es nula, justamente aquellas que anulan la segunda forma fundamental. Estas direcciones se llaman direcciones asint´oticas y satisfacen las ecuaciones diferenciales e(u′ )2 + 2f u′ v ′ + g(v ′ )2 = 0. Por cada punto de una superficie con curvatura de Gauss negativa pasan dos l´ıneas asint´oticas (salvo reparametrizaciones de las mismas). Ejercicio 244. Demu´estrese que si una curva es l´ınea de curvatura o asint´otica, entonces tambi´en lo son las reparametrizaciones de la misma. ¿Por qu´e? Ejercicio 245. Calc´ ulense los puntos umbilicales y las l´ıneas de curvatura de la superficie x(u, v) = (u, v, u2 + v 2 ). Ejercicio 246. Demu´estrese que las curvas coordenadas de la superficie       u+v u−v u+v (u−v)/2 (u−v)/2 x(u, v) = e cos ,e sin , 2 2 2 son l´ıneas asint´ oticas. Ejercicio 247. Demu´estrese que una geod´esica que es l´ınea asint´ otica Ejercicio 248. Est´ udiese si los paralelos y meridianos del toro de revoluci´ on son l´ıneas de curvatura o l´ıneas asint´ oticas. Ejercicio 249. Demu´estrese que una condici´ on necesaria y suficiente para que las curvas coordenadas de una parametrizaci´on sean l´ıneas de curvatura en un entorno de un punto no umbilical es que F = f = 0. Ejercicio 250. Dada la superficie  x(u, v) = u − u3 /3 + uv 2 , v − v 3 /3 + vu2 , u2 − v 2 ,

calc´ ulense las curvaturas principales en un punto arbitario. Demu´estrese que las l´ıneas de curvatura son las curvas coordenadas y que las l´ıneas asint´ oticas est´ an dadas por u + v = cte., u − v = cte.. Ejercicio 251. Calc´ ulense las l´ıneas asint´ oticas del hiperboloide x2 +y 2 −z 2 = 1.

´ GEOMETRICA ´ 6.4. UNA INTERPRETACION DE LA CURVATURA DE GAUSS109

6.4.

Una interpretaci´ on geom´ etrica de la curvatura de Gauss

En esta secci´on estudiaremos la definici´on original de curvatura dada por Gauss. Sea R sea un disco cerrado de radio r centrado en un punto (u0 , v0 ) y contenido en el dominio de la superficie parametrizada x : U ⊂ R2 → R3 . Consideremos la normal unitaria N : U → R3 asociada a la parametrizaci´on, es decir, xu × xv N= . kxu × xv k

Como kNk ≡ 1, la imagen N(R) ⊂ S 2 , siendo S 2 la esfera unitaria. La definici´on original de curvatura de una superficie en un punto que propuso Gauss es la siguiente: Area(N(R)) l´ım . (6.3) r→0 Area(x(R))

Ejercicio 252. Calc´ ulese directamente el valor de este l´ımite cuando la superficie es a) un plano, b) una esfera de radio ρ, c) un cilindro. Tenemos que Area(N(R)) =

Z

R

kNu × Nv kdudv.

Calculemos Nu × Nv . Puesto que hN, Ni = 0, tenemos que hNu , Ni = hNv , Ni = 0, es decir, Nu y Nv son tangentes a la superficie. Por tanto se pueden escribir en t´erminos de xu y xv : Nu = a11 xu + a21 xv Nv = a12 xu + a22 xv para ciertas funciones aij = aij (u, v) que ahora vamos a determinar. Si escribimos las u ´ltimas ecuaciones de forma matricial, tenemos que   Nu Nv = xu xv W, (6.4)

siendo

  a11 a12 W= . a21 a22

CAP´ITULO 6. LAS FORMAS FUNDAMENTALES

110

Usando las propiedades del producto vectorial, vemos que Nu × Nv = (a11 a22 − a12 a21 ) xu × xv = det W · xu × xv . Calculemos la matriz W. Para ello, escribamos  de forma algo  la ecuaci´on (6.4) m´as compacta denotando dN = Nu Nv , dx = xu xv , con lo que dN = dx · W.

(6.5)

Como dx⊤ dx = I, y dx⊤ dN = −II,

multiplicando por la izquierda por dx⊤ ambos miembros de (6.5) obtenemos que dx⊤ dN = −II = dx⊤ dx · W = I · W, y, por tanto, W = −I−1 II. Recordemos que la matriz −W apareci´o anteriormente cuando calcul´abamos las curvaturas principales. La matriz W se llama matriz de Weingarten. Ahora podemos calcular det W = det(−I−1 II) =

det II = K, det I

y por tanto Area(N(R)) =

Z

R

kNu × Nv kdudv =

Z

R

|K|kxu × xv kdudv.

Entonces el l´ımite (6.3) se puede escribir como R |K|kxu × xv kdudv Area(N(R)) = l´ım RR l´ım r→0 r→0 Area(x(R)) kxu × xv kdudv R |K(u∗ , v ∗ )|kxu × xv k(u∗ , v ∗ )Area(R) = l´ım , r→0 kxu × xv k(u∗∗ , v ∗∗ )Area(R) donde (u∗ , v ∗ ), (u∗∗ , v ∗∗ ) ∈ R son ciertos puntos cuya existencia est´a asegurada gracias al teorema del valor medio del c´alculo integral. Por tanto, obtenemos que Area(N(R)) |K(u0 , v0 )| = l´ım . r→0 Area(x(R))

6.5. EL TEOREMA EGREGIO DE GAUSS

111

Es posible dar una interpretaci´on al signo de la curvatura de Gauss. Para ello consideremos una parametrizaci´on cr del borde del disco R positivamente orientada en el plano, es decir, que gira en sentido antihorario. Diremos que la aplicaci´on normal N : U → S 2 preserva la orientaci´on en (u0 , v0 ) cuando la curva N ◦ cr en la esfera S 2 gira alrededor del punto N(x(u0 , v0 )) en el mismo sentido que la curva x ◦ cr gira en la superficie alrededor del punto x(u0 , v0 ). Para poder comparar ambos sentidos de giro, observemos que el plano tangente a la esfera en el punto N(u0 , v0 ) es el mismo que el de la superficie en el punto de par´ametros (u0 , v0 ). Ejercicio 253. Demu´estrese esta u ´ltima afirmaci´on. Podemos escoger como lado para observar la superficie aquel hacia el que se˜ nala N(u0 , v0 ) y, desde esta posici´on, decidir si ambas curvas giran en el mismo sentido. En el caso de que sea as´ı, asignaremos un signo positivo al ´area de Area(N(R)) y signo negativo en caso contrario. Pues bien, se puede demostrar que el que ambas orientaciones coincidan en entornos arbitrariamente peque˜ nos del punto es equivalente a que la curvatura de Gauss sea positiva, con lo que resulta Area(N(R)) . r→0 Area(x(R))

K(u0 , v0 ) = l´ım

Ejercicio 254. Parametr´ıcese un entorno del polo norte de la esfera x2 + y 2 + z 2 = 1 y consid´erese la circunferencia cr (θ) = (r cos θ, r sin θ). Dib´ ujese la curva N ◦ cr y dec´ıdase si tiene o no la misma orientaci´ on que la curva x ◦ cr . ¿Cu´al es el signo de la curvatura de Gauss? H´agase el mismo ejercicio para la superficie z = xy en el punto (0, 0, 0).

6.5.

El teorema egregio de Gauss

Tanto la definici´on de curvatura que hemos dado en funci´on de las curvaturas principales como la definici´on original de Gauss parecen depender del ambiente en el que se encuentra la superficie, ya que ambas hacen uso del vector normal a la superficie. El teorema egregio de Gauss demuestra un resultado sorprendente: la curvatura de Gauss depende en realidad s´olo de la primera forma fundamental, es decir, es un objeto perteneciente a la geometr´ıa intr´ınseca de la superficie. Para demostrar este resultado, fij´emonos

112

CAP´ITULO 6. LAS FORMAS FUNDAMENTALES

en el numerador de la curvatura de Gauss: eg − f 2 = hxuu , Nihxvv , Ni − hxuv , Ni2 = (hxuu , xu × xv ihxvv , xu × xv i − hxuv , xu × xv i)/(EG − F 2 ) = (det(xuu , xu , xv ) det(xvv , xu , xv ) − det(xuv , xu , xv )2 )/(EG − F 2 ), (6.6) y el numerador de esta u ´ltima expresi´on se puede escribir como  ⊤  ⊤ xuv xuu   ⊤ ⊤   det xuv x⊤  xv det xu det xvv xu xv − det x⊤ u u xv xv  ⊤    ⊤ ⊤ ⊤ ⊤ xuu xvv x⊤ x x x x x x x x x u v uv u v uu uu uv uv uv  − det  x⊤ . x⊤ x⊤ x⊤ x⊤ = det  x⊤ u xvv u xu u xv u xuv u xu u xv ⊤ ⊤ ⊤ ⊤ ⊤ ⊤ xv xvv xv xu xv xu xv xuv xv xu xv xu  ⊤    ⊤ ⊤ ⊤ xuu xvv x⊤ x⊤ uu xu xuu xv uv xuv xuv xu xuv xv E F  − det  x⊤ E F  = det  x⊤ u xvv u xuv ⊤ ⊤ xv xvv F G xv xuv F G

(6.7)

Notemos ahora que

1 x⊤ uu xu = Eu 2 1 ⊤ xvv xv = Gv 2 ⊤ ⊤ x⊤ x = (x uu v u xv )u − xu xvu 1 = Fu − Ev 2 ⊤ ⊤ xu xvv = (xu xv )v − x⊤ uv xv 1 = F v − Gu 2 y que 1 x⊤ uv xu = Ev 2 1 x⊤ uv xv = Gu . 2 De esta forma, todos los coeficientes de las matrices que aparecen en la ecuaci´on (6.6) se escriben en t´erminos de la primera forma fundamental salvo ⊤ x⊤ uu xvv y xuv xuv . Pero notemos que estos dos coeficientes tienen el mismo

113

6.5. EL TEOREMA EGREGIO DE GAUSS

menor complementario, de forma que podemos escribir     ⊤ ⊤ ⊤ ⊤ 0 x⊤ xuu xvv − x⊤ uv xu xuv xv uv xuv xuu xu xuu xv E F  x⊤ E F −det x⊤ (eg−f 2 )(EG−F 2 ) = det  u xuv u xvv ⊤ ⊤ xv xuv F G xv xvv F G (6.8) y observar que la diferencia ⊤ ⊤ ⊤ x⊤ uu xvv − xuv xuv = (xu xvv )u − xu xvvu

⊤ − (x⊤ u xuv )v − xu xuvv



⊤ = (x⊤ u xvv )u − (xu xuv )v 1 1 1 = (Fv − Gu )u − ( Ev )v = − (Evv − 2Fuv + Guu ). 2 2 2

Por tanto K=

=

eg − f 2 2 EG 1−F − (Evv − 2Fuv + Guu ) 2 Fv − 12 Gu 1 G 2 v

1 E 2 u

Fu − 21 Ev 0 − 1 Ev E F 12 Gu F G 2 (EG − F 2 )2

1 E 2 v

E F

1 G 2 u



(6.9) F G

Por tanto hemos visto que para calcular la curvatura de Gauss es necesario u ´nicamente el conocimiento de la primera forma fundamental, es decir, de medidas intr´ınsecas sobre la superficie. Ejercicio 255. Demu´estrese que en el caso de que las l´ıneas coordenadas sean ortogonales la curvatura de Gauss viene dada mediante la expresi´ on      1 Ev Gu √ K = −√ . − √ EG EG v EG u

Sponsor Documents

Or use your account on DocShare.tips

Hide

Forgot your password?

Or register your new account on DocShare.tips

Hide

Lost your password? Please enter your email address. You will receive a link to create a new password.

Back to log-in

Close