Integral múltiple

De Wikipedia, la enciclopedia libre
Saltar a: navegación, búsqueda

Una integral múltiple es un tipo de integral definida aplicada a funciones de más de una variable real, por ejemplo, f(x,y) ó f(x,y,z).

La doble integral como el volumen bajo una superficie. La región rectangular abajo de la figura es el dominio de integración, mientras que la superficie es la gráfica de la función de dos variables de la integral.

Introducción[editar]

De la misma manera en que la integral de una función positiva f(x) de una variable definida en un intervalo puede interpretarse cómo el área entre la gráfica de la función y el eje x en ese intervalo, la doble integral de una función positiva f(x,y) de dos variables, definida en una región del plano xy, se puede interpretar como el volumen entre la superficie definida por la función y el plano xy en ese intervalo. Al realizar una "integral triple" de una función f(x,y,z) definida en una región del espacio xyz, el resultado es un hipervolumen, sin embargo es bueno notar que si f(x,y,z)=1 el resultado se puede interpretar como el volumen de la región de integración. Para integrales de órdenes superiores, el resultado geométrico corresponde a hipervolúmenes de dimensiones cada vez superiores.

La manera más usual de representar una integral múltiple es anidando signos de integración en el orden inverso al orden de ejecución (el de más a la izquierda es el último en ser calculado), seguido de la función y los diferenciales en orden de ejecución. El dominio de integración se representa sobre cada signo de integral, o a menudo es abreviado por una letra en el signo de integral de más a la derecha:


\iint \ldots \int_\mathbf{D}\;f(x_1,x_2,\ldots,x_n) \;\mathbf{d}x_1 \mathbf{d}x_2\!\ldots\mathbf{d}x_n

Es importante destacar que no es posible calcular la función primitiva o antiderivada de una función de más de una variable por lo que las integrales múltiples indefinidas no existen.

Definición[editar]

Una forma relativamente sencilla de definir las integrales múltiples es mediante su representación geométrica como la magnitud del espacio entre el objeto definido por la ecuación x_{n+1} = f(x_1,...,x_n) y una región T en el espacio definido por los ejes de las variables independientes de la función f (si T es una región cerrada y acotada y f está definida en ésta). Por ejemplo, si n=2, el volumen situado entre la superficie definida por x_3 = f(x_1,x_2) y una región T en el plano x_1x_2 es igual a alguna integral doble, si es que, como se mencionó, f está definida en T.

T puede dividirse en una partición interior \Delta formada por m subregiones rectangulares sin solapamiento que estén completamente contenidas en T. La norma ||\Delta|| de esta partición está dada por la diagonal más larga en las m subregiones.

Si se toma un punto (x_{1i}, x_{2i},..., x_{ni}) que esté contenido dentro de la subregión con dimensiones \Delta x_{1i}\Delta x_{2i} ...\Delta x_{ni} para cada una de las m subregiones de la partición, se puede construir un espacio con una magnitud aproximada a la del espacio entre el objeto definido por x_{n+1} = f(x_1,...,x_n) y la subregión i. Este espacio tendrá una magnitud de:

f(x_{1i}, x_{2i},\ldots, x_{ni})\Delta A_{i} = f(x_{1i}, x_{2i},\ldots, x_{ni})\Delta x_{1i}\Delta x_{2i} \ldots \Delta x_{ni}

Entonces se puede aproximar la magnitud del espacio entero situado entre el objeto definido por la ecuación x_{n+1} = f(x_1,...,x_n) y la regiónTmediante la suma de Riemann de las magnitudes de los m espacios correspondientes a cada una de las subregiones:

\sum_{i=1}^m f(x_{1i}, x_{2i},\ldots, x_{ni})\Delta x_{1i}\Delta x_{2i} \ldots \Delta x_{ni}

Esta aproximación mejora a medida que el número m de subregiones se hace mayor. Esto sugiere que se podría obtener la magnitud exacta tomando el límite. Al aumentar el número de subregiones disminuirá la norma de la partición:

\lim_{m \to \infty}\sum_{i=1}^m f(x_{1i}, x_{2i},\ldots, x_{ni})\Delta x_{1i}\Delta x_{2i} \ldots \Delta x_{ni} = \lim_{\|\Delta\| \to 0}\sum_{i=1}^m f(x_{1i}, x_{2i},\ldots, x_{ni})\Delta x_{1i}\Delta x_{2i} \ldots \Delta x_{ni}

El significado riguroso de éste último límite es que el límite es igual L si y sólo si para todo \varepsilon > 0 existe un \delta > 0 tal que

\left | L -\sum_{i=1}^m f(x_{1i}, x_{2i},\ldots, x_{ni})\Delta x_{1i}\Delta x_{2i} \ldots \Delta x_{ni} \right \vert  < \varepsilon

