Cambio de variable
Contents
4.4. Cambio de variable#
Antes de empezar esta sección, te recomendamos que pongas al día tus conocimientos para funciones de una variable. Puedes consultar:
En este mismo Jupyter Book, la sección 4.1., Sección Revisión una variable: cálculo de primitivas.
Wikipedia: https://es.wikipedia.org/wiki/Métodos_de_integración#Integración_por_cambio_de_variable
…
4.4.1. Matriz jacobiana y cambio de variable#
El método de cambio de variable esencialmente revierte la regla de la cadena para derivadas. En otras palabras, nos ayuda a integrar composición de funciones. A menudo, este método es útil para reducir la complejidad del cálculo de algunas integrales. Hay muchos casos en los que con las variables iniciales no sabremos resolver la integral, ya sea por la función a integrar, o por la complejidad de la expresión de los intervalos de integración. En estos casos, buscaremos cambios de variables que nos faciliten la resolución de la integral.
Theorem (Cambio de variable en integrales dobles)
Sea \(R\) una región en el plano \(XY\) y sea \(S\) una región en el plano \(UV\). Sean \(\mathbf{F}: S \subset \mathbb{R}^2 \rightarrow \mathbb{R}^2\) y \(G: R \subset \mathbb{R}^2 \rightarrow \mathbb{R}\) dos funciones tales que \(\mathbf{F}\) es inyectiva y \(\mathbf{F}(u,v) = (F_1(u,v),F_2(u,v)) \subset R\), donde \(F_1\) y \(F_2\) tienen primeras derivadas parciales continuas. Si \(G\) es continua en \(R\) y \(\det \mathrm{J}\mathbf{F}(u,v) \neq 0\) en \(S\), entonces
donde \(\mathrm{J}\mathbf{F}\) es la matriz jacobiana de \(\mathbf{F}\).
Example
Sea \(R\) la región delimitada por las líneas
![../../_images/5.4_Ejemplo_1.png](../../_images/5.4_Ejemplo_1.png)
Evalúa la integral
Solución:
Como es complejo establecer los intervalos de integración sobre \(R\), vamos a definir una región \(S\) más sencilla. Para ello, consideramos el siguiente cambio de variable:
Ahora sí, fácilmente podemos deducir los intervalos de integración de las nuevas variables a partir de la nueva región \(S = \Big{\{}(u,v) \in \mathbb{R}^2 \Big{/} 1 \le u \le 4, -4 \le v \le 0\Big{\}}\). Además, despejando \(x\) e \(y\) a partir del cambio de variable propuesto, obtenemos la función \(\mathbf{F}: S \subset \mathbb{R}^2 \rightarrow \mathbb{R}^2\) tal que
Como
por el resultado teórico anterior, obtenemos
Podemos usar Sympy
para que nos confirme lo bien que trabajamos:
import sympy as sp
x, y, u, v = sp.symbols('x y u v', real=True) # define las variables simbólicas x, y, u, v
# Definimos las funciones F y G como matrices
F = sp.Matrix([ 1/3*(2*u+v), 1/3*(u-v) ])
G = sp.Matrix([ 9*x*y ])
# Definimos la nueva función a integrar
GoF_expr = G.subs(x,1/3*(2*u+v)).subs(y,1/3*(u-v))
# Calculamos el determinante de la matriz jacobiana asociada a F
det_jac_F = sp.det( F.jacobian([u,v]) )
display(det_jac_F)
# Calculamos la integral sobre la nueva región de integración
res = sp.integrate(GoF_expr*abs(det_jac_F), (v, -4, 0), (u, 1, 4))
display('Valor de la integral: ')
display(sp.simplify(res[0]))
'Valor de la integral: '
4.4.2. Integrales dobles en coordenadas polares#
En \(\mathbb{R}^2\), si la función a integrar o la región de integración tienen alguna clase de simetría radial, a menudo resulta conveniente aplicar la transformación de coordenadas cartesianas a coordenadas polares.
Theorem (Cambio de variable a la forma polar)
Sea \(R\) una región en el plano \(XY\) y sea \(S\) una región en el plano \(r\theta\). Sean \(\mathbf{F}: S \subset \mathbb{R}^2 \rightarrow \mathbb{R}^2\) y \(G: R \subset \mathbb{R}^2 \rightarrow \mathbb{R}\) dos funciones tales que \(\mathbf{F}(r,\theta) = (r\cos(\theta),r\sin(\theta)) \subset R\). Si \(S = \{(r,\theta) \in \mathbb{R}^+ \times [0,2\pi] | \alpha \le \theta \le \beta, 0 \le g_1(\theta) \le r \le g_2(\theta)\}\), donde \(g_1\) y \(g_2\) son continuas en \([\alpha,\beta]\) y \(G\) es continua en \(R\), entonces
Fíjate que, en este caso, como la aplicación que define el cambio de variable es \(\mathbf{F}(r,\theta) = (r\cos(\theta),r\sin(\theta)) \subset R\), su matriz jacobiana será
y, por tanto, el valor absoluto del determinante de esta matriz será
Example
Sea \(R\) la región anular que se encuentra entre los dos círculos
![../../_images/5.4_Ejemplo_2.png](../../_images/5.4_Ejemplo_2.png)
Evalúa la integral
Solución:
Utilizando el cambio de variable a la forma polar, tenemos \(\mathbf{F}: S \subset \mathbb{R}^2 \rightarrow \mathbb{R}^2\)
donde los intervalos de integración de las nuevas variables son:
Por el resultado teórico anterior, obtenemos
Vamos a hacer que trabaje un poco nuestro oráculo de confianza. Con la ayuda de Sympy
, vamos a comprobar que los cálculos realizados en el ejemplo anterior son correctos:
import sympy as sp
x, y, r, th = sp.symbols('x y r th', real=True) # define las variables simbólicas x, y, r, th
# Definimos las funciones F y G como matrices
F = sp.Matrix([ r*sp.cos(th), r*sp.sin(th) ])
G = sp.Matrix([ x**2 + y ])
# Definimos la nueva función a integrar
GoF_expr = G.subs(x,r*sp.cos(th)).subs(y,r*sp.sin(th))
# Calculamos el determinante de la matriz jacobiana asociada a F
det_jac_F = sp.det( F.jacobian([r,th]) )
# Calculamos la integral sobre la nueva región de integración
res = sp.integrate(GoF_expr*abs(det_jac_F), (r, 1, sp.sqrt(5)), (th, 0, 2*sp.pi))
# res = sp.integrate(GoF_expr[0]*abs(det_jac_F), (r, 1, sp.sqrt(5)), (th, 0, 2*sp.pi))
display('Valor de la integral: ')
display(sp.simplify(res[0]))
# display(sp.simplify(res))
'Valor de la integral: '