AC 5 Volumen de La Hiperesfera

Descargar como pdf o txt
Descargar como pdf o txt
Está en la página 1de 25

Volúmenes de una hiperesfera en M dimensiones

Test
Tecnológico de Monterrey

Martes, 24 de Septiembre

Introducción al problema
En la Mecánica Estadı́stica frecuentemente necesitamos conocer el volumen de una hiperes-
fera en M dimensiones. En este problema lo vamos a calcular. A continuación, se describen las
instrucciones:

Parte 1: El volumen de una hiperesfera en M dimensiones se define como:


Z
VM = dx1 . . . dxM (1)
x21 +x22 +···+x2M ≤R2

Parte 2: Mostrar que VM = RM CM , donde CM es el volumen de la (hiper)esfera unitaria en


M dimensiones.

Parte 3: Para calcular vamos a integrar una función que depende solamente del radio ge-
neralizado en M dimensiones, no de las coordenadas angulares generalizadas. Una función
conveniente (¿por qué?) para este propósito es la siguiente:

f (x1 , . . . , xM ) = exp −(x21 + · · · + x2M )



(2)

Mostrar que:
Z ∞
f (x1 , . . . , xM ) dx1 . . . dxM = π M/2 (3)
−∞

Parte 4: Expresar la integral (3) del inciso (3) en coordenadas esféricas generalizadas. Su-
gerencia: Ya que VM = RM CM , un diferencial de volumen correspondiente a un cascarón
ubicado entre los radios R y R + dR se puede expresar como dVM = ∂V ∂R dR.
M

Parte 5: Evaluar la integral del inciso (4). El resultado dependerá de CM .

Parte 6: Igualar las expresiones de los incisos (3) y (5) para calcular CM . Usar la función Γ.

Indicaciones adicionales
Tu solución debe seguir todos y cada uno de los pasos antes mencionados. Debes demostrar que
sabes hacer cálculos matemáticos, como estas integrales sencillas. Si requieres hacer algún cambio
de variable, indı́calo claramente y describe con palabras tus procedimientos, es decir, argumenta
por qué haces tal o cual paso.
Imagina que le estás explicando a alguien de tercer semestre tu procedimiento, pero de forma
escrita. NO SE ACEPTA QUE HAGAS TODOS LOS DESARROLLOS MATEMÁTICOS Y AL
FINAL O AL INICIO DESCRIBAS LOS PASOS. ESO NO ES PROFESIONAL, ESO NO ES DE
UN ALUMNO DE 5º SEMESTRE.

1
Hiperesferas en M dimensiones Test

Parte 1: Definición del volumen de una hiperesfera en M dimensio-


nes
Para comenzar con este problema, nos centramos en la definición del volumen de una hiper-
esfera en M dimensiones. Recordemos que el volumen de una hiperesfera es el análogo en varias
dimensiones de la conocida fórmula del volumen de una esfera en 3D. En este caso, una hiperesfera
en M dimensiones se define por el conjunto de puntos (x1 , x2 , . . . , xM ) tales que:

x21 + x22 + · · · + x2M ≤ R2


El volumen VM de la hiperesfera está dado por la siguiente integral múltiple sobre las coorde-
nadas cartesianas x1 , x2 , . . . , xM :
Z
VM = dx1 . . . dxM
x21 +x22 +···+x2M ≤R2

La expresión antes mostrada es una generalización de la idea clásica que tenemos sobre el
volumen de una esfera en 3 dimensiones pero extrapolado a M dimensiones. La integral se toma
sobre todas las variables cartesianas que satisfacen la condición de estar dentro de una región
esférica de radio R.

Interpretación fı́sica
Este es un concepto importante donde trabajaremos con configuraciones en espacios de muchas
dimensiones. Un ejemplo tı́pico es el espacio de fase de un sistema con muchas partı́culas, donde
cada partı́cula añade tres dimensiones adicionales a dicho espacio, una idea que estudiamos en
cuarto semestre, en la clase teórica de Mecánica Clásica.
En M dimensiones, la suma de los cuadrados de las coordenadas x21 + x22 + · · · + x2M ≤ R2
representa una superficie de igual radio en ese espacio, similar a cómo en 3D, la condición x2 +
y 2 + z 2 = R2 define la superficie de una esfera.

Motivación matemática
En nuestro intento de ver el problema desde una vista más matemática utilizamos la geometrı́a
multivariada y, por tanto, análisis de sistemas complejos. El volumen de la hiperesfera nos permitirá
calcular cantidades relacionadas con el espacio disponible en configuraciones de alta dimensionali-
dad.
La integral múltiple que define el volumen en coordenadas cartesianas puede ser complicada
de evaluar directamente. Sin embargo, durante el desarrollo de cada uno de los puntos propues-
tos por esta actividad, la expresión (1) será transformada a coordenadas esféricas generalizadas,
permitiéndonos un manejo más ’natural’ en los cálculos.

2
Hiperesferas en M dimensiones Test

Parte 2: Demostración de VM = RM CM
Nuestro objetivo es demostrar que el volumen VM de una hiperesfera de radio R en M dimen-
siones se puede expresar como:

VM = RM CM
donde CM es el volumen de la hiperesfera unitaria en M dimensiones (es decir, una hiperesfera
de radio 1).

Entendiendo la relación entre el radio y el volumen en dimensiones superiores


En el espacio tridimensional, el volumen V3 de una esfera de radio R está dado por:
4
V3 = πR3
3
Observamos que el volumen es proporcional al cubo del radio. Este hecho, aunque no contamos
con una muestra de volúmenes, nos empieza a dar pistas sobre que, en general, el volumen de
una hiperesfera en M dimensiones deberı́a ser proporcional a RM . Nuestra meta es formalizar y
demostrar esta intuición.

