ISSN en trmite

Fascculo 3 Notas de
matemtica

Norberto Fava
Guillermo Keilhauer

Integrales sobre la esfera Sn-1

Departamento de Matemtica
Facultad de Ciencias Exactas y Naturales
Universidad de Buenos Aires
2014
Notas de matemtica

Fascculo 3

Comit Editorial:

Carlos Cabrelli (Director)

Departamento de Matemtica, FCEyN, Universidad de Buenos Aires
E-mail: cabrelli@[Link]

Gabriela Jernimo
Departamento de Matemtica, FCEyN, Universidad de Buenos Aires
E-mail: jeronimo@[Link]

Claudia Lederman
Departamento de Matemtica, FCEyN, Universidad de Buenos Aires
E-mail: clederma@[Link]

Leandro Vendramin
Departamento de Matemtica, FCEyN, Universidad de Buenos Aires.
E-mail: lvendramin@[Link]

ISSN en trmite

Derechos reservados
2014 Departamento de Matemtica, Facultad de Ciencias Exactas y Naturales,

Universidad de Buenos Aires.

Departamento de Matemtica
Facultad de Ciencias Exactas y Naturales
Universidad de Buenos Aires
Ciudad Universitaria Pabelln I
(1428) Ciudad de Buenos Aires
Argentina.
[Link]
e-mail. secre@[Link]
tel/fax: (+54-11)-4576-3335
 Integrales sobre la esfera S n1
N. Fava y G. Keilhauer

Resumen
Las coordenadas esfericas y las integrales sobre la esfera juegan un papel destacado en
muchas cuestiones del An alisis en los Espacios Eucldeos. Otros enfoques del mismo
tema se hallar
an entre los tems de la bibliografa.

1. Coordenadas esf
ericas. Sera u
til comenzar por un estudio cuidadoso de las co-
ordenadas esfericas en el espacio eucldeo de dimension n. Nos referimos a la aplicacion
x = T (r, 1 , . . . , n1 ) de una parte de Rn en s mismo, definida por las formulas

x1 = r cos 1
x2 = r sen 1 cos 2
x3 = r sen 1 sen 2 cos 3
.. ..
(1) . .
xn1 = r sen 1 sen 2 sen n2 cos n1
xn = r sen 1 sen 2 sen n2 sen n1 ,

donde r > 0, 0 < i < para 1 6 i 6 n 2 y 0 < n1 < 2.

Las restricciones sobre los argumentos son para lograr que T sea un difeomorfismo entre
dos subconjuntos abiertos de Rn , en la forma que precisamos a continuacion:

El punto = (1 , . . . , n1 ) vara en el intervalo L = (0, )n2 (0, 2) de Rn1 , de modo

que (r, ) vara en el conjunto P = (0, ) L.

De las ecuaciones anteriores se sigue facilmente que T (P ) G, donde G = {x : xn 6=

0 o xn1 < 0}. En efecto, si xn = 0 entonces n1 = y por consiguiente xn1 < 0.

Por induccion sobre n se probaran las siguientes afirmaciones:

1. Si x = T (r, ), entonces r2 = x21 + x22 + + x2n ;

2. T aplica biyectivamente P sobre G;

1
 3. El jacobiano de T es J = rn1 senn2 1 senn3 2 sen n2 .

Antes de probarlas hagamos dos observaciones: (i ) el complemento de G consta de los

puntos x de Rn que satisfacen las relaciones xn = 0, xn1 > 0; (ii ) el jacobiano de T , que
escribimos abreviadamente en la forma

(2) J = rn1 g(),

donde g() = senn2 1 senn3 2 sen n2 , es positivo en P . Luego, la aplicacion inversa

T 1 es de clase C en G.

Para probar las afirmaciones anteriores expresamos T en la forma T = T2 T1 , donde T2

y T1 se definen por medio de los siguientes sistemas de ecuaciones:

T2 T1

x1 = x1 x1 = r cos 1
x2 = r1 cos 2 r1 = r sen 1
x3 = r1 sen 2 cos 3 2 = 2
x4 = r1 sen 2 sen 3 cos 4 3 = 3
.. .. .. .. ..
. . . . .
xn1 = r1 sen 2 sin n2 cos n1 n2 = n2
xn = r1 sen 2 sen n2 sen n1 n1 = n1

Notemos que T2 act

ua como la identidad en la primera coordenada y transforma las
un unas ecuaciones analogas a las de T aunque en el espacio de dimension n 1;
restantes seg
en tanto que T1 transforma las dos primeras coordenadas seg
un las ecuaciones polares del
plano y act
ua como la identidad en las restantes.

Suponiendo n > 3 y denotando por P 0 el conjunto analogo de P en el espacio Rn1 ,

tendremos:
T T
1
P R P 0
2
G.

Supongamos ahora que las tres afirmaciones son verdaderas en el espacio eucldeo de
dimension n 1. Entonces tendremos:

r2 = x21 + r12 = x21 + x22 + + x2n .