para toda partición \Delta de la regiónT(que satisfaga ||\Delta|| < \delta), y para todas las elecciones posibles de (x_{1i}, x_{2i},..., x_{ni}) en la iésima subregión. Esto conduce a la definición formal de una integral múltiple:

Sifestá definida en una región cerrada y acotadaTdel definido por los ejes de las variables independientes de f, la integral defsobreTestá dada por:
 \iint \ldots \int_{T} \,f(x_1,\ldots,x_n) dx_1 \ldots \,dx_n = \lim_{\|\Delta\| \to 0}\sum_{i=1}^m f(x_{1i}, x_{2i},\ldots, x_{ni})\Delta x_{1i}\Delta x_{2i} \ldots \Delta x_{ni}
siempre que el límite exista. Si el límite existe se dice quefes integrable con respecto a T.

Propiedades[editar]

Las integrales múltiples comparten muchas de las propiedades de las integrales simples. Si f y g son funciones continuas en una región cerrada y acotada D en un espacio Rn y c una constante con respecto a todas las variables involucradas entonces se puede demostrar que:

1.

\iint \ldots \int_D\;c f(x_1,x_2,\ldots,x_n) \;dx_1 dx_2 \ldots dx_n 
= c \iint \ldots \int_\mathbf{D}\;f(x_1,x_2,\ldots,x_n)
 \;dx_1 dx_2 \ldots dx_n

2.

\iint \ldots \int_{D} \,[f(x_1,x_2,\ldots,x_n) \pm g(x_1,x_2,\ldots,x_n)]
 dx_1\,dx_2 \ldots \,dx_n =
\iint \ldots \int_{D} \,f(x_1,x_2,\ldots,x_n) dx_1\,dx_2 
\ldots \,dx_n \pm \iint \ldots \int_{D} \,g(x_1,x_2,\ldots,x_n) dx_1\,
dx_2 \ldots \,dx_n

3.

Si \quad f(x_1,x_2,\ldots,x_n) \geq 0, entonces:
\iint \ldots \int_{D} \,f(x_1,x_2,\ldots,x_n) dx_1\,dx_2
 \ldots \,dx_n \geq 0

4.

Si \quad f(x_1,x_2,\ldots,x_n) \geq g(x_1,x_2,\ldots,x_n), entonces:
\iint \ldots \int_{D} \,f(x_1,x_2,\ldots,x_n) dx_1\,dx_2 \ldots
 \,dx_n \geq \iint \ldots \int_{D} \,g(x_1,x_2,\ldots,x_n) dx_1\,dx_2 
\ldots \,dx_n

5.

Sea D la unión entre dos regiones, D1 y D2, que no solapan entre sí, entonces:
\iint \ldots \int_{D} \,f(x_1,x_2,\ldots,x_n) dx_1\,
dx_2 \ldots \,dx_n =
\iint \ldots \int_{D_1} \,f(x_1,x_2,\ldots,x_n) dx_1
\,dx_2 \ldots \,dx_n + \iint \ldots \int_{D_2} \,f(x_1,x_2,\ldots,x_n)
 dx_1\,dx_2 \ldots \,dx_n


Integrales múltiples e Integrales iteradas[editar]

Las integrales múltiples están estrechamente relacionadas con las integrales iteradas, mismas que son necesarias para resolver las integrales múltiples. La diferencia entre integrales múltiples e iteradas consiste en que una se refiere al concepto matemático de integral (aplicado a varias variables) y otra al procedimiento por el cual se resuelve la integral múltiple. Si la expresión

\int_a^b\int_c^d f(x,y)\,dy\,dx

se refiere a una integral iterada, la parte externa

\int_a^b \cdots \, dx

es la integral con respecto a x de la función de x:

g(x)=\int_c^d f(x,y)\,dy.

Una integral doble, en cambio está definida con respecto a un área en el plano xy. La integral doble existe si y sólo si las dos integrales iteradas existen y son iguales. Es decir, si la integral doble existe, entonces es igual a la integral iterada, sin importar si el orden de integración es dydx ó dxdy, y por lo general uno la calcula calculando una sola de estas. Sin embargo, a veces las dos integrales iteradas existen sin ser iguales y en este caso no existe la integral doble, ya que se tiene:

\int_a^b\int_c^d f(x,y)\,dy\,dx \neq \int_c^d\int_a^b f(x,y)\,dx\,dy.

De una manera más formal, el Teorema de Fubini afirma que

\int_{A\times B} |f(x,y)|\,d(x,y)<\infty,

Esto es, si la integral es absolutamente convergente, entonces la integral doble es igual a la integral iterada.

\int_{A\times B} f(x,y)\,d(x,y)=\int_A\left(\int_B f(x,y)\,dy\right)\,dx=\int_B\left(\int_A f(x,y)\,dx\right)\,dy.

