Universidad de Córdoba
Departamento de Matemáticas y Estadística
Pregrado en Matemáticas
Teorema de Cambio de Variable en Rn
Teoremas y sus demostraciones
Alumnos
César Rico Garcés
Yerlis Enith Pérez Reyes
Mauricio Negrette Berrocal
Jeimy María Acevedo Fontalvo
María Mercedes Furnieles Calderon
Sebastián Enrique Rodríguez Muñoz
Profesor:
Abraham Arenas Tawil
Fecha de entrega: 19 de noviembre de 2024.
Doy mi palabra de que he realizado esta actividad con integridad académica.
Parte I
Teoremas propuestos y resueltos
1
Capítulo 1
Cambio de Variables
Teorema 1.1(Teorema del Cambio de Variable) Sea φ una función estrictamente
creciente y continua tal que
φ([A, B]) = [a, b]
Suponga que α es una función monótona creciente en [a, b] y f ∈ R(α) en [a, b]. Definamos las
funciones g y β en [A, B] por:
β(y) = α(φ(y))
g(y) = f (φ(y))
Entonces g ∈ R(β) en [A, B] y
Z B Z b
g(y)dβ(y) = f (x)dα(x).
A a
Teorema 1.2(Teorema Fundamental del Cálculo) Sea f ∈ R en [a, b]. Si existe una
función F diferenciable en [a, b] tal que
F ′ (x) = f (x), x ∈ [a, b]
entonces Z b
f (x)dx = F (b) − F (a)
a
Sea φ una función continuamente diferenciable que asigna [α, β] a [a, b]. Para cada función
continua f en [a, b], tenemos la siguiente fórmula de cambio de variables:
Z β Z φ(α)
′
f (φ(y))φ (y) dy = f (x) dx. (2,1) (1.1)
α φ(β)
La fórmula proviene de una aplicación directa del Teorema Fundamental del Cálculo. Sea F (x)
una función primitiva de f , es decir, F ′ = f . Considere la función compuesta g(y) = F (φ(y)).
Por la regla de la cadena,
2
g ′ (y) = F ′ (φ(y))φ′ (y) = f (φ(y))φ′ (y).
Según el teorema fundamental del cálculo,
Z β Z β
′
g(β) − g(α) = g (y) dy = f (φ(y))φ′ (y) dy.
α α
Por otro lado,
Z φ(α)
g(β) − g(α) = F (φ(β)) − F (φ(α)) = f (x) dx.
φ(β)
Por tanto, la fórmula es válida.
Cuando φ asigna [α, β] biyectivamente a [a, b], φ es estrictamente creciente con φ(α) = a,
φ(β) = b o es estrictamente decreciente con φ(α) = b, φ(β) = a. En el primer caso φ′ no es
negativo o en el segundo caso no es positivo. Entonces (2.1) se convierte en la fórmula
Z β Z b
f (φ(y))|φ′ (y)| dy = f (x) dx. (2,2) (1.2)
α a
En las dos primeras secciones ampliaremos (1.2) a una dimensión superior. En las dos últimas
secciones consideramos una extensión de (1.1).
La fórmula del cambio de variables
Sean D1 y D2 dos regiones en Rn . (Aquí nos ocupamos principalmente de n = 2,3.) Una
función biyectiva de D1 a D2 se llama C 1 -difeomorfismo si y su inverso ambos son
continuamente diferenciables.
Para una función diferenciable Φ de D1 a Rn , su matriz jacobiana ∇Φ está dada por
(∂Φi /∂xj ) , i, j = 1, 2, · · · , n, es decir,
∂Φ1 ∂Φ1
··· ···
∂x1 ∂xn
∂Φ2 ∂Φ2
··· ···
∂x1 ∂xn
··· ··· ··· ···
∂Φn ∂Φn
··· ···
∂x1 ∂xn
El determinante de la matriz jacobiana se llama jacobiano de Φ. Se denotará por JΦ .
3
Según el teorema de la función inversa, una función C 1 de una región D en Rn a Rn que es
uno- a uno y cuyo jacobiano siempre es diferente de cero, establece un C 1 -diffeomorfismo entre
D y su imagen Φ(D). Este hecho se utilizará implícita y frecuentemente a continuación.
Teorema 1.3(Fórmula de cambio de variables) Sea Φ un difeomorfismo C 1 de D1 a D.
Para cualquier función continua f en D,
Z Z
f (x)dx = f (Φ(y)) |JΦ (y)| dy (2.3)
D D1
Aquí dx y dy se refieren a la integración sobre una región n-dimensional. Para n = 2, en
nuestra notación habitual, esta fórmula se lee como,
ZZ ZZ
f (x, y)dA(x, y) = f (Φ(u, v)) |JΦ (u, v)| dA(u, v)
D D1
y para n = 3,
ZZZ ZZZ
f (x, y, z)dV (x, y, z) = f (Φ(u, v, w)) |JΦ (u, v, w)| dV (u, v, w).
Ω Ω1
Las fórmulas de integración para las coordenadas polares, las coordenadas cilíndricas y las
coordenadas esféricas son casos especiales de este teorema.
Coordenadas polares
En el caso de las coordenadas polares, tomamos n = 2 y Φ(r, θ) = (r cos θ, r sin θ). Entonces
cos θ −r sen θ
JΦ = det(J) =
sen θ r cos θ
= (cos θ)(r cos θ) − (−r sen θ)(sen θ)
= r cos2 θ + r sen2 θ
= r(cos2 θ + sen2 θ) = r > 0,
por lo que la fórmula (2.3) se convierte en
ZZ ZZ
f (x, y)dA(x, y) = f (r cos θ, r sin θ)rdA(r, θ)
D D1
4
Ejemplo de Cambio de Variables en Coordenadas Polares
Supongamos que queremos integrar la función f (x, y) = x2 + y 2 sobre el círculo de radio 1
centrado en el origen.
Figura 1.1: Ejemplo 1
En coordenadas cartesianas, la integral sería:
ZZ
(x2 + y 2 ) dx dy
x2 +y 2 ≤1
Para cambiar a coordenadas polares:
La región de integración en coordenadas polares es 0 < r ≤ 1 y 0 ≤ θ < 2π. 2. La función
f (x, y) en coordenadas polares es f (r, θ) = r2 .
La integral se convierte en:
ZZ ZZ
2
r · r dr dθ = r3 dr dθ
0≤r≤1, 0≤θ<2π 0≤r≤1, 0≤θ<2π
Evaluamos la integral:
5
Z 2π Z 1 Z 1
dθ r3 dr = 2π r3 dr
0 0 0
La integral interior es:
" #1
Z 1
3 r4 1
r dr = =
0 4 0
4
Por lo tanto, la integral total es:
1 π
2π · =
4 2
Este es el valor de la integral de x2 + y 2 sobre el círculo unitario usando el cambio de variables
a coordenadas polares.
Coordenadas cilíndricas
En el caso de las coordenadas cilíndricas, tomamos n = 3 y Φ(r, θ, z) = (r cos θ, r sin θ, z).
Entonces
cos θ −r sen θ 0
JΦ = det(J) = sen θ r cos θ 0
0 0 1
= (cos θ)(r cos θ) − (−r sen θ)(sen θ)
= r cos2 θ + r sen2 θ
= r(cos2 θ + sen2 θ) = r > 0,
y (2.3) se convierte en
ZZZ ZZ
f (x, y, z)dV = f (r cos θ, r sin θ, z)rdV (r, θ, z)
Ω Ω1
Ejemplo de Cambio de Variables en Coordenadas Cilindricas
Integrar f (x, y, z) = x2 + y 2 + z 2 sobre un cilindro de radio 1 y altura 2, centrado en el eje z.
Queremos integrar la función f (x, y, z) = x2 + y 2 + z 2 sobre un cilindro de radio 1 y altura 2,
centrado en el eje z. En coordenadas cartesianas, la región de integración es:
0 < x2 + y 2 ≤ 1 y −1≤z ≤1
En coordenadas cilíndricas, tenemos:
6
Figura 1.2: Ejemplo 2
x = r cos θ
y = r sin θ
z=z
La función f (x, y, z) = x2 + y 2 + z 2 se convierte en:
f (r, θ, z) = (r cos θ)2 + (r sin θ)2 + z 2 = r2 (cos2 θ + sin2 θ) + z 2 = r2 + z 2
En coordenadas cilíndricas, los límites son:
0<r≤1
0 ≤ θ < 2π
−1 ≤ z ≤ 1
Como calculamos previamente, el determinante del jacobiano es r. Por lo tanto, el elemento de
volumen dV en coordenadas cilíndricas es r dr dθ dz.
7
La integral original se convierte en:
ZZZ ZZZ
(x2 + y 2 + z 2 ) dx dy dz = (r2 + z 2 ) · r dr dθ dz
V V′
Evaluamos la integral en coordenadas cilíndricas:
ZZZ Z 1 Z 2π Z 1
2 2
(r + z ) · r dr dθ dz = dz dθ (r3 + rz 2 ) dr
V′ −1 0 0
Descomponemos la integral:
Z 1 Z 2π Z 1 Z 1
3 2
dz dθ r dr + z r dr
−1 0 0 0
Evaluamos las integrales internas:
" #1
Z 1
3 r4 1
r dr = =
0 4 0
4
" #1
Z 1
r2 1
r dr = =
0 2 0
2
Así, tenemos:
Z 1 Z 2π
1 1
dz dθ + z2 ·
−1 0 4 2
Evaluamos las integrales restantes:
Z 2π
dθ = 2π
0
!
Z 1
1 z2 Z 1
1 Z 1 2
z
+ dz = dz + dz
−1 4 2 −1 4 −1 2
Para la primera integral:
Z 1
1 1 1 1 1
dz = [z]1−1 = (1 − (−1)) = · 2 =
−1 4 4 4 4 2
Para la segunda integral:
Z 1
z2 1Z 1 2
dz = z dz
−1 2 2 −1
8
" #1
Z 1
2 z3 1 1 1 2
z dz = = (1 − (−1)3 ) = (1 − (−1)) = · 2 =
−1 3 −1
3 3 3 3
1 2 1
· =
2 3 3
Sumando estas dos partes:
1 1 3 2 5
+ = + =
2 3 6 6 6
Finalmente, multiplicamos por 2π:
5 10π 5π
2π · = =
6 6 3
Por lo tanto, la integral de x2 + y 2 + z 2 sobre el cilindro es:
ZZZ
5π
(x2 + y 2 + z 2 ) dx dy dz =
V 3
Coordenadas esféricas
En el caso de las coordenadas esféricas, tomamos n = 3 y
Φ(ρ, φ, θ) = (ρ sen φ cos θ, ρ sen φ sen θ, ρ cos φ), φ ∈ (0, π), θ ∈ [0, 2π)
Entonces
sen φ cos θ ρ cos φ −ρ sen φ sen θ
JΦ = det(J) = sen φ sen θ ρ cos φ sen θ ρ sen φ cos θ
cos φ −ρ sen φ 0
= sen φ cos θ(−(ρ sen φ cos θ · −ρ sen φ)) − ρ cos φ cos θ(−(ρ sen φ cos θ · cos φ))
+ (−ρ sen φ sen θ)(sen φ sen θ · −ρ sen φ − ρ cos φ sen θ · cos φ)
= sen φ cos θ(ρ2 sen2 φ cos θ) + ρ2 cos2 φ cos2 θ sen φ + ρ2 sen φ sen2 φ sen2 θ
+ ρ2 sen2 θ sen φ cos2 φ
= ρ2 sen φ(sen2 φ cos2 θ + cos2 φ cos2 θ + sen2 φ sen2 θ + sen2 θ cos2 φ)
= ρ2 sen φ(sen2 φ(sen2 θ + cos2 θ) + cos2 φ(sen2 θ + cos2 θ))
= ρ2 sen φ(sen2 θ + cos2 θ)(sen2 φ + cos2 φ)
= ρ2 sen φ > 0,
y (2.3) se convierte en
ZZZ ZZZ
f (x, y, z)dV = f (ρ sin φ cos θ, ρ sin φ sin θ, ρ cos φ)ρ2 sin φdV (ρ, φ, θ).
Ω Ω1
9
Ejemplo de Cambio de Variable en coordenadas esfericas
Supongamos que queremos calcular la integral de la función f (x, y, z) = x2 + y 2 + z 2 sobre la
esfera de radio R.
Figura 1.3: Ejemplo 3
En coordenadas esféricas, las transformaciones son:
x = ρ sin φ cos θ
y = ρ sin φ sin θ
z = ρ cos φ
La función f (x, y, z) en coordenadas esféricas es:
f (ρ, φ, θ) = (ρ sin φ cos θ)2 + (ρ sin φ sin θ)2 + (ρ cos φ)2
= ρ2 sin2 φ(cos2 θ + sin2 θ) + ρ2 cos2 φ
= ρ2 (sin2 φ + cos2 φ)
= ρ2
En coordenadas esféricas:
0<ρ≤R
10
0<φ<π
0 ≤ θ < 2π
Estos límites cubren toda la esfera de radio R.
El jacobiano para el cambio de coordenadas esféricas a cartesianas es ρ2 sin φ.
La integral original se convierte en:
ZZZ ZZZ
2 2 2
(x + y + z ) dx dy dz = ρ2 · ρ2 sin φ dρ dφ dθ
V V ′
Donde V ′ representa la región en coordenadas esféricas.
Integramos primero respecto a ρ, luego φ, y finalmente θ:
1. Integración respecto a ρ:
" #R
Z R
ρ5 R5
ρ4 dρ = =
0 5 0
5
2. Integración respecto a φ: Z π
sin φ dφ = [− cos φ]π0 = 2
0
3. Integración respecto a θ: Z 2π
dθ = 2π
0
Multiplicamos estas integrales:
!
ZZZ
R5 4πR5
ρ2 dρ dφ dθ = · 2 · 2π =
V′ 5 5
Por lo tanto, la integral de x2 + y 2 + z 2 sobre la esfera de radio R en coordenadas esféricas es
4πR5
5
.
Prueba de la fórmula
En esta sección presentamos una prueba más detallada de (2.3). Seguiremos el tratamiento del
profesor Tom Wan. Se pueden encontrar más discusiones en el libro de M. Spivak: Calculus on
Manifolds. Para prepararnos para ello, necesitamos las siguientes proposiciones.
Proposición 2.2. Sea Φ1 : D2 → D1 y Φ2 : D3 → D2 sean dos C 1 -funciones. Entonces Φ =
Φ1 ◦ Φ2 satisface
∇Φ = ∇Φ1 · ∇Φ2
11
(producto matricial) y
JΦ = JΦ1 JΦ2
Prueba. Sean x = Φ1 (y) y y = Φ2 (z) de modo que x = Φ(z). Tenemos
Φi (z) = (Φ1 ◦ Φ2 )i (z) = Φ1i (Φ21 (z), · · · , Φ2n (z))
Por la regla de la cadena,
∂Φi X ∂Φ1i ∂Φ2k
(z) = (y) (z)
∂zj k ∂yk ∂zj
que es precisamente el producto matricial
∇Φ(z) = ∇Φ1 (y) · ∇Φ2 (z)
La segunda fórmula se deriva de la propiedad del determinante: det AB = det A det B. En
efecto;
JΦ = det(∇Φ(z)) = det(∇Φ1 (y) · ∇Φ2 (z))
= det(∇Φ1 (y)) det(∇Φ2 (y))
= JΦ1 JΦ2
Proposición 2.3. Sea Φ : D1 → D un difeomorfismo C 1 . Entonces
JΦ−1 JΦ = 1
En particular, JΦ ̸= 0 en D1 .
Prueba. Tenemos Φ−1 (Φ(x)) = x. Por la Proposición 2.2 y usando el hecho de que la matriz
jacobiana de la función identidad es la matriz identidad, ∇Φ−1 · ∇Φ es igual a la matriz
identidad, y la fórmula sigue tomando determinante de ambos lados.
Nótese que
I(x) = (Φ−1 ◦ Φ)(x).
Luego, por la Proposición 2.2 tenemos que
∇I(x) = ∇Φ−1 (x) · ∇Φ(x) ⇔ I(x) = ∇Φ−1 (x) · ∇Φ(x)
de donde
det(I(x)) = det(∇Φ−1 (x) · ∇Φ(x)) ⇔ 1 = det(∇Φ−1 (x)) det(∇Φ(x))
⇔ 1 = JΦ−1 JΦ
12
Así, JΦ−1 JΦ = 1 como se quería probar.
Proposición 2.4. Sean Φ1 : D2 → D1 y Φ2 : D3 → D2 dos C 1 -difeomorfismos y
Φ = Φ1 ◦ Φ2 : D3 → D1 . Supongamos que (1.3) es válido para Φ1 y Φ2 . Entonces también es
válido para Φ.
Prueba. Sean f y g continuos en D1 y D2 respectivamente. Por suposición, tenemos
Z Z
f (x)dx = f (Φ1 (y)) |JΦ1 | (y)dy
D1 D2
y
Z Z
g(y)dy = g (Φ2 (z)) |JΦ2 | (z)dz
D2 D3
Como f (Φ1 (y)) |JΦ1 | (y) es continua en D2 , tomándolo como g, tenemos
Z Z
f (x)dx = f (Φ1 (y)) |JΦ1 | (y)dy
D1 D2
Z
= f (Φ1 (Φ2 (z))) |JΦ1 (Φ2 (z))| |JΦ2 | (z)dz
D3
Z
= f (Φ(z)) |JΦ | (z)dz. (Proposición 2.2)
D3
Ahora limitémonos a n = 2.
Proposición 2.5. La fórmula de cambio de variables (1.3) se cumple en los dos casos
siguientes:
(a) Φ es un C 1 -difeomorfismo de la forma Φ(u, v) = (φ(u, v), v) en D; y
(b) Φ es un C 1 -difeomorfismo de la forma Φ(u, v) = (u, ψ(u, v)).
Prueba. Probamos (a) mientras que (b) se puede demostrar de manera similar. Tomaremos D
13
como un rectángulo [a, b] × [c, d]. En primer lugar,
∂Φ1 ∂Φ1
∂x1 ∂x2
JΦ =
∂Φ2 ∂Φ2
∂x1 ∂x2
∂φ ∂φ
∂u ∂v
=
∂v ∂v
∂u ∂v
∂φ ∂φ
= ∂u ∂v
0 1
∂φ
=
∂u
Según la Proposición 2.3, ya sea ∂φ/∂u > 0 o ∂φ/∂u < 0 en D. Supongamos que es lo
primero. Bajo Φ, la línea vertical (u, v), donde u ∈ [a, b] es fija, se asigna a (φ(u, v), v). Esta es
la gráfica de φ(u, ·) sobre [c, d]. Dado que para cada v fijo, φu > 0, φ (u1 , v) < φ (u2 , v) para
u1 < u2 , la imagen de D bajo φ tiene la forma:
{(x, y) : φ(a, y) ≤ x ≤ φ(b, y), y = v ∈ [c, d]}
Según el teorema de Fubini,
Teorema 1.4(Teorema de Fubini) Sean Q = A × B donde A y B son rectángulos de Rk y
Rn respectivamente, y f : Q → R una función acotada, escrita como f (x, y) para x ∈ A y
y ∈ B. Para cada x ∈ A se definen las integrales
Z Z
I(x) = f (x, y) dy, I(x) = f (x, y) dy
y∈B y∈B
Si f es integrable sobre Q, entonces I y I son integrables sobre A y además
Z Z Z
f= I(x) dx = I(x) dx
Q x∈A x∈A
es decir
Z Z Z Z Z
f= f (x, y) dy dx = f (x, y) dx dy
Q x∈A y∈B y∈B x∈A
14
ZZ Z d Z φ(b,y)
f (x, y)dA = f (x, y)dxdy
Φ(D) c φ(a,y)
Z dZ b ∂φ
= f (φ(u, y), y)
dudy
c a ∂u
Z dZ b
∂φ
= f (φ(u, v), v) dudv
c
ZZ a ∂u
= f (Φ(u, v)) |JΦ | (u, v)dA(u, v)
D
Tenga en cuenta que en la segunda línea, hemos utilizado el cambio de variables x = φ(u, y)
para que dx = φu du.
Cuando ∂φ/∂u < 0, Φ(D) se convierte
{(x, y) : φ(b, y) ≤ x ≤ φ(a, y), y = v ∈ [c, d]}
De manera similar a lo anterior, tenemos
ZZ Z d Z φ(a,y)
f (x, y)dA = f (x, y)dxdy
Φ(D) c φ(b,y)
Z dZ a ∂φ
= f (φ(u, y), y)
dudy
c b ∂u
Z dZ b
∂φ
= f (φ(u, v), v) dudv
c a ∂u
ZZ
= f (Φ(u, v)) |JΦ | (u, v)dA(u, v)
D
Ahora demostramos la fórmula de cambio de variables (1.3) para un Φ general. Por simplicidad,
solo consideraremos n = 2 y tomaremos D1 sea un rectángulo. El caso general es esencialmente
el mismo. En primer lugar, dado que Φ = (φ1 , φ2 ) es un difeomorfismo C 1 , su jacobiano
∂φ1 ∂φ2 ∂φ1 ∂φ2
JΦ = −
∂u ∂v ∂v ∂u
nunca se anula. Por lo tanto, ∂φ1 /∂u o ∂φ1 /∂v no es cero en ningún punto. Podemos arreglar
una partición en D tan fina que cada subrectángulo Rij pertenezca a A o B donde
( )
∂φ1
A = Rij : > 0 in Rij
∂u
y
15
( )
∂φ1
B = Rij : < 0 in Rij
∂v
Usando
ZZ
f (Φ(u, v)) |JΦ | (u, v)dA
D
X ZZ X ZZ
= f (Φ(u, v)) |JΦ | (u, v)dA + f (Φ(u, v)) |JΦ | (u, v)dA
Rij ∈A Rij Rij ∈B Rij
vemos que basta con establecer la fórmula bajo el supuesto adicional R ∈ A o R ∈ B. (Hemos
escrito R para Rij por simplicidad).
Consideremos solo R ∈ A. (El otro caso se puede manejar de manera similar). Consideramos
las funciones Φ1 (u, v) = (φ1 (u, v), v) y Φ2 (s, t) = (s, h(s, t)) donde h(s, t) = φ2 Φ−1
1 (s, t) .
Dado que JΦ1 = ∂φ1 /∂u ̸= 0, Φ1 es un C 1. -diffeomorfismo de R en su imagen En particular, la
función inversa Φ−1
1 existe ahora.
Φ2 (Φ1 (u, v)) = Φ2 (φ1 (u, v), v)
= (φ1 (u, v), h (φ1 (u, v), v))
= (φ1 (u, v), h (Φ1 (u, v))
= (φ1 (u, v), φ2 (u, v))
= Φ(u, v)
Por las Proposiciones 2.4 y 2.5, vemos que (2.3) se cumple para Φ.
16
Parte II
Conclusiones Generales
17
El teorema del cambio de variable en Rn es una herramienta fundamental para la evaluación de
integrales múltiples, ya que permite simplificar el proceso de integración mediante la
transformación tanto del dominio como de la función integranda. Este procedimiento se realiza
utilizando una función de transformación adecuada y el determinante jacobiano, que ajusta los
cambios en la escala y la orientación introducidos por la transformación. Al expresar la integral
en términos de las nuevas variables, no solo se facilita el cálculo, sino que también se obtiene
una mayor flexibilidad para abordar problemas complejos. Tiene aplicaciones amplias y
profundas en matemáticas y física, donde se utiliza para analizar sistemas en coordenadas más
convenientes, como las coordenadas polares, cilíndricas o esféricas, optimizando la resolución
de problemas en áreas como la mecánica.
18