2
 T es biyectiva porque lo son T1 y T2 . Es claro que 0 < 1 < por ser r1 > 0.

Denotando por Ji el jacobiano de Ti (i = 1, 2), tendremos:

(x1 , r1 , 2 , , n1 ) (x1 , r1 )
J1 = = = r;
(r, 1 , 2 , , n1 ) (r, 1 )

y en virtud de la hipotesis inductiva,

(x1 , x2 , , xn1 , xn ) (x2 , , xn1 , xn )

J2 = =
(x1 , r1 , 2 , , n1 ) (r1 , 2 , , n1 )

= r1n2 senn3 2 senn4 3 sen n2 .

La demostracion concluye teniendo en cuenta que J = J2 J1 y reemplazando r1 por su

valor r sen 1 .

2. Elemento de
area de una hipersuperficie. Una hipersuperficie H se define local-
mente por medio de una aplicacion de clase C

: U Rn ,

donde U es un conjunto abierto de Rn1 . Escribimos x = (u) y representamos los vectores

de Rn en forma de matriz columna:

x1 1 (u)
. .
x= . .
. = . .
xn n (u)

Suponemos que la derivada 0 (u) tiene rango n 1 en cada punto de U . Bajo esta
hipotesis, decimos que es una inmersion y el par (U, ) se llama carta local o sistema local
de coordenadas de H.

Los vectores
x
vk (u) = = (k = 1, 2, . . . , n 1)
uk uk
son linealmente independientes y generan el subespacio tangente Hx en el punto x = (u)
de H.

3
 Denotando por N (u) un vector unitario en el complemento ortogonal de Hx , el elemento
de area de H se define por medio de la formula

(3) d = g(u) du,

donde:

(4) g(u) = | det(N (u), v1 (u), . . . , vn1 (u))|.

La justificacion heurstica reside en el hecho de que g(u) representa el volumen del par-
aleleppedo engendrado por N (u) y los vectores tangentes v1 (u) . . . , vn1 (u) que generan
su base. Dicho volumen es igual al area de la base debido a que el vector N (u) es ortogonal
al subespacio tangente y tiene norma 1. La justificacion se refuerza con la ecuacion
n1
X
(u + h) (u) = hk vk (u) + o(|h|),
k=1

en la que suponemos 0 6 hk 6 k (1 6 k 6 n 1). Esa ecuacion da idea de como se

transforma, aproximadamente, el intervalo de Rn1 de lados [uk , uk + k ].

La integral de una funcion f (x) sobre H se define por medio de la formula

Z Z
f d = f ((u)) g(u) du.
H U

Otra forma u
til de expresar la funcion g es la siguiente:

v
u n1
uX
(5) g(u) = t J 2 (u), k
k=1

donde
(1 , . . . , k , . . . , n )
Jk (u) =
(u1 , u2 , . . . , un1 )
y el smbolo k indica la supresion de k .

Para demostrar (5) basta estudiar la forma lineal

fu (x) = det(x, v1 (u), . . . , vn1 (u)),

cuyo n
ucleo es el subespacio tangente H(u) . Denotando por ek los vectores de la base
canonica de Rn , tendremos:
n
! n
X X
fu (x) = fu xk ek = xk fu (ek ) = hx, (u)i,
k=1 k=1

4
donde (u) es el vector de componentes fu (ek ) = (1)k1 Jk (u), 1 6 k 6 n. Puesto que
(u) es ortogonal a cada uno de los vectores vk (u) existe un escalar = (u) tal que
(u) = N (u)) y por consiguiente,

g(u) = |fu (N (u))| = |hN (u), N (u)i| = || = |(u)|,

lo que demuestra (5).

En general no basta una sola carta local para describir globalmente una hipersuperficie.
Un estudio mas general del tema requiere los conceptos de variedad diferenciable y particion
de la unidad.

Con ayuda de las expresiones anteriores se demuestra que el area tiene una significacion
geometrica intrnseca; es decir, que es invariante tanto bajo cambio de coordenadas locales
como bajo tranformaciones ortogonales del espacio. Ahora estamos mejor preparados para
estudiar el caso particular de la esfera unitaria.

area de la esfera unitaria Sn-1 . La aplicacion

3. Elemento de

x0 = T (1, ) ( L),

donde T es la transformacion del primer parrafo y vara en el intervalo L, abarca la esfera

unitaria S n1 con excepcion de un arco de meridiano conjunto de area nula, interseccion
de la esfera con un semi-hiperplano.

Conviene escribir T en la forma x = rx0 . Los vectores:

x0
vi () = (i = 1, 2, . . . , n 1),
i
que son linealmente independientes por lo que se vera enseguida, generan el subespacio tan-
gente en el punto x0 . Por otro lado, x0 es un vector unitario ortogonal a dicho subespacio.
En efecto, por derivacion con respecto a i en la ecuacion hx0 , x0 i = 1, se obtiene:
x0
hx0 , i = 0 (1 6 i 6 n 1).
i
El jacobiano de la transformacion T , que es positivo y hemos escrito en la forma rn1 g(),
puede expresarse tambien del modo siguiente:
x0 x0 x0
   