Demostración matemática detallada


Partimos de la definición del volumen VM de una hiperesfera de radio R en M dimensiones:
Z
VM = dx1 dx2 . . . dxM (4)
x21 +x22 +···+x2M ≤R2

Esta integral representa el volumen del conjunto de puntos (x1 , x2 , . . . , xM ) que satisfacen x21 +
x22 + · · · + x2M ≤ R2 .

Cambio de variables
Para relacionar VM con CM , realizaremos un cambio de variables que normalice la hiperesfera
de radio R a una hiperesfera unitaria. Definimos:
xi
yi =para i = 1, 2, . . . , M (5)
R
Lo que logramos con este cambio de variables es esencialmente una transformación de escala
que ’encoge’ la hiperesfera de radio R a una hiperesfera de radio 1, esto al eliminar dividiendo por
el valor del radio R a toda componente correspondiente a cada una de las M dimensiones. Puedes
observar como diversos autores realizan este mismo procedimiento como en [5, 11].

Dominio de integración transformado


Con este cambio, el dominio de integración se transforma de la siguiente manera:

x21 + x22 + · · · + x2M ≤ R2 =⇒ R2 y12 + y22 + · · · + yM


2
≤ R2 =⇒ y12 + y22 + · · · + yM
2

≤1

Por lo tanto, el dominio de integración ahora es la hiperesfera unitaria en M dimensiones.

3
Hiperesferas en M dimensiones Test

Cálculo detallado del jacobiano del cambio de variables


A ti, joven ficticio de tercer semestre al que le estoy explicando mi tarea, te preguntarás
por qué y quién es Jacob?, permı́teme contarte un poco al respecto.
Cuando realizamos un cambio de variables en una integral múltiple, necesitamos ajustar el
elemento de volumen para asegurar que la integral en las nuevas variables represente el mismo valor
que en las variables originales. Este ajuste se hace mediante el jacobiano J de la transformación,
que es el determinante de la matriz de derivadas parciales de las variables originales respecto a las
nuevas variables.

Definición del jacobiano Para un cambio de variables de M dimensiones, donde las variables
originales xi son funciones de las nuevas variables yj , el jacobiano J se define como:
 
∂(x1 , x2 , . . . , xM ) ∂xi
J= = det (6)
∂(y1 , y2 , . . . , yM ) ∂yj
Es decir, J es el determinante de la matriz jacobiana J cuyos elementos son las derivadas
∂xi
parciales ∂y j
.

Expresión de las variables originales en términos de las nuevas variables En nuestro


caso, estamos realizando el cambio de variables:

xi = Ryi para i = 1, 2, . . . , M (7)


donde R es una constante positiva (el radio de la hiperesfera) y yi son las nuevas variables.