Esto ocurre, cuando f es una función acotada y tanto A como B son regiones acotadas también. Esto se entiende fácilmente pensando que si la función o la región del dominio no están acotadas, la integral múltiple no puede existir.

La notación

\int_{[a,b]\times[c,d]} f(x,y)\,dx\,dy

se puede usar si se desea ser enfático al referirse a una integral doble y no a una iterada.

Métodos de integración[editar]

Funciones constantes[editar]

En el caso de funciones constantes, el resultado es trivial: simplemente multiplíquese el valor de la función constante c por la medida del dominio de integración. Si c = 1, y es integrada a través de una región de R2 esto da el área de la región, mientras que si es una región de R3 da el volumen de la región y así sucesivamente.

Por ejemplo:

D = \{ (x,y)| \ 2 \le x \le 4 \ ; \ 3 \le y \le 6 \} y f(x,y) = c\,\!
Integrando f sobre D:
\int_3^6 \int_2^4 \ c \ dx\, dy = c \cdot \mbox{area}(D) = c \cdot (3 \cdot 2) = c \cdot 6.

Uso de simetrías[editar]

En el caso de un dominio en el que exista simetría al menos respecto de uno de los ejes, y donde la función para integrar contenga al menos una función impar con respecto a esa variable, la integral se vuelve nula (ya que la suma de cantidades iguales con signo opuesto es cero). Por ejemplo

Dada f(x,y) = 2 \ \sin(x) - 3 \ y^3 + 5 y que T = x^2 + y^2 \le 1 es el dominio de integración del disco de radio 1 centrado en el origen.
Usando la propiedad lineal de las integrales, la integral se descompone en tres partes:
\iint_T (2 \ \sin(x) - 3 \ y^3 + 5) \ dx \, dy = \iint_T 2 \ \sin(x) \ dx \, dy - \iint_T 3 \ y^3 \ dx \ dy + \iint_T 5 \ dx \ dy

Ya que tanto 2 sin(x) como 3y3 son funciones impares, y existe simetría tanto con respecto al eje x como con respecto al eje y, las primeras dos integrales se nulifican, de tal forma que la integral original es igual únicamente a la tercera.

\iint_T (2 \ \sin(x) - 3 \ y^3 + 5) \ dx \, dy = \iint_T 5 \ dx \, dy = 5 \pi

Cambio de variables[editar]

A menudo, es útil para reducir la complejidad de la integral cambiar una variable por otra que resulte más cómoda, sin embargo esto exige el cambio de la región de integración, además de añadir un factor de corrección al diferencial conocido como determinante jacobiano (en valor absoluto o módulo). El cambio de una variable por otra es en un sentido geométrico, una transformación desde un espacio hasta otro, y es esta transformación la que exige estos ajustes.

Si se utiliza una transformación que siga la relación:

f(y_1,\ldots,y_n) \rightarrow f(y_1(x_1,x_2,\ldots,x_n),\ldots,y_n(x_1,x_2,\ldots,x_n))

Entonces se puede utilizar el jacobiano de la transformación para simplificar la integral

J=\frac{D(y_1,\ldots,y_n)}{D(x_1,\ldots,x_n)}= \begin{vmatrix} \frac{\partial y_1}{\partial x_1} & \cdots & \frac{\partial y_1}{\partial x_n} \\ \vdots & \ddots & \vdots \\ \frac{\partial y_n}{\partial x_1} & \cdots & \frac{\partial y_n}{\partial x_n} \end{vmatrix}

Integrando la función transformada en el dominio de integración correspondiente a las variables x, y multiplicando por el valor absoluto del determinante jacobiano y por la serie de diferenciales, se obtiene una integral múltiple que es igual a la integral original, si es que esta existe.

\iint \ldots \int_{D} \,f(y_1,\ldots,y_n) dy_1 \ldots \,dy_n = \iint \ldots \int_{T} \,f(x_1,\ldots,x_n) |J| dx_1 \ldots \,dx_n

A continuación se dan algunos ejemplos de estas transformaciones.

Coordenadas Polares[editar]

La transformación de coordenadas rectangulares a polares. Se puede notar que el área de la región polar es distinta que la de la región rectangular, lo que justifica la necesidad del jacobiano. También se puede demostrar que si se consiera \rho = (\rho_1 + \rho_2)/2 (el radio medio), el área de la región polar es efectivamente \rho \Delta \rho \Delta \theta.

En un espacio R2, un dominio de integración que tenga una simetría circular es muchas veces suceptible de ser transformado de coordenadas rectangulares a polares, lo que significa que cada punto P (x, y) del dominio de una integral doble tomará su valor correspondiente en coordenadas polares mediante la siguiente transformación:

f(x,y) \rightarrow f(\rho \ \cos \theta,\rho \ \sin \theta )

