Gradiente y Laplaciano en Polares

Comentarios sobre el Calculo de Gradientes y
Laplacianos en otros Sistemas de Coordenadas

Stefano Garcia
Aprendamos mediante un ejemplo heuristico.
Problema P3. Ayudantia 7:
Clacule el gradiente en coordenadas polares de
f(r, ) = ln(r) +
2
Recuerde que polares signica:
x = r cos()
y = r sin()
Ademas, necesitara saber que los vectores bases de las coordenadas polares son:
r =
_
cos()
sin()
_
= cos()
i+sin()
j (1a)
=
_
sin()
cos()
_
=sin()
i+cos()
j (1b)
Primero necesitamos entender la pregunta: "Encuentre el gradiente de f en coordenadas polares". Eso
signica: Expresar el vector

f ya no en la base

i y

j como estamos acostumbrados, sino que en la base
r y

que es la base de las coordenadas polares.
Partamos entonces con lo unico que tenemos: La forma de

f en cartecianas.
f =
_
f
x
,
f
y
_
=
f
x
i +
f
y
j (2)
* Notemos que
f
x
i +
f
y
j es el vector en si (piensenlo como una echa en R

2
), mientras que
_
f
x
,
f
y
_
son las coordenadas del vector

f en la base cartesiana

i y

j.
Notemos ahora que si cambiamos de sistema coordenado (cambiar

i y

j por r y

en este
caso), las coordenadas que deniran al mismo vector

f (el cual no se ha movido de su lugar
en el plano R
2
) cambiaran tambin , pero aun asi, seguira siendo verdad que

f =
f
x
i +
f
y
j,
sin embargo, ahora

f =
_
f
x
,
f
y
_
.
1
** Como ejemplo, imaginen el vector a = (2, 3) en cartesianas, esto es, vectores base

i = (1, 0)
y

j = (0, 1). Ahora piensen que invierten la direccion de eje x, es decir, hacemos un cambio de
coordenadas u =
i y v =

j; pues bien, bajo estas nuevas coordenadas, el vector a se escribira
como (2, 3), pero noten que el vector a no se ha movido a ningun lado durante todo este
rato, es decir, sigue situado en el mismo sitio de R
2
, lo unico que paso, fue que sus coordenadas
cambiaron. Por ende, uno puede escribir a = 2
i + 3
j = 2 u + 3 v ya que se trata del mismo

vector.
Necesitamos hallar entonces, f
x
y f
y
. Como f = f(r, ), nos conviene dejarlas expresadas en terminos
de f
r
y f
. Entonces:
f
x
= f
r
r
x
+ f
x
(3a)
f
y
= f
r
r
y
+ f
y
(3b)
Entonces, ahora debemos hallar r
x
r
y

x

y
.
Manera 1
Usaremos que la matriz Jacobiana de una transformada es la inversa de la matriz Jacobiana de la
transformada inversa. (A veces suele escribirse esto como [JF
1
] = [JF]
1
).
La matriz Jacobiana de la transformacion de (r, ) (x, y) es:
JF =
_
r
x
r
y
x

y
_
(4)
Mientras que la matriz Jacobiana de la transformacion de (x, y) (r, ) es:
JF
1
=
_
x
r
x
y
r
y
_
=
_
cos() r sin()
sin() r cos()
_
[JF
1
]
1
=
1
|JF
1
|
_
r cos() r sin()
sin() cos()
_
=
_
cos() sin()
sin()/r cos()/r
_
(5)
Y por lo que dijimos al comienzo, comparando componente a componente las matrices (4) y (5), somos
capace de establecer:
r
x
= cos() (6a)
r
y
= sin() (6b)
x
= sin()/r (6c)
y
= cos()/r (6d)
2
Manera 2
Otra forma hubiese sido tomar
x = r cos()
y = r sin()
Y derivar abmas respecto a x para obtener el sistema lineal:
1 = r
x
cos() r sin()
x
0 = r
x
sin() + r cos()
x
Luego con Cramer obtenemos:
r
x
=
1 r sin()
0 r cos()
cos() r sin()
sin() r cos()
=
r cos()
r
= cos()
Idem para
x
.
Para obtener r
y
y
y
se procede de manera analoga, solo que ahora partimos derivando respecto de y.
Ahora estamos en condiciones de reemplazar (6) en (3) y obtener:
f
x
= f
r
r
x
+ f
x
=f
x
= f
r
cos() f
sin()
r
(7a)
f
y
= f
r
r
y
+ f
y
=f
y
= f
r
sin() + f
cos()
r
(7b)
Ahora debemos reemplazar las ecuaciones (7) en (2):
f =
_
f
x
,
f
y
_
=
_
f
r
cos() f
sin()
r
, f
r
sin() + f
cos()
r
)
=
_
f
r
cos() f
sin()
r
_
i +
_
f
r
sin() + f
cos()
r
_
j (8)
Lo cual seria la forma de

f en la base

i y

j (cartesiana).
Calculemos entonces para f(r, ) = ln(r) +
2
. Se tendra:
3
f
r
= 1/r
f
= 2
=

f =
_
cos()
r
2
sin()
r
_
i +
_
sin()
r
+ 2
cos()
r
_
j (9)
Entonces (9) es el gradiente de f en coordenadas cartesianas.
Notemos que este debiese de ser el mismo resultado que conseguiriamos si reemplazasemos
x = r cos() y y = r sin() en f(r, ) y luego simplemente hiciesemos