∂xi
Cálculo de las derivadas parciales Necesitamos calcular las derivadas parciales ∂yj para todos
los valores de i y j.
Dado que xi depende únicamente de yi (no hay dependencia de yj para j ̸= i), tenemos:
(

∂xi (Ryi ) = R si i = j
= ∂yj (8)
∂yj 0 si i ̸= j
Esto se debe a que:

Para i = j, derivamos xi = Ryi respecto a yi , obteniendo R.

Para i ̸= j, xi = Ryi no depende de yj , por lo que la derivada parcial es cero.

Construcción de la matriz jacobiana La matriz jacobiana J es una matriz de M × M cuyos


∂xi
elementos son ∂y j
:
 ∂x ∂x1 ∂x1

∂y1
1
∂y2 . . . ∂y M
 ∂x2 ∂x2 ∂x2 
 ∂y1 ∂y2 . . . ∂y 
M 
J= .
 . . . (9)
 .. .. .. .. 
∂xM ∂xM
∂y1 ∂y2 . . . ∂x
∂yM
M

Sustituyendo las derivadas parciales calculadas, obtenemos:

4
Hiperesferas en M dimensiones Test

 
R 0 ... 0
0 R ... 0
J=. . . (10)
 
.. 
 .. .. .. .
0 0 ... R
Esta es una matriz diagonal con R en cada entrada de la diagonal principal y ceros en todas
las demás posiciones.

Cálculo del determinante de la matriz jacobiana El determinante de una matriz diagonal


es igual al producto de los elementos de la diagonal. Por lo tanto:
M   M
Y ∂xi Y
J = det(J) = = R = RM (11)
∂yi
i=1 i=1

Esto se debe a que:

Cada elemento diagonal es R, y hay M elementos diagonales.

No hay contribución al determinante de los elementos fuera de la diagonal, ya que son ceros.

Interpretación El jacobiano J = RM representa el factor de escala volumétrico entre el elemento


de volumen en las variables originales dx1 dx2 . . . dxM y el elemento de volumen en las nuevas
variables dy1 dy2 . . . dyM :

dx1 dx2 . . . dxM = J dy1 dy2 . . . dyM = RM dy1 dy2 . . . dyM (12)
Esto significa que al escalar cada coordenada xi por R, el volumen elemental se escala por
RM , lo cual es consistente con la idea de que el volumen en M dimensiones escala con la M -ésima
potencia del factor de escala lineal. Para leer un poco más acerca del jacobiano puedes recurrir a
los siguientes recursos [1, 3, 8]

Transformación de la integral
La integral original se transforma utilizando el cambio de variables y el jacobiano:

Z
VM = dx1 dx2 . . . dxM (13)
x21 +···+x2M ≤R2
Z
∂(x1 , x2 , . . . , xM )
= dy1 dy2 . . . dyM (14)
y12 +···+yM
2 ≤1 ∂(y1 , y2 , . . . , yM )
Z
= RM dy1 dy2 . . . dyM (15)
y12 +···+yM
2 ≤1
Z
= RM dy1 dy2 . . . dyM (16)
y12 +···+yM
2 ≤1

5
Hiperesferas en M dimensiones Test

Identificación con CM
Reconocemos que la integral:
Z
CM = dy1 dy2 . . . dyM
2 ≤1
y12 +y22 +···+yM

es precisamente el volumen de la hiperesfera unitaria en M dimensiones.


Por lo tanto, podemos escribir:
Z !
VM = RM dy1 dy2 . . . dyM = RM CM (17)
y12 +···+yM
2 ≤1

Interpretación fı́sica y matemática


La relación VM = RM CM indica que el volumen de una hiperesfera en M dimensiones es
proporcional a la M -ésima potencia del radio. Esta relación generalizada nos permite entender
cómo el volumen cambia con el radio en espacios de dimensionalidad arbitraria. Esto es coherente
con nuestra intuición geométrica:

En 1 dimensión (M = 1): Una ’hiperesfera’ es simplemente un segmento de longitud 2R,


y su ”volumen”(longitud) es proporcional a R1 .

En 2 dimensiones (M = 2): La hiperesfera es un cı́rculo, y su área es proporcional a R2 .

En 3 dimensiones (M = 3): La hiperesfera es la esfera habitual, y su volumen es propor-


cional a R3 .

Ejemplo numérico
Para ilustrar esta relación utilizando información con la que ya estamos familiarizados, consi-
deremos el caso de M = 2:
El volumen de la hiperesfera unitaria en 2D (el área de un cı́rculo de radio 1) es:

C2 = π(1)2 = π
Entonces, el volumen de un cı́rculo de radio R es:

V 2 = R 2 C2 = R 2 π
Lo cual coincide con la fórmula conocida para el área de un cı́rculo.
Ahora, para un segundo ejemplo que es fácil de comprobar con ayuda de una tabla básica de
geometrı́a o, para los más generales, la realización de una clásica integral triple en coordenadas
esféricas, consideraremos el caso de M = 3:
El volumen de la hiperesfera unitaria en 3D (el volumen de una esfera de radio 1) es:

C3 =
3
Entonces, el volumen de una esfera de radio R es:

6
Hiperesferas en M dimensiones Test

4π 3
V3 = R3 C2 = R
3
Lo cual coincide con la fórmula conocida para el volumen de una esfera para cualquier radio.
Habiendo realizado todo esto podemos decir que demostramos que el volumen de una hiperesfera
en M dimensiones es VM = RM CM , donde CM es el volumen de la hiperesfera unitaria.

Parte 3: Cálculo detallado de la integral gaussiana multidimensional


Nuestro objetivo es demostrar que:
Z ∞
f (x1 , x2 , . . . , xM ) dx1 dx2 . . . dxM = π M/2
−∞
donde la función f está definida como:

f (x1 , x2 , . . . , xM ) = exp −(x21 + x22 + · · · + x2M )




Justificación de la elección de la función f


Primero responderemos esa cuestión inicial del por qué? es conveniente esta expresión f . Bien,
la elección de la función f es especialmente conveniente por las siguientes razones:

1. Dependencia radial: La función depende únicamente de r2 = x21 + x22 + · · · + x2M , que es el


cuadrado de la distancia desde el origen en M dimensiones. Esto significa que f es una función
radialmente simétrica, lo que simplifica el cambio a coordenadas esféricas generalizadas y
permite separar variables angulares y radiales.

2. Propiedades integrales conocidas: La integral de la función gaussiana unidimensional


2 √
e−x es bien conocida y resulta en π. Extender este resultado a múltiples dimensiones nos
permitirá calcular la integral multidimensional de manera más sencilla.

3. Facilidad para manejar lı́mites infinitos: La función f decrece rápidamente a cero cuando
xi → ±∞, lo que asegura la convergencia de la integral en el intervalo de −∞ a ∞.

Cálculo de la integral multidimensional


Queremos calcular la integral:
Z ∞
exp −(x21 + x22 + · · · + x2M ) dx1 dx2 . . . dxM

I= (18)
−∞
Procederemos paso a paso, justificando cada uno.

7
Hiperesferas en M dimensiones Test

Paso 1: Separación de variables


Observamos que el exponente en el integrando es una suma de términos que dependen única-
mente de cada variable xi :

−(x21 + x22 + · · · + x2M ) = −x21 − x22 − · · · − x2M


Debido a que las variables xi están separadas en el exponente y los lı́mites de integración son
los mismos para cada variable (−∞ a ∞), podemos expresar la integral multidimensional como
un producto de integrales unidimensionales independientes y posteriormente utilizar una expresión
más elegante empleando el sı́mbolo de producto o como también es conocido, usando el sı́mbolo de
la ’Piatoria’:

Z ∞ Z ∞ Z ∞
2 2 2
I= e−x1 dx1 e−x2 dx2 · · · e−xM dxM (19)
−∞ −∞ −∞
M Z ∞ 
−x2i
Y
= e dxi (20)
i=1 −∞

Paso 2: Cálculo de la integral gaussiana unidimensional


La integral gaussiana unidimensional, es decir, la integral correspondiente a la dimensión x1 es:
Z ∞
2
I1 = e−x dx (21)
−∞
Aunque este resultado es conocido, lo demostraremos detalladamente para asegurar una com-
prensión completa.

Demostración de I1 = π Consideramos el cuadrado de I1 :
Z ∞  Z ∞  Z ∞Z ∞
2 −x2 −y 2 2 −y 2
I1 = e dx e dy = e−x dx dy (22)
−∞ −∞ −∞ −∞
Aquı́ hemos introducido una nueva variable y para representar la segunda integral.

Cambio a coordenadas polares Para evaluar esta integral doble, cambiamos a coordenadas
polares en el plano (x, y):

x = r cos θ (23)
y = r sin θ (24)
∂(x, y)
dx dy = dr dθ = r dr dθ (25)
∂(r, θ)

El jacobiano del cambio de variables es r.

Lı́mites de integración Como es bien sabido, en el cambio de coordenadas a polares para el


plano bidimensional, la variable r va de 0 a ∞. Y, por otro lado, θ va de 0 a 2π, que corresponde
a una ’vuelta’ completa en el plano.

8
Hiperesferas en M dimensiones Test

Integral en coordenadas polares La integral se convierte en:


Z 2π Z ∞
2
I12 = e−r r dr dθ (26)
0 0

Evaluación de la integral radial Calculamos primero la integral respecto a r:


Z ∞   
−r2 1 −r2 ∞ 1 1
e r dr = − e = 0− − = (27)
0 2 0 2 2
2
La exponencial e−r tiende a cero cuando r → ∞, y es igual a 1 cuando r = 0.

Evaluación de la integral angular La integral respecto a θ es:


Z 2π
dθ = 2π (28)
0

Resultado total Por lo tanto:


 
1
I12 = 2π = π =⇒ I1 = π (29)
2

Paso 3: Regreso a la integral multidimensional


Sabemos ahora que:
Z ∞ √
2
e−xi dxi = π para cada i = 1, 2, . . . , M (30)
−∞
Por lo tanto, la integral multidimensional es:
M
Y √ √ M
I= π= π = π M/2 (31)
i=1

Interpretación geométrica
La integral I representa el volumen bajo la superficie definida por f (x1 , x2 , . . . , xM ) en todo el
espacio M -dimensional. El resultado π M/2 es una extensión natural del resultado unidimensional a
múltiples dimensiones.

Aplicación en coordenadas esféricas


Aunque en este caso no fue necesario debido a que la función f utilizada resulto en ser una
expresión conocida, es instructivo ver cómo se podrı́a calcular esta integral utilizando coordenadas
esféricas generalizadas, esto aunque podrı́a ser tedioso en primera instancia, resultara ser útil en
partes posteriores del problema. Para leer más acerca de las coordenadas generalizadas puedes
recurrir a las siguientes fuentes [1–3, 6, 8, 10].

9
Hiperesferas en M dimensiones Test

Paso 1: Cambio a coordenadas esféricas generalizadas


En M dimensiones, las coordenadas esféricas se definen por:

x1 = r cos ϕ1 (32)
x2 = r sin ϕ1 cos ϕ2 (33)
x3 = r sin ϕ1 sin ϕ2 cos ϕ3 (34)
..
. (35)
xM −1 = r sin ϕ1 sin ϕ2 . . . sin ϕM −2 cos ϕM −1 (36)
xM = r sin ϕ1 sin ϕ2 . . . sin ϕM −2 sin ϕM −1 (37)
El elemento de volumen en coordenadas esféricas es:
M −1
!
Y
dx1 dx2 . . . dxM = rM −1 sinM −k−1 ϕk dr dϕ1 dϕ2 . . . dϕM −1 (38)
k=1

Paso 2: Reescritura de la integral


La integral se transforma en:
M −1
! !
Z ∞ Z
−r2
Y
I= e rM −1 dr sinM −k−1 ϕk dϕ1 dϕ2 . . . dϕM −1 (39)
r=0 ángulos k=1

Paso 3: Separación de las integrales radial y angular


Definimos:
M −1
Z !
Y
SM = sinM −k−1 ϕk dϕ1 dϕ2 . . . dϕM −1 (40)
ángulos k=1
Este término representa el ’área’ de la superficie de una hiperesfera unitaria en M dimensiones.
La integral total es entonces:
Z ∞
2
I = SM e−r rM −1 dr (41)
0

Paso 4: Cálculo de la integral radial


La integral radial es:
Z ∞
2
R= e−r rM −1 dr (42)
0
Realizamos el cambio de variable u = r2 =⇒ du = 2rdr:

1
rM −1 dr = u(M −1)/2−1 du (43)
2Z
1 ∞ −u (M −1)/2
=⇒ R = e u du (44)
2 0

10
Hiperesferas en M dimensiones Test

Reconocemos que esta es la definición de la función Gamma:


Z ∞
Γ(z) = e−u uz−1 du (45)
0
Por lo tanto:
 
1 M
R= Γ (46)
2 2

Paso 5: Cálculo de SM
Sabemos que el volumen de la hiperesfera unitaria en M dimensiones es:

2π M/2
Vunidad = (47)
MΓ M

2
Y la superficie de la hiperesfera unitaria es:
1
SM = M Vunidad = 2π M/2 M
 (48)
Γ 2

Paso 6: Resultado final


Entonces, la integral total es:

I = SM × R (49)
!  
1 1 M
= 2π M/2 M
 Γ (50)
Γ 2
2 2
= π M/2 (51)
Lo cual coincide con el resultado obtenido anteriormente mediante la separación de variables.

Parte 4: Expresión de la integral en coordenadas esféricas genera-


lizadas
Nuestro objetivo es expresar la integral:
M
!
Z ∞ X
I= exp − x2i dx1 dx2 . . . dxM (52)
−∞ i=1
en coordenadas esféricas generalizadas. Además, utilizaremos la sugerencia de que, dado que
VM = RM CM , el diferencial de volumen correspondiente a un cascarón esférico entre los radios R
y R + dR se puede expresar como:

∂VM
dVM = dR = M RM −1 CM dR (53)
∂R

11
Hiperesferas en M dimensiones Test

Cambio a coordenadas esféricas generalizadas


Recuperando del desarrollo de la parte 3, sabemos que en un espacio de M dimensiones, las
coordenadas esféricas generalizadas (r, θ1 , θ2 , . . . , θM −1 ) se definen de la siguiente manera:

x1 = r cos θ1 (54)
x2 = r sin θ1 cos θ2 (55)
x3 = r sin θ1 sin θ2 cos θ3 (56)
..
. (57)
xM −1 = r sin θ1 sin θ2 . . . sin θM −2 cos θM −1 (58)
xM = r sin θ1 sin θ2 . . . sin θM −2 sin θM −1 (59)

Interpretación de las coordenadas


r: Distancia radial desde el origen hasta el punto (x1 , x2 , . . . , xM ).
θ1 , θ2 , . . . , θM −2 : Ángulos que varı́an entre 0 y π, análogos al ángulo polar en tres dimensiones.
θM −1 : Ángulo que varı́a entre 0 y 2π, análogo al ángulo azimutal en tres dimensiones.
Estas coordenadas permiten describir cualquier punto en M dimensiones mediante un radio y
M − 1 ángulos, generalizando las coordenadas esféricas conocidas en tres dimensiones. [1–3, 6, 8].

Cálculo del elemento de volumen


El elemento de volumen en coordenadas esféricas generalizadas se obtiene calculando el deter-
minante del jacobiano del cambio de variables:

∂(x1 , x2 , . . . , xM )
dV = dr dθ1 dθ2 . . . dθM −1
∂(r, θ1 , θ2 , . . . , θM −1 )
Este determinante resulta en:
M −2
!
θk (sin θM −1 )0 dr dθ1 dθ2 . . . dθM −1
Y
M −1 M −k−1
dV = r sin (60)
k=1
Notemos que (sin θM −1 )0 = 1, ya que en el último término no aparece la función seno elevado
a ninguna potencia.

Justificación del elemento de volumen El elemento de volumen en coordenadas esféricas


generalizadas se construye multiplicando:
El factor radial rM −1 dr, que representa la expansión del volumen en la dirección radial en M
dimensiones.
Q −1 M −k−1
Los factores angulares M k=1 sin θk dθk , que representan cómo el volumen se distribuye
en las direcciones angulares.
Este elemento generaliza el conocido elemento de volumen en coordenadas esféricas tridimen-
sionales:

dV = r2 sin θ dr dθ dϕ

12
Hiperesferas en M dimensiones Test

Reescritura de la integral en coordenadas esféricas


Dado que la función en el integrando depende únicamente de r, ya que:
M
X
x2i = r2
i=1

podemos reescribir la integral I como:


Z ∞ Z
−r2 M −1
I= e r dr dΩM (61)
r=0 ΩM
R
donde ΩM dΩM representa la integración sobre todas las variables angulares, y dΩM es el
elemento de ángulo sólido en M dimensiones:
M −2
!
Y
dΩM = sinM −k−1 θk dθ1 dθ2 . . . dθM −1 (62)
k=1

Separación de variables La integral total se separa en:

Integral radial:
Z ∞
2
Ir = e−r rM −1 dr (63)
0

Integral angular:
Z
SM = dΩM (64)
ΩM

Por lo tanto, podemos escribir:

I = Ir · S M (65)

Cálculo de la integral angular SM


La integral angular SM es el área de la superficie de una hiperesfera unitaria en M dimensiones.
Existe una expresión conocida para SM :

2π M/2
SM = (66)
Γ M

2

Justificación Esta expresión se obtiene mediante la evaluación de la integral múltiple de los


factores seno elevados a potencias especı́ficas, que aparecen debido a la geometrı́a de la hiperesfera
en M dimensiones. La deducción detallada de SM se encuentra en referencias clásicas de análisis
matemático multidimensional. Puedes observar algunos ejemplos en las referencias adjuntadas a
este párrafo. [4, 8, 9]. Al igual que, puedes leer más acerca de la función gamma en los siguientes
recursos [1, 3, 8].

13
Hiperesferas en M dimensiones Test

Relación con el diferencial de volumen dVM


Recordemos que el volumen de una hiperesfera de radio R es:

VM = RM CM (67)
El diferencial de volumen correspondiente a un cascarón esférico entre los radios R y R + dR
es:

dVM
dVM = dR = M RM −1 CM dR (68)
dR

Comparación con el elemento de volumen El elemento de volumen obtenido en la ecuación


(60) también puede interpretarse como el volumen de un cascarón esférico infinitesimal:

dV = rM −1 dr · dΩM (69)
Al integrar dΩM sobre todas las variables angulares, obtenemos SM :
Z
dΩM = SM (70)
ΩM

Por lo tanto, el elemento de volumen total es:

dV = rM −1 dr · SM (71)

Expresión de dVM en términos de dV


Dado que ambos representan el mismo diferencial de volumen en un cascarón esférico, podemos
igualar:

dVM = dV =⇒ M RM −1 CM dR = rM −1 dr · SM (72)
Observando que R = r y dR = dr, obtenemos:

M CM = SM (73)
Esta relación nos permite expresar CM en términos de SM :

SM
CM = (74)
M

Expresión final de la integral


Sustituyendo Ir de la ecuación (63) y SM de la ecuación (66) en la ecuación (65), obtenemos:

2π M/2
Z 
2
I= e−r rM −1 dr · (75)
Γ M

0 2
SM
Pero recordando que CM = M , podemos escribir:

I = Ir · M C M (76)

14
Hiperesferas en M dimensiones Test

Evaluación de la integral radial Ir


Para calcular Ir , realizamos el cambio de variable:

du
u = r2 =⇒ du = 2rdr =⇒ rdr = (77)
2
Entonces, rM −1 dr = rM −2 (rdr) = 12 u(M −2)/2 du
Sustituyendo en Ir :

Z ∞    du 
−u (M −2)/2
Ir = e u (78)
u=0 2
Z ∞
1
= e−u u(M/2)−1 du (79)
2 0
 
1 M
= Γ (80)
2 2

Donde hemos utilizado la definición de la función Gamma:


Z ∞
Γ(z) = e−u uz−1 du (81)
0

Expresión completa de I
Sustituyendo Ir en la ecuación (76):


 
1 M
I= Γ · M CM (82)
2 2
 
M M
= CM Γ (83)
2 2

Conclusión de la Parte 4
Hemos logrado expresar la integral I en términos de CM , M y la función Gamma:
 
M M
I= CM Γ (84)
2 2
Este resultado es crucial, ya que nos permite relacionar la integral calculada en la Parte 3 con
CM , lo que nos llevará a determinar CM explı́citamente en las siguientes dos secciones faltantes
correspondientes a la parte cuatro y cinco. Una vez más, te recuerdo que puedes leer más acerca
de la función gamma en los siguientes recursos [1, 3, 8].

15
Hiperesferas en M dimensiones Test

Parte 5: Evaluación de la integral y dependencia de CM


En la Parte 4, obtuvimos la expresión de la integral I en términos de CM , M y la función
Gamma:
 
M M
I= CM Γ (85)
2 2
Por otro lado, en la Parte 3, habı́amos encontrado que:

I = π M/2 (86)
Nuestro objetivo en esta sección es igualar estas dos expresiones de I para resolver CM en
términos de π, M y la función Gamma.

Igualación de las expresiones de I


Igualamos las ecuaciones (85) y (86):
 
M M
CM Γ = π M/2 (87)
2 2

Despeje de CM
Despejamos CM de la ecuación (87):

2 π M/2
CM = (88)
MΓ M

2

Esta expresión nos proporciona el valor de CM en términos de constantes y funciones conocidas.

Interpretación del resultado


La ecuación (88) es una fórmula general para calcular el volumen de la hiperesfera unitaria en
M dimensiones:

2π M/2
CM = (89)
MΓ M

2
Este resultado es fundamental, ya que establece una conexión directa entre el volumen de la
hiperesfera unitaria y la función Gamma, permitiendo calcular CM para cualquier valor de M .

Verificación en dimensiones especı́ficas


Para confirmar la validez de nuestra fórmula, evaluaremos CM para valores especı́ficos de M .

16
Hiperesferas en M dimensiones Test

Caso M = 1
1
 √
La función Gamma tiene el valor Γ 2 = π. Sustituyendo en la ecuación (89):

2π 1/2
C1 = (90)
1 · Γ 21


2 π
= √ (91)
1· π
=2 (92)

Esto es consistente con el hecho de que el ’volumen’ de una hiperesfera unitaria en una dimensión
(un intervalo de longitud) es 2.

Caso M = 2
La función Gamma tiene el valor Γ(1) = 1. Sustituyendo en la ecuación (89):

2π 1
C2 = (93)
2 · Γ(1)

= (94)
2·1
=π (95)

Este resultado coincide con el área de un cı́rculo unitario en dos dimensiones.

Caso M = 3

3
 π
La función Gamma tiene el valor Γ 2 = 2 . Sustituyendo en la ecuación (89):

2π 3/2
C3 = (96)
3 · Γ 23


2π 3/2
= √  (97)
3 · 2π
2π 3/2 · 2
= √ (98)
3 π
4π 3/2
= √ (99)
3 π
4π (3/2)−(1/2)
= (100)
3
4π 1
= (101)
3

= (102)
3
Este es el volumen de una esfera unitaria en tres dimensiones, lo que confirma la validez de
nuestra fórmula.

17
Hiperesferas en M dimensiones Test

Con esta corroboracion contra valores que ya conocemos, podemos decir que ha quedado de-
mostrado que al igualar las dos expresiones de I, podemos obtener una fórmula general para CM
en términos de π, M y la función Gamma.

Parte 6: Cálculo de CM usando la función Gamma


En esta sección, consolidaremos nuestros resultados y profundizaremos en la importancia de la
función Gamma en el cálculo de CM .

Fórmula general para CM


Como resultado final, tenemos:

2π M/2
CM = (103)
MΓ M

2
Esta fórmula nos permite calcular el volumen de la hiperesfera unitaria en cualquier dimensión
M.

Importancia de la función Gamma


La función Gamma Γ(z) es una extensión del factorial para números reales y complejos, definida
por:
Z ∞
Γ(z) = e−t tz−1 dt para Re(z) > 0 (104)
0

Propiedades relevantes:

Γ(n) = (n − 1)! para n entero positivo.


 √
Γ 12 = π.

Γ(z + 1) = zΓ(z).

Estas propiedades son esenciales para simplificar y calcular CM en dimensiones especı́ficas.

Ejemplo de cálculo para M = 4


Calculemos C4 utilizando la fórmula general.
La función Gamma para M = 4 es:
 
4
Γ = Γ(2) = 1! (105)
2
Sustituyendo en la ecuación (103):

18
Hiperesferas en M dimensiones Test

2π 4/2
C4 = (106)
4 · Γ(2)
2π 2
= (107)
4·1
π2
= (108)
2
π2
El volumen de una hiperesfera unitaria en cuatro dimensiones es 2 .

Comprobación de coherencia inicial sobre el desarrollo


A lo largo de este trabajo, hemos:

Definido el volumen de una hiperesfera en M dimensiones y demostrado que VM = RM CM .

Calculado una integral gaussiana multidimensional y encontrado su valor en términos de


π M/2 .

Expresado la integral en coordenadas esféricas generalizadas y relacionado I con CM y la


función Gamma.

Igualado las expresiones de I obtenidas por dos métodos diferentes para resolver CM .

Verificado que nuestros resultados son coherentes con casos conocidos en dimensiones bajas.

Con esta primera recapitulación y habiendo utilizado una vez nuestra expresión para la cuarta
dimensión hemos completado con éxito el cálculo de CM utilizando la función Gamma.

Parte 7: Recapitulación y análisis adicional


Recapitulación de los pasos teóricos
A lo largo de este trabajo, hemos abordado el problema de calcular el volumen de una hiperesfera
en M dimensiones, siguiendo los siguientes pasos:

1. Definición del volumen de una hiperesfera en M dimensiones: Iniciamos definiendo el


volumen VM mediante una integral múltiple sobre el dominio que satisface x21 +x22 +· · ·+x2M ≤
R2 .

2. Demostración de que VM = RM CM : Mediante un cambio de variables, establecimos que


el volumen de una hiperesfera de radio R es proporcional a RM multiplicado por CM , el
volumen de la hiperesfera unitaria.

3. Cálculo de la integral gaussiana multidimensional: Calculamos la integral gaussiana


para M dimensiones y el resultado fue I = π M/2 .

19
Hiperesferas en M dimensiones Test

4. Expresión de la integral en coordenadas esféricas generalizadas: Reescribimos la


integral en coordenadas esféricas y relacionamos el elemento de volumen con el diferencial
dVM de un cascarón esférico.

5. Evaluación de la integral y dependencia de CM : Al evaluar la integral radial y angular,


encontramos una expresión de I en términos de CM y la función Gamma.

6. Cálculo de CM usando la función Gamma: Igualando las expresiones de I, obtuvimos


2π M/2
la fórmula general CM = .
MΓ M 2

Cálculos para dimensiones de 1 a 4 y comparación con ecuaciones conocidas


Calcularemos CM para M = 1, 2, 3, 4 y compararemos los resultados con las fórmulas conocidas.

Caso M = 1


2π 1/2 2 π
C1 =  = √ =2
1 · Γ 12 π
V1 = R1 C1 = 2R

Esto coincide con la longitud de un segmento de −R a R, lo cual es coherente con nuestra


intuición en una dimensión.

Caso M = 2

2π 1 2π
C2 = = =π
2 · Γ(1) 2·1
V2 = R2 C2 = πR2

Este es el área de un cı́rculo de radio R, lo cual coincide con la fórmula conocida.

Caso M = 3

2π 3/2 2π 3/2 4π 1 4π
C3 = 3
 =  √  = =
3·Γ 2 3· 2 π 3 3

4
V3 = R3 C3 = πR3
3
Este es el volumen de una esfera en tres dimensiones, coincidiendo con la fórmula clásica.

20
Hiperesferas en M dimensiones Test

Caso M = 4

2π 2 2π 2 π2
C4 = = =
4 · Γ(2) 4·1 2
2
π 4
V4 = R4 C4 = R
2
Este es el volumen de una hiperesfera en cuatro dimensiones. Aunque menos familiar, este
resultado es consistente con las fórmulas obtenidas.

Evolución del volumen de una hiperesfera en función de M


A continuación, presentamos una gráfica que muestra cómo varı́a el volumen de una hiperesfera
unitaria (R = 1) al aumentar el número de dimensiones M .

Volumen
5
Volulemes mmáximos
Volumen de la hiperesfera unitaria

0
5 10 15 20
Dimensiones (M)

Figura 1: Evolución del volumen de una hiperesfera unitaria (R = 1) en función del número de
dimensiones M . En este gráfico, las tres lı́neas verticales resaltan tres valores, estos pertenecen a
los tres volúmenes más altos para las hiperesferas donde, se puede observar claramente que el valor
del volumen máximo de una hiperesfera se da en la dimensión cinco, el hecho de que conforme
aumentan las dimensiones el volumen decae tendiendo a cero será discutido más adelante. Gráfica
dada por elaboración propia.

Código en Julia para el cálculo y la gráfica


Para calcular los volúmenes y generar la gráfica, utilizamos el siguiente código en Julia:

using SpecialFunctions
using Plots
pgfplotsx()

# Función para calcular C_M

21
Hiperesferas en M dimensiones Test

function volumen_hiperesfera(M)
return (2 * pi^(M / 2)) / (M * gamma(M / 2))
end

# Rango de dimensiones
M_values = 1:20
volumes = [volumen_hiperesfera(M) for M in M_values]

# Imprimir tabla de volúmenes


println("Dimensión\tVolumen")
for (M, V) in zip(M_values, volumes)
println("$M\t\t$V")
end

# Graficar el volumen vs dimensiones


plot(M_values, volumes, marker=:o, xlabel="Dimensiones (M)",
ylabel="Volumen de la hiperesfera unitaria", color=:purple,
label="Volumen", legend=:topright, grid=false)

# Agregar lı́neas verticales punteadas en los tres valores más altos


vline!([4, 5, 5], linestyle=:dash, color=:black, label="Volulemes mmáximos")

savefig("volumen_hiperesfera.pdf")

Explicación del código:

Importamos las librerı́as SpecialFunctions para acceder a la función Gamma y Plots para
generar la gráfica. La función pgfplotsx() se utiliza para que las gráficas sean compatibles
con LATEX.

Definimos la función volumen hiperesfera(M) que calcula CM usando la fórmula obtenida.

Calculamos los volúmenes para dimensiones de 1 a 20 y almacenamos los resultados en la


lista volumes.

Imprimimos una tabla en la ventana de comandos que muestra la dimensión y el volumen


correspondiente.

Generamos la gráfica del volumen en función de M y la guardamos en formato pdf.

Salida de la tabla de volúmenes


Al ejecutar el código anterior, se mostrará en la ventana de comandos una tabla como la si-
guiente:

Dimensión Volumen
1 2.0
2 3.141592653589793
3 4.188790204786391
4 4.934802200544679

22
Hiperesferas en M dimensiones Test

5 5.263789013914324
6 5.167712780049969
7 4.724765970331401
8 4.058712126416768
9 3.2985089027387065
10 2.550164039877345
11 1.8841038793898997
12 1.3352627688545893
13 0.9106287547832835
14 0.5992645293207919
15 0.3814432808233041
16 0.2353306303588932
17 0.14098110691713998
18 0.08214588661112819
19 0.04662210680255007
20 0.025806891390014044

Análisis del comportamiento del volumen al aumentar M


Observamos que el volumen de la hiperesfera unitaria aumenta inicialmente con el número de
dimensiones, alcanza un máximo y luego disminuye tendiendo a cero para dimensiones muy altas.

Explicación del fenómeno: Este comportamiento puede parecer contraintuitivo al principio,


pero tiene una explicación matemática y geométrica:

Competencia entre factores crecientes y decrecientes: La fórmula del volumen

2π M/2
CM =
MΓ M

2

implica dos factores contrapuestos al aumentar M :

1. El término π M/2 crece exponencialmente con M .


2. El término Γ M

2 crece aún más rápido que exponencialmente (superexponencial), ya
que la función Gamma generaliza el factorial.

Máximo del volumen: Debido a esta competencia, el volumen alcanza un máximo alrededor
de M = 5 y luego disminuye.

Disminución
 en altas dimensiones: En dimensiones muy altas, el término en el denomi-
nador Γ M
2 domina, haciendo que el volumen tienda a cero.

Intuición geométrica: En espacios de alta dimensionalidad, el ’espacio’ se vuelve cada vez


más ’hueco’. La mayor parte del volumen de un hipercubo unitario (el análogo multidimensional
de un cubo de lado 1) se concentra lejos del origen y cerca de las esquinas. La hiperesfera, al
estar centrada en el origen, ocupa una fracción cada vez menor del volumen total disponible en el
hipercubo circundante.

23
Hiperesferas en M dimensiones Test

Analogı́a: Imaginemos una esfera dentro de un cubo en tres dimensiones. El volumen de la esfera
es menor que el del cubo, pero ocupa una porción significativa. Al aumentar las dimensiones, la
esfera (hiperesfera) ocupa una porción cada vez más pequeña del hipercubo, debido a la creciente
’lejanı́a’ de las esquinas en el espacio multidimensional.
Estas interpretaciones tomaron base en el video divulgativo [7].

24
Hiperesferas en M dimensiones Test

Referencias
[1] George B. Arfken and Hans J. Weber. Mathematical Methods for Physicists. Academic Press,
7th edition, 2012.

[2] L. E. Blumenson. A derivation of n-dimensional spherical coordinates. The American Mathe-


matical Monthly, 67(1):63–66, 1960.

[3] Mary L. Boas. Mathematical Methods in the Physical Sciences. John Wiley & Sons, 3rd edition,
2006.

[4] Wikipedia contributors. N-sphere — wikipedia, the free encyclopedia. https://en.


wikipedia.org/wiki/N-sphere, 2024. Accessed: 2024-09-22.

[5] Howard E. Haber. The volume of a hypersphere. https://scipp.ucsc.edu/~haber/


archives/physics116A10/volume_10.pdf, 2010.

[6] Thomas C. Hull. Volume of an n-ball. Rose-Hulman Undergraduate Mathematics Journal,


10(2):3–7, 2009.

[7] Mates Mike. Cómo ver una esfera en 5d. YouTube, 2018. Accessed: 2024-09-23.

[8] K. F. Riley, M. P. Hobson, and S. J. Bence. Mathematical Methods for Physics and Engineering.
Cambridge University Press, 3rd edition, 2006.

[9] L. M. Rozman. Surface area and volume of the n-dimensional sphere. https://www.phys.
uconn.edu/~rozman/Courses/P2400_17S/downloads/nsphere.pdf, 2017.

[10] Eric W. Weisstein. Curvilinear coordinates. https://mathworld.wolfram.com/


CurvilinearCoordinates.html, 2023. Accessed: 2024-09-23.

[11] Maxime Wolf. The volume of the n-ball. https://www.maximewolf.com/blog/2024/


The-Volume-of-the-nball/, 2024.

25

También podría gustarte