Por ejemplo:

Si la función es f(x,y) = x + y\,\!
aplicando la transformación se obtiene la función fácilmente integrable con respecto a \phi y a \rho.
f(\rho, \phi) = \rho \cos \phi + \rho \sin \phi = \rho \ (\cos \phi + \sin \phi ).

Se pueden obtener funciones incluso más simples:

Si la función es f(x,y) = x^2 + y^2\,\!
Uno tiene:
f(\rho, \theta) = \rho^2 (\cos^2 \theta + \sin^2 \theta) = \rho^2\,\!

Si aplica la identidad trigonométrica pitagórica de senos y cosenos.

El determinante jacobiano de la transformación es:

\frac{\partial (x,y)}{\partial (\rho, \theta)} = 
\begin{vmatrix}
\cos \theta & - \rho \sin \theta \\
\sin \theta & \rho \cos \theta 
\end{vmatrix} = \rho

El cual se obtiene insertando las derivadas parciales de x = ρ cos(θ), y = ρ sin(θ) en la primera columna con respecto a ρ y en la segunda con respecto a \theta.

Por lo tanto, una vez transformada la función, y multiplicada por su determinante jacobiano, ésta es igual a la integral original:

\iint_D f(x,y) \ dx\, dy = \iint_T f(\rho \cos \theta, \rho \sin \theta) \rho \, d \rho\, d \theta.

Coordenadas Esféricas[editar]

Gráfica de las coordenadas esféricas.

Cuando existe simetría esférica en un dominio en R3, es posible utilizar una transformación hacia coordenadas esféricas para simplificar una integral triple. La función es transformada por la relación:

f(x,y,z) \longrightarrow f(\rho \cos \theta \sin \phi, \rho \sin \theta \sin \phi, \rho \cos \phi)\,\!

El determinate jacobiano de la transformación es el siguiente:

\frac{\partial (x,y,z)}{\partial (\rho, \theta, \phi)} = 
\begin{vmatrix}
\frac{\partial (x)}{\partial (\rho)} & \frac{\partial (x)}{\partial (\theta)} & \frac{\partial (x)}{\partial (\phi)}\\
\frac{\partial (y)}{\partial (\rho)} & \frac{\partial (y)}{\partial (\theta)} & \frac{\partial (y)}{\partial (\phi)}\\
\frac{\partial (z)}{\partial (\rho)} & \frac{\partial (z)}{\partial (\theta)} & \frac{\partial (z)}{\partial (\phi)}
\end{vmatrix}
= 
\begin{vmatrix}
\cos \theta \sin \phi & - \rho \sin \theta \sin \phi  &  \rho \cos \theta \cos \phi\\
\sin \theta \sin \phi &   \rho \cos \theta \sin \phi & \rho \sin \theta \cos \phi \\
\cos \theta & 0 & - \rho \sin \phi
\end{vmatrix} = - \rho^2 \sin \phi

Tomando el valor absoluto del determinante se obtiene el factor que se debe añadir a la integral.

Por lo tanto los diferenciales dx dy dz se transforman en ρ2 sin(φ) dρ dθ dφ.

Finalmente se obtiene la fórmula de integración:

\iiint_D f(x,y,z) \, dx\, dy\, dz = \iiint_T f(\rho \cos \theta \sin \phi, \rho \sin \theta \sin \phi, \rho \cos \phi) \rho^2 \sin \phi \, d\rho\, d\theta\, d\phi.

Coordenadas Cilíndricas[editar]

Gráfica de las Coordenadas Cilíndricas (Se muestra el ángulo θ como φ).

El uso de coordenadas cilíndricas para transformar una integral triple, es conveniente especialmente cuando el dominio de integración presenta simetría alrededor del eje z. La función se transforma mediante la siguiente relación.

f(x,y,z) \rightarrow f(\rho \ \cos \theta,\rho \ \sin \theta, z)

El determinate jacobiano de la transformación es el siguiente:

\frac{\partial (x,y,z)}{\partial (\rho, \theta,z)} = 
\begin{vmatrix}
\cos \theta & \sin \theta & 0 \\
- \rho \sin \theta & \rho \cos \theta & 0 \\
0 & 0 & 1 \end{vmatrix} = \rho

Por lo tanto, se puede derivar la siguiente fórmula de integración:

\iiint_D f(x,y,z) \, dx\, dy\, dz = \iiint_T f(\rho \cos \theta, \rho \sin \theta, z) \rho \, d\rho\, d\theta\, dz.

Véase también[editar]

Referencias[editar]

  • Roland E. Larson, Robert P. Hosteler, Bruce H. Edwards (1999). «Integración Múltiple». Cálculo Volumen 2. México D.F.: McGrawHill. ISBN 970-10-2756-6.