C ÁLCULO D IFERENCIAL EN Rn Francisco Javier Pérez González Departamento de Análisis Matemático Universidad de Granada Septiembre de 2012
I
Licencia. Este texto se distribuye bajo una licencia Creative Commons en virtud de la cual se permite: Copiar, distribuir y comunicar públicamente la obra. obr a. Hacer obras derivadas. Bajo las condiciones siguientes:
Reconocimiento. Debe reconocer los créditos de la obra de la manera especificada por el autor o BY:
el licenciador (pero no de una manera que sugiera que tiene su apoyo o apoyan el uso que hace de su obra).
$ \
C
No comercial. No puede utilizar esta obra para fines comerciales. transfor ma esta obra, o genera una obra derivada, Compartir bajo la misma licencia. licencia. Si altera o transforma sólo puede distribuir la obra generada bajo una licencia idéntica a ésta.
Universidad de Granada Dpto. de Análisis Matemático
Prof. Javier Pérez Cálculo Diferencial en Rn
´Indice general
1. Estructura Estructura euclídea euclídea y topologí topologíaa de Rn
1
1.1. .1. Pro Producto esc escalar lar y norma euc euclí líddeos eos . . . . . . . . . . . . . . . . . . . . . . .
1
1.2. .2. Esp Espacio cios normad mados y esp espacio acioss métr métric icoos . . . . . . . . . . . . . . . . . . . . .
4
1.3. .3. Topolog logía de un espa spacio cio métr métric icoo . . . . . . . . . . . . . . . . . . . . . . . . .
6
1.4. Continuidad . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 12 1.5. Límite funcional . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 20 1.6. 1.6. Contin Continuid uidad ad y lím límite itess de campo camposs escala escalares res y vecto vectoria riales les . . . . . . . . . . . . . 21
2. Derivadas Derivadas parciales y extremos extremos relativos relativos de campos escalares escalares
30
2.1. 2.1. Deri Derivvadas adas parc parcia iale les. s. Vecto ectorr grad gradie ient ntee . . . . . . . . . . . . . . . . . . . . . . . 31 2.1.1. 2.1.1. Inter Interpre pretac tació iónn geo geomét métric ricaa de las deri deriva vada dass parc parcial iales es . . . . . . . . . . . 32 2.1.2 .1.2.. Cam Campos escal scalar ares es dife iferen renciab iables les . . . . . . . . . . . . . . . . . . . . . 33 2.2. .2. Rec Rectas tas tan tangente entess y plan lanos tan tangent entes . . . . . . . . . . . . . . . . . . . . . . . 38 2.3. .3. Der Derivadas parcia rciale less de orden sup superio rior . . . . . . . . . . . . . . . . . . . . . . 42 2.4. 2.4. Teore eorema ma de Taylo aylorr. Extr Extrem emos os rela relati tivvos . . . . . . . . . . . . . . . . . . . . . . 46
3. Derivación Derivación de campos vectoriales vectoriales
54
3.1. 3.1. Deri Derivvada ada de un camp campoo vect vector oria ial. l. Matr Matriz iz jaco jacobi bian anaa . . . . . . . . . . . . . . . . 54 3.1.1. 3.1.1. El espacio espacio normado normado L.Rn ; Rm / . . . . . . . . . . . . . . . . . . . . . 55 3.1.2. 3.1.2. Regla Regla de la cadena. cadena. Deriva Derivadas das parciale parcialess de funcion funciones es compu compuesta estass . . . 59 3.1.3. 3.1.3. Teorema eorema del del valo valorr medio medio para para campo camposs vecto vectoria riales les . . . . . . . . . . . . 64
II
Índice general
III
3.2. Teorema de la función inversa . . . . . . . . . . . . . . . . . . . . . . . . . . 65 3.3. Teorema de la función implícita . . . . . . . . . . . . . . . . . . . . . . . . . 71 3.4. Variedades diferenciables en Rn . . . . . . . . . . . . . . . . . . . . . . . . . 81 3.5. Espacios tangente y normal . . . . . . . . . . . . . . . . . . . . . . . . . . . . 87 3.6. Extremos condicionados . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 90 3.6.1. Cálculo de extremos absolutos . . . . . . . . . . . . . . . . . . . . . . 98
Universidad de Granada Dpto. de Análisis Matemático
Prof. Javier Pérez Cálculo Diferencial en Rn
Cap´ıtulo
1
Estructura eucl´ıdea y topolog´ıa de Rn
1.1. Producto escalar y norma euclídeos Como sabes, Rn es un espacio vectorial en el que suele destacarse la llamada base canónica formada por los vectores fe1 ; e2 ; : : : ; en g donde e k es el vector cuyas componentes son todas nulas excepto la que ocupa el lugar k que es igual a 1.
1.1 Definición. Dados dos vectores x D .x1 ; x2 ; : : : ; xn /, y D .y1 ; y2 ; : : : ; yn / se define su producto escalar por: n
˝ ˇ ˛D X xy
xk y k
D x1y1 C x2y2 C C xnyn
kD1
Este producto escalar se llama producto escalar euclídeo. Observa que el producto escalar de dos vectores no es un vector sino un número real. La notación x y es frecuentemente usada en los libros de Física para representar el producto escalar de los vectores x e y . .
Propiedades del producto escalar. 1. x x > 0, y x x
˝ ˇˇ ˛
0.
˝ ˇˇ ˛ D ” D ˝ ˇˇ ˛ D ˝ ˇˇ ˛ 2 ˝ C ˇˇ ˛ D ˝ ˇˇ ˛ C ˝ ˇˇ ˛ 0
x
2. Simetría. x y
y x
3. Linealidad . ˛ x
ˇy z
para todos x ; y Rn . ˛ xz
ˇ y z para todos ˛; ˇ R y para todos x ; y ; z Rn .
1
2
2
Producto escalar y norma euclídeos
2
Dichas propiedades del producto escalar se deducen fácilmente de su definición.
1.2 Definición. La norma euclídea de un vector x D .x1 ; x2 ; : : : ; xn / se define por
p
n
q ˇ X k k D ˝ ˇ ˛ D x
xx
2
xk2
kD1
Desigualdad de Cauchy-Schwarz. Para todos x ; y
2 Rn se verifica que
ˇˇ˝ ˇˇ ˛ˇˇ
xy 6 x
k k2ky k2
ˇ˝ ˇ ˛ˇ D k k k k
Además, supuesto que x e y no son nulos, la igualdad x y x 2 y 2 equivale a que exista un número R tal que x y (es decir, los vectores x e y están en una misma recta que pasa por el origen).
2
D
Propiedades de la norma euclídea. 1.
kxk2 > 0, y kxk2 D 0 ” x D 0. 2. Homogeneidad . kxk2 D jjkxk2 para todo x 2 Rn y todo 2 R. 3. Desigualdad triangular. Para todos x ; y 2 Rn se verifica que kx C y k2 6 kxk2 C ky k2 Además, supuesto que x e y no son nulos, la igualdad k x C y k2 D kxk2 C ky k2 equivale a que hay un número > 0 tal que x D y (es decir, los vectores x e y están en una misma semirrecta que pasa por el origen).
1.3 Definición. Se dice que los vectores x e y son ortogonales, y escribimos x ? y , cuando su producto escalar es cero. Se dice que un vector x es ortogonal a un conjunto de vectores E Rn cuando x es ortogonal a todo vector en E . Un conjunto de vectores no nulos que son mutuamente ortogonales se dice que es un conjunto ortogonal de vectores; si, además, los vectores tienen todos norma 1 se dice que es un conjunto ortonormal de vectores. Una base vectorial que también es un conjunto ortogonal (ortonormal) se llama una base ortogonal (ortonormal). Si x e y son vectores no nulos, el vector
Q
.x/ y
se llama proyección ortogonal de x sobre y .
˝ ˇˇˇ ˛ D˝ ˇ ˛ xy y y
y
El vector x y .x/ es ortogonal a y . En particular, si y es un vector unitario (de norma 1) entonces el vector x x y y es ortogonal a y .
Q ˝ˇ˛ ˇ
Universidad de Granada Dpto. de Análisis Matemático
Prof. Javier Pérez Cálculo Diferencial en Rn
Producto escalar y norma euclídeos
3
Ejercicios propuestos 1. Particulariza las definiciones y propiedades anteriores para el caso n D 1, es decir, para R. 2. Prueba la desigualdad de Cauchy-Schwarz. Sugerencia. Comprueba que la ecuación x y x y D 0, en la que es un número real arbitrario y x e y son vectores que se suponen fijos, es un trinomio de segundo grado en la variable . Ten en cuenta que dicho trinomio toma siempre valores mayores o iguales que cero lo que proporciona información sobre su discriminante.
˝
ˇ
˛
3. Prueba la desigualdad triangular. Sugerencia. Una estrategia para probar desigualdades entre normas euclídeas es elevar al 2 cuadrado. La desigualdad k x C y k22 6 kxk2 C ky k2 es equivalente a la desigualdad triangular pero es fácil de probar desarrollando k x C y k22 D x C y x C y y usando la desigualdad de Cauchy-Schwarz.
˝ ˇ ˛
4. Teorema de Pitágoras . Prueba que los vectores x e y son ortogonales si, y solo si
kx C y k22 D kxk22 C ky k22 5. Prueba que el vector x
Q
y
.x/ es ortogonal a y .
1.4 Definición. La distancia euclídea en Rn es la aplicación d2 W Rn Rn ! RCo definida por: d2 .x; y / D kx y k2
.x; y Rn /
2
La distancia euclídea entre los vectores x e y es el número d2 .x; y /. Las siguientes propiedades de la distancia euclídea se deducen fácilmente de las de la norma euclídea.
Propiedades de la distancia euclídea. 1. d2 .x; y / > 0, y d2 .x; y / D 0 ” x D y . 2. Simetría. d2 .x; y / D d2 .y ; x/ para todos x ; y 2 Rn . 3. Homogeneidad . d2 .x; y / D jj d2 .x; y / para todos x ; y 2 Rn y todo 2 R. 4. Desigualdad triangular. d2 .x; y / 6 d2 .x; z/ C d2 .z; y / para todos x ; y ; z 2 Rn .
Universidad de Granada Dpto. de Análisis Matemático
Prof. Javier Pérez Cálculo Diferencial en Rn
Espacios normados y espacios métricos
4
1.2. Espacios normados y espacios métricos Las propiedades de la norma y de la distancia euclídeas en Rn se pueden abstraer dando lugar a los conceptos de espacio normado y espacio métrico.
1.5 Definición. Sea X un espacio vectorial real. Una norma sobre X es una aplicación k k W X ! RCo que verifica las siguientes propiedades: 1. kxk D 0 ” x D 0. 2. Homogeneidad . kxk D jjkxk para todo 2 R y todo x 2 X . 3. Desigualdad triangular . kx C y k 6 kxk C ky k para todos x ; y 2 X . El par ordenado .X ; k k/ se llama un espacio normado. Naturalmente, sobre un mismo espacio vectorial pueden considerarse distintas normas, cada una de ellas da lugar a un espacio normado diferente. Para tener en cuenta este hecho se dice que un espacio normado es un par ordenado .X ; k k/ formado por un espacio vectorial real X y una norma. No obstante, con frecuencia se dice simplemente “sea X un espacio normado” y se sobreentiende que X es un espacio vectorial real en el que está definida una norma concreta.
1.6 Ejemplos. En Rn suelen considerarse, además de la norma euclídea, la norma de la suma, k k1 , y la norma del máximo, k k1, definidas para todo x 2 Rn por: n
X kk D j x
1
xk
kD1
j
kxk1 D maK x fjxk j W 1 6 k 6 ng En el espacio vectorial, B .A/, de todas las funciones reales acotadas definidas en un con junto no vacío A R, se define la norma uniforme dada para toda f 2 B .A/ por:
kf k1 D sup fjf .t /j W t 2 Ag : En el espacio vectorial, C .Œa; b/, de todas las funciones reales continuas definidas en un intervalo Œa; b, se define la norma integral de orden 1 dada para todo f 2 C .Œa; b/ por: b
kf k1
D jf .t /j d t a
1.7 Definición. Sea E un conjunto cualquiera no vacío. Una distancia en E es una aplicación d W E E ! RCo que verifica las siguientes propiedades: 1. d.x; y/ D 0 ” x D y . 2. Simetría. d .x; y/ D d.y; x/ para todos x ; y 2 E . 3. Desigualdad triangular . d.x; y/ 6 d.x; z/ C d.z; y/ para todos x ; y; z 2 E . Universidad de Granada Dpto. de Análisis Matemático
Prof. Javier Pérez Cálculo Diferencial en Rn
Espacios normados y espacios métricos
5
El par ordenado .E; d/ se llama un espacio métrico. Los elementos de un espacio métrico suelen llamarse puntos de dicho espacio métrico. Naturalmente, sobre un mismo conjunto pueden considerarse distintas distancias, cada una de ellas da lugar a un espacio métrico diferente. Para tener en cuenta este hecho se dice que un espacio métrico es un par ordenado .E; d/ formado por un conjunto no vacío y una distancia. No obstante, con frecuencia se dice simplemente “sea E un espacio métrico” y se sobreentiende que en E está definida una distancia concreta. Dado un espacio normado, .X ; k k/, la aplicación d W X X ! RCo dada por: d.x; y / D kx y k
2
.x; y X /
es una distancia en X que se llama distancia asociada a la norma. Todo espacio normado se considera siempre como espacio métrico con la distancia asociada a su norma.
Ejercicios propuestos 6. Prueba que k k 1 y k k 1 son normas en desigualdades:
Rn
y que para todo x 2 Rn se verifican las
kxk1 6 kxk2 6 kxk1 6 nkxk1 7. Prueba que k k1 es una norma en B .A/. 8. Sea n 2 N y S .n/ D f1; 2; : : : ; ng. Prueba que los espacios vectoriales B .S .n// y Rn son isomorfos. ¿Qué relación hay entre la norma uniforme y la del máximo en dichos espacios? 9. Prueba que k k1 es una norma en C .Œa; b/. 10. Sea .E; d/ un espacio métrico. Prueba la desigualdad:
jd.x; y/ d.y; z/j 6 d.x; z/
2
.x; y; z E/
(1.1)
11. Sea .E; d/ un espacio métrico y sean x1 ; x2 ; : : : ; xn puntos de E . Prueba que: n1
d.x1 ; xn / 6
X
d.xj ; xj C1 /
j D1
12. Prueba que la distancia asociada a una norma verifica las propiedades de la definición 1.7. Comprueba que dicha distancia es invariante por traslaciones y es homogénea. 13. Prueba que en todo espacio vectorial real se pueden definir normas y en todo conjunto se pueden definir distancias.
Universidad de Granada Dpto. de Análisis Matemático
Prof. Javier Pérez Cálculo Diferencial en Rn
Topología de un espacio métrico
6
1.3. Topología de un espacio métrico Sea .E; d/ un espacio métrico. Dados un punto a 2 E y un número r > 0 definimos el conjunto (al que por ahora no pondremos nombre): B.a; r /
D fx 2 E W d.x; a/ < r g
Observa que el conjunto así definido depende claramente de la distancia d por lo que una notación más apropiada sería Bd .a; r /, pero dicha notación es incómoda y solamente se usa cuando se consideran varias distancias diferentes en un mismo contexto.
1.8 Definición. Se dice que un conjunto A E es abierto en el espacio métrico .E; d/ si para cada punto x 2 A hay un número rx > 0 tal que B.x; rx / A. Por convenio, el conjunto vacío, Ø, se considera abierto. 1.9 Proposición. Sean .E; d/ un espacio métrico, a 2 E y r > 0. Se verifica que el conjunto B.a; r / es abierto. Dicho conjunto se llama bola abierta de centro a y radio r . 1.10 Proposición (Propiedades de los conjuntos abiertos de un espacio métrico) . En todo es pacio métrico .E; d/ se verifica que: 1. Los conjuntos E y Ø son abiertos. 2. La unión de cualquier colección de conjuntos abiertos es un conjunto abierto. 3. La intersección de una colección finita de conjuntos abiertos es un conjunto abierto. Las propiedades anteriores se expresan diciendo que si T es la clase de todos los conjuntos abiertos de un espacio métrico .E; d/ , entonces T es una topología en E . Se dice que dicha topología está asociada a la distancia d.
1.11 Definición. Se dice que un conjunto F E es cerrado en el espacio métrico .E; d/ si su complementario E nF es abierto en dicho espacio métrico. 1.12 Proposición (Propiedades de los conjuntos cerrados de un espacio métrico) . En todo espacio métrico .E; d/ se verifica que: 1. Los conjuntos E y Ø son cerrados. 2. La intersección de cualquier colección de conjuntos cerrados es un conjunto cerrado. 3. La unión de una colección finita de conjuntos cerrados es un conjunto cerrado.
1.13 Proposición. Sean .E; d/ un espacio métrico, a 2 E y r > 0. Se verifica que el conjunto B.a; r /
D fx 2 E W d.x; a/ 6 r g
es cerrado. Dicho conjunto se llama bola cerrada de centro a radio r .
Recuerda que todo espacio normado se considera automáticamente como espacio métrico con la distancia asociada a la norma, por lo que las definiciones y resultados anteriores tienen sentido también para espacios normados. En particular, la topología de un espacio normado es la topología asociada a la distancia asociada a la norma del mismo. En R se considera siempre la topología y la distancia asociadas al valor absoluto, es decir, R se considera siempre como el espacio normado . R; j j/. Universidad de Granada Dpto. de Análisis Matemático
Prof. Javier Pérez Cálculo Diferencial en Rn
Topología de un espacio métrico
7
Ejercicios propuestos 14. Prueba que si d es una distancia en E la función W E E ! R dada por: .x; y/
también es una distancia en E .
y/ D 1 Cd.x; d.x; y/
Sugerencia. Prueba que la aplicación t 7!
t 1
.x; y E/
2
es creciente en RCo .
C t
15. Prueba que un conjunto A E es abierto en el espacio métrico .E; d/ si, y sólo si, A es unión de bolas abiertas. 16. Describe las bolas abiertas y cerradas en los espacios normados .R; j j/, .R2 ; k .R2 ; k k1 /, . R2 ; k k1 / y . B .Œ0; 1/; k k1 /.
k2 /,
17. Da un ejemplo de una colección de conjuntos abiertos (resp. cerrados) en . R; j j/ cuya intersección (resp. unión) no sea un conjunto abierto (resp. cerrado). 18. Prueba que todo intervalo abierto (resp. cerrado) en rrado) en . R; j j/.
R
es un conjunto abierto (resp. ce-
19. Sean .Ei ; di / i D 1; 2; : : : ; n espacios métricos. En el conjunto producto cartesiano de los E i , E D E1 E2 En , definimos: ..x1 ; x2 ; : : : ; xn /; .y1 ; y2 ; : : : ; yn //
D maK x fdi .xi ; yi / W 1 6 i 6g
a) Prueba que .E; / es un espacio métrico. Se dice que .E; / es el espacio métrico producto de los espacios métricos .Ei ; di / i D 1; 2; : : : ; n. b) Describe las bolas abiertas en .E; /.
1.14 Definición. Sea .E; d/ un espacio métrico y A E .
Decimos que un punto x 2 E es adherente al conjunto A si toda bola abierta centrada en x tiene puntos de A. El conjunto de todos los puntos adherentes a A se llama la adherencia de A y se representa por A. Decimos que un punto x 2 E es un punto de acumulación del conjunto A si toda bola abierta centrada en x tiene puntos de A distintos de x . El conjunto de todos los puntos de acumulación de A se llama la acumulación de A y se representa por A 0 . El conjunto de todos los puntos adherentes a A y a E=A se llama la frontera de A y se representa por Fr.A/. Decimos que un punto x 2 A es un punto interior al conjunto A si hay alguna bola abierta centrada en x contenida en A. El conjunto de todos los puntos interiores de A se llama el ı
interior de A y se representa por A. Decimos que un punto x 2 A es un punto aislado en A si hay alguna bola abierta B.x; r / tal que B.x; r / \ A D fx g. Universidad de Granada Dpto. de Análisis Matemático
Prof. Javier Pérez Cálculo Diferencial en Rn
Topología de un espacio métrico
8
Ejercicios propuestos 20. Sea .E; d/ un espacio métrico y A E . Prueba que: 1. La adherencia de A es el más pequeño conjunto cerrado que contiene a A. Por tanto A es cerrado si, y sólo si A D A. 2. El interior de A es el más grande conjunto abierto contenido en A . 3. La frontera de A son los puntos adherentes a A que no son interiores de A . 4. La adherencia de A es la unión de A con la frontera de A .
21. Prueba que los puntos adherentes de A que no son de acumulación son puntos aislados de A . 22. Da ejemplos de conjuntos en . R; j j/ y en . R2 ; k k2 / que no sean abiertos ni cerrados. 23. Prueba que el supremo y el ínfimo de un conjunto no vacío y acotado de números reales son adherentes a dicho conjunto. 24. Calcula la adherencia en R de los conjuntos Q y R=Q. Conceptos topológicos En todo espacio métrico conviven dos estructuras, la propia de espacio métrico, .E; d/, y la de espacio topológico, .E; T /, donde T es la topología asociada a la distancia d. Sucede que en los espacios métricos es cómodo definir algunos conceptos usando la distancia, pero a veces dichos conceptos pueden formularse también únicamente en términos de abiertos, sin referencia ninguna a la distancia original. Los conceptos que pueden formularse únicamente en términos de abiertos se llaman topológicos. Es importante saber cuando un concepto es topológico porque puede ocurrir que distancias distintas sobre un mismo conjunto E den lugar a la misma topología y, por tanto, las propiedades que solamente dependen de la topología son las mismas para ambas distancias. Todos los conceptos introducidos en la definición 1.14 dependen aparentemente de la distancia porque en todos ellos se utilizan bolas abiertas. Pero es fácil comprobar que de hecho son conceptos topológicos. Para ello basta con que compruebes que en dichas definiciones puede sustituirse la expresión “bola abierta centrada en x” por “conjunto abierto que contiene a x ” sin que ello afecte para nada a los conceptos allí definidos.
1.15 Definición. Se dice que dos distancias d y sobre un mismo conjunto E son equivalentes si ambas definen la misma topología en E . Se dice que dos normas sobre un mismo espacio vectorial son equivalentes si sus distancias asociadas son equivalentes.
1.16 Proposición. Sean k k y jjj jjj dos normas sobre un espacio vectorial X . Equivalen las afirmaciones: a)
k k y jjj jjj son normas equivalentes.
Universidad de Granada Dpto. de Análisis Matemático
Prof. Javier Pérez Cálculo Diferencial en Rn
Topología de un espacio métrico
9
b) Existen números m > 0 , M > 0 verificándose que: m x 6 x
k k jjj jjj 6 M kx k
.x X /
2
Convergencia en espacios métricos 1.17 Definición. Se dice que una sucesión f xng de puntos de un espacio métrico .E; d/ es convergente si hay un elemento x 2 E tal que fd.xn ; x/g ! 0, es decir, para todo " > 0 existe un m" 2 N tal que para todo n > m" se verifica que d .xn ; x/ < ". Es fácil comprobar que si hay algún elemento x 2 E que verifica la condición de convergencia anterior dicho elemento x es único. Tal elemento se llama límite de la sucesión fxn g y escribimos lKım fxn g D x o, simplemente, fxn g ! x . n!1
Observa que fxn g ! x quiere decir que para toda bola B .x; "/ existe un m " 2 N tal que para todo n > m" se verifica que x n 2 B.x; "/.
Ejercicios propuestos 25. Prueba que el concepto de sucesión convergente es topológico. Por tanto, distancias equivalentes tienen las mismas sucesiones convergentes. 26. Para cada n 2 N sea f n W Œ0; 1 ! R la función dada por f .x/ D x n. Prueba que la sucesión ff n g converge a la función nula f D 0 en el espacio . C .Œ0; 1/; k k1 / pero no es convergente en . C .Œ0; 1/; k k1 /. 27. Sea fxm g una sucesión de puntos de . Rn ; k k2 /. Representaremos por x m.k/ la coordenada k ésima del vector xm. Prueba que fxmg ! x si, y sólo si, fxm .k/g ! x.k/ para 1 6 k 6 n. Esto es, la convergencia en . Rn ; k k2 / equivale a la convergencia por coordenadas. 28. Sea fxmg una sucesión de puntos de . Rn ; k k2 /. Prueba que fxmg ! x si, y sólo si, se verifican las siguientes dos condiciones: a) xm y ! x y 8y 2 Rn . b) fkxm k2 g ! kxk2 .
˚˝ ˇ ˛ ˝ ˇ ˛
1.18 Proposición. Sean .E; d/ un espacio métrico, A E y x 2 E. Equivalen las siguientes afirmaciones:
a) x es un punto adherente a A. b) x es límite de una sucesión de puntos de A.
La condición de Cauchy para una sucesión de números reales puede formularse para sucesiones de puntos de un espacio métrico. Universidad de Granada Dpto. de Análisis Matemático
Prof. Javier Pérez Cálculo Diferencial en Rn
Topología de un espacio métrico
10
1.19 Definición (Sucesión de Cauchy). Se dice que una sucesión fxn g de puntos de un espacio métrico .E; d/ cumple la condición de Cauchy o que es una sucesión de Cauchy, si para todo " > 0 existe un m " 2 N tal que para todos p > m" , q > m" se verifica que d .xp ; xq / < ".
Se dice que .E; d/ es un espacio métrico completo si toda sucesión de puntos de E que verifica la condición de Cauchy es convergente. Un espacio normado que es completo como espacio métrico se llama un espacio de Banach. El concepto de sucesión de Cauchy no es topológico, es decir, dos distancias equivalentes pueden no tener las mismas sucesiones de Cauchy.
1.20 Definición. Sea .E; d/ un espacio métrico y A un subconjunto no vacío de E . Si el conjunto de números reales:
D fd.x; y/ W x 2 A; y 2 Ag
C
está mayorado, se dice que A está acotado, en cuyo caso el número sup .C / se llama diámetro de A . Se dice que una sucesión fxn g está acotada si el conjunto fxn W n 2 N g está acotado.
Ejercicios propuestos 29. Prueba que toda sucesión convergente es de Cauchy. 30. Prueba que toda sucesión de Cauchy está acotada. 31. Prueba que si una sucesión de Cauchy tiene una sucesión parcial convergente entonces es convergente. 32. Sea .X ; k k/ un espacio normado. Prueba que un conjunto A X está acotado si, y sólo si, existe un número M > 0 tal que kx k 6 M para todo x 2 A. 33. Prueba que . Rn ; k k2 / es un espacio de Banach. 34. Prueba que el espacio . B .A/; k k1 / es completo. 35. Para cada n 2 N sea f n W Œ0; 2 ! R la función dada por: f .x/
D
x n ; 0 6 x 6 1; 1; 1 6 x 6 2.
Prueba que la sucesión f f n g es de Cauchy en el espacio .C .Œ0; 2/; k convergente.
k1 / pero no es
36. Sea I R un intervalo. Prueba que el conjunto C .I / de las funciones reales continuas en I es un conjunto cerrado en el espacio . B .I /; k k1 /.
Universidad de Granada Dpto. de Análisis Matemático
Prof. Javier Pérez Cálculo Diferencial en Rn
Topología de un espacio métrico
11
37. Sea ` 1 el espacio normado de las sucesiones acotadas de números reales con la norma uniforme, es decir, ` 1 D .B .N /; k k1 /. Para cada n 2 N sea ı n W N ! R la sucesión definida por 1; si k D n; ın .k/ D 0; si k ¤ n. Prueba que la sucesión fın g no tiene ninguna sucesión parcial convergente.
1.21 Teorema ( Bolzano – Weierstrass en .Rn ; k k2 /). Toda sucesión acotada de puntos de .Rn ;
k k2/ tiene alguna sucesión parcial convergente en .Rn; k k2/.
1.22 Teorema (Hausdorff ). Todas las normas en un espacio vectorial de dimensión finita son equivalentes.
Como consecuencia de este teorema en Rn hay una única topología que procede de una norma, dicha topología se llama la topología de la norma y es la que consideraremos siempre en Rn . En dicha topología los abiertos son uniones de bolas abiertas para alguna norma. Puesto que el concepto de sucesión convergente es topológico, deducimos, como consecuencia del teorema de Hausdorff, que las sucesiones convergentes en Rn son las mismas para todas las normas.
1.23 Corolario. En todo espacio normado de dimensión finita la convergencia de una sucesión de puntos equivale a la convergencia por coordenadas.
Otra consecuencia del teorema anterior y de la proposición 1.16 es que las sucesiones de Cauchy y las sucesiones acotadas en un espacio normado de dimensión finita son las mismas para todas las normas. Hemos visto en el ejercicio 33 que . Rn ; k k2 / es un espacio de Banach, es decir en dicho espacio las sucesiones convergentes coinciden con las de Cauchy. Deducimos que en Rn con cualquier norma las sucesiones de Cauchy coinciden con las sucesiones convergentes. A partir de aquí es fácil probar los siguientes resultados.
1.24 Teorema (Complitud). Todo espacio normado de dimensión finita es un espacio de Banach.
1.25 Teorema ( Bolzano–Weierstrass ). Toda sucesión acotada de puntos de un espacio normado de dimensión finita tiene alguna sucesión parcial convergente.
1.26 Definición. Un subconjunto K de un espacio métrico se dice que es compacto si toda sucesión de puntos de K tiene alguna sucesión parcial que converge a un punto de K . Puesto que el concepto de sucesión convergente es topológico, la compacidad es una propiedad topológica.
1.27 Proposición. Todo conjunto compacto de un espacio métrico es cerrado y acotado. Como consecuencia de lo visto en el ejercicio 37, se tiene que la bola cerrada B .0; 1/ en el espacio `1 no es un conjunto compacto aunque, evidentemente, dicho conjunto es cerrado y acotado.
1.28 Teorema (Caracterización de la compacidad en dimensión finita ). Los conjuntos com pactos en un espacio normado de dimensión finita son los conjuntos cerrados y acotados.
Universidad de Granada Dpto. de Análisis Matemático
Prof. Javier Pérez Cálculo Diferencial en Rn
Continuidad
12
Ejercicios propuestos 38. Sea K un subconjunto compacto no vacío de Rn . Sea r > 0 un número fijo y definamos A D x2K B.x; r / donde B .x; r / indica la bola cerrada de centro x y radio r para una norma que se supone fijada en Rn . Prueba que A es compacto.
S
1.4. Continuidad 1.29 Definición. Sean .E; d/, .F ; / espacios métricos, A E y f W A ! F . Se dice que f es continua en un punto ˛ 2 A si para todo " > 0 existe ı > 0 tal que para todo x 2 A con d.x; ˛/ < ı se verifica que .f .x/;f .˛// < ". Si f es continua en todos los puntos de un conjunto B A se dice que f es continua en B . Si f es continua en A se dice simplemente que f es continua.
1.30 Proposición (Caracterización secuencial de la continuidad ). Sean .E; d/ , .F ; / espacios métricos, A
E , f W A ! F y ˛ 2 A. Equivalen las siguientes afirmaciones:
a) f es continua en ˛. b) Para toda sucesión xn
f g ! ˛ con xn 2 A se verifica que ff .xn/g ! f.˛/.
Como el concepto de sucesión convergente es topológico, del resultado anterior se sigue que la continuidad es una propiedad topológica.
Ejercicios propuestos 39. Enuncia apropiadamente y prueba los siguientes resultados: a) La composición de funciones continuas es continua. b) Cualquier restricción de una función continua es continua. 40. Prueba que toda función es continua en los puntos aislados de su conjunto de definición. 41. Prueba que en todo espacio normado la norma es una aplicación continua. 42. Prueba que la imagen de un compacto por una función continua es un compacto. 43. Prueba que todo compacto en
R tiene máximo y mínimo.
44. Propiedad de compacidad. Sea f W K ! R una función continua definida en un con junto compacto K de un espacio métrico. Entonces f alcanza en K un máximo y un mínimo absolutos. Universidad de Granada Dpto. de Análisis Matemático
Prof. Javier Pérez Cálculo Diferencial en Rn
Continuidad
13
45. Sean .E; d/ un espacio métrico, f i W E ! E , i D 1; 2; : : : ; n funciones continuas. Consideremos el espacio métrico producto .E n ; / y definamos F W E n ! E n por F .x1 ; x2 ; : : : ; xn /
D .f 1.x1/; f 2.x2/ ; : : : ; fn .xn//
Prueba que F es continua.
46. Sea F W Rn
! Rn una función continua tal que para todos x; y 2 Rn se verifica que kx y k 6 kF.x/ F.y /k. Prueba que la imagen de F es un conjunto cerrado en Rn.
1.31 Proposición (Carácter local de la continuidad ). Sean .E; d/ , .F ; / espacios métricos, A E , f A F y ˛ A. Supongamos que hay alguna bola abierta B.˛; r / tal que la B.˛; r / A es continua en ˛ . Entonces f es continua restricción f jB de f al conjunto B en ˛.
W !
2
D
\
A es un conjunto abierto y f jU es continua, entonces f es continua
En particular, si U en U .
1.32 Proposición ( Caracterización topológica de la continuidad ). Sean .E; d/ y .F ; / es pacios métricos y f E
W ! F . Equivalen las siguientes afirmaciones:
a) f es continua en E. b) La imagen inversa por f de todo abierto de F es un abierto de E. c) La imagen inversa por f de todo cerrado de F es un cerrado de E .
1.33 Corolario. Sea .E; d/ un espacio métrico y f WE ! R una aplicación continua. Entonces
para todo t R se verifica que:
2
f 2 W g f 2 W g b) Los conjuntos fx 2 E W f .x/ 6 t g , fx 2 E W f .x/ > t g y fx 2 E W f .x/ D t g son cerrados a) Los conjuntos x E f .x/ < t y x E f .x/ > t son abiertos en E.
en E .
1.34 Corolario. Sean d y dos distancias en un conjunto no vacío E. Equivalen las afirmaciones siguientes: a) d y son distancias equivalentes. b) Las aplicaciones identidad entre los espacios métricos .E; d/ y .E; / son continuas. c) Los espacios métricos .E; d/ y .E; / tienen las mismas sucesiones convergentes.
Sea .E; d/ un espacio métrico, A un subconjunto no vacío de E y notemos djA la restricción de d a A A; es decir, djA W A A ! RCo es la aplicación definida por: djA .x; y/ D d.x; y/
8.x; y/ 2 A A
Es inmediato comprobar que djA es una distancia en A que se llama distancia inducida en A por la distancia d. Por tanto .A; djA / es un espacio métrico. Se dice que .A; djA / es un subespacio métrico de .E; d/. Universidad de Granada Dpto. de Análisis Matemático
Prof. Javier Pérez Cálculo Diferencial en Rn
Continuidad
14
No suele distinguirse entre la restricción d jA de la distancia d al conjunto A y la distancia d, y ambas se representan con la misma letra d, aunque no hay que olvidar que .A; d/ y .E; d/ son espacios métricos diferentes. Naturalmente, los conceptos que se han definido antes para un espacio métrico .E; d/ conservan su significado cuando nos referimos a un subespacio métrico. Por ejemplo, la expresión “A es un subespacio completo de E ” quiere decir que A, considerado como subespacio métrico de .E; d/, es un espacio métrico completo: toda sucesión de Cauchy de puntos de A converge a un punto de A . Es frecuente, sin embargo, no usar la palabra “subespacio” sino “subconjunto” con el convenio de que todo subconjunto no vacío de un espacio métrico .E; d/ se considera automáticamente como subespacio métrico.
Ejercicios propuestos 47. Da ejemplos de subconjuntos de
R que no sean
completos.
48. Prueba que todo conjunto compacto de un espacio métrico es completo. 49. Prueba que todo conjunto completo de un espacio métrico es cerrado. 50. Prueba que todo conjunto cerrado de un espacio métrico completo es completo. 51. Sea f W R ! R una función continua e inyectiva. Definamos W R R ! RCo por .x; y/
D jf .x/ f .y/j
8.x; y/ 2 R R
Prueba que: a) es una distancia en R equivalente a la usual. b) Si la función f está mayorada o minorada entonces el espacio métrico .R; / no es completo.
52. En el intervalo E D0; 1 se define la distancia: .x; y/
ˇD ˇ ˇ 1
1
x
y
ˇˇ ˇ
8.x; y/ 2 E E
Prueba que .E; / es un espacio métrico completo y que la distancia es equivalente a la distancia usual en E . Como sabemos, todo espacio métrico se considera como espacio topológico con la topología asociada a la distancia, por tanto, todo subconjunto no vacío A de un espacio métrico .E; d/ también se considera como espacio topológico con la topología que tiene como subespacio métrico .A; d/. Los abiertos en esta topología se llaman abiertos relativos y los cerrados se llaman cerrados relativos de A . Teniendo en cuenta que las bolas abiertas, BA .a; r /, de un punto ˛ 2 A en .A; d/ son de la forma: BA .a; r /
D fx 2 A W d .x; a/ < r g D fx 2 E W d .x; a/ < r g \ A D B.a; r / \ A
es fácil describir los abiertos y los cerrados relativos de A . Universidad de Granada Dpto. de Análisis Matemático
Prof. Javier Pérez Cálculo Diferencial en Rn
Continuidad
15
1.35 Proposición. Sea A un subconjunto no vacío de un espacio métrico E . Se verifica que: a) Los abiertos relativos de A son las intersecciones de los abiertos de E con A. b) Los cerrados relativos de A son las intersecciones de los cerrados de E con A.
Consideremos ahora una función f definida en un subconjunto no vacío A de un espacio métrico .E; d/ con valores en otro espacio métrico .F ; /. La continuidad de f en un punto ˛ 2 A quiere decir que para todo " > 0 existe un ı > 0 tal que para todo x 2 A con d.x; ˛/ < ı se verifica que .f.x/; f .˛// < ". Observa que aquí solamente intervienen el conjunto A , la distancia d trabajando en A y el espacio .F ; /. Es decir, para definir la continuidad de f en un punto ˛ 2 A podemos olvidarnos del espacio .E; d/ y quedarnos con el subespacio .A; d/. En particular, la continuidad de f en A quiere decir que la aplicación f del espacio métrico .A; d/ en el espacio métrico .F ; / es continua. Como consecuencia de las proposiciones 1.32 y 1.35, deducimos el siguiente resultado.
1.36 Proposición. Sean E y F espacios métricos, Ø ¤ A E y f W A ! F . Equivalen las siguientes afirmaciones:
a) f es continua en A. b) La imagen inversa por f de todo abierto de F es un abierto relativo de A. c) La imagen inversa por f de todo cerrado de F es un cerrado relativo de A.
Continuidad uniforme Piensa un par de minutos antes de responder a la siguiente pregunta. Supongamos que f es una función real continua en un intervalo I . ¿Es cierto que si tomamos valores x ; y 2 I muy próximos entre sí los correspondientes valores de la función f .x/ ; f .y/ también están muy próximos entre sí? Si tu respuesta ha sido afirmativa, como suele ser, te equivocas. Considera la función continua f W0; 1 ! R dada por f .x/ D 1=x . Los puntos 10 10 y 1020 están muy próximos entre sí: 10 10 1020 < 1010 , pero f .1010 / D 1010 y f .1020 / D 1020 están muy distantes entre sí. No hay nada extraño en este comportamiento. A cualquier función continua cuya gráfica tenga una asíntota vertical le pasa lo mismo: hay puntos muy próximos entre sí en los que la función toma valores muy distantes entre sí. Pero también hay funciones continuas y acotadas que se comportan de forma parecida. Considera la función continua gW0; 1 ! R dada por g.x/Dsen.1=x/. Es una función acotada: el mayor valor que toma es 1 y el menor valor que toma es 1, de hecho se tiene que g.0; 1/ D 1 Œ1; 1. Sea n un número natural. Los puntos xn D 2nC=2 e yn D 2n 1=2 están, para n suficientemente grande, muy próximos entre sí; de hecho f xn yn g ! 0. Los valores que toma en ellos la función g.xn / D 1 y g.yn / D 1 distan entre sí 2 unidades (que es la máxima distancia que puede haber entre valores tomados por esta función). Si lo piensas un poco, te darás cuenta de que en ambos ejemplos este comportamiento se debe a que las funciones f y g “oscilan mucho” en intervalos arbitrariamente pequeños. Conviene precisar la idea de “oscilación en un intervalo”. Se define la oscilación de una función real f en un intervalo J contenido en el dominio de Universidad de Granada Dpto. de Análisis Matemático
Prof. Javier Pérez Cálculo Diferencial en Rn
Continuidad
16
definición de f como: !.f; J /
D
sup f .J / Kınf f .J /; si f .J / está acotado I C1; si f .J / no está acotado:
En otros términos: la oscilación de f en J es la longitud del intervalo más pequeño que contiene a f .J /. Para la función f .x/ D 1=x se tiene que ! f; Œ1=2n; 1=n D n y ! f; 0; 1=n DC1. Para la función g.x/ D sen.1=x/ tenemos que ! g; Œ1=.2n C =2/; 1=.2n =2/ D 2. Estas funciones tienen una oscilación “grande” en intervalos arbitrariamente pequeños. En algunas circunstancias interesa poder controlar el tamaño de la oscilación de una función de manera que dicha oscilación sea menor que una cierta cantidad fijada, " > 0, en cualquier intervalo de longitud menor que un cierto número ı > 0 . Las funciones para las que esto puede hacerse cualquiera sea " > 0 , se llaman uniformemente continuas.
Se dice que una función f W I ! R es uniformemente continua en un intervalo I si para todo " > 0 es posible encontrar un ı > 0 de manera que siempre que J sea un intervalo contenido en I de longitud menor que ı , se verifica que la oscilación de f en J es menor o igual que " . Teniendo en cuenta que !.f; J /6"”jf .x/ f .y/j6" para todos x; y 2 J , la definición dada puede expresarse de forma equivalente como sigue. Una función f es uniformemente continua en un intervalo I si para todo " > 0 es posible encontrar un ı > 0 de manera que siempre que x , y sean puntos de I con jx y j 6 ı , se verifica que jf .x/ f .y j 6 ". Simbólicamente: C
C
8" 2 R 9 ı 2 R
W jxx;yy2j 6I ı �jf .x/ f .y/j 6 "
Esta definición puede generalizarse fácilmente para funciones entre espacios métricos. Llegamos así a la siguiente definición.
1.37 Definición. Sean .E; d/, .F ; / espacios métricos, A un subconjunto no vacío de E y f W A ! F una función definida en A . Se dice que f es uniformemente continua en A si para todo " > 0 es posible encontrar un ı > 0 de manera que siempre que x , y sean puntos de A con d.x; y/ 6 ı , se verifica que .f .x/ ; f .y// 6 ". Simbólicamente:
8" 2 RC 9 ı 2 RC W
d.x; y/ 6 ı x; y A
2
�
.f.x/ ; f .y// 6 "
(1.2)
1.38 Observaciones. El concepto de “continuidad uniforme” es un concepto global: depende del comportamiento de la función en todo un conjunto. No tiene sentido decir que una función es uniformemente continua en un punto: la continuidad uniforme no es un concepto local.
Es muy interesante comparar las definiciones de continuidad puntual 1.29 y de continuidad uniforme 1.37. Resulta evidente que la continuidad uniforme en A implica la continuidad en todo punto de A: toda función uniformemente continua en un conjunto es continua en dicho conjunto.
Universidad de Granada Dpto. de Análisis Matemático
Prof. Javier Pérez Cálculo Diferencial en Rn
Continuidad
17
En general, no es cierto que una función continua en un conjunto A sea uniformemente continua en A como lo prueban los ejemplos dados al principio. Pero hay una situación particular en la que dicha afirmación sí es cierta. Este es el contenido del siguiente teorema. Se trata de un resultado importante en el que pueden destacarse aportaciones de varios matemáticos. Dirichlet ya lo incluyó en sus lecciones de 1862 y en 1872 Heine dio una primera demostración del mismo. Posteriormente Weierstrass, Borel y Lebesgue generalizaron el resultado inicial.
1.39 Teorema (Heine). Toda función continua en un conjunto compacto K de un espacio métrico con valores en un espacio métrico es uniformemente continua en K .
Por los resultados vistos hasta ahora, queda claro que la compacidad es una propiedad que tiene importantes consecuencias de muy variado tipo. El siguiente resultado nos dice que dicha propiedad debe usarse con cuidado en espacios normados de dimensión infinita. Se dice que un espacio normado es localmente compacto si las bolas cerradas son compactos. Teniendo en cuenta que en un espacio normado dos bolas cerradas cualesquiera de radios positivos son homeomorfas, resulta que un espacio normado es localmente compacto si, y sólo si, hay alguna bola cerrada de radio positivo que sea compacta.
1.40 Teorema (Riesz, 1918). Un espacio normado es de dimensión finita si, y sólo si, es localmente compacto.
Ejercicios propuestos 53. Justifica con un ejemplo que la imagen de una sucesión de Cauchy por una función continua puede no ser una sucesión de Cauchy. 54. Prueba que la imagen de una sucesión de Cauchy por una función uniformemente continua es una sucesión de Cauchy. 55. Prueba que un espacio normado es localmente compacto si, y sólo si, hay alguna bola cerrada de radio positivo que sea compacta. 56. Sea `1 el espacio normado de las sucesiones acotadas de números reales con la norma uniforme, es decir, `1 D .B .N /; k k1 /. Supongamos que A B.0; 1/ es tal que B.0; 1/ B.x; 1=2/. Prueba que A es un conjunto infinito.
[
x2A
Algunas clases de funciones continuas 1.41 Definición. Sean .E; d/, .F ; / espacios métricos, A un subconjunto no vacío de E y f W A ! F . Se dice que f es lipchiciana si existe M > 0 tal que .f.x/ ; f .y// 6 M d.x; y/
8x; y 2 A:
Si M < 1 se dice que f es una aplicación contractiva. Cualquier M > 0 que satisfaga la desigualdad anterior se llama una constante de Lipschitz para f . Universidad de Granada Dpto. de Análisis Matemático
Prof. Javier Pérez Cálculo Diferencial en Rn
Continuidad
18
Es claro que las funciones lipchicianas son uniformemente continuas.
1.42 Proposición. Sean .E; k k/ , .F ; jjj jjj/ espacios normados y T W E ! F una aplicación lineal. Equivalen las siguientes afirmaciones: a) T es continua. b) T es continua en 0. c) Existe M > 0 tal que
jjjT .x/jjj 6 M kxk para todo x 2 E.
d) T es lipchiciana.
Ejercicios propuestos 57. Prueba que la aplicación x 7! kx k es lipchiciana con constante de Lipschitz 1. 58. Sea .E; d/ un espacio métrico y Ø ¤ A E . Sea dist. ; A/ W E definida por: 8x 2 E dist.x; A/ D Kınf fd.x; a/ W a 2 Ag .
! R la aplicación
Prueba que dicha aplicación es lipchiciana con constante de Lipschitz 1. Prueba también que x 2 A si, y sólo si, dist .x; A/ D 0.
59. Sea A Rn un conjunto abierto no vacío y distinto de
D 2
K p
x A
Rn . Para cada p n
1
W kxk 6 p; dist.x; R nA/ > n
2 N se define:
Prueba que: ı
a) K p es compacto y K p K pC1. 1
b)
[
pD1
K p
D A.
c) Si K A es compacto, existe un p 2 N tal que K K p .
60. Sean F Rn un conjunto cerrado y K Rn un conjunto compacto no vacíos y disjuntos F \ K D Ø. Sea k k una norma en Rn . Prueba que hay puntos a 2 F y b 2 K tales que ka bk D Kınf fkx y k W x 2 F ; y 2 K g Sugerencia. Considera la función x 7! dist.x; F /.
61. Sean F un subespacio finito dimensional de un espacio normado .E; k k/ y a 2 E . Prueba que existe b 2 F tal que: ka bk D dist.a; F /
Universidad de Granada Dpto. de Análisis Matemático
Prof. Javier Pérez Cálculo Diferencial en Rn
Continuidad
19
Deduce que dados n 2 N y una función f 2 C .Œ0; 1/, existe una función polinómica h de grado menor o igual que n tal que: 1
1
jf .x/ h.x/j dx 6 jf .x/ p.x/j dx 0
0
cualquiera sea la función polinómica p de grado menor o igual que n .
1.43 Proposición. Toda aplicación lineal de un espacio normado de dimensión finita en otro espacio normado es continua.
Reciben el nombre de campos escalares las funciones definidas en subconjuntos de Rn que toman valores en R. Un campo escalar es, por tanto, una función real que depende de n variables. Un campo escalar de una variables es, simplemente, una función real de variable real; un campo escalar de dos variables es una función definida en un subconjunto del plano que toma valores reales; un campo escalar de tres variables es una función definida en un subconjunto del espacio que toma valores reales. Los campos escalares de una o dos variables se pueden visualizar por medio de sus representaciones gráficas que son, respectivamente, curvas en el plano o superficies en el espacio. No es posible visualizar campos escalares de tres o más variables porque sus gráficas están en espacios de dimensión mayor o igual que 4. Naturalmente, los campos escalares se pueden sumar, multiplicar y dividir (donde no se anule el denominador) al igual que lo hacemos con las funciones reales. Es fácil probar, usando la caracterización secuencial de la continuidad, que dichas operaciones conservan la continuidad. Ejemplos de campos escalares continuos lo proporcionan las proyecciones sobre los ejes coordenados. Representaremos por k la aplicación k W Rn ! R que a cada vector x D .x1 ; x2 ; : : : ; xn / 2 Rn hace corresponder su coordenada k -ésima en la base canónica: k .x/
D k .x1; x2 ; : : : ; xn/ D xk
Dicha aplicación se llama proyección coordenada k -ésima. Es una aplicación lineal. Los campos escalares más sencillos son las funciones polinómicas de varias variables. Dichas funciones se obtienen como sumas de productos de proyecciones sobre los ejes coordenados y son, por tanto, continuas. Para n D 3 las proyecciones sobre los ejes coordenados son 1 .x; y; z/
D x;
2 .x; y; z/
D y;
3 .x; y; z/
Dz Un producto de estas funciones es una función de la forma f .x; y; z/ D x my p z q donde m; p; q son números naturales o nulos. Las funciones polinómicas en tres variables son combinaciones lineales de este tipo de funciones. Universidad de Granada Dpto. de Análisis Matemático
Prof. Javier Pérez Cálculo Diferencial en Rn
Límite funcional
20
Las funciones racionales de n variables son las funciones de la forma R.x1 ; x2 ; : : : ; xn /
P .x1 ; x2 ; : : : ; xn / D Q.x ;x ;:::;x / 1
n
2
Donde P .x1; x2 ; : : : ; xn / y Q.x1 ; x2 ; : : : ; xn / son funciones polinómicas de n variables. El dominio natural de definición de una función racional es el conjunto de puntos donde no se anula el denominador D fx 2 Rn W Q.x/ ¤ 0g. Las funciones racionales son continuas en su dominio natural de definición.
Componiendo funciones continuas reales de una variable con funciones polinómicas y racionales en varias variables obtenemos muchísimos ejemplos de campos escalares continuos. Aquí tienes unos pocos. f .x; y/
D sen.xy/;
f .x; y/
2
D log.1 C x C y
2
/;
C xy 2 C xz 2 f .x; y; z/ D 2 C arctg.xyz/ 1
1.5. Límite funcional El concepto de límite para funciones reales de una variable real, que ya conoces de primer curso, se generaliza con facilidad para funciones entre espacios métricos.
1.44 Definición. Sean .E; d/, .F ; / espacios métricos, A E , f W A ! F y ˛ 2 A 0 . Se dice que f tiene límite en ˛ si hay un elemento ` 2 F con la propiedad de que para todo " > 0 existe un ı > 0 tal que para todo x 2 A con 0 < d.x; ˛/ < ı se verifica que .f.x/;`/ < ". Es fácil comprobar que solamente puede haber un elemento ` 2 F que cumpla la condición anterior. Se dice que ` es el límite de f en el punto ˛ y escribimos lKım f .x/ D `. x!˛
Con el abuso usual de notación, podemos escribir: lKım f .x/ D `
” 82
" RC
9ı 2 RC W
0 < d.x; ˛/ < ı x A
�
.f.x/;`/ < " (1.3) 2 Una notación mejor para el límite sería x!˛ l Kım f .x/ D ` para indicar que solamente se consix2A deran valores de la variable x 2 A. Esta precisión es importante porque la existencia y el valor x!˛
de un límite pueden depender del conjunto A donde se considera definida la función. No obstante, no usaremos esta notación porque es algo incómoda pero no debes olvidar que la variable x siempre toma valores en el conjunto donde se supone que está definida la función.
Observa que en la definición se supone que ˛ 2 A 0 es un punto de acumulación de A . Esto asegura que para todo ı > 0 , hay puntos x 2 A que cumplen la condición 0 < d.x; ˛/ < ı. Por 2 A, en cuyo caso f puede no estar definida en ˛ . otra parte, puede ocurrir que ˛ 6 La condición 0 < d .x; ˛ / < ı es una forma de escribir x ¤ ˛ y d.x; ˛ / < ". Es decir, en la condición de límite (1.3) no interviene para nada, en el caso de que f esté definida en ˛ , el valor que f pueda tener en ˛ . Podemos escribir (1.3) usando bolas abiertas: lKım f .x/ D `
x!˛
” 8
2 \ A �f .x/ 2 B.`; "/ B .`;"/ 9Bd .˛;ı/ W ¤
Universidad de Granada Dpto. de Análisis Matemático
x Bd .˛;ı/ x ˛
Prof. Javier Pérez Cálculo Diferencial en Rn
Continuidad y límites de campos escalares y vectoriales
21
Puedes comprobar que si en lo anterior sustituyes “bolas abiertas centradas en ˛ o en ` ” por “abiertos que contengan a ˛ o a `” la condición que se obtiene es equivalente. Es decir, el concepto de límite de una función en un punto es un concepto topológico. La siguiente caracterización del límite funcional mediante sucesiones es una generalización directa de la que se estudia en primer curso.
1.45 Proposición (Caracterización secuencial del límite funcional) . Sean .E; d/ , .F ; / espa-
E , f W A ! F y ˛ 2 A 0. Equivalen las siguientes afirmaciones: a) lKım f .x/ D `. x!˛ b) Para toda sucesión fxn g de puntos de A distintos de ˛ , x n 2 Anf˛ g , que converge a ˛ , fxng !d ˛ , se verifica que ff .xn/g ! `. cios métricos, A
La relación entre límite funcional y continuidad es la que ya conoces de primer curso.
1.46 Proposición. Sean .E; d/ , .F ; / espacios métricos, A E , f W A ! F y ˛ 2 A \ A 0 .
Equivalen las siguientes afirmaciones: a) f es continua en ˛. b) lım f .x/
K
x!˛
D f.˛/.
1.47 Proposición ( Carácter local del límite funcional ). Sean .E; d/ , .F ; / espacios métricos, A E , f A F y ˛ A 0 . Supongamos que hay alguna bola abierta B.˛; r / tal que la restricción f jB de f al conjunto B B.˛; r / A tiene límite en ˛ igual a `. Entonces f también tiene límite en ˛ igual a `.
W !
2
D
\
1.6. Continuidad y límites de campos escalares y vectoriales Las funciones que más nos interesan en este curso son los campos escalares y las funciones vectoriales. Un campo escalar de n variables es una función definida en un conjunto no vacío A Rn que toma valores reales. Una función vectorial de n variables es una función definida en un conjunto no vacío A Rn que toma valores vectoriales en un espacio Rm con m > 2. Por comodidad en la expresión, llamaremos a veces campos vectoriales a las funciones vectoriales aunque dicho nombre suele reservarse para aquellas funciones vectoriales definidas en un subconjunto de un espacio vectorial y que toman valores en dicho espacio vectorial.
1.48 Definición (Componentes de una función vectorial) . Sea A Rn y F W A ! Rm una función vectorial. Para cada x 2 A se tiene que F .x/ es un vector de Rm . Para 1 6 k 6 m definamos f k D k ı F W A ! R donde k es la proyección coordenada k -ésima. Las funciones f k así definidas son campos escalares que se llaman componentes de F . Para cada x 2 A se tiene que m
F.x/ D .f 1.x/; f 2 .x/ ; : : : ; fm .x// m
donde fuk W 1 6 k 6 mg es la base canónica de R .
X D
f k .x/uk
kD1
Escribiremos F D .f 1 ; f 2 ; : : : ; fm / para indicar que F es un campo vectorial cuyos campos escalares componentes son los f k ( 1 6 k 6 m). Universidad de Granada Dpto. de Análisis Matemático
Prof. Javier Pérez Cálculo Diferencial en Rn
Continuidad y límites de campos escalares y vectoriales
22
Ejercicios propuestos
62. Sea F D .f 1; f 2 ; : : : ; fm / W A ! Rm un campo vectorial de n variables. Prueba que:
1. F es continuo en un punto ˛ 2 A si, y sólo si, todos los campos escalares f k son continuos en ˛. 2. lKım F.x/ D ˇ ” lKım f k .x/ D ˇk .1 6 k 6 m/. x!˛
x!˛
Teniendo en cuenta este ejercicio, para estudiar límites o continuidad de campos vectoriales es suficiente hacerlo para campos escalares. Campos escalares definidos en un mismo conjunto se pueden sumar, multiplicar y dividir. Así mismo, campos vectoriales con igual número de funciones componentes pueden sumarse, y también podemos multiplicar un campo vectorial por un campo escalar. Todo ello en el supuesto de que los campos considerados estén definidos en un mismo conjunto. También pueden componerse campos vectoriales entre sí o con campos escalares. Todas estas operaciones se comportan respecto a la continuidad y al límite funcional en la forma que cabe esperar. Tiene interés definir el concepto de campo escalar divergente en un punto que generaliza lo que ya conoces de primer curso. 1.49 Definición. Sea f W A ! R donde A Rn y ˛ 2 A 0 . Se dice que f es positivamente divergente en ˛ si para todo M > 0 existe ı > 0 tal que para todo x 2 A con 0 < kx ˛k < ı se verifica que jf .x/j > M . Simbólicamente escribimos lKım f .x/ D C1. x!˛
Se dice que f es negativamente divergente en ˛ si f es positivamente divergente en ˛ , en cuyo caso escribimos lKım f .x/ D 1. x!˛
Se dice que f es divergente en ˛ si jf j es positivamente divergente en ˛, en cuyo caso escribimos lKım f .x/ D 1. x!˛
En la definición anterior k k es cualquier norma en Rn , pues sabemos que todas ellas son topológicamente equivalentes por lo que podemos usar la que queramos para estudiar propiedades topológicas como son límites o continuidad. En lo que sigue vamos a considerar el caso más sencillo de un campo escalar de dos variables: f W A ! R donde Ø ¤ A R2 . Supongamos que ˛ 2 A 0 y queremos estudiar el límite lKım f .x/. En R2 tenemos que x D .x; y/, ˛ D .a; b/, por lo que el límite anterior se escribe x!˛ f .x; y/ (y no debes olvidar el paréntesis que indica vector: usualmente en la forma lKım .x;y/!.a;b/
.x; y/ y .a; b/, porque no tiene sentido escribir
lKım f .x; y/). Siempre que te encuentres en
x;y!a;b
esta situación el resultado de intentar evaluar la función f en .a; b/ va a ser una indeterminación, es decir, una expresión del tipo 0=0,1=1 o 11 y otras que ya conoces de primer curso. Veamos un ejemplo. sen.x 2 C y 2 / 1.50 Ejemplo. Sea f .x; y/ D . Esta función es un campo escalar definido (y 2 2 2
x y .0; 0/ . Si tratas de evaluar f .x; y/ en .0; 0/ obtienes la indeter-
C
g continuo ¿por qué?) en R nf minación 0=0. En este caso, como en otros parecidos, hay una forma sencilla de proceder. Basta Universidad de Granada Dpto. de Análisis Matemático
Prof. Javier Pérez Cálculo Diferencial en Rn
Continuidad y límites de campos escalares y vectoriales
23
sen.t / donde t D x 2 C y 2 . Dicho de otra t sen.t / forma, el campo escalar f es la composición de la función g.t / D con la función polinódarse cuenta de que nuestra función es de la forma
t
mica h.x; y/ D x 2 C y 2 . Llegado aquí ya debes intuir lo que vale el límite
l Kım
.x;y/!.0;0/
f .x; y/:
cuando .x; y/ está muy próximo a .0; 0/ el número t D x 2 C y 2 está muy próximo también a sen.t / 0 por lo que está muy próximo a 1. t
Justifiquemos que
lKım
D 1. Para ello definamos g.0/ D 1 , con lo que la sen t D 1). Ahora, función g resulta continua en 0 (porque ya conoces el famoso límite l Kım t !0 t está claro que definiendo f .0; 0/ D 1 la función f así extendida es continua en .0; 0/ porque es composición de funciones continuas: la función h.x; y/ D x 2 C y 2 que es continua y h.0; 0/ D 0 f .x; y/ D 1. y la función g que es continua y g.0/ D 1. Luego lKım .x;y/!.0;0/ .x;y/!.0;0/
f .x; y/
Otra forma de proceder, más cómoda en este caso, es trabajar con sucesiones. Para ello sea f.xn; yn /g ! .0; 0/ con .xn; yn/ ¤ .0; 0/. Pongamos zn D xn2 C yn2. Claramente, fzng ! 0 y zn ¤ 0 por lo que: sen .xn2 C yn2 / sen .zn / h.xn ; yn /
D
C yn2 D
xn2
zn
! 1:
Como esto es cierto para toda sucesión en las condiciones anteriores, deducimos, por la propof .x; y/ D 1. sición 1.45, que lKım .x;y/!.0;0/
Ejercicios propuestos 63. Estudia la continuidad de los siguientes campos escalares y el límite en el punto ˛ indicado en cada caso: 1 a) f .x; y/ D .1 C xy/ xy definido en RC RC , ˛ D .0; 0/. log .1 C x 2 C y 2 / b) f .x; y/ D definido en R2 nf.0; 0/g, ˛ D .0; 0/. 2 2 x
2
c) f .x; y/ D
e x ey x2
2
Cy 1 definido en R2 nf.0; 0/g, ˛ D .0; 0/.
C y2 cos .x 2 C y 2 / 1 d) f .x; y/ D definido en R2 nf.0; 0/g, ˛ D .0; 0/. .x 2 C y 2 /2 sen x arctg y e) f .x; y/ D definido en A D .x; y/ 2 R2 W xy ¤ 0 , ˛ D .0; 0/. xy x2 C y2 f) f .x; y/ D definido en R2 nf.0; 0/g, ˛ D .0; 0/. x 2 C y 2 C 1 1 1 cos.x C y/ g) f .x; y/ D 2 definido en RC RC , ˛ D .0; 0/. 2 x C y C 2xy
˚
p
Universidad de Granada Dpto. de Análisis Matemático
Prof. Javier Pérez Cálculo Diferencial en Rn
Continuidad y límites de campos escalares y vectoriales
24
Para no tener que indicar en cada caso el conjunto donde está definido un campo escalar, te recuerdo que si un campo escalar viene dado por medio de funciones elementales su dominio natural de definición es el conjunto más grande en el que las operaciones que definen a dicho campo tienen sentido para valores reales de las variables.
Límites direccionales y límite de un campo escalar a lo largo de una curva En la recta real solamente podemos aproximarnos a un punto ˛ 2 R por dos direcciones: por valores mayores y por valores menores que ˛ . Como sabes esto da lugar a los límites laterales en ˛ . Muy diferente situación tenemos si ˛ 2 R2 pues podemos acercarnos a ˛ por infinitas direcciones, a saber, todas las de la forma ˛ C t u donde u es un vector unitario en R2 y t 2 R. Naturalmente, si existe el límite xlK! ım˛ f .x/ D ` entonces para todo conjunto B A con ˛ 2 B 0 x2A
también se verifica que xlK! ım˛ f .x/ D `. En particular, si B D f˛ C t u W t 2 Rg \ A, entonces se x2B
tiene que:
lKım f .x/ D lKım f .˛ C t u/
x!˛ x2B
t !0
Este límite se llama límite de f en ˛ según la dirección dada por el vector u . Según acabamos de decir, si existe xlK! ım˛ f .x/ D ` entonces también debe ser lKım f .˛ C t u/ D ` para todo vector unitario u .
t !0
x2A
Podemos generalizar considerando en vez de una recta que pasa por ˛ otro tipo de curvas que pasan por ˛, por ejemplo, parábolas γ.t / D .t ; t 2 / o γ.t / D .t 2 ; t /. En general, si γ.t / D .x.t /; y.t // es una curva que pasa por ˛ , γ.0/ D ˛, el límite lKım f .γ.t // se llama límite de f t !0 γ en ˛ a lo largo de la curva . Si existe el límite de f en ˛ entonces existe el límite de f a lo largo de toda curva que pasa por ˛ y todos ellos coinciden con el valor del límite. Estos resultados pueden ser útiles para probar que un cierto límite no existe: para ello basta encontrar que no existe el límite a lo largo de una determinada curva o dirección, o bien probar que los límites a lo largo de dos curvas o direcciones distintas son diferentes. Puesto que un vector unitario en R2 con la norma euclídea es de la forma u D .cos #; sen #/, y ˛ D .a; b/, suelen escribirse los límites direccionales en la forma: lKım f .˛ C t u/ D lKım f .a C t cos #; b C t sen # /
t !0
t !0
Es frecuente escribir una recta que pasa por .a; b/ en la forma γ.t / D .a C t ; b C t / donde 2 R es un parámetro fijo – la pendiente de la recta – en cuyo caso los límites direccionales en .a; b/ vienen dados por: lKım f .a C t ; b C t / t !0
Cuando .a; b/ D .0; 0/ las rectas por el origen suelen representarse en la forma γ .t / D .t ; t / o γ.t / D .t ; t / donde es un parámetro fijo. Con frecuencia es más cómodo estudiar los límites en .0; 0/. Esto puede hacerse siempre mediante una traslación pues: lKım
.x;y/!.a;b/
Universidad de Granada Dpto. de Análisis Matemático
f .x; y/
f .x C a; y C b/ D .x;y/!.0;0/ lKım
Prof. Javier Pérez Cálculo Diferencial en Rn
Continuidad y límites de campos escalares y vectoriales
25
Hay que observar que la existencia de todos los límites direccionales en un punto siendo, además, todos ellos iguales, no garantiza la existencia del límite en dicho punto. y
1.51 Ejemplo. Sea f .x; y/ D sen .x 2 C y 2 / definido en A D .x; y/ 2 R2 ; x ¤ 0 . Estux diemos el límite de f en ˛ D .0; 0/.
˚
Consideremos rectas por el origen γ.t / D .t ; t /. Tenemos f .γ.t // D sen..1 C 2 /t 2 /, por lo que lKım f .t ; t / D 0. Por tanto, los límites según todas las direcciones que pasan por t !0 el origen (excepto la dirección del eje de ordenadas γ.t / D .0; t / en la que no está definida f ) existen y son todos iguales a 0. Observa que si γ.t / D .x.t /; y.t // es una curva con γ.0; 0/ D 0, se tiene que: f .γ.t //
y.t / D x.t / sen.x.t /2 C y.t /2 /
Como las curvas son funciones continuas l Kım sen.x.t /2 C y.t /2 / D sen.0/ D 0. Por lo que para t !0
toda curva en estas condiciones que verifique que el cociente lKım f .γ.t // D 0.
y.t /
x.t /
esté acotado se tendrá que
t !0
A la vista de lo anterior, podemos buscar una curva γ.t / D .x.t /; y.t // con γ.0; 0/ D 0 tal y.t /
D 1 y que sen.x.t /2 C y.t /2 / tienda a 0 más lentamente. Como sabemos que t !0 x.t / sen.t / t para t ! 0, la idea es tomar una curva .t / D .x.t /; y.t // que verifique: y.t / 1 jt j 2 D D x.t / t x.t / t 2 1 C t 4 � jt j y.t / D x.t /2 C y.t /2 D t 1 C t 4 que lKım
γ
9> = >>;
8ˆˆ < ˆˆ:
s s
Para dicha curva (que pasa por .0; 0/) tenemos que: f .γ.t //
D f .x.t /; y.t // D t 12 sen .jt j/ D j1t j senj.t jjt j/ � lt !0 Kım f . .t // D C1 γ
El campo escalar f diverge positivamente a lo largo de dicha curva, por lo que dicho campo no tiene límite en .0; 0/.
Límites iterados Para estudiar la existencia del límite de un campo escalar de dos variables suelen usarse los límites iterados. Consideremos un campo escalar f W A ! R donde A es de la forma A D I J nf.0; 0/g siendo I , J intervalos abiertos en R que contienen al 0. Esta situación no es restrictiva dado el carácter local del límite. Supongamos que para todo x 2 I nf0g existe el límite lKım f .x; y/. Tal límite dependerá en general del valor fijado de x 2 I nf0g. Definimos: y!0
W nf0g ! R
G I
Universidad de Granada Dpto. de Análisis Matemático
G.x/
D y!0 lKım f .x; y/
(1.4)
Prof. Javier Pérez Cálculo Diferencial en Rn
Continuidad y límites de campos escalares y vectoriales
26
Si la función G tiene límite en 0 dicho límite se llama un límite iterado de f en .0; 0/ y suele representarse en la forma: lKım G.x/ D lKım
x!0
x!0
K
lım f .x; y/
y!0
(1.5)
Análogamente, si para todo y 2 J nf0g existe el límite lKım f .x; y/, podemos definir la función: x!0
W nf0g ! R
H J
H .y/
D x!0 lKım f .x; y/
(1.6)
Si la función H tiene límite en 0 dicho límite se llama un límite iterado de f en .0; 0/ y suele representarse en la forma: lKım H .y/ D lKım
y!0
y!0
K
lım f .x; y/
x!0
(1.7)
Los límites iterados pueden no existir incluso aunque exista el límite de f en .0; 0/.
1.52 Proposición. En la situación que estamos considerando, supongamos que se verifica alguna de las siguientes afirmaciones: Los dos límites iterados (1.5) y (1.7) existen y son distintos. La función G definida en (1.4) existe pero no existe el límite de G en 0. La función H definida en (1.6) existe pero no existe el límite de H en 0. Entonces no existe el límite de f en .0; 0/. Si existe solamente uno de los límites iterados o bien ambos y coinciden, entonces el límite de f en .0; 0/ puede no existir, pero si existe es igual al valor de dicho límite iterado.
Observa que los límites iterados son límites de funciones de una variable por lo que suelen ser fáciles de calcular.
Ejercicios propuestos 64. Estudia la existencia del límite en .0; 0/ para cada uno de los siguientes campos escalares definidos en R2 nf.0; 0/g. a) f .x; y/ D d) f .x; y/ D
x2
y
p C x2
y2
xy
x2
C y2
y2 b) f .x; y/ D 2 x C y2 .x C y/2 e) f .x; y/ D 2 x C y2
c) f .x; y/ D
x 2y
x4
C y2 x sen y C y sen x f) f .x; y/ D x2 C y2
En los ejemplos de campos escalares de dos variables que hemos considerado aparece con frecuencia la expresión x 2 C y 2 . En los campos escalares de tres variables aparece con frecuencia la expresión x 2 C y 2 C z 2 . Esto no es un capricho matemático sino que procede de la Física: Universidad de Granada Dpto. de Análisis Matemático
Prof. Javier Pérez Cálculo Diferencial en Rn
Continuidad y límites de campos escalares y vectoriales
27
los potenciales newtonianos y los potenciales eléctricos varían en razón inversa a la distancia euclídea y las fuerzas de atracción gravitatoria y eléctrica varían en razón inversa al cuadrado de la distancia euclídea. La expresión x 2 C y 2 se simplifica usando coordenadas polares, es decir, escribiendo .x; y/ D . cos #; sen #/ con lo cual x 2 C y 2 D 2.
1.53 Proposición. Sea f W A ! R donde A D B2 ..0; 0/; r /nf.0; 0/g y B2 ..0; 0/; r / es la bola euclídea abierta en R2 centrada en el origen de radio r . Definamos la función:
W
! R
g 0; r Œ Œ0; 2
dada por g.;#/
D f . cos #; sen #/
Entonces equivalen las siguientes afirmaciones: a)
lKım
.x;y/!.0;0/
f .x; y/
D `.
20; ıŒ se verifica que sup fjf . cos #; sen # / `j W 0 6 # 6 2 g 6 ":
b) Para todo " > 0 existe ı > 0 tal que para todo
En el estudio de límites conviene tener presente la desigualdad jxy j 6 12 .x 2 C y 2 / que se deduce de 0 6 .jx j jy j/2 .
Ejercicios propuestos 65. Estudia la existencia del límite en .0; 0/, y calcula su valor cuando exista, para cada uno de los siguientes campos escalares definidos en R2 nf.0; 0/g. x 2y x3 y3 log .1 C x 2 C y 2 / f .x; y/ f .x; y/ D D a) f .x; y/ D 2 b) c) 2 2 2
Cy sen y log.1 C x 2/ d) f .x; y/ D x2 C y2 x
x
Cy
e) f .x; y/ D
p
jx j C jy j
x sen y
y sen x .x 2 C y 2 / 3 2
66. Estudia la existencia del límite en .0; 0/, y calcula su valor cuando exista, para cada uno de los siguientes campos escalares definidos en RC RC . x ey Cy ex 1 1 sen x C sen y a) f .x; y/D b) f .x; y/D c) f .x; y/D 2 C 2 sen.xy/ x
Cy
x
Cy
x
y
Conjuntos conexos Queremos definir una clase de conjuntos en los espacios métricos que desempeñe un papel análogo al de los intervalos en la recta real. Los intervalos pueden caracterizarse fácilmente con el teorema de Bolzano.
1.54 Proposición. Sea Ø ¤ C R. Equivalen las afirmaciones: a) C es un intervalo. b) Para toda función continua f C
W ! R se verifica que f .C / es un intervalo.
Universidad de Granada Dpto. de Análisis Matemático
Prof. Javier Pérez Cálculo Diferencial en Rn
Continuidad y límites de campos escalares y vectoriales
28
La implicación a/ ) b/ es el famoso teorema de Bolzano. La otra implicación es un sencillo ejercicio que debes hacer.
1.55 Definición. Sea E un espacio métrico y Ø ¤ C E se dice que C es un conjunto conexo si para toda función continua f W C ! R se verifica que f .C / es un intervalo. Los conexos de R son, en virtud de la proposición anterior, los intervalos. Como la continuidad es una propiedad topológica, el concepto de conjunto conexo es topológico. Seguidamente lo vamos a caracterizar de una forma más cómoda, pues la definición que hemos dado hace intervenir a todas las funciones continuas f W C ! R y eso la hace incómoda de usar en casos concretos. Es muy fácil caracterizar los conjuntos que no son conexos. Un conjunto C no es conexo si, y sólo si, existe alguna función continua f W C ! R tal que f .C / no es un intervalo. Que f .C / 2 f .C / tal que 1; zŒ\f .C / ¤ Ø y no sea un intervalo equivale a que existe un z 2 R con z 6 1 z; C1Œ\f .C / ¤ Ø. Pero entonces se tiene que C D f . 1; zŒ/ [ f 1 .z; C1Œ/. Puesto que f es continua, sabemos, por la proposición 1.36 que los conjuntos A D f 1 . 1; zŒ/ y B D f 1 .z; C1Œ/ son abiertos relativos de C . Observa que A ¤ Ø y B ¤ Ø, A \ B D Ø y C D A [ B .
Hemos probado así que si C no es conexo entonces existen abiertos relativos de C disjuntos y no vacíos A y B tales que C D A [ B . Se dice que el par fA; B g es una partición no trivial de C por abiertos relativos. Recíprocamente, sea C un conjunto no vacío de un espacio métrico y supongamos que hay una partición no trivial de C por abiertos relativos fA; B g. Entones es fácil comprobar usando la proposición 1.36 que la función f W C ! R dada por: f .x/
1; si x 1; si x
D
2 A; 2 B.
es continua. Y como su imagen f .C / D f1; 1g no es un intervalo concluimos que C no es conexo. Hemos probado que C no es conexo equivale a que exista una partición no trivial de C por abiertos relativos. Dicho de otra forma:
1.56 Proposición. Un conjunto Ø ¤ C
E de un espacio métrico E es conexo si, y sólo si, la única partición de C por abiertos relativos es la trivial. Es decir, si A y B son abiertos relativos de C con A \ B D Ø y C D A [ B entonces el par fA; B g ha de ser la partición trivial fA; B g D fØ; C g. Al igual que vimos con la compacidad, la conexión se conserva por funciones continuas.
1.57 Proposición. Sean E y F espacios métrico, C un subconjunto conexo de E y f W C ! F una función continua. Entonces se verifica que f .C / es conexo.
Ejercicios propuestos 67. Prueba la proposición 1.57. Universidad de Granada Dpto. de Análisis Matemático
Prof. Javier Pérez Cálculo Diferencial en Rn
Continuidad y límites de campos escalares y vectoriales
29
68. Prueba que la unión de conexos con intersección no vacía es conexo. Los conjuntos abiertos y conexos se llaman dominios. Los dominios de Rn tienen una caracterización muy útil.
1.58 Proposición. Sea Rn un conjunto abierto. Equivalen las siguientes afirmaciones: a) es conexo. b) Todo par de puntos de puede unirse por medio de una curva sin salirse de . Es decir, para cada par de puntos a ; b existe una función continua γ Œ0; 1 Rn , tal que γ.0/ a , γ.1/ b y para todo t Œ0; 1 se tiene que γ.t / .
D
2
2
2
W
!
D
Intuitivamente, un dominio es un conjunto abierto de un solo trozo.
Universidad de Granada Dpto. de Análisis Matemático
Prof. Javier Pérez Cálculo Diferencial en Rn
Cap´ıtulo
2
Derivadas parciales y extremos relativos de campos escalares
Curvas en Rn Como ya debes saber, una curva en Rn es una aplicación continua γ W Œa; b ! Rn . El punto γ.a/ se llama origen y el punto γ.b/ extremo de la curva. Naturalmente, como γ.t / 2 Rn podremos expresarlo por medio de sus componentes en la base canónica que serán funciones de t . γ.t /
D . 1.t /; γ
γ2 .t / ; : : : ; γn .t //
Las funciones γk .t / se llaman funciones componentes de γ. Se dice que γ es derivable en un punto t cuando todas sus funciones componentes son derivables en dicho punto, en cuyo la derivada de γ en t es, por definición, el vector γ 0 .t /
D . 10.t /; 20.t / ; : : : ; n0.t // 0 .t 0 / ¤ 0, se define la recta tangente a γ
γ
γ
γ en el punto a (aunque Dado un punto a D γ.t 0 / tal que γ es más apropiado decir en el punto t 0 ) como la recta de ecuación paramétrica a C t γ 0 .t 0/, es decir, la recta que pasa por a con vector de dirección γ 0 .t 0 /.
Cuando se interpreta γ.t / como la función de trayectoria de un móvil, entonces su velocidad en un instante t es el vector γ 0 .t / y su rapidez es kγ 0 .t /k2 . La distancia que recorre dicho móvil entre dos instantes t D a y t D b viene dada por b
k
γ 0 .t / 2
a
k dt :
30
Derivadas parciales. Vector gradiente
31
También consideraremos en lo que sigue funciones γ W I ! Rn definidas y continuas en un intervalo I R no necesariamente cerrado y acotado. Las seguiremos llamando, por abuso de lenguaje, curvas.
2.1. Derivadas parciales. Vector gradiente Hemos visto en el capítulo 1 que los conceptos de continuidad y límite para funciones reales de una variable se generalizan fácilmente para campos escalares de varias variables. No ocurre lo mismo con el concepto de derivabilidad el cual no puede generalizarse de forma inmediata. La razón es que el concepto de derivabilidad hace intervenir la división de números reales, pues una derivada es un límite de cocientes incrementales, y en Rn no podemos dividir por vectores, es decir, la estructura algebraica de Rn no permite generalizar algo parecido a un “cociente incremental”. Si f es un campo escalar de dos o más variables, la expresión
f .a/ xa
f .x/
no tiene ningún sentido.
Otra diferencia importante es que en la recta real, R, solamente podemos acercarnos a un punto de ella a través de la propia recta, mientras que en Rn para n > 2 hay muchísimas más posibilidades de acercarse a un punto dado; por ejemplo, podemos acercarnos a través de cualquier curva que pase por dicho punto. Surge así una primera idea que consiste en acercarse a un punto dado a través de una recta dada. Parece que esta situación es más parecida a lo que conocemos para funciones reales de una variable. Las cuestiones que vamos a estudiar son de naturaleza algebraica y topológica, no dependen para nada de la norma que se use; por eso, en todo lo que sigue, representaremos con k k una norma en Rn . Si para ti es más cómodo, puedes considerar que siempre que aparezca una norma en Rn se trata de la norma euclídea. Cuando en algún resultado concreto sea preciso considerar la norma euclídea la representaremos en la forma usual por k k2 . Una dirección o un vector unitario en Rn es un vector de norma euclídea 1.
Dados un punto a 2 Rn y una dirección u, la recta que pasa por a con dirección u es la imagen de la aplicación γ W Rn ! R dada por γ.t / D a C t u, es decir, es el conjunto de puntos fa C t u W t 2 Rg. 2.1 Definición. Sea f un campo escalar definido en un conjunto abierto Rn , sea a 2 y u una dirección. Se define la derivada de f en a en la dirección u como el límite Du f .a/
D t l!Kım0 f .a C t ut / f .a/
(2.1)
supuesto, claro está, que dicho límite exista. La derivada direccional de un campo escalar f en un punto a en la dirección del vector ek de la base canónica, se llama derivada parcial (de primer orden) de f en a respecto a la variable k -ésima. Está definida por: Dek f .a/
D lKım
f .a
t !0
D
lKım
xk !ak
C t ek / f .a/ D lKım
f .a1 ; : : : ; ak
C t ; : : : ; a / f .a ; : : : ; a
t t 0 f .a1 ; : : : ; xk ; : : : ; an / f .a1 ; : : : ; ak ; : : : ; an / xk ak
Universidad de Granada Dpto. de Análisis Matemático
!
n
1
k ; : : : ; an /
t
(2.2)
Prof. Javier Pérez Cálculo Diferencial en Rn
Interpretación geométrica de las derivadas parciales y se representa con los símbolos D k f .a/ y
@f @xk
32
.a/.
Observa que las derivadas que acabamos de definir son derivadas de funciones reales de una variable real pues para calcular la derivada de un campo escalar f en un punto a en la dirección u lo que se hace es derivar en t D 0 la función t 7! f .a C t u/ que es una función real de una variable real. Observa que la segunda igualdad de (2.2) nos dice que, para calcular la derivada parcial Dk f .a/ , lo que se hace es derivar f respecto a la variable k-ésima considerando fijas las demás variables. Por eso se llaman derivadas parciales.
2.1.1. Interpretación geométrica de las derivadas parciales Es importante que entiendas el significado de las derivadas parciales de una función en un punto. Para poder visualizarlo vamos a considerar un campo escalar f de dos variables definido en un abierto R2 . Fijemos un punto .a; b/ 2 . Las derivadas parciales de f en .a; b/ son, por definición f .x; b/ f .a; b/ D t l!Kım0 f .a C t ; b/t f .a; b/ D x!a lKım xa f .a; b C t / f .a; b/ f .a; y/ f .a; b/ D lKım D2 f .a; b/ D lKım t y b t !0 y!b Es decir, lo que hacemos es derivar las funciones parciales x 7! f .x; b/ y y 7! f .a; y/ en los puntos x D a e y D b respectivamente. La gráfica de f , es decir, el conjunto S D f.x; y; f .x; y// W .x; y/ 2 g es una superficie D1 f .a; b/
en R3 . Las funciones
γ1 .x/
D .x; b; f .x; b//;
γ2 .y/
D .a; y; f .a; y//
son curvas contenidas en dicha superficie que pasan por el punto .a; b/. Dichas curvas se obtienen cortando la superficie S por los planos y D b y x D a respectivamente. Los vectores tangentes a dichas curvas en los puntos γ1 .a/ y γ2 .b/ son, respectivamente γ1 0 .a/
D .1; 0; D1 f .a; b//;
γ2 0 .b/
D .0; 1; D2 f .a; b//
En la figura (2.1) se ha representado la gráfica de f y las curvas obtenidas cortándola por los planos x D a e y D b junto a sus vectores tangentes en el punto .a; b/
Cuando un campo escalar f tiene derivadas parciales en todos los puntos de un conjunto abierto Rn , podemos definir las funciones derivadas parciales (de primer orden) de f , Dk f W ! R que a cada punto x 2 E hace corresponder el número Dk f .x/. Dichas funciones son también campos escalares.
2.2 Definición. Sea f un campo escalar. Se define el vector gradiente de f en un punto a como el vector
r f .a/ D
D1 f .a/; D2 f .a/ ; : : : ; Dn f .a/
supuesto, claro está, que dichas derivadas parciales existan. Universidad de Granada Dpto. de Análisis Matemático
Prof. Javier Pérez Cálculo Diferencial en Rn
Campos escalares diferenciables
33
Figura 2.1. Derivadas parciales
Supongamos que f es una función real de una variable real. La derivabilidad de f en un punto a 2 R se expresa por
f .a/ f .x/ f .a/ f 0 .a/.x a/ 0 D f .a/ ” x!a D0 lKım lKım x!a xa xa Recuerda que la recta de ecuación cartesiana y D f .a/ C f 0 .a/.x a/ es la recta tangente a f .x/
la gráfica de f en el punto .a; f .a//.
Si ahora f es un campo escalar definido en un abierto Rn , cuyo vector gradiente r f .a/ está definido en un punto a 2 , podemos considerar el hiperplano en RnC1 de ecuación cartesiana xnC1 D f .a/ C r f .a/ x a . Este hiperplano pasa por el punto .a; f .a// 2 RnC1 y es la generalización natural de la recta tangente a la gráfica de una función. Observa el parecido formal entre las expresiones
˝ ˇ ˛
y
D f .a/ C f 0 .a/.x a/;
xnC1
˝ ˇˇ ˛
D f .a/ C r f .a/ x a
Ambas representan hiperplanos (un hiperplano en R2 es una recta) y la segunda se deduce de la primera sustituyendo la derivada por el vector gradiente y el producto usual de números reales por el producto escalar de vectores. Esto nos lleva a la siguiente definición.
2.1.2. Campos escalares diferenciables 2.3 Definición. Sea f un campo escalar definido en un abierto Rn y sea a 2 . Supongamos que está definido el vector gradiente r f .a/. Se dice que f es diferenciable o derivable en a si se verifica que lKım
f .x/
x!a
Definamos R.x; a/
Universidad de Granada Dpto. de Análisis Matemático
˝ ˇ ˛ ˇ ˝r ˇ ˛
f .a/ r f .a/ x a D 0 kx ak
D f .x/ f .ak/x akf .a/ x
(2.3)
a
Prof. Javier Pérez Cálculo Diferencial en Rn
Campos escalares diferenciables
34
La igualdad (2.3) dice que lKım R.x; a/ D 0. Con lo que, otra forma equivalente de escribir la x!a igualdad (2.3) es la siguiente. f .x/
donde
˝ ˇ ˛
D f .a/ C r f .a/ x a C R.x; a/kx ak
lKım R.x; a/ D 0 (2.4)
x!a
La siguiente proposición expresa lo dicho en la definición anterior de una forma más abstracta. Recuerda que si ˛ W Rn ! R es una aplicación lineal entonces para todo x D n
.x1 ; x2 ; : : : ; xn /
X D
xj ej en Rn se tiene que:
j D1 n
˛. x/
X D
xj ˛. ej /
j D1
donde O˛ D ˛. e1 /;˛.e2 / ; : : : ; ˛ .en / .
˝ ˇˇ ˛
D ˛O x
2.4 Proposición. Sea f un campo escalar definido en un abierto Equivalen las siguientes afirmaciones:
Rn y sea a
2 .
a) f es diferenciable en a. b) Existe una aplicación lineal ˛ Rn
! R tal que f .x/ f .a/ ˛. x a/ D0 lKım ! kx ak W
x
a
(2.5)
Además, caso de que se cumpla b) se tiene que ˛ es la aplicación lineal dada por ˛. x/ f .a/ x .
˝ ˇ˛
D r
2.5 Definición. Sea f un campo escalar diferenciable en un punto a. El hiperplano en de ecuación cartesiana xnC1
RnC1
˝ ˇ ˛
D f .a/ C r f .a/ x a
se llama hiperplano tangente a f en a o hiperplano tangente a la gráfica de f en el punto .a; f .a//.
2.6 Proposición. Sea f un campo escalar diferenciable en un punto a y sea u una dirección en Rn . Entonces se verifica que
˝ ˇ˛ D C D C ˝r ˝ ˇ ˇ ˛˛C C D C r C C D ˝r ˇˇ ˛ C K C Du f .a/
D r f .a/ u
Demostración . En la igualdad (2.4) pongamos x f .a
C t u/
Deducimos que lKım
t !0
f .a
C t u/ f .a/ t
Universidad de Granada Dpto. de Análisis Matemático
f .a/ f .a/
t
a
t u con lo que obtenemos
f .a/ t u
R.a
f .a/ u
R.a
f .a/ u
lım R.a
t !0
k k t u; a/jt j
t u; a/ t u
t u; a/
jt j D r f .a/ u : t
˝ ˇˇ ˛
Prof. Javier Pérez Cálculo Diferencial en Rn
Campos escalares diferenciables
35
2.7 Corolario. Sea f un campo escalar diferenciable en un punto a con vector gradiente no nulo en a. a) La dirección en la que la derivada direccional de f en a es máxima es la dirección dada f .a/ por el gradiente, es decir, la dirección u . f .a/ 2
D kr r k
b) La dirección en la que la derivada direccional de f en a es mínima es la dirección f .a/ opuesta a la dada por el gradiente, es decir, la dirección v . f .a/ 2
D kr r k
Demostración . Las afirmaciones hechas son consecuencia de la proposición anterior y de la desigualdad de Cauchy–Schwarz, pues para toda dirección w se tiene que
˝ ˇˇ ˛
j r f .a/ w j 6 kr f .a/k2 kw k2 D kr f .a/k2 Y la igualdad se da si, y solo si, hay un número 2 R tal que w D r f .a/. Tomando normas en esta igualdad se deduce que jj D 1=kr f .a/k2 , es decir las direcciones w que hacen máximo jD f .a/j D j r f .a/ w j son u D kr r ff..aa//k y v D kr r ff..aa//k . 2 2 w
˝ ˇ˛ ˇ ˇ r D r ˇ kr k D kr k ˝r ˇr ˛ D kr
Para la primera se tiene que Du f .a/
f .a/
f .a/
f .a/
1
f .a/
2
f .a/
f .a/
f .a/
2
k2
que es el valor máximo que puede tener una derivada direccional. Análogamente, para la segunda se tiene que Dv f .a/
Dkr f .a/k2
que es el valor mínimo que puede tener una derivada direccional.
El resultado anterior nos dice que el vector gradiente en un punto señala la dirección en la que el campo tiene máximo crecimiento en dicho punto. Mientras que en la dirección opuesta a la del vector gradiente en un punto el campo tiene máximo decrecimiento.
2.8 Proposición. Sean f un campo escalar definido en un abierto Rn y γ W I ! Rn una
curva en Rn que toma valores en . Supongamos que γ es derivable en un punto t 0 y que f es diferenciable en el punto a γ.t 0 / . Entonces se verifica que la función h.t / f .γ.t // es derivable en t 0 y su derivada viene dada por:
D
2
D
n
h 0 .t 0 /
ˇ X D ˝r ˇ ˛ D f .a/
γ 0 .t 0 /
Dk f .a/γk0 .t 0 /
(2.6)
kD1
Si suponemos que f es diferenciable en y que γ es derivable en I entonces para todo t se tiene que:
2 I
n
h 0 .t /
ˇ ˛DX D ˝r f .γ.t //
γ 0 .t /
Dk f .γ.t //γk0 .t /
(2.7)
kD1
Universidad de Granada Dpto. de Análisis Matemático
Prof. Javier Pérez Cálculo Diferencial en Rn
Campos escalares diferenciables
36
Demostración . Se tiene que h.t /
˝
ˇˇ ˇˇ ˇ
˛
h.t 0/ D f . .t // f . .t 0// D r f .a/ .t / .t 0 / C R. .t /; .t 0//k .t / .t 0/k Dividiendo por t t 0 tenemos k .t / .t 0/k h.t / h.t 0 f . .t // f . .t 0 // .t / .t 0 / f .a/ R. .t /; .t 0 // D D r C t t 0 t t 0 t t 0 t t 0 .t / .t 0 / D 0.t 0/ se deduce que Teniendo en cuenta que lKım t !t t t 0 h.t / h.t 0 / D r f .a/ 0.t 0/ lKım t !t t t 0 γ
γ
γ
γ
γ
γ
γ
γ
γ
γ
γ
γ
γ
γ
γ
γ
γ
γ
γ
0
˝ ˇˇ ˛ γ
0
como queríamos demostrar.
2.9 Proposición. Sean f un campo escalar diferenciable en un abierto Rn y g W I ! R
una función derivable en un intervalo I . Supongamos que f ./ I y sea h g f Entonces el campo escalar h es diferenciable en y su gradiente viene dado por:
D ı W ! R.
r h.x/ D g 0.f .x//r f .x/
8x 2
(2.8)
Demostración . Sean a 2 y b D f .a/ 2 I , fijos en lo que sigue. Como g es derivable en b la función: '.t /
D g.t /t g.b/ g 0 .b/ t 2 I nfb g ; '.b/ D 0 b
es continua. Sustituyendo en la igualdad: g.t /
g.b/ D .g 0.b/ C '.t //.t b/
D f .x/, b D f .a/ y usando la igualdad (2.4), tenemos: g.f .x// g.f .a// D g 0 .b/ C '.f.x// .f .x/ f .a//D D g 0.b/ r f .a/ x a C g 0.b/ C '.f.x// kx akR.x; a/ C '.f.x// r f .a/ x a
t
˝ ˇˇ ˛
Por tanto: h.x/
0
˝ ˇˇ ˛
h.a/ g .b/ r f .a/ x a D kx ak
Podemos suponer que la norma, desigualdad de Cauchy-Schwarz:
ˇˇ ˇ
'.f.x//
Universidad de Granada Dpto. de Análisis Matemático
˝ ˇˇ ˛
g 0 .b/
˝r ˇˇ ˛
C '.f.x// R.x; a/ C '.f .x// fk.xa/ xak
a
k k, que hemos fijado en Rn es la euclídea para aplicar la
˝r ˇ ˛ ˇˇ k k ˇ f .a/ x x
a
a
j
jkr f .a/k
6 '.f.x//
Prof. Javier Pérez Cálculo Diferencial en Rn
Campos escalares diferenciables
37
Y teniendo en cuenta que lKım R.x; a/ D 0 y lKım '.f.x// D 0, deducimos que x!a
x!a
h.a/ g 0 .b/ r f .a/ x a D0 lKım ! kx ak Lo que prueba que h es diferenciable en a y que r h.a/ D g 0 .f .a//r f .a/.
˝ ˇˇ ˛
h.x/
x
a
Que un campo escalar tenga derivadas parciales en un punto es una propiedad muy débil. xy ; f .0; 0/ D 0 tiene derivadas parciales nulas Por ejemplo, el campo escalar f .x; y/ D 2 2 x
en .0; 0/ pero no es continuo en dicho punto.
Cy
La propiedad de ser diferenciable es mucho más fuerte que tener derivadas parciales.
2.10 Proposición. Un campo escalar diferenciable en un punto es continuo en dicho punto. Para campos escalares diferenciables hay un análogo al teorema del valor medio para funciones derivables reales. El segmento que une dos puntos u; v 2 Rn es el conjunto Œu; v Df u C t .v u/ W 0 6 t 6 1g.
2.11 Teorema (Del valor medio para campos escalares) . Sea f un campo escalar diferenciable
Rn. Sean a; b 2 y supongamos que el segmento Œa; b . Entonces 2 Œ a; b tal que se verifica la igualdad: f .b/ f .a/ D r f .c/ b a
en un abierto existe un punto c
˝ ˇ ˛
De donde se deduce:
jf .b/ f .a/j 6 kb ak2 sup fkr f .x/k2 W x 2 Œa; bg El siguiente resultado proporciona una condición suficiente de diferenciabilidad muy útil.
2.12 Teorema (Condición suficiente de diferenciabilidad). Un campo escalar que tiene derivadas parciales continuas en un conjunto abierto es diferenciable en todo punto de dicho conjunto.
La proposición 2.9 puede usarse para estudiar la diferenciabilidad de campos escalares. Con frecuencia un campo escalar f puede expresarse en la forma f .x/ D g.p.x// donde p es una función polinómica en n variables. Como las funciones polinómicas en n variables tienen derivadas parciales continuas, son diferenciables y, por la proposición antes citada, deducimos que el campo escalar f .x/ D g.p.x// es diferenciable en todo punto x tal que g sea derivable en p.x/. sen t 2.13 Ejemplo. La función g.t / D , para t ¤ 0 y g.0/ D 1 es derivable en R. Por tanto, el t sen .x 2 C y 2 / campo escalar f .x; y/ D y f .0; 0/ D 1 es diferenciable en R2 . 2 2 x
Universidad de Granada Dpto. de Análisis Matemático
Cy
Prof. Javier Pérez Cálculo Diferencial en Rn
Rectas tangentes y planos tangentes
38
En la práctica suele suponerse que los campos escalares tienen derivadas parciales continuas. Esta hipótesis garantiza que son diferenciables y es suficiente para justificar la mayoría de los resultados que siguen. Es sabido que una función derivable en un intervalo con derivada nula es constante. Para campos escalares hay un resultado análogo. Observa la hipótesis de que el campo esté definido en un dominio.
2.14 Proposición. Un campo escalar definido en un dominio con derivadas parciales nulas en todo punto del mismo es constante.
Ejercicios propuestos 69. Sean f; g W ! R campos escalares diferenciables en un abierto . Prueba que los campos escalares f C g, f g son diferenciables en y calcula sus vectores gradientes. Si, además, g.x/ ¤ 0 para todo x 2 entonces prueba que el campo escalar diferenciable en y calcula su vector gradiente.
f g
es
70. Sean f i , 1 6 i 6 n funciones reales derivables en un intervalo abierto I . Definamos f W I n ! R por: n
f .x1 ; x2 ; : : : ; xn /
X D
j D1
f j .xj /
8x D .x1; x2 ; : : : ; xn/ 2 I n
Prueba que el campo escalar f es diferenciable en I n .
71. Sean f; g W ! R campos escalares definidos en un abierto Rn . Sea a 2 y supongamos que f es diferenciable en a, f .a/ D 0 y g es continua en a. Prueba que f g es diferenciable en a.
2.2. Rectas tangentes y planos tangentes En esta sección vamos calcular rectas y planos tangentes a curvas y superficies considerando las distintas formas en que éstas pueden venir dadas. Mi propósito es esencialmente práctico, a saber, que entiendas la forma de proceder en cada caso; por lo que no me preocupo de justificar con detalle todo lo que digo.
Curvas en el plano Una curva en el plano puede venir dada de tres formas:
a) Como la gráfica de una función y D f .x/ donde x 2 I siendo I un intervalo de R.
Universidad de Granada Dpto. de Análisis Matemático
D f.x; f .x// W x 2 I g Prof. Javier Pérez Cálculo Diferencial en Rn
Rectas tangentes y planos tangentes
39
t / D .x .t /; t /; y .t // //. b) Por medio de ecuaciones paramétricas γ.t /
D
γ.I / I /
t /; y .t // // W t 2 I g D f.x .t /;
c) De forma implícita como el conjunto de puntos g.x ; y / D 0 donde se anula una función diferenciable diferenciable de dos variables. variables.
n D
2
2 W g .x ; y /
.x ; y / R
o D 0
Suele usarse la siguiente terminología. Si h.x ; y / es un campo escalar diferenciable, las curvas de ecuación implícita h.x ; y / D c o, lo que es igual h.x ; y / c D 0 , donde c es una constante, se llaman curvas de nivel. Dichas curvas se obtienen cortando la gráfica de h con planos de la forma z D c . Estas curvas son las que ves representadas en los mapas topográficos. g .x ; y / D f .x / y ) y también es Observa que a) es un caso particular de c) (basta considerar considerar g. un caso particular de b) (basta considerar γ.x / D .x ; f .x //).
La tangente en un punto de viene dada en cada caso como sigue.
a0 ) La tangente en un punto .a . a; b / D .a; f .a// 2 es la recta de ecuación cartesiana /.x a/ b D f 0 .a/.x
y
El vector .1 . 1; f 0 .a// es tangente a en el punto .a . a; b / y el vector .f 0 .a/; 1/ es ortogonal . a; b /. a en el punto .a
b0 ) La tangente en un punto γ .t 0 / D .a; b/ 2 es la recta de ecuaciones paramétricas .x ; y /
D
γ .t 0 /
C t
γ0 .t 0 /
D .a; b/ C t .t .x 0.t 0/; y 0.t 0//
El vector γ0 .t 0 / D .x 0 .t 0 /; y 0 .t 0 // es tangente a en .a . a; b /. . a; b / 2 es la recta de ecuación implícita c0 ) La tangente en un punto .a
˝r ˇˇ
g .a ; b / .x
˛
a; y b/ D 0
Se supone que r g.a; b/ ¤ 0 pues en otro caso, la tangente en .a; b / no está definida. El . a; b /. vector gradiente r g.a; b/ es ortogonal a en el punto .a Estas últimas afirmaciones requieren alguna justificación. Para ello, supongamos que co . a; b /. Es decir, hay nocemos una representación paramétrica local de en torno al punto .a ˛. t // D .˛1 .t /; t /; ˛2 .t // t // 2 que pasa por el punto .a; b / y que es una curva de la forma ˛.t 1 derivable . Pongamos ˛.t ˛. t 0 / D .a; b /. Por lo visto en b0 ), sabemos que la tangente a en .a; b / es la recta que pasa por el punto .a; b / con vector de dirección ˛ 0 .t 0 /. Pongamos h.t / t / D g .˛.t .˛. t // //. En virtud de la igualdad ( 2.7 t / D r g .˛.t .˛. t // // ˛ 0 .t / t / . Pero 2.7), ), tenemos que h 0 .t / 0 0 h.t / t / D 0, por lo que h .t / t / D r g .˛.t .˛. t // // ˛ .t / t / D 0. Resulta así que el vector r g .˛.t .˛. t // // es 0 ortogonal al vector tangente ˛ .t / t /. En particular, el vector r g .a; b / es ortogonal al vector
˝
1 El teorema de la función implícita,
ˇˇ ˛
˝
ˇ ˛
que se verá más adelante, garantiza la existencia de dicha curva siempre que
el vector gradiente r g .a; b/ ¤ 0.
Universidad de Granada Dpto. de Análisis Matemático
Prof. Javier Pérez Cálculo Diferencial en Rn
Rectas tangentes y planos tangentes
40
˛ 0 .t 0 / tangente a en .a; b /. Concluimos que la recta que pasa por .a; b / y tiene como vector ortogonal g.a; b / es la recta tangente a en .a . a; b /, pero dicha recta es justamente la recta de ecuación cartesiana g.a; b/ .x a; y b/ 0.
r
˝r ˇ
˛
D
De lo antes visto, merece la pena destacar la siguiente propiedad. El vector gradiente g .x ; y / de un campo escalar es ortogonal en todo punto .x ; y / (en el que g .x ; y / 0) a la curva de nivel que pasa por dicho punto.
r
¤
r
Superficies en R3 Una superficie S en el espacio R3 puede venir dada de tres formas: función y D f .x ; y / donde .x . x ; y / 2 A siendo A un conjunto de R2 . a) Como la gráfica de una función
D D f.x ; y; f .x ; y// W .x ; y / 2 Ag
S
t / D .x .s; t /; t /; y .s; t /; /; z .s; t // t // donde .s . s; t / t / 2 A R2 . b) Por ecuaciones paramétricas γ.s; t /
D D
S
γ.A/
D f.x .s; t /; t /; y .s; t /; t /; z .s; t // t // W .s; t / t / 2 Ag
g .x ; y ; z / D 0 donde se anula una función c) De forma implícita como el conjunto de puntos g. diferenciable diferenciable de tres variables. variables.
n D D
S
3
2 W g .x ; y ; z /
.x ; y ; z / R
o D 0
g .x ; y ; z / D f .x ; y / z ) y tambi Observa que a) es un caso particular de c) (basta considerar g. también én γ t / D .s; t ; f .s; t // t //). El plano tangente en un es un caso particular de b) (basta considerar .s; t / punto de S viene dado en cada caso como sigue. . a; b; c / D .a; b; f .a; b // 2 S es a0 ) El plano tangente en un punto .a es el plano de ecuación cartesiana z
Los vectores 1; 0;
@f /.x a/ C .a; b /.y /.y b / f .a; b/ D @@f x .a; b/.x @y
@f @f @x
.a; b / y 0; 1; @f @x
@f @f @y
.a; b / son tangentes a S en .a; b; c / y el vector en .a
.a; b /;
@f @f @y
.a; b /; 1
. a; b; c /. es ortogonal a S en el punto .a
b0 ) El plano tangente en un punto γ .s0 ; t 0 / D .a; b; c / 2 S es es el plano de ecuaciones paramétricas .x ; y ; z /
Donde
@γ @s
D
.s0 ; t 0 /
Universidad de Granada Dpto. de Análisis Matemático
γ .s0 ; t 0 /
D
C s @@s .s0; t 0/ C t @@t .s . s0 ; t 0 / γ
γ
@x @y @y @z @ z .s .s0 ; t 0 /; .s .s0 ; t 0 /; .s0 ; t 0 / @s @s @s
Prof. Javier Pérez Cálculo Diferencial en Rn
Rectas tangentes y planos tangentes
41
y
D
@γ .s . s0 ; t 0 / @t
@x
@y @y @z @ z .s0 ; t 0 /; .s . s0 ; t 0 /; .s . s0 ; t 0 / @t @t @t
Dichos vectores son tangentes a S en .a . a; b; c /.
. a; b; c / 2 S es c0 ) El plano tangente en un punto .a es el plano de ecuación implícita
ˇˇ
˝r
g .a ; b ; c / .x
˛
a; y b; z c / D 0
. a; b; c / no está Se supone que r g .a; b; c / ¤ 0 pues en otro caso, el plano tangente a S en .a definido. El vector gradiente r g.a; b; c / es ortogonal a S en el punto .a . a; b; c /. g .x ; y ; z / es un campo escalar, las superficies de ecuación implícita g. g .x ; y ; z / D c o, lo Si g. que es igual g.x ; y ; z / c D 0, donde c es una constante, se llaman superficies de nivel (cuando el campo se interpreta como un potencial se llaman superficies equipotenciales). De lo dicho en c0 ), se sigue que el vector gradiente r g.x ; y ; z / es ortogonal en todo punto .x ; y ; z / (en el que r g .x ; y ; z / ¤ 0) a la superficie de nivel que pasa por dicho punto.
Curvas en R3 Una curva en el espacio puede venir dada de dos formas.
a) Como intersección de dos superficies S 1 y S 2 .
t / D .x .t /; t /; y .t /; t /; z .t // t // donde t 2 I b) Por medio de ecuaciones paramétricas γ.t / intervalo.
D
γ.I / I /
R e I es
un
t /; y .t /; t /; z .t // t // W t 2 I g D f.x .t /;
La tangente en un punto de viene dada en cada caso como sigue. . a; b; c / 2 es la recta intersección de los planos tangentes a S 1 y a0 ) La tangente en un punto .a .a; b; c /. Por ejemplo, si las superficies vienen dadas por sus ecuaciones a S 2 en .a ecuaciones implícitas. S 1 S 2
D D
.x ; y ; z / R3 f .x ; y ; z / .x ; y ; z / R3 g .x ; y ; z /
˚˚
2 W 2 W
Dn ˇˇ ˇˇ
D0 D0
3
2 W g .x ; y ; z / D f .x ; y ; z /
.x ; y ; z / R
Entonces, las ecuaciones implícitas de la recta tangente son
8< ˝r : ˝r
o D 0
˛ ˛
c / D 0 g .a; b; c / .x a; y b; z c / D 0 Donde se supone que los vectores gradiente r f .a; b; c /; r g.a; b; c / son linealmente inf .a; b; c / .x
a; y
b; z
dependientes pues, en otro caso, la recta tangente a la curva en .a . a; b; c / no está definida.
b0 ) La tangente en un punto γ .t 0 / D .a; b; c / 2 es la recta de ecuaciones paramétricas .x ; y ; z /
D
γ .t 0 /
C t
γ0 .t 0 /
D .a; b; c / C t .t .x 0.t 0/; y 0.t 0/; z 0.t 0//
. a; b; c /. El vector γ0 .t 0 / D .x 0 .t 0 /; y 0 .t 0 /; z 0 .t 0 // es tangente a en .a
Universidad de Granada Dpto. de Análisis Matemático
Prof. Javier Pérez Cálculo Diferencial en Rn
Derivadas parciales de orden superior
42
2.3. Derivadas parciales de orden superior Supongamos un campo escalar f con derivadas parciales Dk f en un abierto Rn . Las funciones Dk f son también campos escalares de n variables que podemos, cuando se dejen, volver a derivar parcialmente en puntos de . Obtenemos de esta forma las derivadas parciales de segundo orden de f , es decir las funciones Dj k f D Dj .Dk f /, que se representan también de las formas: Dj k f .x/
@2 f
D @x @x j
.x/;
Djj f .x/
k
@2f
D @x 2 .x/ j
De manera análoga se definen las derivadas parciales de tercer orden de f como las derivadas parciales de las derivadas parciales de segundo orden de f , Dijk f D D i .Dj .Dk f //, y se representan también de las formas: 3
Dijk f .x/
@ f D @x @x @x i
j
3
I
.x/ Djjj f .x/ k
3
f D @@xf 3 .x/I Djj k f .x/ D @x@2@x j
j
.x/ k
Las derivadas parciales de cuarto orden son las derivadas parciales de las derivadas parciales de tercer orden y así sucesivamente. Es natural preguntarse si el orden en que se realizan las derivadas debe ser o no tenido en cuenta. Afortunadamente, en la mayoría de los casos podemos olvidarlo porque se verifica el siguiente utilísimo resultado.
2.15 Teorema (Schwarz). Las derivadas parciales de orden menor o igual que k de un campo escalar con derivadas parciales de orden k continuas solamente dependen del número de veces que se deriva parcialmente respecto de cada variable, pero el orden en que se realicen dichas derivaciones no afecta para nada al resultado final.
2.16 Definición. Se dice que un campo escalar f es de clase C k en un abierto si f tiene derivadas parciales de orden k continuas en , en cuyo caso escribimos f 2 C k ./. Se dice que f es de clase C 1 en si tiene derivadas parciales continuas de todos órdenes en .
Ejercicios propuestos
Como para calcular derivadas parciales de una función de varias variables se consideran fijas todas las variables menos aquella respecto a la que se deriva, calcular derivadas parciales es lo mismo que derivar funciones de una variable. Solamente debes tener cuidado para darte cuenta qué tipo de función es la que tienes que derivar porque ello puede depender de la variable respecto de la que derivas. Por ejemplo, la función f .x; y/ D x y cuando fijas y (para derivar respecto a x ) es una función potencia (la variable está en la base y el exponente está fijo) y cuando fijas x (para derivar respecto a y ) es una función exponencial (la variable está en el exponente y la base está fija). Te recuerdo que es muy frecuente, sobre todo en libros de Física e ingenierías diversas, representar las funciones por letras. Así, lo que los matemáticos solemos escribir f .x; y/ D cos.xy/ C xy 2 , para indicar que f es una función de dos variables x e y cuyo valor en el punto .x; y/ viene dado por cos .xy/ C xy 2 , suele expresarse de forma menos Universidad de Granada Dpto. de Análisis Matemático
Prof. Javier Pérez Cálculo Diferencial en Rn
Derivadas parciales de orden superior
43
precisa en la forma z D cos.xy/ C xy 2 , cuyo significado es exactamente el mismo que el anterior cambiando f por z . Naturalmente, en vez de z puede usarse cualquier otro símbolo que sea distinto de x e y . Tienes que acostumbrarte a esta notación y entender cuándo una letra representa una variable y cuándo representa una función.
73. Calcula las derivadas parciales de primer orden de los campos escalares: (a) f .x; y/Dx 2y Cz 2 x Cy sen.xz / (b) z D.x 2Cy 3 / ex y (c) w Dx ez Cz ey Cxyz . 74. Calcula las derivadas parciales de primer y segundo orden del campo escalar: f .x; y; z/
D 1 C yxy2 C z 2
75. Calcula las derivadas parciales de primer y segundo orden de los campos escalares: (a) z D sen cos exy
(b) w D log 4 C arctg.x=y/
(c) u D tg .xy/z
(d) v D arctg z x y
76. Estudia la existencia de las derivadas en el origen según las distintas direcciones del campo escalar f W R2 ! R definido por: x.x 2
y 2/ C f .x; y/ D x2 y2
si jx j ¤ jy j;
f .x; y/
D0
si jx j D jy j
¿Es f continuo en .0; 0/? Te recuerdo que una dirección viene dada por un vector de norma euclídea 1. Si a y b son puntos de Rn la dirección del punto a hacia el punto b viene dada por el vector
p C
77. Calcula la derivada direccional de f .x; y/ D log.1 C la dirección hacia el origen. 78. Calcula la derivada direccional de z.x; y/ D arctg dirección hacia el punto .2; 1/.
x2
xy
x2
C y2
a . kb ak b
y 2 / en el punto .1; 2/ en
en el punto .1; 1/ en la
79. Calcula valores de a , b y c para que la derivada direccional de la función f .x; y; z/
D axy 2 C byz C cz 2x 3
en el punto .1; 2; 1/ tenga un valor máximo igual a 64 en la dirección del eje OZ.
80. Considera la curva dada por las ecuaciones paramétricas x .t /Det C cos t , y.t /Det C sen t . Calcula la ecuación de la recta tangente en el punto .x.0/; y.0//. Universidad de Granada Dpto. de Análisis Matemático
Prof. Javier Pérez Cálculo Diferencial en Rn
Derivadas parciales de orden superior
44
81. Calcula, para los siguientes campos escalares, el vector normal en P 0 a la curva de nivel que pasa por dicho punto. 1. f .x; y/ D arctg 2. f .x; y/ D
!
y
p C 1
x2
sen.x C y/ 2 C cos.x y/
C
y2
P 0
D .1; 1/.
D .=2; =4/. xy en el punto .1; 1/ en la dirección 82. Calcula la derivada de h.x; y/ D 1 C log.1 C x 2 y 2 / P 0
dada por el vector ortogonal (de norma 1) en el punto .1; 1/ a la curva de nivel del campo f .x; y/ D x y 3 C x 3 y que pasa por dicho punto.
83. Calcula la ecuación de la recta tangente y de la recta normal a la elipse de ecuación x2 a2
y2
C b2 D 1
en un punto .u; v/ de la misma.
84. Calcula las ecuaciones del plano tangente y de la recta normal a cada una de las siguientes superficies en el punto P o indicado. z2
2x 2 2y2 12 D 0; z log.x 2 C y 2 / D 0; x 2 C y 2 C z 3 2x C 4y C 3z C 1 D 0; 4 x 2 4z 2 D y; z.xy 1/ .x C y/ D 0; z C ez C2x C 2y x 2 y 2 3 D 0;
P o .1; 1; 4/
I
P o .1; 0; 0/ P o .3; 4; 3/
I
P o .0; 0; 1/
I p P .1; 1 C e; 1/ P o .1; 2; 3/ o
85. Halla la ecuación de la tangente a la curva dada como intersección del elipsoide x 2 C 4y 2 C 2z 2 D 27 y el hiperboloide x 2 C y 2 2z 2 D 11 en el punto .3; 2; 1/. 86. Calcula la ecuación de la recta tangente a la curva definida por la intersección de las superficies z D x y , x 2 C y 2 2z D 4 en el punto .3; 1; 3/. Comprueba el resultado expresando la curva por sus ecuaciones paramétricas. 87. Calcula la ecuación de la recta tangente a la curva definida por la intersección de las superficies 4xz D .x C z/y , 3z 2 C y D 5x en el punto .1; 2; 1/. 88. Estudia la diferenciabilidad de los siguientes campos escalares a) f .x; y/ D
e x
2
Cy 2
x2
1 ; f .0; 0/ D 1
C y2
Universidad de Granada Dpto. de Análisis Matemático
b) f .x; y/ D
log .1 C x 2 C y 2 / x2
C y2
; f .0; 0/
D1
Prof. Javier Pérez Cálculo Diferencial en Rn
Derivadas parciales de orden superior
45
89. Estudia si los siguientes campos escalares definidos por f .0; 0/ D 0 y para .x; y/ ¤ .0; 0/ por: a) f .x; y/ D
sen.x 2y/
p C x2
y2
; b) f .x; y/
y xy 1 cos x cos y D sen x sen ; c/ f .x; y/ D x2 C y2 x2 C y2
p
son de clase C 1 en R2 .
p
90. Estudia si los siguientes campos escalares definidos por f .0; 0/ D 0 y para .x; y/ ¤ .0; 0/ por: a) f .x; y/ D xy
x 2 x2
y Cy
2
; b) f .x; y/ 2
y xy D y arctgxx C xy arc tg y ; c) f .x; y/ D arctgxx sen Cy 2
2
2
2
son de clase C 1 en R2 . Calcula en cada caso D 12 f .0; 0/ y D 21 f .0; 0/ e indica si son de clase C 2 en R2 .
91. Estudia la derivabilidad de las normas x 7! kxk2 , x 7! kxk1 y x 7! kxk1 y calcula sus vectores gradiente. 92. Justifica que un campo escalar de clase
k C
también es de clase C q para q 6 k .
93. Justifica que la suma, el producto y el cociente de campos escalares de clase C k también son campos escalares de clase C k . Deduce que toda función racional de n variables es de clase C 1 en su dominio natural de definición. 94. Generaliza las proposiciones 2.8 y 2.9 para campos de clase C k . 95. Un campo escalar f W Rn ! R se dice que es homogéneo de grado n 2 N , si para todo x 2 Rn y para todo t 2 R se verifica que f .t x/ D t n f .x/. Supuesto que f es homogéneo de grado n y derivable, prueba que: n
a)
X
xj Dj f .x/
j D1
8x D .x1; x2; : : : ; xn/ 2 Rn .
D nf .x/
b) Si n > 2, las derivadas parciales Dj f son funciones homogéneas de grado n 1. c) Si n D 1 entonces f es una forma lineal.
96. Prueba que si f es de clase C 2 y es homogéneo de grado 2 se verifica que: f .x/
1
D2
n
X
8x D .x1; x2; : : : ; xn/ 2 Rn
Dijf .0/xi xj
i;j D1
97. Sea f W ! R un campo escalar de clase C 1 . Supongamos que es un conjunto convexo y que 0 2 . Prueba que para todo x D .x1 ; x2 ; : : : ; xn / 2 se verifica que: n
f .x/
Universidad de Granada Dpto. de Análisis Matemático
D f .0/
X C
j D1
1
x D f .t x/ dt j
j
0
Prof. Javier Pérez Cálculo Diferencial en Rn
Teorema de Taylor. Extremos relativos
46 @f
2.17 Observación. ¿Qué notación es mejor para las derivadas parciales? ¿ o D1 f ?. Quizás @x un ejemplo te ayude a decidir: consideremos f .x; y/ D x ey . Tenemos que D1 f .x; y/ D ey y @f
.x; y/ ey . Si ahora permutamos x por y en estas igualdades obtenemos D 1 f .y; x/ ex @x @f y .y; x/ ex . La igualdad D1 f .y; x/ ex es correcta porque, efectivamente, e x es el valor @y de la derivada parcial de f .x; y/ respecto a la primera variable evaluada en el punto .y; x/. La @f @f .y; x/ ex no es correcta porque .x; y/ x ey y si la evaluamos en .y; x/ el igualdad @y @y resultado es y ex .
D D
D
D
D
D
Este ejemplo pone de manifiesto que es mucho mejor la notación de subíndices D 1 f , D2 f para las derivadas parciales que las notaciones
@f @x
y
@f @y
. La razón es que la notación D1 no
tiene nunca ambigüedad: indica derivación respecto a la primera variable y no depende para @f
nada de la forma en que representemos dicha variable, mientras que la notación nos dice @x que dicha variable la hemos representado con la letra x . Pero con frecuencia hay que cambiar las letras con las que representamos las variables.
2.4. Teorema de Taylor. Extremos relativos 2.18 Definición. Sea f un campo escalar definido en un conjunto E Rn . Se dice que f tiene un máximo relativo (resp. mínimo relativo) en un punto a 2 E , si a es un punto interior de E y existe un número r > 0 tal que B.a; r / E y f .x/ 6 f .a/ (resp. f .a/ 6 f .x/) para todo x 2 B.a; r /. Cuando las desigualdades anteriores se verifican de forma estricta para todo x 2 B.a; r / con x ¤ a se dice que el máximo o el mínimo relativo es estricto. Los puntos en los que f tiene un máximo o un mínimo relativos se llaman extremos relativos de f .
2.19 Proposición (Condición necesaria de extremo relativo). Sea f un campo escalar definido en un conjunto E Rn y supongamos que f tiene un extremo relativo en un punto a E y además que el vector gradiente de f en a está definido. Entonces se verifica que f .a/ 0. Es decir, las derivadas parciales de primer orden de f en a son todas nulas.
r
2 D
Demostración . Supongamos que f tiene un máximo relativo en a ysea r > 0 tal que B.a; r / E y f .x/ 6 f .a/ para todo x 2 B.a; r /. Definamos ' W r; r Œ! R por '.t / D f .a C t ek /. La función ' está definida en el intervalo r; r Œ pues para todo t 2 r; r Œ se tiene que ka C t ek ak D jt j < r por lo que a C t ek 2 B.a; r / E. Además, para todo t 2 r; r Œ se tiene que '.t / D f .a C t ek / 6 f .a/ D '.0/. Luego ' tiene en t D 0 un máximo relativo. Además como, por hipótesis, existe Dk f .a/, tenemos que ' es derivable en t D 0. Luego ' 0 .0/ D 0, pero ' 0 .0/ D Dk f .a/.
2.20 Definición. Los puntos donde se anula el gradiente de un campo escalar f se llaman puntos críticos o puntos estacionarios de f . Los puntos críticos de un campo escalar que no son extremos relativos se llaman puntos de silla.
Universidad de Granada Dpto. de Análisis Matemático
Prof. Javier Pérez Cálculo Diferencial en Rn
Teorema de Taylor. Extremos relativos
47
Si f es un campo escalar diferenciable, en los puntos críticos el hiperplano tangente es “horizontal”. La condición necesaria de extremo relativo no es suficiente. Por ejemplo, el campo escalar f .x; y/ D x 2 y 2 tiene un punto crítico en .0; 0/, pero no tiene extremo relativo en dicho punto pues en toda bola centrada en .0; 0/ toma valores positivos y negativos. Al igual que para funciones de una variable, la derivada segunda proporciona una condición suficiente de extremo relativo, para campos escalares de varias variables las derivadas parciales de segundo orden nos van a permitir dar una condición suficiente de extremo relativo. Necesitaremos para ello generalizar el teorema de Taylor para campos escalares de n variables. Para expresar dicho teorema de una forma sencilla es muy conveniente introducir una notación apropiada. Si f es un campo escalar y a; x 2 Rn definimos el operador: dk .f; a; x/ D .x1 D1 C x2 D2 C C xn Dn /kf .a/
(2.9)
Donde se entiende que la potencia k -ésima se desarrolla de forma simbólica. Un par de ejemplos aclarará esto.
0X D@
1 A
n
2
2
d .f; a; x/ D .x1 D1 C x2 D2 C C xn Dn / f .a/ 3
0X 1 A D@ n
3
d .f; a; x/
Observa que d1 .f; a; x/
xi xj Di Dj f .a/
i;j D1
1 A
n
xj Dj f .a/
j D1
0X D@
i;j ;kD1
f .a/ x .
xi xj Dijf .a/
i;j D1
n
xi xj xk Di Dj Dk f .a/
D ˝r ˇˇ ˛
n
X D
X D
xi xj xk Dijkf .a/
i;j ;kD1
2.21 Teorema (Taylor). Sea f un campo escalar de clase C kC1 en un abierto . Sean a 2 , x Rn tales que el segmento Œa; a x está contenido en . Entonces existe un punto c Œa; a x tal que se verifica la igualdad:
2
C
C
k
f .a
C x/ D f .a/
ˇ X C ˝r ˇ ˛ C f .a/ x
j D2
1 j!
dj .f; a; x/ C
1
C 1/!
.k
2
dkC1 .f; c; x/ (2.10)
Demostración . Tenemos que Œa; a C x D f a C t x W 0 6 t 6 1g . Como el conjunto ft 2 R W a C t x 2 g es abierto (es la imagen inversa de por la función continua t 7! a C t x) y contiene al intervalo Œ0; 1, hay un intervalo abierto I Œ0; 1 tal que a C t x 2 para todo t 2 I . Definamos γ W I ! Rn por γ.t / D a C t x. Consideremos la función h W I ! R dada por h.t / D f .γ.t //. Por la proposición 2.8 dicha función es derivable y su derivada viene dada por: n 0
h .t /
ˇˇ ˛ D ˝r ˇˇ ˛ D X D ˝r f .γ.t //
γ0
.t /
f .γ.t // x
j D1
n
Dj f .γ.t //xj
X D
ı
.Dj f
j D1
γ/.t /xj
Puesto que las funciones Dj f tienen derivadas parciales continuas son diferenciables. Podemos calcular la derivada de la función Dj f ı γ volviendo a usar la igualdad (2.9) como hemos Universidad de Granada Dpto. de Análisis Matemático
Prof. Javier Pérez Cálculo Diferencial en Rn
Teorema de Taylor. Extremos relativos
48
hecho antes pero con f sustituido por Dj f . Tenemos que: n
γ 0 .t /
ˇˇ ˛ D X
ı D ˝r
Dj f
.Dj f /. γ.t // x
n
h .t /
n
X X D j D1
Dk .Dj f /.γ.t //xk
kD1
Por tanto: 00
n
X D
Dkj f .γ.t //xk
kD1
n
! X D
Dkj f .γ.t //xk xj
kD1
Dkjf .γ.t //xk xj
j;kD1
D d2 .f;
γ.t /; x/
Análogamente se comprueba que h000 .t / D d3 .f; γ.t /; x/. En general, se tiene que h.j / .t / D dj .f; γ.t /; x/. Como las derivadas de orden k tienen derivadas parciales continuas son diferenciables, esto nos dice que h es k C 1 veces derivable en I con derivada de orden k C 1 continua en I . Podemos aplicar el teorema de Taylor a la función h para obtener que hay algún 2 Œ0; 1 tal que: k 0
h.1/
D h.0/ C h .0/
X C
j D2
1 j!
h.j / .0/
C .k C1 1/! h.kC1/./
Como h.1/ D f .a C x/, h.0/ D f .a/, h.j / .0/ D dj .f; a; x/ y h.kC1/ ./ D dkC1 .f; a C x; x/, la igualdad anterior es la misma (2.10) con c D a C x 2 Œa; x.
2.22 Teorema ( Taylor–Young). Sea f un campo escalar definido en un abierto
Rn y supongamos que f tiene derivadas parciales de segundo orden continuas en . Sea a y r > 0 tal que B.a; r / . Entonces para todo x con x < r se tiene que
kk
n
f .a
C x/ D f .a/
X C
Dk f .a/xk
kD1
2
n
1
X
C2
Dj k f .a/xj xk
j;kD1
Ckxk2'.x/
con lım '. x/
K
0
x!
D0
(2.11)
Demostración . Para cada x con kxk < r se tiene que a C t x 2 B.a; r / para 0 6 t 6 1, por lo que Œa; a C x B.a; r / . Por el teorema de Taylor hay algún punto c 2 Œa; a C x tal que n
f .a
C x/ D f .a/
X C
1
C 2
Dk f .a/xk
kD1
n
X
Dj k f .c/xj xk
j;kD1
Lógicamente, el punto c dependerá de x y será de la forma c D a C x donde es un número en Œ0; 1 que dependerá de x . Tenemos así que: n
C D
f .a x/ f .a/
X C
Dk f .a/xk
kD1
1
C2
Definamos '. x/
1
D 2kx k2
Universidad de Granada Dpto. de Análisis Matemático
n
X
Dj k f .a/xj xk
j;kD1
1
C2
n
X
Dj k f .c/ Dj k f .a/ xj xk
j;kD1
n
X
j;kD1
Dj k f .c/
Dj k f .a/ xj xk Prof. Javier Pérez Cálculo Diferencial en Rn
Teorema de Taylor. Extremos relativos
49
Con lo que, evidentemente, tenemos que: n
f .a
C x/ D f .a/
X C
1
C 2
Dk f .a/xk
kD1
n
X
C kxk2'.x/
Dj k f .a/xj xk
j;kD1
Podemos suponer que la norma que usamos es la euclídea, con lo que jxj j 6 kxk y tenemos: n
1
j'.x/j 6 2
X ˇˇ
Dj k f .c/
j;kD1
Dj k f .a/
n
ˇˇ j jj j X ˇˇ kk xi xj x
6
2
1
2
Dj k f .c/
j;kD1
Dj k f .a/
ˇˇ
Cuando x ! 0 se tiene que c D a C x ! a, por lo que, en virtud de la continuidad de las derivadas parciales de segundo orden, concluimos que l Kım '. x/ D 0.
0
x!
2.23 Definición. Sea f un campo escalar de n variables que tiene derivadas parciales de segundo orden continuas en y sea a 2 . La matriz n n H .f; a/
se llama matriz hessiana de f en a.
D
Dij f .a/
16i;j 6n
Observa que la matriz hessiana es simétrica porque Dij f .a/ D Dj i f .a/. En consecuencia, dicha matriz define una forma cuadrática, que representaremos por Q.f; a/, que viene dada para todo x D .x1 ; x2 ; : : : ; xn / 2 Rn por n
Q.f; a/.x/
t
D x H .f; a/ x .
.
X D
Dj k f .a/xj xk
j;kD1
donde el punto “ ” indica producto matricial y xt es el vector columna x. Con esta notación podemos escribir la igualdad (2.11) en la forma .
f .a
C x/ D f .a/ C r f .a/ x C 12 Q.f; a/.x/ C kxk2 '.x/
˝ ˇˇ ˛
donde
lKım '. x/ D 0 (2.12)
x!
0
Si suponemos que a es un punto crítico de f podemos escribir f .a
C x/ D f .a/ C 21 Q.f; a/.x/ C kxk2'.x/
donde
lKım '. x/ D 0 0
x!
(2.13)
De donde se sigue que f .a
C x/ f .a/ D 1 Q.f; a/.x/ C '.x/ kxk2 2kxk2
donde
lKım '. x/ D 0 0
x!
Teniendo en cuenta que las formas cuadráticas son polinomios homogéneos de grado 2, es decir, Q.f; a/.x/
la igualdad f .a
D 2Q.f; a/.x/, se tiene que 2kx1 k2 Q.f; a/.x/ D 12 Q.f; a/.x=kxk/. Resulta así
C x/ f .a/ D 1 Q.f; a/.x=kxk/ C '.x/ 2 kxk2
Universidad de Granada Dpto. de Análisis Matemático
donde
lKım '. x/ D 0 0
x!
(2.14)
Prof. Javier Pérez Cálculo Diferencial en Rn
Teorema de Taylor. Extremos relativos
50 n i;j D1
2.24 Definición. Una forma cuadrática Q.x/ D
P
˛ij xi xj se llama:
Positiva definida si Q.x/ > 0 para todo x 2 Rn con x ¤ 0. Semidefinida positiva si Q.x/ > 0 para todo x 2 Rn .
Positiva negativa si Q.x/ < 0 para todo x 2 Rn con x ¤ 0. Semidefinida negativa si Q.x/ 6 0 para todo x 2 Rn .
No definida o Indefinida si toma valores positivos y negativos.
2.25 Teorema. Sea f un campo escalar definido en un abierto Rn y supongamos que f
tiene derivadas parciales de segundo orden continuas en . Sea a un punto crítico de f y sea Q.f; a/ la forma cuadrática asociada a la matriz hessiana de f en a.
2
n
Q.f; a/.x/
t
D x H .f; a/ x .
.
X D
Dj k f .a/xj xk
j;kD1
a) Si la forma cuadrática Q.f; a/ es definida positiva entonces f tiene en a un mínimo relativo estricto. b) Si la forma cuadrática Q.f; a/ es definida negativa entonces f tiene en a un máximo relativo estricto. c) Si f tiene un máximo relativo en a entonces la forma cuadrática Q.f; a/ es semidefinida negativa. d) Si f tiene un mínimo relativo en a entonces la forma cuadrática Q.f; a/ es semidefinida positiva. e) Si la forma cuadrática Q.f; a/ es no definida entonces f tiene un punto de silla en a.
Demostración . Como Q.f; a/ es una función polinómica y, por tanto, continua, y la esfera unidad de Rn , S .0; 1/ D fu 2 Rn W kuk D 1g, es un conjunto compacto, dicha función alcanza un mínimo valor y un máximo valor en S .0; 1/. Sea m
D mKın fQ.f; a/.u/ W kuk D 1g ;
D maK x fQ.f; a/.u/ W kuk D 1g
M
a) Supongamos que Q.f; a/ es definida positiva. Entonces se tiene que m > 0 . y, por la igualdad (2.14), tenemos que
C x/ f .a/ D 1 Q.f; a/.x=kxk/ C '.x/ > m C '.x/ donde lKım '.x/ D 0 !0 2 2 kxk2 La condición lKım '. x/ D 0 garantiza la existencia de un número s > 0 tal que j'. x/j < m=4 !0 siempre que 0 < k xk < s . En consecuencia, si en la desigualdad anterior suponemos que 0 < kxk < s , se tiene f .a C x/ f .a/ m m m m C D x > '. / > > 0 2 2 4 4 kxk2 Deducimos que f .a C x / f .a/ > 0 para todo x con 0 < k xk < s . O, lo que es igual, f .z/ f .a/ > 0 para todo z tal que 0 < kz ak < s . Lo que prueba que f tiene en a un f .a
x
x
mínimo relativo estricto.
Universidad de Granada Dpto. de Análisis Matemático
Prof. Javier Pérez Cálculo Diferencial en Rn
Teorema de Taylor. Extremos relativos
51
Los demás puntos se prueban de forma parecida.
Para poder usar el resultado anterior hay que saber clasificar una forma cuadrática. Hay varios procedimientos sencillos para ello. Los dos que siguen a continuación son los que me parecen más cómodos.
Clasificación de formas cuadráticas Sean A D aij
16i;j 6n
una matriz simétrica de números reales y n
QA .x/
Dx Ax .
.
t
X D
aij xi xj
(2.15)
i;j D1
la forma cuadrática definida por A. Los valores propios de A son las raíces del polinomio característico p ./, que se define como el determinante de la matriz A I : p./
ˇ
ˇ
D A I
Es sabido que, en la situación que estamos considerando, las raíces de dicho polinomio son todas reales. Sean j .1 6 j 6 n/ los valores propios de A. Se demuestra que hay una base B D fu1; u2 ; : : : ; un g en Rn tal que para todo vector x 2 Rn se tiene que n
QA .x/
X D
j xj2
j D1
donde .x1 ; x2 ; : : : ; xn / son los coordenadas del vector x en la base B . De aquí se siguen los siguientes criterios.
La forma cuadrática Q A es definida positiva si, y sólo si, todos los valores propios de son positivos.
A
La forma cuadrática QA es definida negativa si, y sólo si, todos los valores propios de son negativos.
A
La cuadrática Q A es no definida si, y sólo si, A tiene valores propios positivos y negativos.
La forma cuadrática Q A es semidefinida positiva si, y sólo si, todos los valores propios de A son mayores o iguales que 0. La forma cuadrática Q A es semidefinida negativa si, y sólo si, todos los valores propios de A son menores o iguales que 0.
Para aplicar estos criterios no es preciso calcular los valores propios de A sino solamente saber cuántos de ellos son positivos, negativos o nulos. Afortunadamente, hay un criterio que nos proporciona esta información sin más que observar los coeficientes del polinomio característico.
Universidad de Granada Dpto. de Análisis Matemático
Prof. Javier Pérez Cálculo Diferencial en Rn
Teorema de Taylor. Extremos relativos
52
2.26 Proposición (Regla de los signos de Descartes). Sea f .x/
D anx n C an1x n1 C C a1 x C a0
un polinomio con coeficientes reales y cuyas raíces son todas números reales. Se verifica entonces que: a) El número de raíces positivas de f (contando multiplicidades) es igual al número de cambios de signo en la sucesión .an ; an1 ; : : : ; a1 ; a0 / de los coeficientes de f . b) El número de raíces negativas de f (contando multiplicidades) es igual al número de cambios de signo en la sucesión .. 1/n an ; . 1/n1 an1 ; : : : ; a1 ; a0 / de los coeficientes de f . x/.
Para contar los cambios de signo en la sucesión de coeficientes se saltan los coeficientes nulos. Por ejemplo, si f .x/ D 2x 6 C x 5 x 3 C x 2 5, la sucesión de coeficientes de f es .2; 1; 0; 1; 1; 0; 1/ cuyo número de cambios de signo es 3.
2.27 Corolario. Sea p./ el polinomio característico de la matriz hessiana de f en a. Entonces. Si p./ tiene grado n , todos sus coeficientes son distintos de cero y tienen igual signo, se verifica que f tiene un máximo relativo estricto en a .
Si p./ tiene grado n , todos sus coeficientes son distintos de cero y van alternando su signo, se verifica que f tiene un mínimo relativo estricto en a .
Otro criterio para estudiar el carácter de la forma cuadrática (2.15) es el criterio de Sylvester se basa en la sucesión de signos de los menores principales de la matriz A. El menor principal de orden k es el determinante k D ai;j 16i;j 6k . Se verifica que:
ˇ ˇ
La forma cuadrática es definida positiva si, y sólo si, k > 0 para 1 6 k 6 n. La forma cuadrática es definida negativa si, y sólo si, . 1/k k > 0 para 1 6 k 6 n.
Observa que cuando la dimensión n es par, si el determinante de la matriz entonces la forma es no definida.
A
es negativo
Podemos particularizar este criterio para el caso de dos dimensiones. Sea A R2 un conjunto abierto y sea f un campo escalar definido en A que tiene derivadas parciales de segundo orden continuas. Supongamos que .a; b/ 2 A es un punto crítico de f y sea H .f;.a; b//
0 D B@
@2 f .a; b/ @x 2 @2 f .a; b/ @x@y
@2 f
.a; b/ @x@y @2 f .a; b/ @y 2
1 CA
la matriz hessiana de f en .a; b/ y notemos det H .f;.a; b// su determinante. @2 f Si detH .f;.a; b// > 0 y 2 .a; b/ > 0 entonces f tiene en .a; b/ un mínimo relativo @x
estricto.
Universidad de Granada Dpto. de Análisis Matemático
Prof. Javier Pérez Cálculo Diferencial en Rn
Teorema de Taylor. Extremos relativos
53
@2 f
Si detH .f;.a; b// > 0 y 2 .a; b/ < 0 entonces f tiene en .a; b/ un máximo relativo @x estricto. Si detH .f;.a; b// < 0 entonces f no tiene extremo relativo en .a; b/. Se dice que .a; b/ es un punto de silla de f . Cuando det H .f;.a; b// D 0 el conocimiento de la matriz hessiana no permite decidir si hay o no hay extremo relativo en .a; b/. Cuando esto sucede puede ser interesante estudiar el comportamiento de las curvas f .a; t C b/ y f .a C t ; b/. Si alguna de dichas curvas no tiene extremo relativo o tienen extremos relativos de distinta naturaleza en t D 0, podemos concluir que en .a; b/ no hay extremo relativo de f .
Ejercicios propuestos 98. Generaliza el teorema de Taylor – Young para funciones de clase
k C .
99. Determinar los extremos relativos de las funciones: f .x; y/ D 2xy 2x 3 y xy 2 C x 3 y 2 I
f .x; y/ D 2x 3 6xy 2 C 3x 2 C 3y2 I
f .x; y/ D x 3 C 3xy 2 15x 12yI
f .x; y/ D x 3 C y 3 xy 2 x C 16I
f .x; y/ D xy C 2x 3 y xy 2 2x 3 y 2 I
f .x; y/ D cos.x/ cos.y/
f .x; y/ D 2x C y C x 2 C xy C y 3 I
f .x; y/ D x 2 y 2 x 2 y 2 I
f .x; y/ D x log y x
f .x; y/ D 2x 4 C y 4 4x 2 2y 2 I f .x; y; z/ D x 3 C y 3 C z 3 xyz 3x 3y 3zI
f .x; y/ D xy.1 x y/I f .x; y; z/ D x 2 C y 2 C z 2 C 3xyz x y zI
f .x; y; z/ D .x C y C z/ e1
x 2 y 2 z 2
100. Trazar un plano que pase por el punto .1; 2; 3/ y que forme con los ejes coordenados un tetraedro de volumen mínimo (el volumen del tetraedro es un tercio del área de la base por la altura). 101. Recta de mínimos cuadrados. Dados n puntos .xi ; yi / 2 R2 , determinar los números ˛ n
y ˇ para que la cantidad
X yi
iD1
2
˛ xi
ˇ
sea mínima. n
102. Dados m puntos a i 2 R , calcular el valor mínimo de la función f .x/ n
Universidad de Granada Dpto. de Análisis Matemático
X D k x
iD1
ai
k22.
Prof. Javier Pérez Cálculo Diferencial en Rn
Cap´ıtulo
3
´ de campos vectoriales Derivacion
3.1. Derivada de un campo vectorial. Matriz jacobiana Una función vectorial es cualquier función que toma valores en un espacio vectorial de dimensión mayor mayor que 1. Las curvas en el plano o en el espacio son funciones vectoriales de una variable. Ahora nos interesa considerar funciones vectoriales de varias variables. La siguiente definición expresa la idea básica del cálculo diferencial: aproximar localmente una función por una aplicación lineal.
3.1 Definición. Sea F W ! Rm , donde Rn es un abierto, y sea a 2 . Se dice que F es diferenciable o derivable en a si hay una aplicación lineal T W Rn ! Rm tal que: lKım
x!
0
F.x/ F.a/ T.x a/
kx ak
D0
(3.1)
Dicha aplicación lineal, si existe, es única; se llama diferencial o derivada de F en el punto a y se representa por D F.a/. La siguiente proposición pone las cosas mucho más claras y, de paso, probamos la unicidad de la aplicación lineal T que verifica (3.1).
3.2 Proposición. Sea F D .f 1; f 2 ; : : : ; fm / W !
Rm , donde Rn es un abierto, una función vectorial de n variables y m componentes y sea a . Equivalen las siguientes afirmaciones:
2
a) F es diferenciable en a. componentes de F son diferenciables en a . b) Los campos escalares f 1 ; f 2 ; : : : ; fm
54
El espacio normado L.Rn ; Rm /
55
Supuesto que se verifican a) y b), la matriz de la aplicación lineal D F.a/ en las bases canónicas de Rn y Rm es la matriz cuyas filas son los vectores gradiente f i .a/ , esto es la matriz de m filas y n columnas Dj f i .a/ 16i 6m . Dicha matriz se llama matriz jacobiana de
F en a y se representará por J F.a/.
r
16j 6n
La aplicación lineal D F.a/ Rn
W
! Rm está definida para todo x 2 Rn por D F.a/.x/ D J F .a/ xt .
donde “ ” indica producto matricial y x t es el vector columna x . .
Observa que en el caso particular de que m D 1, es decir, F D f es un campo escalar, la aplicación derivada es una forma lineal sobre Rn que viene dada por:
D r f .a/ x 8x 2 Rn En el caso en que n D 1, es decir, F D es una curva en Rm , la aplicación derivada es una Df .a/.x/
˝ ˇ˛
γ
m
aplicación lineal de R en R que viene dada por:
D t 0.a/ 8t 2 R Y en el caso en que n D m D 1, esto es F D h es una función real de variable real, la aplicación D γ.a/.t /
γ
derivada de h en a es la aplicación lineal de R en R dada por: Dh.a/.t /
D t h 0 .a/
8t 2 R
La idea básica del cálculo diferencial es aproximar (localmente) una función por su derivada, que es una aplicación lineal. Puesto que la aplicación derivada es una buena aproximación local de la función, cabe esperar que la función herede localmente algunas de las propiedades de su aplicación derivada. Las propiedades de una aplicación lineal son muy fáciles de estudiar y lo que se hace en casi todos los resultados del cálculo diferencial es trasladar localmente las propiedades de la aplicación derivada a la función. Un ejemplo aclarará esto. Considera una función real h derivable en un intervalo I ; la aplicación derivada de h en cada punto a 2 I es la aplicación lineal de R en R Dh.a/.t / D t h 0 .a/. Si para a 2 I es h 0 .a/ > 0 dicha aplicación es creciente, pero entonces también la función h es creciente en I . Este ejemplo puedes llevarlo a otras situaciones más complicadas: si la aplicación derivada es una biyección lineal ¿será la función localmente inversible? Claro está, la respuesta que se obtiene a este tipo de preguntas es siempre de tipo local porque la derivada es una buena aproximación local y solamente proporciona información local sobre el comportamiento de la función.
3.1.1. El espacio normado L.Rn ; Rm / Antes de seguir conviene que consideremos el conjunto L.Rn ; Rm / de todas las aplicaciones lineales de Rn en Rm . Con las operaciones usuales de suma y producto por un número real, dicho conjunto es un espacio vectorial real de dimensión nm. En dicho espacio se pueden considerar muchas normas, por ejemplo, podemos identificar una aplicación lineal T con su matriz
Universidad de Granada Dpto. de Análisis Matemático
Prof. Javier Pérez Cálculo Diferencial en Rn
El espacio normado L.Rn ; Rm /
56
en las bases canónicas de Rn y Rm , T D t kj
16k 6m, 16j 6n
seguidamente dicha matriz la identifi-
camos con un vector de Rnm y definimos la norma euclídea de T como la norma euclídea de dicho vector:
pX X m
kTk2 D
n
a2kj
(3.2)
kD1 j D1
De esta forma identificamos totalmente el espacio .L.Rn ; Rm /; k k2 / con .Rnm ; k k2 /. Observa que si representamos por tk la fila k-ésima y por tj la columna j-ésima de la matriz T se tiene que: m
n
X X kk D k k D k k 2 2
T
2 k 2
tj
t
2 2
j D1
kD1
La norma que acabamos de definir en el espacio L.Rn ; Rm / tiene algunas propiedades interesantes. Notemos con k k2 la norma euclídea en Rn y en Rm . Sea x 2 Rn . Tenemos que:
T.x/ D T xt D t1 x ; t2 x ; : : : ; tm x
˝ ˇˇ ˛ ˝ ˇˇ ˛ ˝ ˇˇ ˛ X X ˇ ˝ ˛ ˇ kD k k k k Dk k k k .
Por lo que, usando la desigualdad de Cauchy–Schwarz, tenemos: m
kT.x/ Es decir:
m
2
2 2
tk x 6
kD1
tk
2 2
x
2 2
T
2 2
x
2 2
kD1
kT.x/k2 6 kTk2 kxk2
Además, si S 2 L.Rm ; Rq / y su matriz en las bases canónicas es S
(3.3)
D sik
16i 6q , 16k 6m
se tiene
que la matriz asociada a la aplicación H D S ı T 2 L.Rn ; Rq / es la matriz producto H D S T . Como el elemento h ij de la fila i, columna j de H viene dado por
D E ˇ D si tj
hij
Tenemos que: q
q
n
n
2
q
n
D E X X X X X X ˇ ˇ k ı k D k kk kD X ! 0@X 1A D k k k k Dk k k k S T
2 2
h2ij
iD1 j D1
j
si t
6
2 2
tj
2 2
iD1 j D1
n
2 i 2
tj
s
iD1
si
6
iD1 j D1
q
Por tanto:
.
2 2
S
2 2
T
2 2
j D1
kS ı Tk2 6 kSk2 kTk2
(3.4)
La ventaja que tiene la norma euclídea en L.Rn ; Rm / es que es fácil de calcular. Sin embargo, las desigualdades (3.3) y (3.4) dependen esencialmente de que las normas consideradas en los espacios Rn y Rm sean las normas euclídeas. Podemos proceder de otra forma más abstracta, y por tanto más cómoda en teoría, para definir normas en L.Rn ; Rm / que sigan verificando las desigualdades anteriores aunque las Universidad de Granada Dpto. de Análisis Matemático
Prof. Javier Pérez Cálculo Diferencial en Rn
El espacio normado L.Rn ; Rm /
57
normas fijadas en Rn y Rm no sean las euclídeas. Supongamos, pues, que es una norma en Rn y una norma en Rm . Sea T 2 L.Rn ; Rm /. Como la esfera unidad en .Rn ; /, es decir el conjunto fx 2 Rn W .x/ D 1g, es compacto y la aplicación x 7! T.x/ es continua por ser composición de las aplicaciones continuas y T , por la propiedad de compacidad dicha aplicación debe alcanzar un máximo absoluto en la esfera unidad de . Rn ; /. Definimos:
˛. T/
D maK x
8T 2 L.Rn; Rm /
T.x/ .x/ D 1
˚ W
(3.5)
Es muy fácil comprobar que la aplicación T 7! ˛. T/ así definida es una norma en L.Rn ; Rm /. Dicha norma se llama la norma de operadores en L.Rn ; Rm / asociada a las normas y . Dicha norma verifica que:
T.x/ 6 ˛. T/.x/
8x 2 Rn
(3.6)
Además, si S 2 L.Rm ; Rq /, en Rq tenemos fijada una norma , y representamos por ˇ la norma de operadores en L.Rm ; Rq / asociada a las normas y , y por γ la norma de operadores en n q L.R ; R / asociada a las normas y , se verifica que: γ.
S ı T/ 6 ˇ.S/˛. T/
(3.7)
Demasiadas normas ¿verdad? Una notación en exceso detallista acaba siendo confusa, por eso en la práctica todas las normas en los espacios que estamos considerando se representan igual con la notación usual para normas k k. Si te fijas un poco esto no puede dar lugar a confusión. Consideremos, pues, que en Rn y en Rm se han fijado normas que representamos con el mismo símbolo k k. Dichas normas inducen una norma de operadores que también representaremos con el mismo símbolo k k. La desigualdad (3.6) se expresa ahora:
kT.x/k 6 kTkkxk
8x 2 Rn
(3.8)
Observa que aunque en esta desigualdad hay tres normas en tres espacios distintos no cabe confusión, pues kxk es la norma fijada en Rn , kT.x/k es la norma fijada en Rm y kTk es la norma de operadores correspondiente, y no cabe otra interpretación. Análogamente, la desigualdad (3.7) se expresa:
kS ı Tk 6 kSkkTk
(3.9)
Y no hay posibilidad de confusión porque cada norma está en su espacio de operadores correspondiente. En lo que sigue tú eliges: o bien consideras que todas las normas son euclídeas o bien entiendes que en cada caso se han fijado normas en los espacios Rn que inducen las correspondientes normas de operadores. Es más cómodo y más general adoptar este último punto de vista. Todos los conceptos propios de espacios normados tienen perfecto sentido en el espacio n m L.R ; R /. De hecho, como tal espacio es de dimensión finita en él son válidos los teoremas 1.22, 1.24, 1.25 y 1.28. Por supuesto, tiene perfecto sentido considerar la convergencia de sucesiones en L.Rn ; Rm /. Representaremos por GL.Rn / el conjunto de todas las aplicaciones lineales biyectivas (isomorfismos) de Rn en Rn . Notaremos por L.Rn / el espacio vectorial de los operadores lineales de Rn en Rn . Universidad de Granada Dpto. de Análisis Matemático
Prof. Javier Pérez Cálculo Diferencial en Rn
El espacio normado L.Rn ; Rm /
58
3.3 Proposición. Sea T 2 L.Rn /. Equivalen las siguientes afirmaciones: a) T GL.Rn /.
2
b) T es inyectivo. c) T es sobreyectivo.
3.4 Proposición. a) GL.Rn / es un conjunto abierto en L.Rn /. Concretamente: Para todo T GL.Rn / se verifica que:
2
2 S
n
L.R
/
1
W kT Sk < kT1k
En particular, si I es la identidad en Rn entonces B.I; 1/ b) La aplicación T
GL.Rn /
GL.Rn /.
7! T1 de GL.Rn/ en GL.Rn/ es continua.
Ejercicios propuestos 103. Sea T 2 L.Rn; Rm /. ¿Cuál es la derivada de T en un punto a 2 Rn ? 104. Prueba que una función vectorial diferenciable en un punto a es continua en a. 105. Prueba que la definición dada por la igualdad (3.5) es una norma en L.Rn ; Rm / y que se verifican las desigualdades (3.6) y (3.7). 106. Supongamos que fSn g ! S en el espacio L.Rm; Rq / y que fTn g ! T en el espacio n m n q L.R ; R /. Prueba que fSn ı Tn g ! S ı T en el espacio L.R ; R /. 107. Sea a 2 Rn y ˛ W Rn ! R la aplicación lineal dada por ˛.x/ D a x . Calcula k˛ k cuando en Rn se consideran las normas k k2 , k k1 y k k1 .
˝ ˇˇ ˛
108. Sea T W R2 ! R2 la aplicación lineal de .R2 ; k k1 / en sí mismo, cuya matriz en las bases canónicas es
a b c d
. Prueba que
kTk D maK x fjaj C jc j; jbj C jd jg 109. Sea F W ! Rm , donde Rn es abierto, una función diferenciable en un punto a 2 . Prueba que para cada " > 0 existe ı > 0 tal que para todos x ; y 2 B.a; ı/ se verifica que
kF.x/ F.y / DF.a/.x y /k 6 "kx y k: Si D F.a/ es inyectiva, entonces existe > 0 tal que F es inyectiva en B.a; /.
Universidad de Granada Dpto. de Análisis Matemático
Prof. Javier Pérez Cálculo Diferencial en Rn
Regla de la cadena. Derivadas parciales de funciones compuestas
59
Si F es diferenciable en a, definiendo: R.a; x/
La igualdad lKım
0
x!
D F .x/ F.ak/xDaFk.a/.x a/
F.x/ F.a/ D F.a/.x a/
kx ak
D0
(3.10)
puede escribirse de forma equivalente como:
F.x/ D F.a/ C D F.a/.x a/ C kx akR.a; x/ con
lKım R.a; x/ D 0 (3.11)
x!a
3.5 Definición. Se dice que un campo vectorial es de clase C k (k 2 N [f1g) en un abierto si lo son sus campos escalares componentes. 3.6 Proposición. Sea F D .f 1 ; f 2 ; : : : ; fn / W ! Rm donde es un abierto en Rn . Equivalen las afirmaciones:
a) F es de clase C 1 en . b) F es diferenciable en y la aplicación x
k C
7! D F.x/ de en L.Rn; Rm / es continua.
Además, como consecuencia de lo visto en el ejercicio 92, todo campo vectorial de clase también es de clase C q para q 6 k .
3.1.2. Regla de la cadena. Derivadas parciales de funciones compuestas 3.7 Teorema ( Regla de la cadena ). Sean F W B ! Rm , B Rn , y G W A ! Rn , A Rq , funciones vectoriales tales que G.A/ B de manera que la composición H D F ı G W A ! Rm está definida. Supongamos que G es diferenciable en un punto a 2 A y que F es diferenciable en el punto G .a/ 2 B. Entonces se verifica que la función compuesta H es diferenciable en a , y su diferencial viene dada como la composición de las respectivas diferenciales: D H.a/
D DF.G.a// ı DG.a/
(3.12)
Además, si F y G son de clase C k entonces H también es de clase C k .
Observa que la composición tiene sentido pues D G.a/ W Rq ! Rn y D F..G.a// W Rn ! Rm , por lo que la composición es una aplicación lineal de Rq a Rm , como debe ser, pues H es una función vectorial de q variables y m componentes. La expresión de la igualdad (3.12) por medio de matrices jacobianas es:
D J F.G.a// J G.a/ (3.13) Poniendo H D .h1 ; h2 ; : : : ; hm /, F D .f 1 ; f 2 ; : : : ; fm /, G D .g1 ; g2 ; : : : ; gq /; notando las variables por x D .x1; x2 ; : : : ; xn / 2 Rn , y D .y1 ; y2 ; : : : ; ym / 2 Rq , y escribiendo b D G.a/, J H .a/
.
tenemos que:
@hi
@yj
.a/
16i 6m 16j 6q
D
Universidad de Granada Dpto. de Análisis Matemático
@f i @xk
.b/
16i 6m 16k 6n
.
@gk @yj
.a/
16k 6n 16j 6q
b
D G.a/
Prof. Javier Pérez Cálculo Diferencial en Rn
Regla de la cadena. Derivadas parciales de funciones compuestas
60
De donde se sigue: @hi @yj
n
.a/
X D
kD1
@f i @xk
.b/
@gk @yj
.a/
b
D G.a/
.1 6 i 6 m; 1 6 j 6 q/
(3.14)
Esta igualdad constituye la regla de la cadena para derivadas parciales y es importante que aprendas a aplicarla y que entiendas lo que dice. Voy a intentar facilitarte las cosas. Primero, lo más frecuente es que F sea un campo escalar. Supongamos, pues, que en lo anterior, F D f es un campo escalar, en cuyo caso h D f ı G también es un campo escalar. La igualdad (3.14) queda ahora: @h @yj
n
.a/
X D
kD1
@f @xk
.b/
@gk @yj
.a/
b
D G.a/
.1 6 j 6 q/
(3.15)
En esta igualdad se interpreta que la función G W A ! E Rn lo que hace es un “cambio de variables”. Hablando familiarmente, podemos decir, que las “variables antiguas” de la función f , esto es las x D .x1; x2 ; : : : ; xn / 2 E se han sustituido por “variable nuevas” y D .y1 ; y2 ; : : : ; yq / 2 A y la función f se ha “expresado en estas nuevas variables” dando lugar a la función h . La relación entre unas variables y otras viene dada por:
D gk .y1; y2 ; : : : ; yq /;
xk
16k 6n
(3.16)
De esta manera podemos interpretar la igualdad ( 3.15) en la forma siguiente: Para derivar la función nueva h, respecto a una nueva variable yj , se deriva la función antigua f respecto a cada una de sus variables xk y se multiplica por la derivada de cada una de ellas x k D gk .y1 ; y2 ; : : : ; yq / respecto a la variable yj . Ya se ve que la situación está pidiendo que hagamos algunas simplificaciones que, además, son las que se hacen siempre en la práctica porque, aunque son algo confusas, facilitan mucho los cálculos. Lo primero que se hace es identificar las funciones g k que introducen las nuevas coordenadas con las coordenadas antiguas x k , es decir, vemos las coordenadas antiguas como funciones de las nuevas y esto lo escribimos en la forma siguiente.
D xk .y1; y2 ; : : : ; yq /;
xk
16k6n
(3.17)
Con esta notación, la igualdad (3.15) queda como sigue. @h @yj
n
.a/
X D
kD1
@f @xk
.b/
@xk @yj
.a/
b
D G.a/
.1 6 j 6 q/
(3.18)
Observa el doble papel que desempeña a la derecha de esta igualdad la letra xk ; cuando se deriva respecto de ella representa una variable y cuando ella se deriva respecto de una variable nueva representa una función. La igualdad (3.18) ya es bastante fácil de recordar pero todavía se siguen haciendo en la práctica, sobre en todo en los textos de Física que suelen usar notaciones muy desafortunadas, Universidad de Granada Dpto. de Análisis Matemático
Prof. Javier Pérez Cálculo Diferencial en Rn
Regla de la cadena. Derivadas parciales de funciones compuestas
61
algunas simplificaciones adicionales (y peligrosas). A saber: no se distingue entre la función f y la función h porque, como suele decirse en esos textos aludidos, son “la misma función expresada en distintas variables”. Haciendo la identificación de f con h nos queda lo siguiente. @f @yj
n
.a/
X D
kD1
@f @xk
.b/
@xk @yj
.a/
b
D G.a/
.1 6 j 6 q/
(3.19)
Aquí la letra f desempeña un doble papel: a la izquierda es la función compuesta y a la derecha es la función dada en sus variable iniciales. Todavía suele darse un pasito más que consiste en representar la función f con una letra que suele usarse para representar variables; a saber, la letra z . Esto es frecuente también en textos de Física. Vamos a hacerlo así. @z @yj
n
.a/
X D
kD1
@z @xk
.b/
@xk @yj
.a/
b
D G.a/
.1 6 j 6 q/
(3.20)
Todavía hay algo que podemos simplificar. Habrás observado que siempre indico la relación que hay entre los puntos b y a . Eso es muy importante para entender lo que se hace. Hay que saber dónde se evalúan las derivadas parciales de cada función. Pues bien, eso no se indica jamás en textos de Física. Nunca se indica en dónde se evalúan las derivadas parciales. Así que vamos a suprimirlo. @z @yj
n
X D
kD1
@z @xk
.1 6 j 6 q/
@xk @yj
(3.21)
Debes de familiarizarte con esta igualdad y saber reconocer en ella la igualdad de partida. Y no olvides la forma en que se evalúa esta igualdad. Lo vuelvo a poner. @z @yj
n
.y /
X D
kD1
@z @xk
.G.y //
@xk @yj
.y /
.1 6 j 6 q/
(3.22)
Si tuviéramos que volver a derivar en esta igualdad respecto a una variable yk se derivaría como de costumbre: la derivada de una suma es la suma de las derivadas y para derivar el producto se aplica la regla usual. Pero hay un detalle muy importante y es que la función @z
@z @xk
.G.y // vuelve
con la función G. Por tanto para derivarla @xk hay que aplicarle la misma regla que hemos aplicado para derivar z como función compuesta a ser la función compuesta del campo escalar
y que nos ha llevado a la igualdad anterior. Es por eso que el cálculo de derivadas parciales de segundo orden en funciones compuestas suele ser bastante engorroso y es fácil equivocarse si no se sabe lo que se hace.
3.8 Ejemplo. Vamos a calcular
@z @x
siendo z D u2 C v 5 C 3uv donde u D x 2 C y 2 , v D sen.xy/.
Así es como suelen enunciarse estos ejercicios y debes entender bien el enunciado. Nos están dando una función de las variables .u; v/ a la que llaman z . Esto es la letra z representa una función, a saber, z D u2 C v 5 C 3uv . Nos están dando un cambio de variables por medio de las igualdades u D x 2 C y 2 , v D sen.xy/. Y nos piden calcular Universidad de Granada Dpto. de Análisis Matemático
@z
@x
. Esto último ya nos dice
Prof. Javier Pérez Cálculo Diferencial en Rn
Regla de la cadena. Derivadas parciales de funciones compuestas
62 @z
claramente que debemos ver z como función de x e y , es decir, la letra z en es la función @x que nos dan después de sustituir en ella las nuevas variables, o sea, la función compuesta de z D u2 C v 5 C 3uv con G.x; y/ D .x 2 C y 2 ; sen.xy//. Sabemos que
@z
@z @u @z @v .2u C 3v/2x C .5v 4 C 3u/y cos.xy/ D C D @x @u @x @v @x @z
Si lo dejamos así escrito parece que depende de 4 variables. Pero no es así porque en la @x igualdad anterior las variables son x e y (las nuevas variables) mientras que u y v (las antiguas variables) vienen dadas por u D x 2 C y 2 , v D sen.xy/. Por tanto, es mejor hacer la sustitución, con lo que resulta @z
D .2.x 2 C y 2 / C 3 sen.xy//2x C .5 sen4 .xy/ C 3x 2 C y 2 /y cos.xy/ @x que nos da el valor de la derivada parcial de la función compuesta en un punto .x; y/. En este caso es muy sencillo calcular la función compuesta. Hazlo y comprueba el resultado obtenido.
Ejercicios propuestos Consideremos una función de dos variables x e y , z D z.x; y/, y supongamos que expresamos x e y en función de nuevas variables u y v , lo que indicamos en la forma x D x.u; v/ , y D y.u; v/ . De esta forma la función z es función (función compuesta) de las “variables libres” u y v , a través de las “variables dependientes” x e y . Se trata de calcular las derivadas parciales de z respecto de las nuevas variables u y v . La regla para hacerlo es la siguiente: para derivar una función z
D z.x; y/;
x
D x.u;v/; y D y.u; v/
respecto de una nueva variable, se deriva z respecto de cada una de las antiguas variables y se multiplica por la derivada de cada antigua variable respecto de la nueva variable. Se entiende mejor si lo escribimos simbólicamente @z
@z @x @z @y D C @u @x @u @y @u
En esta igualdad debes darte cuenta de que a la izquierda, como estamos derivando respecto a u, la letra z representa a la función compuesta z Dz.x.u;v/; y.u;v// y la derivada está calculada en un punto .u; v/ . En la parte derecha de la igualdad la letra z representa la función dada z D z.x; y/ y las letras x e y representan variables (cuando se deriva respecto de ellas) y funciones (cuando se derivan respecto de u). Debe entenderse que cuando se sustituye un valor de .u; v/ en la igualdad los valores de x e y deben sustituirse por x D x.u; v/ , y D y.u; v/ . Universidad de Granada Dpto. de Análisis Matemático
Prof. Javier Pérez Cálculo Diferencial en Rn
Regla de la cadena. Derivadas parciales de funciones compuestas
63 @z
111. Sea z D cos.xy/ C ey1 cos x donde x D u2 C v, y D u v 2 . Calcular en el punto @u .u; v/ D .1; 1/. 112. Sea u D .x C y/4 C y 2 .z C x/3 donde x D r s et , y D r s log.1 C t 2 /, z D r 2 s cos t . Calcula
@u @s
cuando r D 2, s D 1, t D 0.
113. Sea z D f .x; y/, y pongamos x D u2 C v 2 , y D u=v. Calcular las derivadas parciales de de z respecto de las nuevas variables u y v en función de las derivadas parciales de z respecto de x e y . 114. Sea u D x 4y C y 2 z 3 C ' .x=y/, donde
Calcular
@u @s
8ˆ< D ˆ: DD
C r s et
x
1
y
r s 2 et
z
r 2 s sen t
cuando r D 2, s D 1, t D 0, sabiendo que ' 0 .3=2/ D 1. @z
115. Sea z D f .x; y/ donde x D s 4 C r 4 , y D 2 r s 2 . Calcula @z @r
.1; 1/
@x
@z
.2; 2/ y
@y
.2; 2/. Siendo
.1; 1/ D 3. D 2 y @z @s
y 116. Prueba que la función F .x; y/ D f . x2y 2 /, donde f es una función real derivable, verifica la igualdad
.x 2
@F 2x y D 0 C y 2 / @F C @x @y
117. Prueba que la función F .u; v/ D f .uv; .u2 v 2 /=2/, donde f W R2 ! R es una función diferenciable, verifica la igualdad 2
.u
Cv
2
/
2
2
2
2
" # C D C @f
@f
@x
@F
@y
@F
@u
@v
118. Sea z D f .x; y/, donde x D cos # , y D sen # . Calcula @z=@ y @z=@# y prueba que 2
2
2
C D C @z
@x
@z
@y
@z
@
1
@z
2
@#
119. Sea g.s; t / D f .s 2 t 2 ; t 2 s 2 /. Prueba la igualdad t
2
@g
@g C D 0. s @s @t
120. Sea u D f .x; y/ donde x D es cos t , y D es sen t . Justifica que 2
@u
2
@u C D e2s 2 2 @x @y
Universidad de Granada Dpto. de Análisis Matemático
2
@u
2
@u C @s 2 @t 2
! Prof. Javier Pérez Cálculo Diferencial en Rn
Teorema del valor medio para campos vectoriales
64
121. Sea z D f .x; y/, donde x D cos # , y D sen # . Prueba que @2z
@2z
@2z
1 @2z
1 @z C D C C @ @x 2 @y 2 @2 2 @# 2
122. Sea z D f .x; y/ donde x D x.u; v/ , y D y.u; v/ . Prueba que @2z @u2
@2 z
D @x 2
2
C @x
@u
2
@2 z @x @y
@2 z
2
C @y
C @y 2
@x@y @u @u
@u
@z @2x @x @u2
@z @2y
C @y @u2
E indica la forma en que se evalúan estas funciones.
123. Sea h.x; y/ D f .u; v/ donde f es de clase C 2 y u D u.x; y/, v D v.x; y/ son funciones de clase C 2 . Prueba que: D11h.x; y/
C D22h.x; y/
D
.D1 u/
2
C .D1v/
2
.D11 f .u;v//
C D22f .u;v///
124. Sean las funciones f .x; y; z/ D .ex Cy 2 ; ez Cy/, g.u; v/ D v 2 C log u para .u; v/ 2 R RC . ¿Qué valor debe tener para que la derivada direccional máxima de g ı f en .0; 0; 0/ sea igual a 1? 125. Sea F W ! Rn un campo vectorial diferenciable e inyectivo en un abierto Rn . Sea W D F./ y supongamos que W es abierto y F 1 es diferenciable en W . Para y 2 W expresa D F1.y / en función de la derivada de F.
3.1.3. Teorema del valor medio para campos vectoriales En el teorema 2.11 vimos una generalización del teorema del valor medio para campos escalares que se expresaba por la igualdad: f .b/
˝ ˇ ˛
f .a/ D r f .c/ b a
(3.23)
No existe una generalización del mismo tipo para funciones vectoriales. Por ejemplo, sea γ .t / D .cos t ; sen t /. Tenemos que γ 0 .t / D . sen t ; cos t / es un vector unitario y, como γ .2 / γ.0/ D .0; 0/, no existe ningún t 2 R que verifique la igualdad γ.2/ γ.0/ D γ 0 .t /2 , es decir
no podemos expresar el incremento de la función igual al valor de la derivada en un punto intermedio por el incremente de la variable. Lo que se generaliza para funciones vectoriales no es una igualdad análoga a (3.23) sino una desigualdad : si tomamos valores absolutos en (3.23)
obtenemos:
jf .b/ f .a/j 6 kr f .c/k2kb ak2
(3.24)
Es esta desigualdad la que se generaliza para campos vectoriales con independencia, además, de las normas que se usen. Recuerda que la principal utilidad del teorema del valor medio es que permite acotar el incremento de una función por el incremento de la variable cuando se conoce una cota de la derivada, por ello nada se pierde con la generalización que vamos a obtener. Probaremos primero un lema que tiene interés por sí mismo. Universidad de Granada Dpto. de Análisis Matemático
Prof. Javier Pérez Cálculo Diferencial en Rn
Teorema de la función inversa
65
3.9 Lema. Sea γ W Œ0; 1 ! R n derivable y supongamos que kγ 0 .t /k 6 M para todo t 2 Œ0; 1.
Entonces se verifica que:
k .1/ γ
γ.0/
k 6 M
3.10 Teorema (Del valor medio para funciones vectoriales) . Sea F una función vectorial diferenciable en un abierto Rn . Supongamos que Œ x; y z Œx; y . Entonces se verifica que:
y que kDf .z/k 6 M para todo
2
kf .y / f .x/k 6 M ky xk 3.11 Corolario. Un campo vectorial derivable en un dominio con derivada nula en todo punto del mismo es constante.
3.2. Teorema de la función inversa Un resultado que debes conocer de primer curso es el siguiente:
Teorema de inversión global para funciones reales. Si f W I ! R es una función derivable en un intervalo I cuya derivada no se anula en ningún punto de I , entonces f es una biyección de I sobre el intervalo J f .I / , y su inversa f 1 J R es derivable en J siendo
D
W !
f 1 0 .y/
1 D f 0 .f 1 .y//
8y 2 J
Recuerda que si f es una función real de una variable real derivable en un punto a, la aplicación derivada de f en a es la aplicación lineal de R en R dada por t 7! f 0 .a/t . Que dicha aplicación lineal sea invertible equivale a que f 0 .a/ ¤ 0. También debes saber que si T 2 L.Rn /, las matrices que representan a T en distintas bases de Rn tienen el mismo determinante, dicho determinante se llama el determinante de la aplicación lineal T y lo representaremos por det .T/. La aplicación lineal T es invertible equivale a que det .T/ ¤ 0.
Si F W ! Rn es un campo vectorial definido en un abierto Rn y diferenciable en un punto a 2 , el determinante de la derivada, det .D F.a//, se llama el jacobiano de F en a. Teniendo en cuenta lo antes dicho, el determinante jacobiano de F en a es el determinante de la matriz jacobiana J F .a/. La derivada D F.a/ es invertible si, y sólo si, det .J F .a// ¤ 0. Pues bien, es fácil dar ejemplos de campos vectoriales de más de una variable que son diferenciables y cuya derivada en todo punto es invertible pero ellos no lo son. Es decir, el anterior teorema de inversión global para funciones reales de una variable real no puede generalizarse literalmente para campos vectoriales diferenciables de varias variables.
3.12 Ejemplo. El campo vectorial de R2 en R2 dado por F .x; y/ D . ex cos y; ex sen y/ es, evidentemente, diferenciable (es de clase C 1 ) en R2 . Su matriz jacobiana en un punto .x; y/ es ex cos y ex sen y J F .x; y/ D x e sen y ex cos y
Universidad de Granada Dpto. de Análisis Matemático
Prof. Javier Pérez Cálculo Diferencial en Rn
Teorema de la función inversa
66
cuyo determinante es igual a e 2x > 0. Sin embargo, dicho campo no es inyectivo en R2 pues F.x; y C 2 / D F.x; y/, por lo que la función F no tiene inversa. Observa, sin embargo, que la función F es inyectiva en R0; 2Œ. Una versión local del teorema de inversión global es la siguiente.
Teorema de inversión local para funciones reales. Si f W I ! R es una función derivable
en un intervalo abierto I cuya derivada no se anula en un punto a I y es continua en dicho I tal que a H y la restricción de f a H es punto. Entonces hay un intervalo abierto H R es derivable en una biyección de H sobre el intervalo J f .H / , y su inversa f 1 J J siendo 1 f 1 0 .y/ y J f 0 .f 1 .y//
D
2
D
2
W !
8 2
Como vamos a ver, este teorema de inversión local sí puede generalizarse para campos vectoriales. Sea, pues, F D .f 1 ; f 2 ; : : : ; fn / W ! Rn un campo vectorial diferenciable en un abierto Rn . Queremos dar condiciones que garanticen que para cada y D .y1 ; y2 ; : : : ; yn / 2 F./ el sistema de ecuaciones: f 1 .x1 ; x2 ; : : : ; xn / f 2 .x1 ; x2 ; : : : ; xn / : :: : :: :: :: :: : :: f n .x1 ; x2 ; : : : ; xn /
D D D
y1 y2
: ::
(3.25)
yn
tiene solución única en x D .x1 ; x2 ; : : : ; xn / 2 . Dicho de otra forma, queremos saber si dicho sistema permite “despejar” x 1 ; x2 ; : : : ; xn en función de y D .y1 ; y2 ; : : : ; yn /. Además, queremos que las soluciones sean funciones diferenciables de y . Así planteado, la pregunta que nos estamos haciendo es si F es inyectiva en . Acabamos de ver con un ejemplo que en general esto no es así: los campos vectoriales no suelen ser funciones inyectivas en dominios muy grandes. Además, las herramientas que vamos a usar proporcionan información sobre el comportamiento local de la función. Por eso, lo primero que hacemos es localizar el problema. Dado a 2 , nos preguntamos si hay conjuntos abiertos U y V Rn con a 2 U tales que para todo y 2 V las ecuaciones 3.25 tienen solución única x 2 U . Naturalmente, las funciones f j pueden ser de muy variada naturaleza y sería una ingenuidad suponer que hay algún tipo de algoritmo que permita resolver dicho sistema. Hay un caso en que sí sabemos resolverlo: cuando las ecuaciones son lineales. Pues bien, en situaciones como ésta lo que se hace es linealizar el problema, es decir, sustituimos la función F por su derivada y trasladamos el problema a dicha aplicación derivada, con lo cual hemos convertido nuestro problema inicial en un problema de álgebra lineal mucho más simple. En nuestro caso, la linealización del problema conduce al sistema de ecuaciones lineales D F.a/.x/ D y : D1 f 1 .a/x1
C D2f 1.a/x2 C C Dnf 1.a/xn D D1 f 2 .a/x1 C D2 f 2 .a/x2 C C Dn f 2 .a/xn D
y1
C D2f n.a/x2 C C Dnf n.a/xn D
yn
: :: :: :: :: : :: :: :
D1 f n .a/x1
Universidad de Granada Dpto. de Análisis Matemático
y2 : ::
Prof. Javier Pérez Cálculo Diferencial en Rn
Teorema de la función inversa
67
Para cada y 2 Rn dicho sistema tiene solución única en x si, y sólo si, el determinante jacobiano det.J F .a// ¤ 0. Un entorno abierto de un punto es un conjunto abierto que contiene a dicho punto.
3.13 Teorema (de la función inversa ). Sea F D .f 1; f 2 ; : : : ; fn / W ! Rn un campo vectorial de clase C 1 en un abierto Rn . Sea a 2 y supongamos que det .J F .a// ¤ 0. Entonces existe un entorno abierto de a , U , tal que F.U / D V es abierto, F es una biyección de U sobre V , la biyección inversa G W V ! U Rn es de clase C 1 en V y para todo y 2 V es D G.y /
D
D F.x/
1
donde x
D G .y /
(3.26)
Si F es de clase C k en entonces G es de clase C k en V .
Demostración . 1 Como es natural, suponemos fijada una norma k k en Rn . Dicha norma puede ser cualquiera pues los conceptos que figuran en el enunciado del teorema son todos de naturaleza algebraico– topológica. Sea n D :.n: : . Consideramos n como espacio métrico con la métrica producto (ver el ejercicio 19). La aplicación ' W n ! R definida por: '. x1 ; x2 ; : : : ; xn /
D det
Dj f i .xi /
16i 6n 16j 6n
es continua porque las derivadas parciales son continuas y el determinante es una suma de n 2 n
productos de la forma
Y ˙
Djif i .xi / donde 1 6 j i 6 n. Como por hipótesis '.a; a; : : : ; a/
¤ 0, iD1 existe una bola abierta en n en la cual ' no se anula. Dicha bola será de la forma B :.n: : B donde B D B.a; r / es una bola abierta en Rn . Para todo x 2 B se tiene que .x; x; : : : ; x/ 2 B :.n: : B por lo que '. x; x; : : : ; x/ D det.J F .x// ¤ 0. Para x ; y 2 B se verifica que Œ x; y B , por lo que podemos aplicar el teorema del valor medio 2.11 a cada componente f i de F para obtener que existe zi 2 Œx; y tal que: f i .y / f i .x/ D r f i .zi / y x i D 1; 2; : : : ; n
ˇˇ ˛
˝
Lo que podemos expresar matricialmente en la forma:
F.y / F.x/ D Dj f i .zi / .y x/t
.
Como . z1 ; z2 ; : : : ; zn / 2 B :.n: : B , tenemos que '.z1; z2 ; : : : ; zn / ¤ 0. Por tanto, si y ¤ x, de la igualdad anterior deducimos que F.y / ¤ F.x/. Hemos probado así que F es inyectiva en B. Pongamos U D B . Hemos probado que det .J F .x// ¤ 0 para todo x 2 U y que F es inyectiva en U . Seguidamente probaremos que la imagen por F de todo conjunto abierto contenido en U es un conjunto abierto en Rn . 1 Esta demostración está tomada del excelente libro Análisis Matemático del profesor Carlos Ivorra Castillo que
puedes bajar de Internet.
Universidad de Granada Dpto. de Análisis Matemático
Prof. Javier Pérez Cálculo Diferencial en Rn
Teorema de la función inversa
68
Sean W U un abierto y v 2 F.W /. Sea c 2 W tal que F .c/ D v . Tomemos s > 0 tal que la bola euclídea cerrada B.c; s/ W . Sea D Fr.B.c; s// D fx 2 Rn W kx ck2 D s g. La aplicación h W ! R dada por h.x/ D kF.x/ v k2 verifica que h.x/ > 0 para todo x 2 y es continua. Por compacidad, h alcanza un mínimo absoluto en un punto x 0 2 . Sea m D h.x0 / > 0 . Probaremos que B.v ; m=2/ F.B.c; s// donde las bolas son bolas euclídeas. Sea y 2 B.v ; m=2/ fijo en lo que sigue. Definimos la función g W B.c; s/ ! R por g.x/ D kF.x/ y k2. Dicha aplicación es continua y, por compacidad, alcanza su mínimo absoluto en un punto z 2 B.c; s/. Como g.c/ D kv y k2 < m=2, tenemos que g.z/ < m=2. Para x 2 se verifica que:
D kF.x/ y k2 > kF.x/ v k2 kv y k2 > m m2 D m2 Luego g.x/ no es el mínimo de g , es decir, x ¤ z. Por tanto z 2 B.c; s/. g.x/
Como g es una función que toma valores positivos, es claro que g 2 alcanza su mínimo en z. Como n
2
g .x/
D kF.x/
X k D y
2 2
.f i .x/
iD1
yi /2
es una función diferenciable 2 , dicho punto debe ser un punto crítico de g 2 . Deducimos que: n 2
Dj g .z/
X D 2
Dj f i .z/.f i .z/
iD1
yi / D 0
1 6 j 6 n
Igualdades que podemos expresar matricialmente como:
y /t D 0 Como det.J F .z// ¤ 0, concluimos que debe ser F.z/ D y . Como y es cualquier punto en B.v ; m=2/ y z 2 B.c; s/, hemos probado que F.B.c; s// B.v ; m=2/ lo que implica que B.v ; m=2/ F.W /. Por tanto F.W / es un conjunto abierto. En particular V D F.U / es un conjunto abierto. La aplicación F es una biyección de U sobre V . Sea G W V ! U Rn la biyección inversa definida por: 8y 2 V G.F.x// D x 8x 2 U ; F.G.y // D y Si O Rn es abierto, W D O \ U U es abierto, y se tiene que G1 .O / D G1 .W / D F.W / J F .z/ .F.z/ .
es abierto. Las proposiciones 1.35 y 1.36 nos dicen que G es continua en V .
Probaremos que G es diferenciable en V . Sea v 2 V y c D G.v / 2 U . Pongamos T D D F.c/ y S D T1. La función ' W ! R definida por '.c/ D 0 y para x 2 U , x ¤ c, por la igualdad:
kF.x/ F.c/ T.x c/k D '.x/kx ck (3.27) es continua en c porque F es diferenciable en c. Pongamos y D F.x/ 2 V donde x 2 U , con lo que G.y / D x. Tenemos: kG.y / G.v / S.y v /k D kS y v T.G.y / G.v // k6 6 kSkky v T.G.y / G.v //k D por (3.27) D D kSk'.G.y //kG.y / G.v /k
2 Para eso trabajamos con la norma euclídea.
Universidad de Granada Dpto. de Análisis Matemático
Prof. Javier Pérez Cálculo Diferencial en Rn
Teorema de la función inversa
69
Por otra parte:
kG.y / G.v /k D kx ck D kS ı T.x c/k 6 kSkkT.x c/k6 6 kSk .kF.x/ F.c/ T.x c/k C kF.x/ F.c/k/ D D kSk .'.x/kx ck C kF.x/ F.c/k/ D D kSk .'.G.y //kG.y / G.v /k C ky v k/ Como la función ' ı G W V ! R es continua en v y '.G.v // D '. c/ D 0, existe ı > 0 tal que 1 '. G.y //kSk 6 siempre que ky v k < ı . Por tanto, para ky v k < ı tenemos: 2 kG.y / G.v /k 6 21 kG.y / G.v /k C kSkky v k Esto es:
kG.y / G.v /k 6 2kSkky v k
Con lo cual resulta finalmente que:
kG.y / G.v / S.y v /k 6 '.G.y //2kSk2ky v k Como, lKım '.G.y //D0, deducimos que G es diferenciable en v y además D G.v /D.D F.c//1 . ! Como esto es válido para todo v 2 V hemos probado que G es diferenciable en V y para todo punto y 2 V es: 8y 2 V D G.y / D .D F.G.y ///1 (3.28) Esta igualdad nos dice que D G D .D F ı G/1. Como la aplicación y 7! G.y / de V en U Rn es continua, la aplicación x 7! D F.x/ de U en GL.Rn / es continua y la aplicación T 7! T1 de GL.Rn / en L.Rn / es continua, deducimos que la aplicación y 7! D G.y / de V en L.Rn / y
v
es continua en V , es decir G es de clase C 1 en V .
Finalmente, probaremos que si F es de clase C k también lo es G . Acabamos de ver que esto es cierto para k D 1. Supongamos que es cierto para k y veamos que en tal caso también lo es para k C 1. Sea F D .f 1; f 2 ; : : : ; fn / de clase C kC1. Eso quiere decir que los campos escalares componentes de F , esto es las funciones f i , tienen derivadas parciales continuas de orden k C 1 y, por tanto, las funciones Dj f i tienen derivadas parciales continuas de orden k . Teniendo en cuenta la igualdad (3.28) y la regla de Cramer para calcular la inversa de una matriz cuadrada, deducimos que las funciones Dj Gi (donde G i , 1 6 i 6 n, son los campos escalares componentes de G ) son funciones racionales de las funciones Dj f i ı G cuyo denominador es siempre el determinante de J F .G.y // que sabemos no se anula para y 2 V . Como, por la hipótesis de inducción, G es de clase C k , las funciones Dj f i ı G son también de clase C k y, siendo las funciones racionales de clase C 1, concluimos, por la regla de la cadena, que las funciones Dj Gi son de clase C k en V , es decir G es de clase C kC1 en V .
Se dice que una aplicación entre espacios métricos es una aplicación abierta si la imagen por dicha aplicación de cualquier abierto es un abierto. Sean U , V abiertos en Rn . Un difeomorfismo de U sobre V es una biyección de U sobre V que es diferenciable en U y cuya inversa es diferenciable en V . Si ambas aplicaciones son de clase C k se dice que el difeomorfismo es de clase C k . Universidad de Granada Dpto. de Análisis Matemático
Prof. Javier Pérez Cálculo Diferencial en Rn
Teorema de la función inversa
70
Sea F W ! Rn un campo vectorial diferenciable en un abierto Rn . Se dice que F es un difeomorfismo local en si para todo punto x 2 hay un entorno abierto U x de dicho punto tal que F.U x/ D V x es un abierto y F es un difeomorfismo de U x sobre V x.
3.14 Corolario. Sea F W ! Rn un campo vectorial de clase C k en un abierto Rn y supongamos que detJ F .x/ ¤ 0 para todo x 2 . Entonces F es un difeomorfismo local de clase k C
en y es una aplicación abierta.
Si, además, F./.
F es inyectivo en entonces es un difeomorfismo de clase
k C
de sobre
Un difeomorfismo F de un abierto U Rn sobre un abierto V Rn suele interpretarse como una función que define nuevas coordenadas en V ; de hecho, a veces se llama a los difeomorfismos cambios de variables o cambios de coordenadas: las coordenadas según F de un punto y 2 V son las coordenadas cartesianas del punto x 2 U tal que F.x/ D y .
3.15 Ejemplo (Coordenadas polares en el plano) . El campo vectorial F W RC ; Œ! R 2 dado por F.;#/ D . cos #; sen #/, es inyectivo en D RC ;Œ. Como F es de clase 1 1 C (es, de hecho, de clase C ) y det J F .;#/ D > 0 , concluimos que F es un difeomorfismo de clase C 1 de sobre F./ D R2 nf.x; 0/ W x 6 0g. Dado .x; y/ 2 F./ el par .;# / 2 tal que F.;#/ D .x; y/ son las coordenadas polares de .x; y/. Una consecuencia interesante del teorema de la función inversa es la siguiente.
3.16 Proposición. Sea F W ! Rm una función de clase C 1 en un abierto Rn . Supongamos que la derivada de F en todo punto de es sobreyectiva. Entonces F es una aplicación abierta.
Demostración . Sean a 2 y fe1 ; e2 ; : : : ; em g la base canónica de Rm . Como D F .a/ W Rn ! Rm es sobreyectiva, existen vectores ui 2 Rn tales que D F .a/.ui / D ei , 1 6 i 6 m. Definamos S W Rm ! Rn por: m
S
m
X ! X D i ei
iD1
i ui
iD1
Claramente, S es lineal y DF .a/ ı S .ei / D e i , 1 6 i 6 m. Por tanto D F .a/ ı S D IRm (la identidad en Rm ). Definamos S a D a ı S donde a W Rn ! Rn es la traslación de vector a: a .x/ D a C x para todo x 2 Rn . Como a D Sa .0/ 2 , el conjunto W a D S1 ./ es un entorno a m abierto de 0 en R . Podemos aplicar ahora el teorema de la función inversa a la función Ga W W a ! Rm definida por G a .x/ D F ı Sa.x/ para todo x 2 W a . Claramente G a es de clase C 1 en W a y D Ga .0/ D D F .a/ ı S D IRm . Obtenemos así que hay un entorno abierto V a W a de 0 tal que Ga es un difeomorfismo de V a sobre el abierto Ga .V a/. .U / \ V a es un entorno Sea U un conjunto abierto y sea a 2 U . El conjunto U a D S1 a abierto de 0 contenido en V a, por lo que Ga U a / es abierto y Ga .0/ D F .a/ 2 Ga U a /. Como: .U / Ga U a/ Ga S1 a
D
ı S / S1.U / D F
.F
a
a
.U / Sa S1 a
F .U /
Se sigue que F .a/ es un punto interior de F .U /. Luego todos los puntos de F .U / son interiores y F .U / es abierto.
Universidad de Granada Dpto. de Análisis Matemático
Prof. Javier Pérez Cálculo Diferencial en Rn
Teorema de la función implícita
71
3.3. Teorema de la función implícita Sea f .x; y/ una función de dos variables con derivadas parciales de primer orden continuas y consideremos la ecuación f .x; y/ D 0 . Las soluciones de dicha ecuación representan una curva en el plano. Bueno, hablando con propiedad pueden representar algo más general que una curva. Para que te convenzas de ello basta que consideres la ecuación
D .x2 C y 2 1/.2.x 1/2 C 3.y 2/2 1/.y x 2/ D 0 la función f se anula en los puntos de la circunferencia x 2 C y 2 D 1, de la parábola y D x 2 y de la elipse 2.x 1/2 C 3.y 2/2 D 1. Por tanto la ecuación f .x; y/ D 0 representa la unión f .x; y/
de todas esas curvas.
Figura 3.1. Conjunto dado por f .x; y/ D 0
Ese conjunto (ver figura (3.1)) no es exactamente una curva pero localmente se parece a una curva. La palabra “localmente” quiere decir que si fijamos un punto .a; b/ tal que f .a; b/ D 0 entonces hay una bola abierta B..a; b/; r / tal que el corte de dicha bola con el conjunto de puntos M D f.x; y/ W f .x; y/ D 0g es una curva. De hecho, no es cierto que la condición anterior se verifique para todos los puntos .a; b/ tales que f .a; b/ D 0. Dicha condición falla en los puntos donde se cortan dos de las curvas cuya unión forma M , pues es claro que en dichos puntos el conjunto M no parece localmente una curva. Pues bien, en dichos puntos se anula el vector gradiente de f y en ellos la recta tangente no está definida. Este ejemplo te ayudará a entender lo que sigue. Volvamos al caso general de una función de dos variables f .x; y/ con derivadas parciales continuas de primer orden. Consideremos ahora la ecuación f .x; y/ D 0 desde otro punto de vista. Intuitivamente, una ecuación es una condición que debe ligar a una de las variables, es decir, que si en la igualdad f .x; y/ D 0 se fija un valor de x entonces el valor de y queda determinado de manera única por dicho valor de x . A veces esto es verdad como en el siguiente ejemplo. Consideremos f .x; y/ D y 3 C y ex C sen x
Fijado un valor de x la ecuación f .x; y/ D 0 es un polinomio de tercer grado en y que tiene una única solución real pues su derivada respecto de y es 3y 2 C ex que no se anula. Es decir, en este caso es cierto que la igualdad y3
Universidad de Granada Dpto. de Análisis Matemático
C y ex C sen x D 0
(3.29)
Prof. Javier Pérez Cálculo Diferencial en Rn
Teorema de la función implícita
72
define de manera única a y como función de x , en el sentido de que fijado un valor de x , hay un único y D ' .x/ que verifica dicha igualdad, esto es, la función '.x/ está definida por la condición: '.x/3 C '.x/ ex C sen x D 0 (3.30) Se dice que la función ' está implícitamente definida por la igualdad (3.29). Puedes calcular con wxMaxima el valor de dicha función y comprobarás que es bastante complicada. El hecho es que la mejor forma de trabajar con la función ' es la igualdad (3.30) que la define. Por ejemplo, si queremos calcular la derivada de ' en un punto basta con que derivemos dicha igualdad para obtener 3' 0 .x/'.x/2
C ' 0.x/ ex C'.x/ ex C cos x D 0
lo que permite calcular ' 0 .x/ en función de ' .x/. En general, no es cierto que una igualdad de la forma f .x; y/ D 0 permita despejar una variable en función de la otra. Para convencerte, considera el primer ejemplo que pusimos. Ni tan siquiera una igualdad tan sencilla como x 2 C y 2 1 D 0 permite despejar una variable como función de la otra pues es claro que para cada valor que fijemos de una variable (comprendido entre -1 y 1) hay dos posibles valores de la otra que verifican dicha igualdad. Relacionemos ahora los dos puntos de vista que hemos considerado. Pongamos
n D
2
2 W f .x; y/
.x; y/ R
o D 0
Si la igualdad f .x; y/ D 0 permitiera despejar y en función de x , es decir, definiera una función y D '.x/ por la condición f .x; y/
D 0 ” y D '.x/
entonces se tendría que (llamando I al intervalo donde está definida ' )
n D
2
2 W f .x; y/
.x; y/ R
o D Df 0
W 2 g
.x; '.x// x I
es decir, el conjunto sería la gráfica de ' , que, como sabemos, es un tipo muy particular de curva. Pero ya hemos visto que el conjunto puede ser una “curva” mucho más general que la gráfica de una función. Pero incluso en este caso, dicha “curva” es localmente, excepto en los puntos donde se anula el gradiente, una gráfica de una función. Las consideraciones anteriores se pueden llevar al caso de una función de tres variables f .x; y; z/ considerando ahora la “superficie” definida por la ecuación f .x; y; z/ D 0. La pregunta ahora es si fijados un valor de x y otro de y queda determinado de manera única un valor de z D '.x; y/ que verifica dicha ecuación. En caso afirmativo tendríamos que la superficie de ecuación f .x; y; z/ D 0 coincidiría con la gráfica de ' . Ya puedes suponer que esto no es cierto en general pues la mayoría de las “superficies” no son gráficas de funciones. Consideremos ahora el caso general en que tenemos una función vectorial F W ! Rm de clase C 1 definida en un abierto RnCm . Estamos interesados en estudiar el conjunto donde Universidad de Granada Dpto. de Análisis Matemático
Prof. Javier Pérez Cálculo Diferencial en Rn
Teorema de la función implícita
73
se anula F, es decir, el conjunto M Dfx 2 W F.x/ D 0g. Dicho conjunto está determinado por el sistema de m ecuaciones: f 1 .x1 ; x2 ; : : : ; xn ; xnC1 ; : : : ; xnCm / f 2 .x1 ; x2 ; : : : ; xn ; xnC1 ; : : : ; xnCm / : : :: : :: : :: : : : :: :: : :: :: : :: : : : :: f m .x1 ; x2 ; : : : ; xn ; ; xnC1 ; : : : ; xnCm /
D D D
0 0
: ::
(3.31)
0
En Mecánica estas ecuaciones se interpretan como “ligaduras” que deben satisfacer las variables. Tenemos en total m ligaduras que, intuitivamente, deben “atar” a m de las variables en el sentido de que el valor de las mismas quede determinado por el valor de las otras n . Por comodidad, consideremos las últimas m variables. La pregunta que nos hacemos es si el sistema de ecuaciones (3.31) permite “despejar” .xnC1 ; xnC2 ; : : : ; xnCm / en función de .x1; x2 ; : : : ; xn /. Para responder a esta pregunta lo primero que hacemos, siguiendo una estrategia ya conocida, es localizar el problema. La siguiente notación es cómoda: si x D .x1; x2 ; : : : ; xn / 2 Rn y y D .y1 ; y2 ; : : : ; ym / 2 Rm , pondremos . x; y / D .x1 ; x2 ; : : : ; xn ; y1 ; y2 ; : : : ; ym / 2 RnCm , es decir, identificamos RnCm con Rn Rm . Supongamos que . a; b/ 2 satisface el sistema de ecuaciones (3.31), es decir F.a; b/ D 0. Nos preguntamos si hay entornos abiertos, A Rn de a y B Rm de b, tales que para todo x 2 A existe un único y 2 B tal que .x; y / sea solución de (3.31). Cuando esto es así, existe una función ' W A ! B que a cada x 2 A hace corresponder el único y 2 B que verifica que es solución de (3.31), es decir, F .x; '. x// D 0. En tal caso, también se tiene que M \ .A B/ D f.x; '. x// W x 2 Ag, es decir, el conjunto M \ .A B/ es la gráfica de la función ' . El siguiente paso es linealizar el problema, esto es, sustituir la función F por su derivada en .a; b/ y trasladar el problema que nos ocupa a dicha derivada. El sistema de ecuaciones lineales que resulta escrito en forma matricial es: J F .a; b/ .x; y /t
D0 Se trata de un sistema de m ecuaciones lineales con nC m incógnitas y lo que queremos saber es .
si dicho sistema permite despejar y en función de x . La respuesta es conocida: basta exigir que el determinante de las últimas m columnas de la matriz J F .a; b/ sea distinto de cero. Llamando J 1 la matriz de las primeras n columnas de J F .a; b/ y J 2 a la matriz de las últimas m columnas tenemos que: J F .a; b/ .x; y /t .
D J 1 xt C J 2 y t D 0 ” y t D .J 21 J 1/ xt .
.
.
.
Después de lo que antecede, el teorema que sigue es lo que se espera.
3.17 Teorema ( de la función implícita). Sea F D .f 1 ; f 2 ; : : : ; fm / W G ! vectorial de clase C 1 definida en un abierto
det DnCj f i .a; b/
Entonces existen un entorno abierto W ' U Rm de clase C 1 tal que:
W !
Rm una función RnCm . Sea . a; b/ tal que F.a; b/ 0 y
2
16i 6m 16j 6m
D
¤0
de .a; b/ , un entorno abierto U de a y una función
f.x; y / 2 W W F.x; y / D 0g D f.x; '.x// W x 2 U g Universidad de Granada Dpto. de Análisis Matemático
Prof. Javier Pérez Cálculo Diferencial en Rn
Teorema de la función implícita
74
Además, para todo x U se verifica que
2
det DnCj f i .x; '.x//
y J ' .x/
D
16i 6m 16j 6m
¤0
1
DnCj f i .x; '.x//
16i 6m 16j 6m
.
Dj f i .x; '.x//
Si F es de clase C k entonces ' es de clase C k .
16i 6m 16j 6n
Demostración . Definimos el campo vectorial H W ! RnCm por:
D .x; F.x; y // 8.x; y / 2 : H es de clase C 1 y su matriz jacobiana en . x; y / 2 viene dada por: J H .x; y /
D
H .x; y /
Inn D1 f 1 .x; y / : :: : :: :: :: :: D1 f m .x; y /
::: : :: :::
0nm Dn f 1 .x; y / DnC1 f 1 .x; y / : :: :: :: :: :: : : :: :: :: :: : :: Dn f m .x; y / DnC1 f m .x; y /
::: : :: :::
DnCm f 1 .x; y / : :: : :: :: :: :: DnCm f m .x; y /
Tenemos que: detJ H .x; y / D det DnCj f i .x; y /
16i 6m 16j 6m
Por lo que det J H .a; b/ ¤ 0. Podemos aplicar el teorema de la función inversa al campo vectorial H para obtener que hay un entorno abierto de H.a; b/D.a; 0/, que podemos tomarlo de la forma U V donde U Rn es un entorno abierto de a y V Rm es un entorno abierto de 0 , y un entorno abierto de . a; b/, W , tal que H es un difeomorfismo de clase C 1 de W sobre U V . Sea G W U V ! W RnCm el difeomorfismo inverso de HjW . Teniendo en cuenta la definición de H, claramente G es de la forma:
G.x; z/ D .x; ˛.x; z//
8.x; z/ 2 U V
donde ˛ W U V ! Rm es una función vectorial de clase
1 C . Definimos '
W U ! Rm por:
D ˛.x; 0/ 8x 2 U : La función ' es de clase C 1 . Para todo x 2 U se tiene que .x; 0/ 2 U V y G.x; 0/ D .x; '.x//. Por lo que .x; 0/ D H.G.x; 0// D H.x; '.x// D .x; F.x; '.x//, luego F .x; '.x// D 0. Observa que para todo x 2 U se tiene que . x; '.x// 2 W . Luego hemos probado que: f.x; y / 2 W W F.x; y / D 0g f.x; '.x// W x 2 U g Supongamos ahora que . x; y / 2 W y F.x; y / D 0. Entonces H.x; y / D .x; 0/ 2 U V , luego x 2 U y, por tanto, . x; '.x// D G.x; 0/ 2 W . Puesto que H.x; '.x// D .x; 0/, y H es inyectiva en W ha de ser y D '.x/. Hemos probado así que: f.x; y / 2 W W F.x; y / D 0g f.x; '.x// W x 2 U g : '.x/
Universidad de Granada Dpto. de Análisis Matemático
Prof. Javier Pérez Cálculo Diferencial en Rn
Teorema de la función implícita
75
Queda probado que:
f.x; y / 2 W W F.x; y / D 0g D f.x; '.x// W x 2 U g En particular, se tiene que ' .a/ D b. Por el teorema de la inversa, tenemos que para todo .x; y / 2 W es detJ H.x; y / ¤ 0 ; en particular, para todo x 2 U es: detJ H .x; '.x// D det DnCj f i .x; '.x// 1 i m ¤ 0
6 6
16j 6m
Para cada i D 1; 2; : : : ; n y para todo x 2 U se verifica que: f i .x; '.x//
D0
Derivando esta identidad respecto a la variable k -ésima por la regla de la cadena, poniendo ' D .'1 ; '2 ; : : : ; 'm /, tenemos que: m
Dk f i .x; '.x//
X C
DnCj f i .x; '.x//Dk 'j .x/
j D1
D0
Para cada valor fijo de x 2 U y de k D 1; 2; : : : ; n éste es un sistema lineal de m ecuaciones y m incógnitas: las D k 'j .x/ para 1 6 j 6 m. El determinante de la matriz del sistema es: det DnCj f i .x; '.x//
16i 6m 16j 6m
¤0
Y, por tanto, dicho sistema tiene solución única que viene dada por:
Dk '1 .x/ :::::: :::::: Dk 'm .x/
˘ D
DnC1 f 1 .x; '.x// : : : : : : :: :: : :: :: ::: : : : :: :: : :: :: ::: DnC1 f m .x; '.x// : : :
DnCm f 1 .x; '.x// : :: : : : :: :: : : : :: : : : :: :: : : DnCm f m .x; '.x//
1
˘ .
Dk f 1 .x; '.x// :::::: :::::: Dk f m .x; '.x//
˘
Tenemos en total n sistemas de ecuaciones, uno para cada valor de k D 1; 2; : : : ; n. Cada uno de ellos tiene como solución la columna k -ésima de la matriz jacobiana J ' .x/. Las soluciones de todos ellos podemos escribirlas matricialmente en la forma: J ' .x/
D
1
DnCj f i .x; '.x//
16i 6m 16j 6m
.
Dj f i .x; '.x//
16i 6m 16j 6n
Finalmente, si la función F es de clase C k esta propiedad se traslada a las funciones H, G, ˛ y, por tanto, a la función implícita ' .
3.18 Observación. Con las notaciones usadas en la demostración del teorema, como W es un entorno abierto de .a; b/, hay entornos abiertos A1 Rn de a y B Rm de b tales que A1 B W . Como ' W U ! Rm es continua y U es abierto, el conjunto A D '1 .B/ \ A1 es un entorno abierto de a. Para todo x 2 A tenemos que '.x/ 2 B por lo que f.x; '.x/ W x 2 Ag A B . Además A B W . Deducimos fácilmente que: .A
B/ \ F1.0/ D f.x; y / 2 A B W F.x; y / D 0g D f.x; '.x/ W x 2 Ag
Universidad de Granada Dpto. de Análisis Matemático
Prof. Javier Pérez Cálculo Diferencial en Rn
Teorema de la función implícita
76
En la práctica el teorema de la función implícita se aplica en la forma que te explico en los siguientes ejemplos.
3.19 Ejemplo. Comprobar que la ecuación xyz
C sen.z 6/ 2.x C y C x 2y 2 / D 0
define a z como función implícita de .x; y/ en un entorno de .1; 1/, con z.1; 1/ D 6. Comprobar que .1; 1/ es un punto crítico de la función z D z.x; y/.
Solución. Pongamos f .x; y; z/ D xyz C sen.z 6/ 2.x C y C x 2 y 2 / que tiene derivadas @f
@f
parciales continuas de todo orden. Tenemos que D xy C cos.z 6/, por lo que .1; 1; 6/ D @z @z 2 ¤ 0. Como, además, f .1; 1; 6/ D 0, el teorema de la función implícita garantiza que hay una función con derivadas parciales continuas, .x; y/ 7! z.x; y/, definida en un entorno, U , de .1; 1/ tal que z .1; 1/ D 6, y f .x; y; z.x; y//
D 0 para todo
2
.x; y/ U :
Derivando esta identidad tenemos que: @f
@f @z @z C D yz 2.1 C 2xy 2 / C .xy C cos.z 6// D 0 @x @z @x @x @f @f @z @z xz 2.1 C 2x 2 y/ C .xy C cos.z 6// D 0 C D @y @z @y @y
Donde las derivadas parciales de la función implícita z D z.x; y/ están calculadas en un punto .x; y/ 2 U y las de f están calculadas en el punto .x; y; z.x; y//. Haciendo x D y D 1,
D z.1; 1/ D 6, en las igualdades anteriores, se obtiene que .1; 1/ es un punto crítico de z D z.x; y/. z
@z
@x
.1; 1/
@z D @y .1; 1/ D 0, esto es,
El ejemplo anterior es todavía demasiado explícito, nos dice muy claramente lo que hay que hacer. Lo más frecuente es que nos encontremos con ejercicios como el siguiente.
3.20 Ejemplo. Sabiendo que y cos.xz /
Calcular
@z @x
C x 3 ez y z C 1 D 0
.x; y/ y particularizar para el punto .x; y/
(3.32)
D .0; 0/.
Solución. En un ejercicio como este lo más fácil es que en la igualdad (3.32) sustituyas mentalmente z D z.x; y/ y la veas como y cos xz .x; y/
C
x 3 ez.x;y/y
z.x; y/ C 1 D 0
(3.33)
es decir, supones que has calculado para valores de x e y dados la solución respecto a z de la igualdad (3.32). Esta solución (que de hecho no es posible expresar de forma explícita, esto es, que no puede calcularse) la representamos por z D z.x; y/ y es la función implícita definida por la igualdad (3.32) (el teorema de la función implícita que es un teorema de existencia garantiza Universidad de Granada Dpto. de Análisis Matemático
Prof. Javier Pérez Cálculo Diferencial en Rn
Teorema de la función implícita
77
que dicha función existe cuando se cumplen las hipótesis del mismo). Ahora derivamos en la igualdad (3.33) respecto a x para obtener
y sen
xz .x; y/
de donde
z.x; y/
@z @x
C
.x; y/
@z x .x; y/ @x
D
C
3x 2 ez.x;y/y
Cx y @z .x; y/ e @x
yz .x; y/ sen xz .x; y/
x 3 y ez.x;y/y
3
z.x;y/y
@z @x .x; y/ D 0
3x 2 ez.x;y/y
x y sen.xz .x; y// 1
Naturalmente, esta igualdad tiene sentido siempre que el denominador de la fracción sea distinto de cero. Puedes comprobar que si llamas f .x; y; z/ D y cos.xz / C x 3 ezy z C 1 entonces la igualdad anterior es precisamente
@f .x; y; z/ @x @f .x; y; z/ @z
calculada en el punto .x; y; z.x; y//. Para .x; y/ D .0; 0/ se tiene que z.0; 0/ viene dado por la ecuación que se obtiene haciendo x D 0 e y D 0 en la igualdad (3.32) de donde se sigue z.0; 0/ D 1. Además @f .0; 0; z.0; 0// @z
.0; 0; 1/ D 1 ¤ 0 D @f @z
Por lo que @z @x
.0; 0/
D 01 D 0
Ejercicios propuestos 126. Particulariza el teorema de la función inversa: a) Para campos escalares de dos y de tres variables. b) Para funciones vectoriales de R3 en R2 y para funciones vectoriales de R4 en R2 . 127. Calcular las derivadas parciales de primer orden de la función z D z.x; y/ definida implícitamente por yz 4 C x 2 z 3 e xyz D0. Particularizar para el punto .x; y/ D .1; 0/. 128. Calcular las derivadas parciales de primer orden de la función z D z.x; y/ definida implícitamente por z 3 C z ex C cos y D 0. 129. Calcular las derivadas parciales de primer orden de la función z D z.x; y/ dada implícitamente por 3x 2y 2 C 2z 2 xy 2zx 3 C 4zy 3 4 D 0, en el punto .2; 1/ siendo z.2; 1/ D 2.
Universidad de Granada Dpto. de Análisis Matemático
Prof. Javier Pérez Cálculo Diferencial en Rn
Teorema de la función implícita
78
130. Supongamos que la igualdad z2
yCz
g.t /dt C h.t /dt D 0
xy
3xCy
donde g y h son funciones reales derivables, define a z como función implícita de x ; y . Calcular las derivadas parciales de primer orden de z D z.x; y/.
131. Sabiendo que g es una función continua con g.0/ D 3, justifica que xy
g.t / dt C x y D 0 2
x 2 1
define a y como función implícita de x en un entorno del punto .1; 0/. Calcula y 0 .1/.
132. Supongamos que la igualdad F .x; y; z/ D 0 determina implícitamente funciones diferenciables x D x.y; z/, y D y.x; z/, z D z.x; y/. Probar que
@x @y @z
@y @z @x
D 1.
133. Calcular la derivada de la función y D y.x/ definida implícitamente por xy
134.
135.
C 3x2 2y 2 2y D 0 Particularizar para x D 1 sabiendo que y.1/ D 1. Calcular la derivada de la función y D y.x/ definida implícitamente por y log.x 2 C y 2 / 2x y D 0 Particularizar para x D 0 sabiendo que y.0/ D 1. Sea F W ! R2 un campo vectorial de clase C 1 en un abierto R2 cuya derivada no se anula en . Supongamos que F satisface las ecuaciones de Cauchy-Riemann en , es decir, que para todo punto .x; y/ 2 se verifican las igualdades: @u @x
.x; y/
@v .x; y/; D @y
@u @y
.x; y/
@v .x; y/ D @x
Prueba que F es un difeomorfismo local en . Si, además, F es inyectiva, prueba que su inversa también satisface las ecuaciones de Cauchy-Riemann.
136. Prueba que existe un campo vectorial F D .f 1 ; f 2 / de clase C 1 definido en un abierto U R2 tal que F.1; 1/ D .1; 1/ y para todo .x; y/ 2 U se verifican las igualdades: x2
C y 2 2f 1.x; y/f 2.x; y/ D 0 x 3 C y 3 C f 1 .x; y/3 f 2 .x; y/3 D 0 Calcula D F.1; 1/ y
@2 f 1 @x 2
Universidad de Granada Dpto. de Análisis Matemático
.1; 1/.
Prof. Javier Pérez Cálculo Diferencial en Rn
Teorema de la función implícita
79
137. Justifica la existencia de una función ' de clase C 1 definida en un entorno U de cero en R con valores en R tal que:
'.x/ ex Cx e'.x/ D0
8x 2 U Calcula el polinomio de Taylor de segundo orden de ' en x D 0. 1
138. Probar que las ecuaciones: xy 5
C yu5 C zv5 D x 5 y C y 5 u C z 5 v D
1 1
definen a u y a v como funciones implícitas de .x; y; z/ en un entorno de .0; 1; 1; 1; 0/. Calcula
@u @y
.0; 1; 1/ y
@2 u @y 2
.0; 1; 1/.
139. Justifica la existencia de abiertos difeomorfos A y B de R2 verificando que .1; 1/ 2 A, .0; 1/ 2 B y que para cada .x; y/ 2 A existe un único .u; v/ 2 B tal que: x eu
Cy ev D 1 C e u ex Cv ey D e
Calcula la matriz jacobiana en .1; 1/ de la aplicación .x; y/ 7! .u.x; y/;v.x; y//. Calcula la matriz jacobiana en .0; 1/ de la aplicación .u; v/ 7! .x.u;v/; y.u;v//.
140. Sea z D z.x; y/ la función dada implícitamente por 3x 2 y 2
Calcula
@2 z @x@y
C 2z2 xy 2zx 3 C 4zy 3 4 D 0:
en el punto .2; 1/ siendo z .2; 1/ D 2.
141. Comprobar que la ecuación xyz
C sen.z 6/ 2.x C y C x 2y 2/ D 0
define a z como función implícita de .x; y/ en un entorno de .1; 1/, con z.1; 1/ D 6 . Comprobar que .1; 1/ es un punto crítico de la función z D z.x; y/ y estudiar si se trata de un máximo relativo, mínimo relativo o punto de silla.
142. Sea z D z.x; y/ la función dada implícitamente por yz 4 C x 2 z 3 exyz D0. Calcular @2 z
@x@y
en el punto .x; y/ D .1; 0/.
143. Calcular y clasificar los puntos críticos de la función z D z.x; y/ definida implícitamente por la igualdad 2x 2 C 2y 2 C z 2 C 8xz z C 8 D 0. 144. Probar que es posible despejar de manera diferenciable u y v en función de x; y; z en las ecuaciones:
C xz u C yv2 D 3 C 2xv u2v2 D 2 @u @v En un entorno de .x; y; z/ D .1; 1; 1/ y .u; v/ D .1; 1/. Calcula .1; 1; 1/ y .1; 1; 1/. @z @z
Universidad de Granada Dpto. de Análisis Matemático
xy 2 u3 xy
Prof. Javier Pérez Cálculo Diferencial en Rn
Teorema de la función implícita
80
145. Prueba que existen dos campos escalares u y v definidos en un entorno abierto U del punto .1; 1/ verificando que u.1; 1/ D v.1; 1/ D .1; 1/ y para todo .x; y/ 2 U :
x 3 y 3 u.x; y/3 v.x; y/3 x 2 y 2 2u.x; y/v.x; y/
C C C
D D
0 0
! R un campo escalar de clase C 2 definido en un entorno abierto V de @2 h D .1; 1/ y definamos h.x; y/ f .u.x; y/;v.x; y//. Calcula .1; 1/ y sabiendo que @x@y D1 f .1; 1/ D D2 f .1; 1/ D D12 f .1; 1/ D 2. Sea f W R3 ! R3 el campo vectorial definido por: f .x; y; z/ D .x xy; xy xyz; xyz/ Sea D .x; y; z/ 2 R3 W xy ¤ 0 . Justifica que es abierto y que f es un difeomorfismo Sea f W V
146.
˚
de sobre f ./. Calcula la matriz jacobiana del difeomorfismo inverso en el punto .0; 0; 1/.
147. Prueba que hay un abierto U R2 con .1; 0/ 2 U y un difeomorfismo f D .f 1; f 2 / de U sobre V D f .U / que verifica f .1; 0/ D .0; 1/ y para todo .x; y/ 2 U :
x 2 f 1 .x; y/ 2yf 2 .x; y/2 xy yf 1 .x; y/2 3f 1 .x; y/f 2 .x; y/ f 2 .x; y/x 2
C
C
C C
D D
0 1
148. Justifica que las igualdades:
C eu cos v D C C eu sen v D
xx u yv xv yu
1 0
definen localmente a u y a v como funciones de x e y en un entorno de .0; 0/ en el cual se verifica que
@u
@v @u @v D y D . @x @y @y @x
149. Sea A D .x; y; z/ 2 R3 W x D y _ y D z _ z D x . Prueba que A es un conjunto cerrado y que la función f W ! R3 dada por:
˚
f .x; y; z/ D .x C y C z; xy C yz C zx ; xyz/ es un difeomorfismo local en D R3 nA. Sea h W R3 ! R dada por h.u; v ; w / D u C ev C sen w . Prueba que existe un abierto U R3 con .0; 1; 0/ 2 U y un campo escalar g W U ! R de clase C 1 en U con g.0; 1; 0/ D 1 C e1 tal que g ı f D h. Calcula r g .0; 1; 0/.
150. Prueba que las igualdades:
Cyu u2 D y cos v C x 2 u2 D x ev
0 1
definen a u y a v como funciones implícitas de x e y en un entorno del punto .2; 1/ siendo u.2; 1/ D 2, v.2; 1/ D 0. Calcula las derivadas parciales de segundo orden de u en .2; 1/. Universidad de Granada Dpto. de Análisis Matemático
Prof. Javier Pérez Cálculo Diferencial en Rn
Variedades diferenciables en Rn
81
151. Sea F W Rn ! Rn un campo vectorial de clase C 1 y supongamos que hay un número < 1 tal que para todo x 2 Rn es kD F.x/k 6 . Sea G W Rn ! Rn el campo vectorial definido para todo x 2 Rn por G.x/ D x C F.x/. Prueba que G es un difeomorfismo de Rn sobre Rn . 152. Sea F W Rn ! Rn un campo vectorial de clase C 1 y supongamos que para todos x; y 2 Rn se verifica que: kx y k 6 kF.x/ F .y /k
Prueba que D F .a/ es invertible en todo punto a 2 Rn , que F .Rn / es cerrado y que F es un difeomorfismo de Rn sobre Rn .
153. Sea f W R ! R una función de clase C 1 en R y supongamos que hay un número ˛ < 1 tal que jf 0 .t /j 6 ˛ para todo t 2 R. Definamos F W R2 ! R2 por F .x; y/
D .x C f .y/; y C f .x//
Prueba que F es un difeomorfismo de R2 sobre R2 .
3.4. Variedades diferenciables en Rn La palabra “variedad” se usa en matemáticas para describir espacios topológicos que localmente son parecidos al espacio Rk con la topología de la norma para algún k 2 N llamado la dimensión de la variedad. Como ya habrás estudiado en la asignatura de Topología I, el concepto de identidad matemática entre espacios topológicos es el de homeomorfismo: dos espacios topológicos son homeomorfos cuando entre ellos existe una biyección continua y con inversa continua. Una tal biyección se llama un homeomorfismo. Recuerda también que todo subcon junto no vacío, M , de un espacio métrico .E; d/ se considera siempre como espacio métrico con la distancia inducida, y también como espacio topológico con la topología dada por dicha distancia en la cual los abiertos son los abiertos relativos: intersecciones de los abiertos del espacio total E con M . Esta topología se llama la topología inducida en M o la topología relativa de M . Si A Rn y B Rk son conjuntos no vacíos, diremos que son homeomorfos si lo son como espacios topológicos con las respectivas topologías relativas, es decir, existe una biyección W A ! B continua cuya inversa también es continua. Una tal biyección, cuando existe, se llama un homeomorfismo de A sobre B .
Ejercicios propuestos 154. Prueba que todas las bolas abiertas en un espacio normado .X ; k k/ son homeomorfas a la bola abierta unidad, B.0; 1/, y que ésta es homeomorfa a X . 155. Prueba que la esfera euclídea unidad en R3 no es homeomorfa a R2 . Universidad de Granada Dpto. de Análisis Matemático
Prof. Javier Pérez Cálculo Diferencial en Rn
Variedades diferenciables en Rn
82
Consideremos, pues, un conjunto no vacío M Rn . Dicho conjunto es un espacio topológico y será localmente parecido al espacio Rk si para todo punto a 2 M existe un abierto relativo V a M que lo contiene que es homeomorfo a Rk . Como toda bola abierta en Rk es homeomorfa a Rk , y la restricción de un homeomorfismo sigue siendo un homeomorfismo, no se pierde nada, y se gana agilidad en la exposición, si en vez de exigir que V a sea homeomorfo a Rk exigimos que sea homeomorfo a un abierto U Rk . Supongamos, pues, que } W U ! V a es un homeomorfismo. Podemos interpretar que U es un mapa “plano” que representa V a y } es la aplicación que lleva cada punto de dicho mapa a su correspondiente en V a . Digamos que lo que hace la aplicación } es “curvar y pegar” el mapa plano U sobre un trozo V a de M . También queremos que nuestras “variedades” tengan un “espacio tangente” en cada punto que tenga la misma dimensión que la variedad. Para ello vamos a exigir que } sea al menos de clase C 1 y que su diferencial tenga rango k en todo punto. Naturalmente, V a será de la forma V a D W \ M donde W es un entorno abierto de a en Rn .
3.21 Definición. Un conjunto no vacío M Rn es una variedad (diferenciable de clase C q , 1 6 q 6 1) de dimensión k (1 6 k 6 n) si para todo punto a 2 M hay un entorno abierto W Rn de a, un abierto U Rk y una aplicación } W U ! Rn (de clase C q , 1 6 q 6 1) tal que rango J } .x/ D k para todo x 2 U y } es un homeomorfismo de U sobre W \ M . Diremos que un par .} ; U / en las condiciones de la definición anterior es una carta local de M en a. El par .} 1 ; }.U // se llama un sistema de coordenadas locales. Las coordenadas de un punto x 2 }.U / son las coordenadas cartesianas de } 1 .x/. Esta terminología no es universal y muchos textos llaman “carta local” a lo que hemos llamado “sistema de coordenadas locales”. Tanto da, pues conocer .}; U / es lo mismo que conocer .} 1 ; }.U //. También se dice que el par .} ; U / o la aplicación } es una parametrización local de M en a. Una parametrización .}; U / tal que la aplicación } 1 sea una proyección en k -coordenadas, es decir, } 1 .x1 ; x2 ; : : : ; xn / D .xi1 ; xi2 ; : : : ; xik / donde 1 6 i1 < i2 < : : : < ik 6 n, se llama una parametrización cartesiana. Las expresiones escritas entre paréntesis en la definición anterior se sobrentenderán en todo lo que sigue: nuestras variedades serán siempre de clase C q con 1 6 q 6 1.
Las variedades de dimensión 1 se suelen llamar curvas regulares y las de dimensión 2 superficies regulares. Las variedades de dimensión n 1 que están en Rn con n > 3 suelen llamarse hipersuperficies. Con el convenio R0 Df0g (espacio vectorial nulo), puede incluirse en la definición anterior el caso k D 0.
Ejercicios propuestos 156. Prueba que las variedades de dimensión n contenidas en Rn son los conjuntos abiertos de Rn . 157. Prueba que las variedades de dimensión 0 son conjuntos de puntos aislados. En lo que sigue, para no caer en los casos extremos considerados en los ejercicios anterioUniversidad de Granada Dpto. de Análisis Matemático
Prof. Javier Pérez Cálculo Diferencial en Rn
Variedades diferenciables en Rn
83
res, se supondrá que 1 6 k < n . En la definición 3.21 las variedades aparecen como subconjuntos de Rn que localmente son imagen homeomorfa de abiertos de Rk mediante homeomorfismos diferenciables de rango k . Podemos dar ejemplos sencillos de variedades.
3.22 Ejemplo. Sea } W U ! Rn una función vectorial de clase C q en un abierto U Rk tal que rango J } .x/ D k para todo x 2 U y } es un homeomorfismo de U sobre }.U /. Entonces }.U / es una variedad de clase C q y dimensión k . Se dice que dicha variedad está globalmente determinada por la parametrización } . En particular, si U Rk es abierto, F W U ! Rnk es una función de clase C q , entonces el conjunto M D f.x; F.x// W x 2 U g Rn es una variedad de clase C q y dimensión k . Observa que dicha variedad es la gráfica de F. El ejemplo anterior nos dice que la gráfica de toda función de clase C q es una variedad de clase C q . Veamos que toda variedad es localmente la gráfica de una función. Seguiremos aquí los convenios de notación que usamos en el teorema de la función implícita: identificaremos Rk Rnk con Rn . Una expresión de la forma . x; y / donde x 2 Rk , y 2 Rnk , se interpreta como un vector de Rn cuyas primeras k componentes son las de x y cuyas últimas n k componentes son las de y . Notaremos por 1 W Rn ! Rk la proyección en las primeras k coordenadas, esto es, 1 .x; y / D x. Notaremos por 2 W Rn ! Rnk la proyección en las últimas n k coordenadas, esto es, 1 .x; y / D y . Observa que estas proyecciones son aplicaciones abiertas.
3.23 Proposición. Sea M Rn . Las siguientes afirmaciones son equivalentes: a) M es una variedad de clase C q y dimensión k .
b) Para todo punto a M existen, previa una reordenación de las variables en Rn , entornos Rn de a , U Rk de 1 .a/ , y una función ' U Rnk de clase C q tal que abiertos W W M .x; '.x// x U .
\ D f
2
W 2 g
W !
Demostración . Probaremos que a/�b/. Sea a 2 M . Por hipótesis hay un entorno abierto W 1 Rn de a, un abierto Rk y una aplicación } W ! Rn de clase C q tal que rango J } .x/ D k para todo x 2 y } es un homeomorfismo de sobre W 1 \ M . Sea x 0 2 tal que }.x0/ D a. Como rango J } .x0 / D k , reordenando las variables si fuera preciso, podemos suponer que las primeras k filas de J } .x0 / D k tienen determinante distinto de cero. Esto es tanto como decir que la aplicación } 1 D 1 ı } W ! Rk , cuya matriz jacobiana está formada precisamente por las primeras k filas de J } .x0 /, verifica que det J }1 .x0 / ¤ 0. El teorema de la función inversa afirma que hay un entorno abierto V de x 0 tal que }1jV W V ! U es un difeomorfismo de clase C q de V sobre el abierto U D }1 .V /. Observa que } 1 .x0 / D 1 .}.x0// D 1 .a/ 2 U . El conjunto }.V / es un entorno abierto de a en W 1 \ M , luego existe un entorno abierto W 2 Rn de a tal que }.V /
D W 2 \ W 1 \ M D W \ M
donde W D W 2 \ W 1 es un entorno abierto de a. Universidad de Granada Dpto. de Análisis Matemático
Prof. Javier Pérez Cálculo Diferencial en Rn
Variedades diferenciables en Rn
84
1 1 W U ! V es el .x// donde } 1jV Sea h W U ! Rn la aplicación definida por h.x/ D }.}1jV difeomorfismo inverso de } 1jV . Tenemos que 1 .U // D }.V / D W \ M h.U / D }.}1jV
Pongamos } 2 D 2 ı } W V ! Rnk , con lo que para todo x 2 V es }.x/ D .}1 .x/; }2 .x//. Para todo x 2 U es: 1 1 1 1 .x// D .}1 .}1jV .x//;}2 .}1jV .x/// D .x; }2 .}1jV .x// h.x/ D }.}1jV 1 Definamos ' W U ! Rnk por ' .x/ D .}2 ı }1jV /.x/. Dicha aplicación es de clase C q y:
h.U / D fh.x/ W x 2 U g D f.x; '.x// W x 2 U g D W \ M
(3.34)
lo que concluye la demostración. La implicación b/ � a/ es un sencillo ejercicio.
Observa que la función h W U ! Rn antes definida, h .x/ D .x; '.x//, es una aplicación de clase C q y rango J h .x/ D k para todo x 2 U . Además h es una biyección continua de U sobre W \ M cuya biyección inversa, .x; '.x// 7! x, es la proyección en las primeras k coordenadas y, por tanto, continua. Esto es, el par .h; U / es una parametrización cartesiana local de M en a. Hemos probado así el siguiente resultado.
3.24 Corolario. Toda variedad tiene parametrizaciones cartesianas en el entorno de cada uno de sus puntos.
Podemos obtener otra interesante consecuencia de la demostración anterior. Pongamos G D 1 W \.U Rnk / que es un entorno abierto de a. Definamos gW G ! Rk por g.z/D}1jV .1 .z//. Observa que para todo z 2 G se tiene que 1 .z/ 2 U , por lo que la aplicación g está bien definida, es de clase C q y g.z/ 2 V . Por la igualdad (3.34) se tiene que W \ M G , lo que implica que W \ M G \ M , pero G \ M D W \ M \ .U Rnk / W \ M . Luego G \ M D W \ M . Como para todo z 2 W \ M se tiene que z D }. x/ donde x 2 V , resulta que: 1 1 .1 .}.x/// D }1jV .}1 .x// D x D } 1 .z/ g.z/ D g.}.x// D }1jV Por tanto las restricciones de g y de } 1 a G \ M D }.V / coinciden. Teniendo en cuenta que los puntos x0 2 V y a D }.x0 / 2 M pueden ser cualquier par de puntos en esa situación, hemos probado el siguiente resultado.
3.25 Proposición. Sea .};/ una parametrización local de una variedad M
Rn de dimensión k y clase C q . Entonces para todo x 2 existen un entorno abierto V de x , un entorno abierto G de }. x/ y una función g W G ! V Rk de clase C q tal que }.V / D G \ M 1 D gjG\M . y } jG\M Podemos probar ahora con facilidad un resultado que es muy importante cuando se estudia la integración en variedades.
3.26 Teorema. Sea M
Rn una variedad de dimensión k y clase C m , .P ; 1/ y . Q; 2 / cartas locales de M tales que P .1 / \ Q.2 / D W ¤ Ø. Sea U p D P 1 .W / y U q D Q1 .W /. Entonces la aplicación Q1 ı P jU es un difeomorfismo de clase C m de U p sobre U q . p
Universidad de Granada Dpto. de Análisis Matemático
Prof. Javier Pérez Cálculo Diferencial en Rn
Variedades diferenciables en Rn
85
Demostración . Observa que U p 1 y U q 2 son abiertos en Rk y la aplicación F D Q1 ı P jU p es una biyección de U p sobre U q . Sea u 2 U p . Entonces P .u/ 2 W Q.2 /, luego P .u/ D Q.x/ donde x 2 2 . Aplicamos la proposición anterior a la carta . Q; 2 / y al punto x para obtener que existen un entorno abierto V de x , un entorno abierto G de Q.x/ y 1 D gjG\M . una función g W G ! V Rk de clase C m tal que Q.V / D G \ M y QjG\M Sea U D P 1 .G \ M / \ U p 1 . Tenemos que u 2 U y U es abierto. Para todo z 2 U
es:
F.z/ D .Q1 ı P jU p /.z/ D Q1 .P .z// D porque P .z/ 2 G \ M D g.P .z// D .g ı P /.z/ Hemos probado que FjU D g ı P jU . La función g ı P jU es composición de funciones de clase m m m C y, por tanto, es de clase C en U , luego F también es de clase C en U . Hemos probado que F es localmente de clase C m en U p y teniendo en cuenta el carácter local de la derivada, deducimos que F es de clase C m en U p . Razonando de forma análoga, se prueba que la función F1 D P 1 ı QjU q es de clase C m en U q .
El teorema de la función implícita nos dice que el conjunto donde se anula un campo vectorial diferenciable de rango máximo es localmente una gráfica funcional y proporciona interesantes ejemplos de variedades.
3.27 Proposición. Sea FW ! Rnk una función vectorial de clase C q en un abierto Rn . Sea
D f.x; y / 2 W F.x; y / D 0g y supongamos que rangoJ F .x; y / D n k para todo . x; y / 2 M . Entonces se verifica que M es una variedad de clase C q y dimensión k . Se dice que la variedad M está determinada por F y que F.x; y / D 0 son las ecuaciones implícitas de la variedad. M
Demostración . Sea .a; b/ 2 M . Entonces F.a; b/ D 0 y rango J F .a; b/ D n k . No es restrictivo suponer que las últimas n k columnas de la matriz J F .a; b/ son linealmente independientes 3 . El teorema de la función implícita nos dice que hay un entorno abierto W de . a; b/, un entorno abierto U Rk de a y una aplicación de clase C q , ' W U ! Rnk tal que:
f.x; y / 2 W W F.x/ D 0g D W \ M D f.x; '.x// W x 2 U g Basta tener ahora en cuenta lo dicho antes del corolario 3.24.
Como el concepto de variedad es local, podemos localizar el resultado anterior y obtenemos otra caracterización de las variedades.
3.28 Proposición. Sea M Rn . Las siguientes afirmaciones son equivalentes: a) M es una variedad de clase C q y dimensión k .
3 Esto siempre puede conseguirse componiendo F con una conveniente permutación de ejes en Rn .
Universidad de Granada Dpto. de Análisis Matemático
Prof. Javier Pérez Cálculo Diferencial en Rn
Variedades diferenciables en Rn
86
b) Para todo punto a M existen, previa una reordenación de las variables en Rn , un entorno abierto W Rn de a y una función F W Rnk de clase C q tales que rango J F .x/ n k para todo x W y
D
2
W
!
2 W \ M D fx 2 W W F.x/ D 0g
Se dice que la función F determina localmente a M en a y también que F .x/ ecuaciones implícitas locales de M en a.
D 0 son las
Demostración . Probaremos que a/�b/. La forma de proceder es exactamente la misma que en la proposición 3.23. Sea a 2 M . Por hipótesis hay un abierto Rk , un entorno abierto W 1 Rn de a y una aplicación } W ! Rn de clase C q tal que rango J } .x/ D k para todo x 2 y } es un homeomorfismo de sobre W 1 \ M . Sea x 0 2 tal que }.x0 / D a. Como rango J } .x0 / D k , reordenando las variables si fuera preciso, podemos suponer que las primeras k filas de J } .x0 / D k tienen determinante distinto de cero. Esto es tanto como decir que la aplicación }1 D 1 ı } W ! Rk , cuya matriz jacobiana está formada precisamente por las primeras k filas de J } .x0 /, verifica que det J }1 .x0 / ¤ 0. El teorema de la función inversa afirma que hay un entorno abierto V de x 0 tal que }1jV W V ! U es un difeomorfismo de clase C q de V sobre el abierto U D }1 .V /. Observa que } 1 .x0 / D 1 .}.x0// D 1 .a/ 2 U . Definimos las funciones:
W ! Rn h W U Rnk ! Rnk
1 D } ı }1jV h D 2 2 ı ' ı 1 Dichas funciones son de clase C q . Como para todo x D .x1 ; x2 ; : : : ; xn / 2 U Rnk es h.x1; x2 ; : : : ; xn / D .xkC1 'kC1 .x1; : : : ; xk / ; : : : ; xn 'n .x1; : : : ; xk // deducimos que rango J h .x/ D n k para todo x 2 U Rnk . Además: h1 .0/ D x 2 Rn W x 2 U Rnk ^ 2 .x/ D 2 ı ' ı 1 .x/ D D x 2 Rn W 1.x/ D 1 .}.u// ^ 2.x/ D .2 ı '/.1.}.u/// para algún u 2 V D D x 2 Rn W 1.x/ D 1 .}.u// ^ 2.x/ D 2 .}.u// para algún u 2 V D D x 2 Rn W x D }.u/ para algún u 2 V D }.V / Como V es abierto y } es un homeomorfismo sobre su imagen, existe un abierto W 2 Rn tal que }.V / D W 2 \ W 1 \ M . Sea W D .U Rnk / \ W 1 \ W 2 y sea F D hjW la restricción de h a W . Tenemos que W es abierto y a 2 W , F W W ! Rnk es de clase C q , rango J F .x/ D n k y para todo x 2 W : W \ M D W \ .M \ W 1 \ W 2 / D W \ }.V / D fx 2 W W h.x/ D 0g D fx 2 W W F.x/ D 0g
' U
n ˚ ˚ ˚
'
o
Que b/�a/ es consecuencia de la proposición 3.27 y del carácter local del concepto de variedad.
Universidad de Granada Dpto. de Análisis Matemático
Prof. Javier Pérez Cálculo Diferencial en Rn
Espacios tangente y normal
87
Hemos visto que una variedad M Rn de dimensión k puede venir dada (localmente) por medio de sus ecuaciones paramétricas, es decir como la imagen de una parametrización, o por medio de sus ecuaciones implícitas, es decir como el conjunto de puntos donde se anula una función vectorial de n variables y n k componentes y de rango máximo. Es instructivo considerar los casos más simples de variedades: los subespacios afines de Rn . Como debes saber, un subespacio afín de dimensión k de Rn puede venir dado por sus ecuaciones implícitas:
8ˆ< ˆ:
C C C
CC CC CC
D D : :: D
a11 x1 a12 x2 a1n xn a21 x1 a22 x2 a2n xn : :: :: :: :: :: :: :: :: :: :: :: :: :: ank1 x1 ank2 x2 ank nxn
c1 c2 : :: :: : cnk
Dichas ecuaciones representan la intersección de n k hiperplanos y deben ser linealmente independientes, es decir, el rango de la matriz de los coeficientes del sistema debe ser igual a n k . Definiendo F W Rn ! Rnk por:
0X D@ n
F.x1; x2 ; : : : ; xn /
a1j xj
j D1
c1;
n
n
X
X
a2j xj
j D1
c2; : : : ;
j D1
ankj xj
cnk
1 A
Se tiene que F es de clase C 1 y su matriz jacobiana, que es la matriz de los coeficientes del sistema, tiene rango máximo n k , y el subespacio afín viene dado globalmente como el conjunto de los ceros de F. Si en la situación anterior suponemos que las últimas n k columnas de la matriz de los coeficientes del sistema son linealmente independientes, podemos resolver dicho sistema y su solución, poniendo A D .aij /16i 6nk , B D .aij / 16i 6nk y c D .c1; c2 ; : : : ; cnk /, está dada 16j 6k kC16j 6n por: .xkC1 ; xkC2 ; : : : ; xn /t D B1ct B1A:.x1 ; x2 ; : : : ; xk /t .
.
Poniendo
}.x1 ; x2 ; : : : ; xk /
D
x1 ; x2 ; : : : ; xk ; B
1 t
c
.
1
B
A .x1; x2 ; : : : ; xk /
.
t
.
se tiene que } W Rk ! Rn es una parametrización que determina globalmente al subespacio afín. Observa que la imagen de }. Rk / viene dada en función de los parámetros x 1; x2 ; : : : ; xk que toman valores reales cualesquiera. A la vista de estos ejemplos, queda claro que el concepto de variedad que hemos definido refleja una vez más la idea de “linealización” básica del cálculo diferencial.
3.5. Espacios tangente y normal En lo que sigue M Rn será una variedad de dimensión k ( 1 6 k < n ), clase C q y a 2 M un punto de M . Se dice que u 2 Rk es un vector tangente a M en a si existe una curva contenida en M , γW ı;ıŒ! M , que pasa por a, γ .0/ D a, y que tiene tangente en a con dirección u , γ 0 .0/ D u. Universidad de Granada Dpto. de Análisis Matemático
Prof. Javier Pérez Cálculo Diferencial en Rn
Espacios tangente y normal
88
Se define el espacio tangente a M en a como el conjunto de todos los vectores tangentes a M en a. Lo notaremos T M .a/ y veremos enseguida que es un subespacio vectorial de Rn . La variedad afín a C T M .a/ se llama variedad afín tangente a M en a (recta tangente si k D 1, plano tangente si k D 2, hiperplano tangente si k D n 1).
? .a/, El complemento ortogonal en Rn del espacio tangente a M en a , que notaremos T M se llama espacio normal a M en a: ? T M .a/
D x 2 Rn W x u D 0; 8u 2 T M .a/
˝ ˇˇ ˛
˚
? La variedad afín a C T M .a/ se llama variedad afín normal a M en a.
Cuando se conoce una parametrización o las ecuaciones implícitas que representan localmente a M en a es muy fácil calcular los espacios tangente y normal.
3.29 Proposición. Sean F .x/ D 0 las ecuaciones implícitas locales de M en a y .};/ una
D a. Entonces se verifica que: T M .a/ D Ker D F.a/ D Im D}.c/
parametrización local de M en a con }. c/
Demostración . Como las matrices jacobianas de F y } tienen respectivamente rango n k y k , los espacios vectoriales Ker D F.a/ e Im D}.c/ tiene dimensión k . Bastará, por ello, probar que: Im D}.c/ T M .a/ Ker D F.a/ Sea h 2 Im D}. c/ y sea v 2 Rk tal que D}.c/.v / D h. Sea ı > 0 tal que el segmento Œc ı v ; c C ı v y definamos γ W Œc ı v ; c C ı v ! M por γ.t / D }.c C t v /. Tenemos que γ .0/ D a y γ 0 .0/ D D}.c/.v / D h. Lo que prueba que h 2 T M .a/. Sea ahora u 2 T M .a/ y sea γW ı;ıŒ! M una curva con γ.0/ D a y γ 0 .0/ D u. Tomando ı suficientemente pequeño para que la curva γ esté contenida en la parte de M determinada por F se tiene que .t / D .F ı γ/.t / D 0 y por tanto 0 .0/ D D F.a/ u D 0, lo que prueba que u 2 Ker D F.a/. .
La proposición anterior tiene una simetría que merece destacar: si representamos localmente M en a de forma implícita por los ceros de una función F, entonces el espacio tangente T M .a/ viene dado por los ceros de D F.a/; y si representamos localmente M en a como la imagen de una parametrización } , entonces el espacio tangente T M .a/ viene dado como la imagen de D}.c/ siendo }.c/ D a.
3.30 Corolario. Sean F.x/ D 0 las ecuaciones implícitas locales de M en a y pongamos F D .F 1 ; F 2 ; : : : ; F nk / , y sea .}; / una parametrización local de M en a con }.c/ D a y
pongamos }
D .}1; }2; : : : ; }n /. Entonces se verifica que:
a) Los vectores fila de la matriz jacobiana F.a/ , es decir, los vectores Fi .a/ , 16i 6n k , ? forman una base de T M .a/ y la variedad afín tangente a M en a viene dada por:
r
a
C T M .a/ D x 2 Rn W r F i .a/ x a D 0; 1 6 i 6 n k
˚
ˇ ˛
˝
Naturalmente, el espacio afín ortogonal a M en a viene dado por: a
C T M ? .a/
8< D: 2
nk
x Rn x
Universidad de Granada Dpto. de Análisis Matemático
W DaC
X iD1
r
i Fi .a/ i
2R 1 6 i 6 n
9= ; k
(3.35)
(3.36)
Prof. Javier Pérez Cálculo Diferencial en Rn
Espacios tangente y normal
89
b) Los vectores columna de la matriz jacobiana J } .c/ , es decir, los vectores: Dj }. c/
D .Dj }1.c/; Dj }2.c/ ; : : : ; Dj }n.c//
1 6 j 6 k
son una base del espacio T M .a/ y la variedad afín ortogonal a M en a viene dada por: a
C T M ? .a/ D x 2 Rn W Dj }.c/ x a D 0; 1 6 j 6 k
˚
ˇˇ ˛
˝
Naturalmente, el espacio afín tangente a M en a viene dado por:
a
C T M .a/
8< D: 2
k
X
x Rn x
W DaC
j Dj }.c/ j
j D1
2R
1 6 j 6 k
(3.37)
9= ;
(3.38)
Observa que (3.35) y (3.37) son, respectivamente, las ecuaciones implícitas de la variedad afín tangente y del espacio afín normal, mientras que ( 3.36) y (3.38) son, respectivamente, las ecuaciones paramétricas del espacio afín normal y de la variedad afín tangente.
Ejercicios propuestos 158. Comprueba que las rectas y planos tangentes a curvas y superficies estudiados en la sección 2.2 son casos particulares de los resultados generales para variedades. 159. Interpreta en el contexto de la teoría de variedades los ejercicios 83 – 87. 160. Prueba que cada uno de los siguientes subconjuntos de R3 es una variedad y calcula el espacio afín tangente y normal en el punto P que se indica.
( D D ˚ ( D n D
x2
a) M
.x; y; z/R3
b) M
.x; y; z/R3 z
c) M
d) M
y2
W a2 C b2 C c 2 D 1 W D x2 y2
.x; y; z/R2
)
z2
;
;
D .a; b c/.
P
D .0; 0; 0/.
P
2
W C C D ^ C o p W C C D D RC x 2
.x; y; z/R3 .
x2
y2
y2
z2
2/2
1
x
z2
1 ;
1
2
P
y2
1
D4
)
; P .0; 0; 1/.
D
.3; 0; 0/.
161. Prueba que M D .x; y; z/ 2 R3 W x 2 C y 2 2z 2 D 0 ^ x C y C z D 1 es una curva regular y calcula su recta tangente en .1; 1; 1/.
˚ D ˚
.x; y; z/ 2 R3 W x 2 C y 2 C z 2 D 1 ^ z D xy es una curva regular y 162. Prueba que M calcula su recta tangente en .a; b; c/ 2 M .
163. Prueba que la esfera euclídea unidad en Rn : S n1
D x 2 Rn W kxk2 D 1
˚
es una variedad de dimensión n 1 y calcula el hiperplano tangente y la recta normal en un punto a 2 S n1 . Universidad de Granada Dpto. de Análisis Matemático
Prof. Javier Pérez Cálculo Diferencial en Rn
Extremos condicionados
90
164. Sean M Rn y N Rm variedades y f W U ! Rm una función de clase abierto U M tal tal que f .M / M / N . Prueba que D f .a/.T /.T M M .a// T N N .f .a//.
1 C
en un
. 0; 1; 1; 0/ del toro: 165. Calcula bases del espacio tangente y del espacio normal en el punto .0
n D D
M
4
.x1 ; x2 ; x3 ; x4 / R
x12
x22
x32
2 W C D 1 ^ C
o D
x42
1
3.6. Extremos Extremos condic condicion ionado adoss En la teoría de extremos relativos se supone que las variables pueden tomar valores en cualquier punto de un conjunto abierto, es decir, pueden “moverse libremente” en dicho conjunto. En muchos, por no decir que en la mayoría, de los problemas reales las variables no tienen tanta libertad y están obligadas a satisfacer ciertas condiciones que en Física suelen llamarse “ligaduras”. Por ejemplo, supongamos que un móvil se mueve en una curva dada por la . x ; y ; z / 2 la energía cinética del móvil vieintersección de dos superficies; para cada punto .x ne dada por una función conocida f .x ; y ; z / y queremos calcular los puntos de la trayectoria donde dicha energía es máxima o mínima. En esta situación las variables x ; y ; z no son libres . x ; y ; z / 2 . Otro ejemplo; supongamos que la tempesino que deben satisfacer la condición .x ratura en un punto .x ; y ; z / de la superficie terrestre viene dada por una función T . T .x ; y ; z / y queremos calcular los puntos de mayor y menor temperatura. Aquí las variables tampoco son libres pues deben verificar una condición de la forma x 2 C y 2 C z 2 D R2 donde R es el radio de la Tierra. Igualmente, en problemas de optimización de costes o beneficios las variables están siempre sometidas a restricciones que dependen de las condiciones de producción o del mercado. Es importante que comprendas la diferencia entre un problema de extremos relativos “libres” y un problema de extremos condicionados. Considera el siguiente ejemplo. 2
2
3.31 Ejemplo. La función f .x ; y / D xy ex Cy tiene un único punto crítico, el origen, que es un punto de silla. Por tanto dicha función no tiene extremos relativos en R2 . Supongamos que imponemos a las variables variables la condición condición x 2 C y 2 D 1 y queremos calcular el máximo valor de f .x ; y / cuando se verifique que x 2 C y 2 D 1. Fíjate en que el problema es completamente distinto. Ahora solamente nos interesan los valores que toma la función f .x ; y / en el conjunto K
n D
2
2
2 Wx Cy
.x ; y / R
2
o D 1
Sabemos que dicho conjunto es un conjunto compacto; además la función f es continua, por tanto podemos asegurar, de entrada, que tiene que haber algún punto .a en el cual la . a; b / 2 K en función f alcanza su mayor valor en K (y tiene que haber otro donde alcance su menor valor . x ; y / 2 K se en K ). Calcular dicho punto es, en este caso, muy sencillop pues para .x se tiene que f .x ; y / D e x y . Como para .x . x ; y / 2 K se se tiene que y D ˙ 1 x 2 y los valores negativos de f no nos interesan porque porque queremos calcular el mayor valor que toma en K , se sigue que
n p
o
max aK x ff .x ; y / W .x ; y / 2 K g D max aK x e x 1 x 2 W 0 6 x 6 1 Universidad de Granada Dpto. de Análisis Matemático
Prof. Javier Pérez Cálculo Diferencial en Rn
Extremos condicionados
91
Nuestro problema se ha convertido en algo que ya sabes hacer: calcular el máximo absoluto absoluto de p convertido 2 la función h. h .x / D e x 1 x para 0 6 x 6 1. Por cierto, por la desigualdad de las medias se tiene que: x
p D q 1
x2
x 2 .1
x2/ 6 21 p
Y la igualdad se da si, y sólo si, x 2 D 1 x 2, esto es en x D 1= 1 = 2 se alcanza el máximo absoluto de f en K . De hecho, tú has resuelto ejercicios de extremos condicionados aunque puede que no seas consciente de ello. Por ejemplo, seguro que alguna vez has resuelto el siguiente ejercicio.
3.32 Ejemplo. Entre todos los rectángulos cuyo perímetro es igual a 16 calcular el que tiene área máxima. Este ejercicio puedes plantearlo como sigue. Sea f .x ; y / D xy la función que da el área de un rectángulo cuyos lados tienen longitudes x e y . Se trata de calcular el máximo de f .x ; y / cuando las variables verifican la condición 2x C 2y D 16. Por tanto, es un problema de extremos condicionados. Seguro que ahora recuerdas algunos otros ejercicios parecidos a este que has hecho sin saber que estabas haciendo problemas de extremos condicionados. condicionados. La razón es clara: la condición que nos dan es tan sencilla que permite despejar una variable en función de la otra, y D 8 x , con lo que nuestra función se convierte en xy D x .8 x / y el problema queda reducido a calcular el mayor valor de x .8 x / cuando 0 6 x 6 8. Por cierto, por la desigualdad de las medias se tiene que:
p
x .8
x/ 6 4
Y la igualdad se da si, y sólo si, x D 8 x , esto es en x D 4 y el rectángulo de área máxima es el cuadrado. cuadrado. Los ejemplos anteriores ponen de manifiesto que los problemas de extremos condicionados de funciones de dos variables en los que puede utilizarse la condición que nos dan para despejar una variable en función de otra, se reducen fácilmente a problemas de extremos de funciones de una variable. Pero supongamos ahora que cambiamos la condición del ejemplo 1
por la siguiente: x
ex Cy C ey C sin.1 C xy/ xy / D 2
La cosa se complica porque ahora es imposible usar la condición impuesta para despejar una variable en función de la otra. Ahora sí tenemos un auténtico problema de extremos condicionados. Lo antes dicho para funciones de dos variables puedes generalizarlo para funciones de tres variables. Por ejemplo el problema de calcular las dimensiones de un ortoedro de volumen igual a 8 para que su superficie lateral sea mínima, puedes plantearlo como sigue: calcular el mínimo de f .x ; y ; z /
Universidad de Granada Dpto. de Análisis Matemático
D 2xy C 2xz C 2yz Prof. Javier Pérez Cálculo Diferencial en Rn
Extremos condicionados
92
(la función que da la superficie lateral de un ortoedro cuyos lados tiene longitudes x , y , z ) con la condición xyz xy z D 8. Se trata de un problema de extremos condicionados, pero la condición =.xy//, con lo que nuestra dada permite despejar una variable en función de las otras dos, z D 8=.xy 16=y C 16= 16=x , función de la que hay que calcular su función queda 2xy C 2xz C 2yz D xy C 16= 0 < y . Hemos convertido así el problema inicial de extremos mínimo absoluto cuando 0 < x , 0 < condicionados en uno de extremos libres porque ahora las variables .x ; y / se mueven en el abierto RC RC . Por cierto, por la desigualdad de las medias se tiene que:
p 3
.xy/. xy /.xz xz /.yz /.yz /
q D 3
.xyz xy z /2
D
p 3
64 6
xy
C xz C yz 3
Y la igualdad se da si, y sólo si, xy D xz D yz lo que implica que x D y D z D 2. El ortoedro buscado es el cubo. Pero si cambiamos la condición anterior por la siguiente x 2 yz 3
C sen.1 C xz / C y eyx D1
o bien, si imponemos dos condiciones como las siguientes: log.1 C x 2 C y 2 / C sin.1 C xz / 1 D 0;
xy z / C x 2 z 2 3 D 0 e1CyCxCz C cos.xyz
entonces no podemos usar esa condición (o condiciones) para despejar una variable (o dos variables) en función de las otras (de la otra). Ocurre que muchos de los ejercicios que se proponen en los libros de texto para ser resueltos utilizando la teoría de extremos condicionados que vamos a desarrollar seguidamente, se resuelven con frecuencia de forma más sencilla usando otro tipo de técnicas o algunas desigualdades como la de las medias geométrica y aritmética o la de Cauchy - Schwarz. Pero aunque su utilidad práctica pueda ser discutida, la teoría de extremos condicionados tiene un gran interés teórico, especialmente en la formulación formulación Lagrangiana de la Mecánica. En términos generales, el problema clásico de extremos condicionados puede enunciarse W ! R de clase C q y una función como sigue: dados un abierto Rn , un campo escalar f W vectorial g W ! Rnk de clase C q verificando verificando que el conjunto: conjunto:
D D fx 2 W g.x/ D 0 ^ rangoDg.x/ D n k g
M
extremos locales de la función f jM , restricción de f no es vacío, se trata de calcular los extremos a M . Es decir, calcular los máximos y mínimos locales de f .x/ cuando las variables, x D .x1 ; x2 ; : : : ; xn /, están obligadas a satisfacer las condiciones dadas por:
g.x/ D 0 ” gi .x1 ; x2 ; : : : ; xn / D 0 1 6 i 6 n k
(3.39)
donde g D .g1 ; g2 ; : : : ; gnk /.
En Mecánica clásica, las ecuaciones (3.39 ( 3.39)) representan las ligaduras (que los físicos llaman holónomas) que obligan a un punto material o a un sistema más complejo a desplazarse en una determinada determinada región del espacio, y el problema problema de optimizar ciertos campos escalares asociados al movimiento del sistema responde responde al problema de extremos condicionados. condicionados. El planteamiento del problema pone de manifiesto que M es una variedad en Rn de dimensión k y que se trata de determinar los extremos de la restricción de f a esa variedad. Universidad de Granada Dpto. de Análisis Matemático
Prof. Javier Pérez Cálculo Diferencial en Rn
Extremos condicionados
93
3.33 Definición. Sea f W ! R, donde Rn y sea M . Se dice que f tiene un máximo (resp. mínimo) local condicionado por M en un punto a 2 M , si hay un número r > 0 tal que para todo x 2 B.a; r / \ M se verifica que f .a/ > f .x/ (resp. f .a/ 6 f .x/). Se dice que f tiene en a un extremo local condicionado por M cuando f tiene en a un máximo o un mínimo local condicionado por M . Cuando las desigualdades anteriores son estrictas para todo x 2 B.x; r / con x ¤ a se dice que el extremo local condicionado correspondiente es estricto. El siguiente resultado es importante porque permite reducir un problema de extremos condicionados a un problema de extremos relativos ordinario.
3.34 Proposición. Sean f W !
R un campo escalar de clase C q en un abierto Rn y M una variedad de clase C q y dimensión k. Sean a M , .} ; U / una parametrización local de M en a y c el vector de coordenadas de a , es decir, }.c/ a. Se verifica entonces que f tiene en a un extremo local (estricto) condicionado por M si, y sólo si, el campo escalar f } U R tiene un extremo relativo (estricto) del mismo tipo en c .
2
D
ı W !
Parece que todo está resuelto, pero no es así porque en la práctica pocas veces se conocen las parametrizaciones locales de una variedad.
3.35 Proposición. Sean f W !
R un campo escalar de clase C q en un abierto Rn , M una variedad de clase C q y dimensión k y a M . Supongamos que f tiene en a un ? .a/. extremo local condicionado por M . Entonces se verifica que f .a/ T M
2
r
2
Demostración . Sea .}; U / una parametrización local de M en a y sea }. c/ D a. En virtud de la proposición anterior, la hipótesis hecha es lo mismo que afirmar que f ı } tiene un extremo relativo en c , en cuyo caso, como es un campo escalar diferenciable en c 2 U y U es abierto, c debe ser un punto crítico de f ı } , es decir, D.f ı }/. c/ D 0, esto es, Df .a/ ı D}.c/ D 0. Por tanto, para todo h 2 Rk es
˝ ˇˇ
˛
ı D}.c/ .h/ D r f .a/ D}.c/.h/ D 0 Como esto vale para todo h 2 Rk , resulta que el vector r f .a/ es ortogonal a D }.c/.Rk / D Df .a/
T M .a/.
En los dos resultados anteriores no se ha supuesto que la variedad M venga dada de alguna forma concreta, sin embargo el resultado siguiente requiere en su enunciado conocer las ecuaciones implícitas locales de M en a, por lo que de aquí en adelante vamos a suponer que la variedad M es la variedad de dimensión k D n m determinada por una función vectorial g W ! Rm de clase C q en un abierto Rn y tal que rango D g.x/ D m < n para todo x 2 : M D fx 2 W g.x/ D 0g Notaremos g D .g1 ; g2 ; : : : ; gm / los campos escalares componentes de g . Observa que esta hipótesis no es realmente restrictiva porque todas las variedades vienen dadas localmente de esta forma y el problema que nos ocupa es justamente un problema local. En lo que sigue f W ! R será un campo escalar de clase C q en . Universidad de Granada Dpto. de Análisis Matemático
Prof. Javier Pérez Cálculo Diferencial en Rn
Extremos condicionados
94
El siguiente resultado proporciona una condición necesaria para que f tenga en a 2 M un extremo local condicionado por M . Para enunciarlo conviene definir la llamada “función de Lagrange” ˇ W Rm ! R dada por: m
ˇ .x; /
Df .x/
X C
x .x1 ; x2 ; : : : ; xn / ; .1 ; 2 ; : : : ; m / Rm
D
i gi .x/;
iD1
2
D
2
(3.40)
Observa que es una función de n C m variables.
3.36 Teorema (de Lagrange). Supongamos que f tiene en a 2 M un extremo local condicio-
2 Rm tal que .a; ˛/ es un punto crítico de ˇ ,
nado por M . Entonces existe un único vector ˛ es decir, ˇ .a; ˛/ 0.
r
D
Demostración . Por el corolario 3.30 sabemos que los vectores r gi .a/, 1 6 i 6 m forman una ? ? .a/. Por la proposición 3.35 sabemos que r f .a/ 2 T M .a/. Por tanto, hay números base de T M m
˛i
2 R, 16i 6m, únicos, tales que r f .a/
X C iD1
˛i gi .a/
r
D 0. Poniendo ˛ D .˛1; ˛2 ; : : : ; ˛m /,
esta igualdad vectorial equivale a las igualdades: m
Dj f .a/
X C
˛i Dj gi .a/
iD1
D 0 ” Dj ˇ .a; ˛/ D 0
1 6 j 6 n
(3.41)
Puesto que DnCj ˇ .a; ˛/ D gj .a/ D 0 para 1 6j 6m porque a 2 M , se sigue que las igualdades (3.41) equivalen a que . a; ˛/ sea un punto crítico de ˇ .
Los números ˛1 ; ˛2 ; : : : ; ˛m en el teorema anterior se llaman los multiplicadores de Lagrange del punto a . Obtenemos así una regla práctica para proceder en un problema de extremos condicionados: se empieza por resolver el “sistema de Lagrange”: m
Dj f .x1 ; x2 ; : : : ; xn /
X C
i Dj gi .x1 ; x2 ; : : : ; xn /
iD1
gj .x1 ; x2 ; : : : ; xn /
D
0
D
0
9> = >;
1 6 j 6 n
Cm
Es un sistema de n C m ecuaciones en general no lineales para el que no puede haber técnicas generales de resolución. De hecho, puede resolverse de forma exacta en muy pocos casos. El teorema de Lagrange, que no es más que una reformulación equivalente de la proposición 3.35, da una condición necesaria para que un punto a 2 M sea extremo local de f condicionado por M . Daremos a continuación una condición suficiente. Suponemos en lo que sigue que nuestras funciones son de clase C q con q > 2. Dado . a; ˛/ 2 Rm , notaremos Qˇ .a; ˛/ la forma cuadrática sobre Rn cuya matriz en la base canónica es Dij ˇ .a; ˛/ 16i 6n . Se trata, pues, de la aplicación Qˇ .a; ˛/ W Rn ! R definida por:
16j 6n
Qˇ .a; ˛/.x/
n
X D
i;j D1
Universidad de Granada Dpto. de Análisis Matemático
Dij ˇ .a; ˛/xi xj
8x 2 Rn
(3.42)
Prof. Javier Pérez Cálculo Diferencial en Rn
Extremos condicionados
95
Consideraremos también una parametrización local .} ; U / de M en a con }.c/ D a. En estas condiciones, sean h D f ı } W U ! R y Qh .c/ la forma cuadrática sobre Rk dada por: k
Qh .c/.u/
X D
8u 2 Rk :
Dij h.c/ui uj
i;j D1
3.37 Proposición. Sea . a; ˛/ un punto crítico de la función de Lagrange. Con las notaciones anteriores, se verifica que: Qh .c/.u/
DQ
ˇ
.a; ˛/.D}.c/.u//
8u 2 Rk
(3.43)
Demostración . Se trata de hacer un cálculo. La mayoría de los textos lo evitan. Vamos a hacerlo. Notaremos por ei , 1 6 i 6 k , los vectores de la base canónica de Rk . Como h es un campo escalar, se verifica que Di h.c/ D r h.c/ ei D Df .c/.ei /. Para campos vectoriales, como } , esta relación también se mantiene si interpretamos que Di }.c/ es la columna i -ésima de la matriz jacobiana J } .c/. Dicha columna es el vector de Rn que se obtiene derivando las funciones componentes de } con respecto a la variable i -ésima que notaremos por Di }.c/. Por tanto D i }.c/ D D}.c/.ei /. Como hay que calcular derivadas parciales de segundo orden conviene calcular las derivadas parciales de primer orden en un punto genérico u 2 U para volver a derivarlas usando la regla de la cadena. Para u 2 U tenemos:
˝ ˇˇ ˛
Di h.u/
˛ ˝
ˇ ˛ ˇˇ ˛ C ˝r
D D.f ı }/. u/.ei / D Df.}.u// ı D}.u/ .ei / D Df.}.u// D}.u/.ei / D D r f .}.u// Di }.u/ D .r f / ı }/. u/ Di }.u/ (3.44)
ˇ
˝
Y volviendo a derivar otra vez como se deriva un producto escalar: Dij h.u/
˝
D Dj .Di h.u// D Dj ..r f / ı }/. u/ Di }.u/
f .}.u// Dj .Di }/. u/ (3.45)
ˇˇ
˛
Calculemos el primero de estos dos sumandos. Observa que r f W x 7! r f .x/ es un campo vectorial cuyas funciones componentes son las derivadas parciales primeras de f , . r f / ı } es el campo vectorial composición de r f y } , por lo que:
r ı }/. u/ D D..r f / ı }/. u/.ej / D D.r f /.}.u// ı D}.u/.ej /D D D.r f /.}.u//.Dj }.u// Sustituyendo esta igualdad en (3.45) y haciendo u D c tenemos: Dij h.c/ D D.r f /. a/.Dj }. c// Di }.c/ C r f .a/ Dj .Di }/. c/ Dj .. f /
ˇ
˝
˛ ˝ ˇ
Como .a; ˛/ es un punto crítico de la función de Lagrange se verifica Sustituyendo esta igualdad en (3.47) tenemos:
˛ X D
(3.46)
(3.47)
m
r f .a/
˛l gl .a/.
r
lD1
m
Dij h.c/
D˝
r
ˇˇ
D. f /. a/.Dj }.c// Di }.c/
Universidad de Granada Dpto. de Análisis Matemático
˛ X ˝r ˇˇ ˛l
lD1
gl .a/ Dj .Di }/. c/
˛
(3.48)
Prof. Javier Pérez Cálculo Diferencial en Rn
Extremos condicionados
96
Teniendo ahora en cuenta que g l ı }.u/ D 0 para todo u 2 U y para 1 6 l 6 m, y usando la igualdad (3.47) con el campo escalar h D f ı } sustituido por g l ı } , obtenemos: 0
ˇ
˝
˛ ˝ ˇ ˇˇ ˛
D Dij .gl ı }/. c/ D D.r gl /.a/.Dj }.c// Di }.c/ C r gl .a/ Dj .Di }/. c/
de donde se sigue que:
˝r ˇˇ
gl .a/ Dj .Di }/. c/
˛ D ˝
˛
D. gl /.a/.Dj }.c// Di }. c/
r
Sustituyendo esta igualdad en (3.48) resulta:
m
Dij h.c/
D ˝ r *
ˇˇ
D. f /.a/.Dj }. c// Di }.c/
D
m
r
D. f /.a/
X C
˛ C *X !
r
˛l D. gl /.a/.Dj }.c// Di }.c/
lD1
r
ˇ
˛l D. gl /.a/ .Dj }.c// Di }.c/
l D1
ˇˇ
+
+
D (3.49)
Pero esta igualdad es justamente lo que queríamos probar pues: m
Dij ˇ .a; ˛/
X C
D Dij f .a/
˛l Dij gl .a/
l D1
Que coincide con el elemento ij de la matriz jacobiana de: m
r
D. f /. a/
X C
r
˛l D. gl /.a/
lD1
(3.50)
Pues observa que la matriz jacobiana de D .r f /.a/ es la matriz hessiana de f en a y que la matriz jacobiana de D.r gl /.a/ es la matriz hessiana de gl en a. Por tanto, la matriz jacobiana de (3.50) es justamente la matriz Dij ˇ .a; ˛/ 16i 6n .
16j 6n
Lo que dice la igualdad (3.49) es que la forma cuadrática asociada a la matriz hessiana de h D f ı } en el punto c es igual a la restricción de la forma cuadrática asociada a la matriz Dij ˇ .a; ˛/ 16i 6n al espacio tangente a M en a. Para verlo más claro todavía, pongamos,
16j 6n
por comodidad de escritura,
A
D
16i 6n . La forma cuadrática asociada a dicha 16j 6n decir, Qˇ .a; ˛/, y viene dada para todo x Rn
Dij ˇ .a; ˛/
matriz, es justamente la definida en (3.42), es por: Qˇ .a; ˛/.x/
2
D x A xt D A xt x .
.
Para todo u D .u1 ; u2 ; : : : ; uk / 2 Rk tenemos que:
˝ ˇˇ ˛ .
k
D}. c/.u/
X D
uj Dj }.c/
j D1
Universidad de Granada Dpto. de Análisis Matemático
Prof. Javier Pérez Cálculo Diferencial en Rn
Extremos condicionados
97
es un vector del espacio tangente T M .a/ y:
* 0@X D X ˝ D
t
k
Qˇ .a; ˛/.D}. c/.u//
uj Dj }.c/
A
.
j D1
k
ui Di }.c/
iD1
k
ui uj
1 ˇX + Aˇ D ˇˇ ˛ D
A
t
.Dj }.c// Di }.c/
.
i;j D1
D Qh.c/.u/
k
por (3.49)
X D
Dij h.c/ui uj
i;j D1
Como queríamos probar.
Naturalmente, para calcular la restricción de la forma cuadrática Qˇ .a; ˛/ al espacio T M .a/, todo lo que se necesita es conocer una base de dicho espacio, pues si fa1 ; a2 ; : : : ; ak g es una tal base, la matriz de la forma cuadrática restringida a T M .a/ es la matriz .mij /kk cuyos elementos están dados por: at i aj
˝ ˇ˛
mij
D
A
1 6 i; j 6 k
.
(3.51)
Como consecuencia de la proposición 3.34, de la proposición anterior y del teorema 2.25, obtenemos condiciones suficientes para que f tenga en a un extremo local condicionado por M .
3.38 Teorema (Condición suficiente de extremo local condicionado) . Con las notaciones y las hipótesis anteriores, supongamos que . a; ˛/ es un punto crítico de la función de Lagrange y Dij ˇ .a; ˛/ 16i 6n . Sea a1 ; a2 ; : : : ; ak una base del espacio T M .a/ y sea pongamos A
D
f
16j 6n
g
.mij /kk la matriz cuyos elementos mij vienen dados por (3.51). Sea cuadrática:
Q
W Rk ! R la forma
k
Q.u/
X D
mij ui uj
i;j D1
Se verifica entonces que: a) Si la forma cuadrática estricto condicionado por M .
Q es
8u 2 Rk
definida positiva entonces f tiene en a un mínimo local
b) Si la forma cuadrática Q es definida negativa entonces f tiene en a un máximo local estricto condicionado por M . c) Si f tiene en a un máximo local condicionado por M entonces la forma cuadrática es semidefinida negativa.
Q
d) Si f tiene en a un mínimo local condicionado por M entonces la forma cuadrática es semidefinida positiva.
Q
e) Si la forma cuadrática condicionado por M .
Q es
no definida entonces f no tiene en a un extremo local
En general, la matriz A D Dij ˇ .a; ˛/
16i 6n 16j 6n
depende de a y de ˛. Te digo esto porque es
importante que calcules las soluciones del sistema de Lagrange en la forma . a; ˛/ agrupando Universidad de Granada Dpto. de Análisis Matemático
Prof. Javier Pérez Cálculo Diferencial en Rn
D
Cálculo de extremos absolutos
98
cada punto a con sus correspondientes multiplicadores de Lagrange. Cada a tiene asociado un único multiplicador ˛ , pero puede ocurrir que diferentes a tengan asociado el mismo multiplicador porque distintos a pueden dar lugar al mismo r f .a/. Los coeficientes m ij dependen en general de a y de ˛ y para cada par .a; ˛/ resulta, en general, una forma cuadrática Q diferente cuyo carácter debes estudiar. Pese a todo, de hecho, en la práctica sucede con frecuencia que A no depende de ˛, en tal caso puedes olvidarte de los valores de ˛ y considerar solamente los de a. La forma de calcular rápidamente la matriz .mij /kk es la siguiente: sea Bnk una matriz cuyas columnas sean los vectores base de T M .a/, entonces:
D B t A B
.mij /kk
.
.
(3.52)
3.6.1. Cálculo de extremos absolutos Cuando la variedad M es un conjunto compacto, la existencia de extremos absolutos para f está garantizada por la propiedad de compacidad de las funciones continuas. En este tipo de problemas lo único que hay que hacer es calcular los puntos críticos de la función de Lagrange y evaluar en ellos la función f , lo que nos proporcionará directamente los puntos donde f alcanza su máximo y su mínimo absolutos en M . Claro está, se supone que hay un número finito de soluciones del sistema de Lagrange. Situación diferente es cuando queremos calcular los extremos absolutos o relativos de una función f en un conjunto compacto K Rn con interior no vacío y cuya frontera M D Fr.K / es una variedad. Por supuesto, la existencia de extremos absolutos está garantizada pero pueden alcanzarse en el interior o en la frontera. Por tanto, lo que hay que hacer es calcular los puntos críticos de f en el interior de K y las soluciones del sistema de Lagrange para f y M . Obtendremos así, en la práctica, un conjunto finito de puntos y todo lo que hay que hacer es evaluar f en cada uno de ellos, lo que nos proporcionará directamente los puntos donde f alcanza su máximo y su mínimo absolutos en K . Dicho sea de paso, con frecuencia dichos extremos absolutos se alcanzan en la frontera de K . Sucede con frecuencia que la frontera de K Rk , k D 2 o k D 3 , no es una variedad aunque se le parezca mucho; por ejemplo, la frontera puede ser un polígono o un ortoedro o una pirámide. En tales casos lo más sencillo es parametrizar las caras del polígono o del ortoedro y reducir el problema de calcular los extremos en cada parte de la frontera de K a un problema de extremos más sencillo en una o dos variables. Claro, también se puede descomponer la frontera de K como unión de variedades disjuntas: las caras o las aristas más los vértices y aplicar la teoría de extremos condicionados a cada una de ellas. No te lo aconsejo.
Ejercicios propuestos 166. Calcula, aplicando la teoría de extremos condicionados, el ortoedro de volumen igual a 8 y superficie lateral mínima. Observación: este ejercicio lo hemos hecho en la introducción usando la desigualdad de las medias. La teoría de extremos condicionados te permitirá probar que el ortoedro Universidad de Granada Dpto. de Análisis Matemático
Prof. Javier Pérez Cálculo Diferencial en Rn
Cálculo de extremos absolutos
99
solución solamente puede ser el cubo, pero dicha teoría no te dice si el mínimo alcanzado es absoluto.
167. Calcula la distancia mínima de un punto P D .˛; ˇ; γ/ al plano Ax C By C C z D D : a) Usando la teoría de extremos condicionados. b) Usando la desigualdad de Cauchy - Schwarz. 168. Calcula el punto de la superficie 4x 2 C y 2 z 2 D 1 que está más próximo al origen de coordendas. 169. Una curva está dada como la intersección del cilindro x 2 C y 2 D 4 con el plano de ecuación 2x C 3y C z D 3. Calcula el punto de la misma que está más próximo al origen. 170. Sea A una matriz simétrica n n. Aplica el método de los multiplicadores de Lagrange para calcular el máximo absoluto de la forma cuadrática asociada a A condicionado por la esfera euclídea unidad en Rn , para probar que dicho máximo se alcanza en un punto de la forma . a; ˛/ tal que A at D ˛ a. Por tanto ˛ es un valor propio de A. .
171. Calcula el máximo absoluto de f .x; y; z/ D x 3y 2 z en el conjunto:
D
M
172.
.x; y; z/ RC o
2 RCo RCo W x C 2y C 3z D 1
˚
Utiliza el método de los multiplicadores de Lagrange para calcular el área del triángulo de mínima área cuyos lados son los segmentos de los ejes coordenados positivos interceptados por la tangente a la elipse: x2 a2
P D.u; v/
y2
C b2 D 1
en un punto P D .u; v/ de la misma situado en el primer cuadrante y dicho segmento de la tangente.
173. Utiliza el método de los multiplicadores de Lagrange para calcular un punto de la elipse de ecuación x2 a2
y2
C b2 D 1
tal que el segmento determinado por la intersección de la tangente a la elipse en dicho punto con los ejes coordenados tenga longitud mínima.
174. Prueba que el máximo absoluto de f .x1 ; x2 ; : : : ; xn / D .x1 x2 xn /2 sometido a la condición k.x1 ; x2 ; : : : ; xn /k2 D r es igual a .r 2=n/n . Deduce que:
q n
x12 x22
Universidad de Granada Dpto. de Análisis Matemático
1
xn2 6 n
n
X iD1
xi2
8.x1; x2 ; : : : ; xn/ 2 Rn Prof. Javier Pérez Cálculo Diferencial en Rn
Cálculo de extremos absolutos
100
y deduce la desigualdad de las medias: n
p x x x 6 1 n 1 2 n
X
n
8.x1; x2 ; : : : ; xn/ 2 .RCo /n
xi
iD1
¿Cuándo se da la igualdad?
175. Utiliza el método de los multiplicadores de Lagrange para calcular el volumen del ortoedro de máximo volumen inscrito en el elipsoide de ecuación x2 a2
y2
z2
C b2 C c 2 D 1
¿Sabes resolver este ejercicio usando la desigualdad de las medias?
176. Prueba que el máximo absoluto de x 2 y 2 sujeto a la condición x 2p p
C
y 2q q
D r2
donde p > 0 , q > 0 verifican que: 1
1 C D 1 p q
se alcanza cuando x 2p D y 2q y es igual a r 2 . Deduce que si x > 0 e y > 0 se verifica que: xy 6
xp p
q
C yq
¿Cuándo se da la igualdad?
177.
Se quiere construir un pentágono colocando un triángulo isósceles sobre un rectángulo, como se muestra en la figura de la derecha. Si el pentágono tiene un perímetro fijo p , determinar las longitudes de los lados del pentágono que maximizan su área.
178. Calcula los puntos de la curva
(
x2
xy C y 2 z 2 D 1 x2 C y2 D 1
que están más próximos al origen de coordenadas.
179. Calcular la mínima distancia del origen a la superficie de ecuación xy 2 z 3 D 2. 180. Calcula el mayor y el menor valor de la función f .x; y; z/ D x 2y C 5z en la esfera x 2 C y 2 C z 2 D 30. Universidad de Granada Dpto. de Análisis Matemático
Prof. Javier Pérez Cálculo Diferencial en Rn
Cálculo de extremos absolutos
101
181. Sea la curva intersección de la esfera x 2 C y 2 C z 2 D 1 y el plano x C y C z D 1. Calcula los puntos de que están más cerca y más lejos del punto .1; 2; 3/. Justifica que los resultados obtenidos son valores máximos y mínimos absolutos. 182. Encontrar un punto P de coordenadas positivas perteneciente al elipsoide de ecuación x2
2
2
y z C C D 1, tal que el plano tangente al elipsoide en P determine con los ejes a2 b2 c2
coordenados un tetraedro de volumen mínimo. El volumen del tetraedro con vértices en el origen y en los puntos .a; 0; 0/, .0; b; 0/ y .0; 0; c/ (donde a; b; c son números positivos), es igual a
1 6
abc .
183. Calcula los extremos absolutos de f .x; y/ D x 2 C y 2 x y en el conjunto K
n D
2
.x; y/ R
2
2 Wx Cy
2
64
o
^ y > 0
184. Calcula el volumen máximo de una caja rectangular situada sobre el plano X Y que tiene un vértice en el origen y el vértice opuesto, de coordenadas positivas, pertenece al paraboloide z C x 2 C 4y2 D 4. 185. Calcula los extremos p absolutos de f .x; y/ D x 2y 3 B k k2 ..0; 0/; 5/. 186. Calcula los extremos absolutos de f .x; y; z/ D x 2 cerrada.
C x 2 C 6y2 en la bola cerrada
C yz en la bola euclídea unidad
187. a) Clasifica los extremos relativos del campo escalar f .x; y/ D x 3 C y 3 xy 2 x C 16. b) Calcula el máximo y el mínimo absolutos de dicho campo escalar en el conjunto: K
n D
2
2
2 Wx Cy
.x; y/ R
2
o
6 2; y > 0 :
188. Se considera la función f W R2 ! R dada por: 3
f .x; y/
3
2
2
D x3 C y3 x2 y2 C 1
a) Calcula y clasifica los puntos críticos de f en la bola abierta B k k2 ..0; 0/; 2/. b) Calcula los extremos absolutos de f en la bola cerrada B k k2 ..0; 0/; 2/.
189. Calcula los puntos del conjunto A D .x; y/ 2 R2 W x 2 C y 2 6 2x; y > 0 donde el campo escalar f .x; y/ D x 2 2xy C y 2 alcanza su máximo y mínimo absolutos.
˚
190. Calcular los extremos absolutos de f .x; y/ D .x 2 C2y 2 / ex
2
y 2
en el disco x 2 C y 2 64.
191. Calcular los valores máximos y mínimos absolutos de f .x; y; z/ D xy 2 z 3 en la bola x 2 C y 2 C z 2 6 1. Universidad de Granada Dpto. de Análisis Matemático
Prof. Javier Pérez Cálculo Diferencial en Rn