f = (f
x
, f
y
) = f
x
i +f
y
j.
Encontremos ahora el gradiente de f en coordenadas polares.
Para ello retomamos (8) y en ella debemos hacer aparecer los vectores base de la ecuacion (1).
f =
_
f
x
,
f
y
_
=
_
f
r
cos() f
sin()
r
_
i +
_
f
r
sin() + f
cos()
r
_
j
= f
r
_
cos()
i + sin()
j
_
. .
r
+
f
r
_
sin()
i + cos()
j
_
. .
Con lo que obtenemos:

f= f
r
r +
f
=
_
f
r
,
f
r
_
Que es el gradiente de f en coordenadas polares.
Calculemos entonces para f(r, ) = ln(r) +
2
. Se tendra:
f
r
= 1/r
f
= 2
=

f =
1
r
r +
2
r
(10)
Claramente el gradiente de f en coordenadas polares es mucho ms bonito que en cartesianas. Esto es
porque f originalmente era una funcion simple en las variables r y mientras que era signicativamente
ms compleja vista desde las valiables x y y.
4
Problema P4. Ayudantia 7:
Hallar el laplaciano en coordenadas polares. Esto es:

2
f =

f f
Lo primero que hay que notar es que el laplaciano de una funcion es un escalar, como podemos
ver

f es un producto punto y como tal, un escalar. Esto implica que toda la discucion
anterior de los vectores base, ect, esta fuera de lugar aqui. En otras palabras, nu usaremos nada
de vectores esta vez.
Tenemos que en cartesianas, el laplaciano se escribe como:
f =

2
f
x
2
+

2
f
y
2
(11)
Pues lo unico que tenemos que hacer entonces, es escribir (10) en terminos de r, y sus deribadas.
De la ecuacion (7) del ejercicio anterior, tenemos ya que:
f
x
= f
r
r
x
+ f
x
(12a)
f
y
= f
r
r
y
+ f
y
(12b)
Derivamos (11a) respecto de x:
=f
xx
= (f
r
r
x
)
x
+ (f
x
)
x
= f
rx
r
x
+ f
r
r
xx
+ f
x
x
+ f
xx
= (f
rr
r
x
+ f
r
x
)r
x
+ f
r
r
xx
+ (f
r
r
x
+ f
x
)
x
+ f
xx
= f
rr
(r
x
)
2
+ f
r
x
r
x
+ f
r
r
xx
+ f
r
r
x
x
+ f
(
x
)
2
+ f
xx
(13a)
Analogamente, tan solo cambiando las x por y, obtenemos:
f
yy
= f
rr
(r
y
)
2
+ f
r
y
r
y
+ f
r
r
yy
+ f
r
r
y
y
+ f
(
y
)
2
+ f
yy
(13b)
Nuevamente, del ejercicio anterior, ecuacion (6), sabemos que:
r
x
= cos() =r
xx
= (cos())
x
= (sin())
_
sin()
r
_
=
sin
2
()
r
(14a)
r
y
= sin() =r
yy
= (sin())
y
= (cos())
_
cos()
r
_
=
cos
2
()
r
(14b)
5
x
=
sin()
r
=
xx
=
_
sin()
r
_
x
=
__
sin()
r
_
r
r
x
+
_
sin()
r
_
x
_
=
__
sin()
r
2
_
cos() +
cos()
r
_
sin()
r
__
= 2
sin() cos()
r
2
(14c)
y
=
cos()
r
=
yy
=
_
cos()
r
_
y
=
__
cos()
r
_
r
r
y
+
_
cos()
r
_
y
_
=
_
cos()
r
2
_
sin()
sin()
r
_
cos()
r
_
= 2
sin() cos()
r
2
(14d)
Y ahora podemos reemplazar los valores de (14) en (13):
f
xx
= f
rr
cos
2
()f
r
sin() cos()
r
+f
r
sin
2
()
r
f
r
sin() cos()
r
+f
_
sin()
r
_
2
+f
sin() cos()
r
2
(15a)
f
yy
= f
rr
sin
2
()+f
r
sin() cos()
r
+f
r
cos
2
()
r
+f
r
sin() cos()
r
+f
_
cos()
r
_
2
f
sin() cos()
r
2
(15b)
Y sumando las ecuaciones en (15) llegamos a:
f = f
xx
+ f
yy
= f
rr
+ 0 +
1
r
f
r
+
1
r
2
f
+ 0
(16)
f = f
rr
+
1
r
f
r
+
1
r
2
f
Podemos probarlo para calcular el laplaciano de f = ln(r) +

2
. Se tendra:
f
r
= 1/r
f
rr
= 1/r
2
f
= 2
=f =
1
r
2
+
1
r
1
r
+
1
r
2
2 = 2/r
2
(17)
6

Gradiente y Laplaciano en Polares

Cargado por

Copyright:

Formatos disponibles

Gradiente y Laplaciano en Polares

Cargado por

Información del documento

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

Gradiente y Laplaciano en Polares

Cargado por

Copyright:

Formatos disponibles

Comentarios sobre el Calculo de Gradientes y

Laplacianos en otros Sistemas de Coordenadas

j es el vector en si (piensenlo como una echa en R

j = 2 u + 3 v ya que se trata del mismo

Con lo que obtenemos:

Podemos probarlo para calcular el laplaciano de f = ln(r) +

También podría gustarte