x x x x 0
J = det , , ,..., = det x , r ,r ,...,r
r 1 2 n1 1 2 n1

0 0
x0
 
n1 0 x x
=r det x , , ,..., ,
1 2 n1

5
lo que demuestra que los vectores vi () son linealmente independientes. El u
ltimo determi-
nante es la funcion g() de la formula (2). Luego, el elemento de area de la superficie S n1
puede expresarse por medio de la formula

dx0 = g() d = senn2 1 senn3 2 sen n2 d1 dn1 ,

y la formula del cambio de variables para coordenadas esfericas puede escribirse del modo
siguiente:
Z Z Z
(6) f (x) dx = dr r n1
f (rx0 ) dx0 ,
n
R 0 S n1

formula que es valida para cualquier funcion f medible no negativa.

4. Invariancia bajo rotaciones. En el caso de la esfera existe un camino directo para

probar que el area de cada subconjunto de la misma es invariante bajo rotaciones.
Consideremos un conjunto boreliano M L y sea E = T (1, M ) la imagen de M en la
esfera. Denotemos por E el sector correspondiente de la bola unitaria, es decir, el conjunto

E = {rx0 : M, 0 < r 6 1}.

Por la formula del cambio de variables, la medida o volumen de dicho sector es

Z Z Z 1 Z
n1 n1 1
| E | = dx = r g() dr d = dr r g() d = A(E).
E (0,1)M 0 M n

Hemos probado que el area de un conjunto E de la esfera se relaciona con el volumen del
sector correspondiente por medio de la formula

A(E) = n| E |,

cuyo segundo miembro es invariante bajo rotaciones.

En lo que sigue haremos uso de la formula

Z Z
p1 t 2
(p) = t e dt = 2 r2p1 er dr.
0 0
2 2 2 2
Si en (6) se pone f (x) = e|x| = ex1 ex2 exn y se aplican el teorema de Fubini y la
formula anterior, se obtiene el area de la superficie esferica unitaria:

2 n/2
A(S n1 ) = .
(n/2)

6
 5. Integral de un polinomio sobre la esfera. Para calcular una integral de la forma
Z
P (x) dx0 ,
S n1

donde P es un polinomio, basta considerar el caso de un monomio x = x1 1 x2 2 xnn .

Suponemos, pues, que los j son enteros no negativos.

j +1
Poniendo j = 2
, dejaremos como importante ejercicio demostrar la formula:

Z


0 si alg
un j es impar

x dx0 =
S n1

2(1 )(2 ) (n )

si todos los j son pares
(1 +2 + +n )

Problemas

1. Exprese la medida (o volumen) de la bola unitaria B = {x : |x| 6 1}.

2. Pruebe que los vectores vi () del parrafo 3 son mutuamente ortogonales.

3. La funcion f (x) se llama radial si existe f0 : R>0 R tal que

f (x) = f0 (|x|). Probar que si f es radial, entonces
Z Z
n1
f (x) dx = A(S ) rn1 f0 (r) dr
Rn 0

4. Exprese en terminos de la funcion gamma la siguiente integral relacionada con el n

ucleo
de Poisson en Rn : Z
dx
n+1
Rn (1 + |x|2 ) 2

5. Pruebe que si A es un subconjunto boreliano de Rn con la propiedad de que para

cualquier vector v de norma 1, el conjunto Av = {t R : tv A} es de medida nula,
entonces A es de medida nula.

6. Pruebe que si M es un subconjunto convexo de Rn , la frontera de M es de medida nula

y M es medible.

7. Siendo A una matriz simetrica positiva definida, muestre que el volumen del elipsoide
E = {x : hAx, xi 6 1} esta dado por la formula
n
2
|E| = n
( 2 + 1) det A

7
8. Para cada permutacion de los primeros n enteros positivos, pongamos x = (x(1) , . . . , x(n) ),
() = T (1, ) ( L). Muestre que las imagenes (L) cubren toda la esfera S n1
y que para dos permutaciones cualesquiera, la aplicacion 1
de L en s
mismo es
de clase C , lo que da a la esfera una estructura de variedad compacta orientable de
clase C .

Bibliografa

J.A. Baker Integration over Spheres and the Divergence Theorem for Balls, Amer.
Math. Monthly 104, Number 1 (1997), 36-47.

L.E. Blumenson, A Derivation of n-dimensional Spherical Coordinates, Amer. Math.

Monthly 67 (1960), 63-66.

N.A. Fava, G. Keilhauer, A. Larotonda, Surface Integrals, Spherical Coordinates and

the Area Element of S n1 , Revista de la U.M.A. 30 (1981-82), 77-84.

G.B. Folland, How to Integrate a Polynomial Over a Sphere, Amer. Math. Monthly
108 (2001) 446-448.

E.M. Stein, G. Weiss, Introduction to Fourier Analysis on Euclidean Spaces, Princeton

University Press, 1971.