13  Integrales de funciones de varias variables

En este capítulo extendemos el concepto de integral de Riemann visto en el capítulo de integrales a funciones de varias variables, en particular en R2 y R3 donde las integrales múltiples tienen multitud de aplicaciones en Ciencias e Ingenierías.

También veremos que el Teorema Fundamental del Cálculo sigue siendo válido y que el cálculo de integrales con respecto a una variable es el proceso inverso del cálculo de la derivada parcial con respecto a esa misma variable.

13.1 Sumas de Riemann

Las sumas de Riemann que se introdujeron en la pueden generalizarse fácilmente a funciones de varias variables, pero primero debemos generalizar el concepto de intervalo cerrado en Rn.

Definición 13.1 (Intervalo n dimensional) Dados n intervalos cerrados de R Ii=[ai,bi] i=1n, su producto cartesiano I=I1××In se conoce como intervalo n dimensional cerrado o hiperrectángulo.

Nota

En R2, el intervalo [ax,bx]×[ay,by] define un rectángulo, mientras que en R3 el intervalo [ax,bx]×[ay,by]×[az,bz] define una caja.

Definición 13.2 (Partición de un intervalo n dimensional) Dado un intervalo n dimensional cerrado I=I1××In, con Ii=[ai,bi] i=1,,n, y una partición Pi={xi,0=ai,xi,1,xi,mi} de cada intervalo Ii, el producto cartesiano P=P1××Pn se conoce como partición n dimensional de I.

El conjunto de todas las particiones n dimensionales de I se denota P(I).

Definición 13.3 (Sumas de Riemann) Dada una función f:IR acotada en el intervalo n dimensional IRn y una partición n dimensional P de I, se define la suma inferior de Riemann de f respecto de P, y se denota s(f,P), como

s(f,P)=i1=1m1in=1mnmi1,,in(x1,i1x1,i11)(xn,inxn,in1),

donde mi1,,in=inf{f(x1,,xn):xj[xj,ij1,xj,ij] j=1,,n}.

Del mismo modo, se define la suma superior de Riemann de f respecto de P, y se denota S(f,P), como

S(f,P)=i1=1m1in=1mnMi1,,in(x1,i1x1,i11)(xn,inxn,in1),

donde Mi1,,in=sup{f(x1,,xn):xj[xj,ij1,xj,ij] j=1,,n}.

Aunque la definición es un poco enrevesada para Rn, si f(x,y) es una función de dos variables acotada en el intervalo I=[a,b]×[c,d]R2, y P=Px×Py es una partición de I con Px={x0=a,x1,,xn=b} y Px={y0=c,y1,,ym=d}, las sumas de Riemann inferior y superior de f respecto de P son

s(f,P)=i=1nj=1mmij(xixi1)(yjyj1),S(f,P)=i=1nj=1mMij(xixi1)(yjyj1),

donde mij y Mij son el ínfimo y el supremo de {f(x,y):(x,y)[xi1,xi]×[yj1,yj]}.

Gráficamente, si f(x,y) es una función positiva en I, la suma inferior se puede interpretar como la suma de los volúmenes de los prismas de base rectangular con lados [xi1,xi] y [yj1,yj] y altura mij.

Suma inferior de Riemann en R2

Ejemplo 13.1 Veamos cómo calcular las sumas de Riemann de la función f(x,y)=x2y en el rectángulo [0,1]×[0,1] para la partición P={0,0.5,1}×{0,0.5,1}

Esta partición descompone el intervalo bidimensional [0,1]×[0,1]$ en 4 rectángulos (en realidad cuadrados) como los que se ven la siguiente figura.

Partición del intervalo [0,1]×[0,1] en 4 subintervalos.

Como f es creciente tanto en x como en y en el intervalo I, alcanzará el mínimo en el extremo inferior izquierdo de cada subintervalo y el máximo en el extremo superior derecho. Por tanto, la suma inferior de Riemann es

s(f,P)=f(0,0)(0.50)(0.50)+f(0,0.5)(0.50)(10.5)+f(0.5,0)(10.5)(0.50)+f(0.5,0.5)(10.5)(10.5)=00.25+00.25+00.25+0.1250.25=0.03125.

Y la suma superior de Riemann es

S(f,P)=f(0.5,0.5)(0.50)(0.50)+f(0.5,1)(0.50)(10.5)+f(1,0.5)(10.5)(0.50)+f(1,1)(10.5)(10.5)=0.1250.25+0.250.25+0.50.25+10.25=0.46875.

13.2 Integral de Riemann múltiple

Cuando la función f(x,y) es positiva en el intervalo I, la suma inferior de Riemann nos da una aproximación por defecto del volumen encerrado entre la superficie de la función y el plano xy en la región definida por I, mientras que la suma superior de Riemann nos da una aproximación por exceso. Al hacer refinamientos de la partición P, cada vez con mayor número de subintervalos, las sumas inferiores crecen y las superiores decrecen, dando aproximaciones cada vez mejores, de manera que en el límite, cuando el número de subintervalos tiende a , obtendremos el volumen real.

Definición 13.4 (Integral inferior y superior de Riemann) Dada una función f:IR acotada en el intervalo n dimensional IRn y una partición n dimensional P de I, se define la integral inferior de Riemann de f en I, como

If=sup{s(f,P):PP(I)}.

Y se define la integral superior de Riemann de f en I, como

If=inf{S(f,P):PP(I)}.

Definición 13.5 Dada una función f:IR acotada en el intervalo n dimensional IRn, se dice que f es integrable Riemann en I si

If=If,

y a este valor se le llama integral de Riemann o integral definida de f en I y se denota por

Ifo bienIf(x1,,xn)dA, donde dA=dx1dxn.

Proposición 13.1 Si f:IR es integrable Riemann en el intervalo n dimensional IRn, entonces

If=limmi1=1min=1mf(x1,i1,,xn,in)ΔA,

donde P es una partición de I en mn subintervalos de igual tamaño y ΔA=(x1,i1x1,i11)(xn,inxn,in1) es el hipervolumen de los subintervalos de la partición.

En el caso de una función de dos variables f(x,y) la proposición anterior establece que si f es integrable Riemann en un intervalo [a,b]×[c,d], entonces

If=limni=1nj=1nf(xi,yj)ΔA,

donde Pn={x0=a,x1,,xn=b}×{y0=c,y1,,yn=d} y ΔA=(xixi1)(yjyj1), es el área de los rectángulos definidos por cada subintervalo de la partición. Es decir, la integral de Riemann aparece al tomar particiones cada vez más refinadas con subintervalos de igual tamaño y hacer la suma de las áreas de los subintervalos multiplicadas por el valor de la función en el extremo superior derecho del subintervalo.

Nota

Cuando se quieren hacer explícitas las variables de las que depende la función f, la integral de Riemann se suele escribir

If(x,y)dxdy

para funciones de dos variables y

If(x,y,z)dxdydz

para funciones de tres variables.

13.3 Propiedades de las Integrales de Riemann múltiples

Las propiedades de la integral de Riemann para funciones de una variable también se pueden generalizar para integrales de Riemann múltiples.

Teorema 13.1 Si f,g:IR son dos funciones integrables Riemann en un intervalo n dimensional I, entonces

  1. f+g es integrable Riemann en I y I(f+g)=If+Ig.

  2. Para cualquier cR, cf es integrable Riemann en I y Icf=cIf.

  3. Si f(x1,,xn)g(x1,,xn) (x1,,xn)I, entonces IfIg.

Prueba. La demostración es similar a la vista para funciones de una variable y se deja como ejercicio.

13.4 Integrales múltiples iteradas

El cálculo de integrales de Riemann tomando límites de sumas de Riemann resulta poco operativo en la práctica, pero afortunadamente, podemos utilizar el teorema fundamental del cálculo para calcular integrales de manera mucho más rápida y sencilla. En esta sección veremos como calcular integrales dobles y triples como integrales iteradas, que se pueden calcular realizando varias integrales simples.

El procedimiento para reducir el cálculo de una integral múltiple al cálculo de varias integrales simples es, en el fondo, el mismo que se utilizó para calcular derivadas parciales, es decir, realizar integrales simples de las funciones parciales que se obtienen al fijar como constantes todas las variables excepto aquella con respecto a la que se integra.

Veamos el procedimiento para una función de dos variables. Supongamos que f(x,y) es integrable Riemann en el intervalo I=[a,b]×[c,d]. Si mantenemos x constante, podemos calcular la integral definida simple con respecto a y en el intervalo [c,d], es decir, cdf(x,y)dy. Para cada valor de x, esta integral tomará un valor distinto, por lo que podemos definir la función

A(x)=cdf(x,y)dy,

que solo depende de x. Si ahora se integra esta función con respecto a x en el intervalo [a,b], se tiene

abA(x)dx=abcdf(x,y)dydx.

Esta integral se conoce como integral doble iterada, primero con respecto a y y después con respecto a x.

De manera similar, podríamos haber integrado primero con respecto a x y luego con respecto a y, y tendríamos esta otra integral doble iterada

cdabf(x,y)dxdy.

El siguiente teorema establece que una integral múltiple se puede calcular mediante integrales iteradas.

Teorema 13.2 (Fubini) Si f(x,y) es una función continua en un intervalo I=[a,b]×[c,d], entonces

If(x,y)dA=abcdf(x,y)dydx=cdabf(x,y)dxdy.

Prueba. Daremos una demostración intuitiva. Si f es positiva, la integral doble es el volumen del sólido encerrado entre la superficie de f y en plano xy en el intervalo I=[a,b]×[c,c]. Si fijamos el valor de x a un valor constante x0[a,b] e intersecamos este sólido con el plano x=x0, se obtiene una región plana cuya área se puede calcular mediante la integral

cdf(x0,y)dy

De esta manera, para cada x[a,b] la función

A(x)=cdf(x,y)dy

nos dará el área de la región plana que se obtiene de la intersección del sólido con el plano que pasa por x y es perpendicular al eje x. Por tanto, si integramos estas áreas con respecto a y en el intervalo [c,b],

abA(x)dx=abcdf(x,y)dydx

obtendremos el volumen del sólido, tal y como se vió en la .

El mismo razonamiento se puede usar cambiando el orden de integración.

Ejemplo 13.2 La integral doble de la función f(x,y)=122xy2 en el intervalo I=[0,1]×[0,3] es

If(x,y)dA=0103122xy2dydx=01[12y2xyy33]03dx=01366x333dx=01276xdx=[27x3x2]01=273=24.

Podríamos haber llegado a este mismo resultado haciendo la otra integral iterada

If(x,y)dA=0301122xy2dxdy=03[12xx2xy2]01dy=03121y2dy=0311y2dy=[11yy33]03=33333=24.

Importante

Aunque cuando la función es integrable Riemann el orden de integración de las integrales iteradas no importa, a menudo, una integral iterada suele ser más fácil de calcular que la otra.

Ejemplo 13.3 Para calcular la integral de la función f(x,y)=xexy en el intervalo I=[0,1]×[0,1], si primero integramos con respecto a y y después con respecto a x se tiene

If(x,y)dA=0101xexydydx=01[exy]01=01ex1dx=[exx]01=e2.

Sin embargo, si primero integramos con respecto a x y luego con respecto a y se tiene una integral mucho más difícil de calcular.

Calculadora de integrales dobles iteradas. Calculadora de sumas de Riemann

En el caso de que f(x,y) pueda factorizarse como producto de funciones que sólo dependen de x o de y, entonces el cálculo de la integral doble se simplifica.

Proposición 13.2 Si f(x,y)=g(x)h(y) es integrable Riemann en el intervalo I=[a,b]×[c,d], entonces

If(x,y)dA=abg(x)dxcdh(y)dy.

Prueba. If(x,y)dA=abcdf(x,y)dydx=abcdg(x)h(y)dydx(1)=abg(x)cdh(y)dydx(2)=cdh(y)dyabg(x)dx. (1) g(x) no depende de y.
(2) cdh(y)dy no depende de x.

Ejemplo 13.4 La integral doble de la función f(x,y)=e2x3y en el intervalo I=[0,1]×[0,2] es

Ie2x3ydxdy=Ie2xe3ydxdy=01e2xdx02e3ydy=[e2x2]01[e3y3]02=(e2212)(e63+13)

13.5 Integrales sobre regiones no regulares

En muchos problemas reales, como por ejemplo en el cálculo de volúmenes, la región de integración no es un integración no es un intervalo n dimensional o hiperrectángulo, sino una región no regular, es decir, con fronteras curvas. En esta sección veremos cómo calcular integrales dobles de funciones de dos variables para algunos tipos de regiones irregulares acotadas del plano xy.

En general, la integral de Riemann de una función f sobre un intervalo n dimensional puede extenderse a una región irregular acotada R simplemente tomando un intervalo n dimensional I que incluya la región, es decir RI y redefiniendo f como

F(x1,,xn)={f(x1,,xn)si (x,1,,xn)R0si (x,1,,xn)IR.

y calculando la integral de Riemann de F sobre I como se ha visto en la sección anterior, es decir,

Rf(x1,,xn)dA=IF(x1,,xn)dA.

Este procedimiento tiene sentido siempre que f sea continua sobre la frontera de la región R, ya que el valor de F en los puntos que están fuera de la región R es 0, y por tanto, no contribuyen la suma de la integral.

En el caso de funciones de dos variables, existen varios tipos de regiones en el plano xy donde la integral de Riemann se puede calcular mediante integrales iteradas donde los límites de integración dependen de funciones que delimitan la región de integración.

Si R es una región del plano xy delimitada por las gráficas de dos funciones g1(x) y g2(x), es decir,

R={(x,y):axb,g1(x)yg2(x)},

la integral de Riemann f sobre esta región se puede calcular mediante la integral iterada

Rf(x,y)dA=abg1(x)g2(x)f(x,y)dydx.

Del mismo modo, si R es una región del plano xy delimitada por las gráficas de dos funciones g1(y) y g2(y), es decir,

R={(x,y):g1(y)xg2(y),ayb},

la integral de Riemann f sobre esta región se puede calcular mediante la integral iterada

Rf(x,y)dA=abg1(y)g2(y)f(x,y)dxdy.

Ejemplo 13.5 Veamos cómo calcular la integral doble de la función f(x,y)=xy sobre la región encerrada entre las curvas y=x e y=x2. La región encerrada entre las curvas se muestra en el siguiente gráfico.

Región comprendida entre las curvas y=x e y=x2.

Esta región puede expresarse como R={(x,y)R2:0x1,x2yx}, por lo que la integral doble de f sobre esta región puede calcularse mediante la integral iterada

01x2xxydydx=01[xy22]x2xdx=01xx2x42dx=01x3x52dx=[x48x612]01=18112=124.

En este caso, la región también puede expresarse como R={(x,y)R2:yxy,0y1}, por lo que la integral doble de f sobre esta región también puede calcularse mediante la integral iterada

01yyxydxdy=01[x22y]yydy=01yy22ydy=01y2y32dx=[y36y48]01=1618=124.

En este tipo de integrales iteradas, hay que integrar siempre primero con respecto a la variable que está acotada por las funciones que definen la región y después con respecto a la variable que está acotada por valores fijos. Cuando interese realizar la integral iterada en el otro orden, porque resulte más sencilla, previamente hay que cambiar la expresión de la región de integración mediante las funciones inversas, siempre y cuando sea posible.

Ejemplo 13.6 Para calcular la integral de la función f(x,y)=ey2 en la región R={(x,y)R2:0x1,0y2x}, como y está acotada por una función de x, habría que hacer la integral iterada

0102xey2dydx.

Sin embargo, la función ey2 no tiene primitiva elemental, por lo que necesariamente hay que realizar la integral iterada en el otro orden, pero para ello previamente hay que expresar la región como R={(x,y)R2:0xy2,0y2}, con lo que la integral doble sobre esta región es

020y/2ey2dxdy=02ey2[x]0y/2dy=02ey2y2dy=[ey24]02=e4+14.

13.6 Integrales dobles en coordenadas polares

Al que igual que se vio en la , para integrar determinadas funciones es preferible hacer un cambio a coordenadas polares. Recordemos que para pasar de coordenadas cartesianas (x,y) a coordenadas polares (r,θ) se utilizan las fórmulas

r=x2+y2θ=arcsin(yx2+y2)=arccos(xx2+y2),

y para pasar de coordenadas polares a cartesianas se utiliza la fórmula

x=rcos(θ)y=rsen(θ).

En coordenadas polares un intervalo bidimensional I=[a,b]×[α,β] define un sector de disco como el que se muestra a continuación.

Intervalo bidimensional en coordenadas polares.

Si la función f(x,y) es continua sobre un intervalo polar I=[a,b]×[α,β] con βα2π, la integral doble de f sobre esta región se puede calcular mediante la integral iterada

If(x,y)dA=αβabf(rcos(θ),rsenθ)rdrdθ.

Ejemplo 13.7 La gráfica de la función f(x,y)=a2x2y2 es una semiesfera de radio uno centrada en el origen de coordenadas. Para calcular su volumen, la región de integración es el círculo de radio uno centrado en el origen de coordenadas, que puede expresarse como R={(x,y)R2:0x1,a2x2ya2x2}. En este caso la integral en coordenadas rectangulares es complicada y merece la pena hacer la integral en coordenadas polares.

If(x,y)dA=02π0af(rcos(θ),rsenθ)rdrdθ=02π0aa2(rcos(θ))2(rsen(θ))2rdrdθ=02π0aa2r2(cos(θ)2+sen(θ)2)rdrdθ=02π0aa2r2rdrdθ=02π[(a2r2)3/23]0adθ=02πa23dθ=[a2θ3]02π=2πa23.

Y por tanto, el volumen de la esfera de radio uno centrada en el origen es 43πa2, que coincide con la fórmula habitual.

Cuando la región de integración es irregular pero puede expresarse de la forma R={(r,θ)R+×R:αθβ,g1(θ)rg2(θ)}, la integral doble en coordenadas polares puede calcularse mediante la integral iterada

Rf(x.y)dA=αβg1(θ)g2(θ)f(rcos(θ),rsen(θ))rdrdθ.

Ejemplo 13.8 Veamos cómo calcular el volumen encerrado entre la gráfica integral del paraboloide f(x,y)=x2+y2 y el plano xy sobre la circunferencia de ecuación x2+(y21)=1. En coordenadas polares, la circunferencia que delimita la región de integración puede escribirse como R={(r,θ)R+×R:0θπ,0r2sen(θ)}. Por tanto, podemos calcular la integral de f sobre la región R mediante la integral iterada

Rf(x,y)dA=0π02sen(θ)f(rcos(θ),rsenθ)rdrdθ=0π02sen(θ)((rcos(θ))2+(rsenθ)2)rdrdθ=0π02sen(θ)r3drdθ=0π[r44]02sen(θ)dθ=0π4sen(θ)4dθ=0π(1cos(2θ))2dθ=0π12cos(2θ)+cos(2θ)2dθ=0π12cos(2θ)+1+cos(4θ)2dθ=[θsen(2θ)+θ2+sen(4θ)8]0π=32π.

13.7 Cálculo de volúmenes

Cuando f(x,y) es positiva en el intervalo de integración I=[a,b]×[c,b], la integral doble If(x,y)dxdy da el volumen del sólido encerrado entre la superficie de f y el plano xy en la región delimitada por I. Ahora bien, si f es negativa en I, entonces el valor de la integral doble es negativo, por lo que habrá que cambiarla de signo para poder interpretarla como un volumen.

La cosa se complica cuando f puede tomar valores positivos y negativos en el intervalo I. En tal caso, para obtener el volumen tendremos que calcular la integral del valor absoluto de la función, o bien descomponer el intervalo de integración I en subintervalos donde la función sea positiva o negativa, integrar la función en cada uno de estos subintervalos y luego sumar el valor absoluto de estas integrales.

Ejemplo 13.9 Veamos cómo calcular el volumen de la función f(x,y)=x2y en el intervalo I=[1,1]×[1,1]. En el intervalo de integración, la función toma valores negativos en el subintervalo I1=[1,1]×[1,0] y valores positivos en el intervalo [1,1]×[0,1], por lo que el volumen encerrado entre la gráfica de la función y el plano xy viene dado por la suma de las integrales

I1f(x,y)dA=1110x2ydyx=11[x2y22]10dx=11x22dx=[x36]11=16+16=13.

y

I1f(x,y)dA=1101x2ydyx=11[x2y22]01dx=11x22dx=[x36]11=16+16=13,

es decir,

Ix2ydA=I1x2ydA+I2x2ydA=13+13=23.

Tip

Obsérvese que también se puede usar la integral doble para calcular el área de una región plana integrando sobre esta región la función f(x,y)=1, ya que cuanto la altura 1, el volumen coincide con el área de la base.

Del mismo modo se puede usar una integral tripe para calcular el volumen de una región tridimensional.

13.8 Cálculo del area de una superficie

Otra aplicación menos intuitiva de la las integrales dobles es el cálculo del área de la superficie de una función sobre una región R. Si f es una función de dos variables f(x,y), con derivadas parciales continuas sobre un intervalo bidimensional I=[a,b]×[c,d], para aproximar el área de la superficie de f sobre I podemos tomar una partición Px={x0=a,x1,,xn=b} del intervalo [a,b] en n subintervalos de igual tamaño, y otra Py={y0=a,y1,,yn=b} del intervalo [c,d], y descomponer el intervalo I en subintervalos rectangulares Iij mediante la partición Pn=Px×Py. De esta manera, todos los subintervalos Iij tendrán el mismo área ΔA=ΔxΔy.

Si para cada subintervalo Iij, podemos calcular el plano tangente a la superficie de f en el punto (xi,yj,f(xi,yj)), el paralelogramo Tij que define el plano tangente sobre el intervalo Iij nos dará una aproximación del área de la superficie de f sobre el subintervalo Iij, y por tanto, la suma de las áreas de estos paralelogramos tangentes será una aproximación de la superficie de f sobre el intervalo I.

Si para cada paralelogramo Tij consideramos los vectores uij y vij con origen en (xi,yj) y con la misma dirección y longitud que los lados del paralelogramo, tal y como se aprecia en la gráfica de más abajo, se tiene que el área del paralelogramo es ΔTij=|uij×vij|, tal y como se vió en la .

Por otro lado, como uij está contenido en la intersección plano tangente a f en (xi,yj,f(xi,yj)), con el plano y=yj, se tiene que su pendiente es fx(xi,yj), por lo que sus coordenadas son (Δx,0,fx(xi,yj)Δx). Del mismo modo, como vij está contenido en la intersección plano tangente, con el plano x=xi, se tiene que su pendiente es fy(xi,yj), por lo que sus coordenadas son (0,Δy,fy(xi,yj)Δy).

Así pues, el producto vectorial de uij y vij es

uij×vij=|ijkΔx0fx(xi,yj)Δx0Δyfy(xi,yj)Δy|=fx(xi,yj)ΔxΔyify(xi,yj)ΔxΔyj+ΔxΔyk=(fx(xi,yj),fy(xi,yj),1)ΔxΔy,

y el área del paralelogramo Tij es

ΔTij=|uij×vij|=fx(xi,yj)2+fy(xi,yj)2+1ΔxΔy.

De este modo, el área de la superficie de f sobre el intervalo I se puede aproximar mediante la suma

i=1nj=1nΔTij=i=1nj=1nfx(xi,yj)2+fy(xi,yj)2+1ΔxΔy.

Tomando particiones cada vez más refinadas, en el límite cuando n se tiene que el área de la superficie de f sobre el intervalo I viene dado por la integral

AI(f)=Ifx(xi,yj)2+fy(xi,yj)2+1dA.

Ejemplo 13.10 La gráfica de la función f(x,y)=a2x2y2 es la semiesfera de radio a centrada en el origen de coordenadas. El área de la superficie de esta semiesfera puede calcularse mediante la integral

AR(f)=Rfx(xi,yj)2+fy(xi,yj)2+1dA

donde R={(x,y)R:axa,a2x2ya2x2}.

Como las derivadas parciales de f valen

fx(x,y)=xa2x2y2,fy(x,y)=ya2x2y2,

la integral anterior se puede calcular mediante la integral iterada

AR(f)=aaa2x2a2x2x2a2x2y2+y2a2x2y2+1dydx.

Esta integral, no obstante, es difícil de calcular, por lo que merece la pena hacer el cambio a coordenadas polares, de manera que la región de integración polar es R={(r,θ)R+×R:0θ2π,0ra} y la integral se convierte en

AR(f)=02π0ar2cos(θ)2a2r2+r2sen(θ)2a2r2+1rdrdθ=02π0aa2a2r2rdrdθ.

Esta última integral es impropia ya que a2a2r2 cuando ra, pero la integral es finita, ya que

limsa0sa2a2r2rdr=limsaa(a2s2+a)=a2,

y siguiendo con el cálculo de la integral anterior se tiene

AR(f)=02πa2dθ=2πa2,

por lo que el área de la esfera de radio a centrada en el origen de coordenadas es 4πa2, que coincide con la fórmula habitual.

13.9 Aplicaciones físicas

13.9.1 Centro de masas

Una aplicación física inmediata de las integrales dobles o triples es el cálculo de la masa de una sólido plano o del espacio con densidad variable en cada punto. Si δ(x,y) es la densidad en cada uno de los puntos de un sólido que ocupa una región R del plano xy, la masa total del sólido puede calcularse mediante la integral

m=Rδ(x,y)dA,

siempre y cuando la función δ sea integrable en la región R.

Ejemplo 13.11 Una placa metálica tiene la forma de región R={(x,y)R2:0x1,0y1x2} y la densidad en cada punto viene dada por la función δ(x,y)=xy en gr/cm2. Entonces su masa es

m=Rδ(x,y)dA=0101x2xydydx=01x[y22]01x2dx=01x(1x2)22dx=01x2x3+x52dx=[x242x48+x612]01=1428+112=112 gr.

También podemos utilizar las integrales dobles para calcular el momento de cada punto de un sólido R del plano xy con respecto a cada uno de los ejes. Tal y como se vio en la , el momento con respecto a cada eje en un punto (x,y) de un sólido plano R viene dado por el producto de la masa en el punto y la distancia del punto al eje, es decir, el momento con respecto al eje x es xm(xy) y el momento con respecto al eje y es ym(x,y), donde m(x,y) es la masa en el punto (x,y). La suma de los momentos en cada uno de los puntos del sólido región R nos dará el momento de todo el sólido. Así pues, si δ(x,y) es la densidad en cada uno de los puntos de un sólido plano R, los momentos de R con respecto a los ejes x e y pueden calcularse, respectivamente, mediante las integrales

Mx=Ryδ(x,y)dAMy=Rxδ(x,y)dA

A partir de los momentos resulta sencillo calcular el centro de masas o centroide (x¯,y¯) del sólido plano, dividiendo cada momento por la masa del sólido, es decir,

x¯=Mym=Rxδ(x,y)dARδ(x,y)dAy¯=Mxm=Ryδ(x,y)dARδ(x,y)dA

Nota

Un sólido plano se comporta como si toda su masa se concentrara en el centroide, de manera que se mantendría en equilibrio horizontal al apoyarlo sobre este punto.

Ejemplo 13.12 El momento con respecto al eje x de la placa metálica del ejemplo anterior es

Mx=Ryδ(x,y)dA=0101x2xy2dydx=01x[y33]01x2dx=01x(1x2)33dx=01x3x3+3x5x73dx=[x26x44+x66x824]01=1614+16124=124 gr cm.

Del mismo modo, el momento con respecto al eje y es

My=Rxδ(x,y)dA=0101x2x2ydydx=01x2[y22]01x2dx=01x2(1x2)22dx=01x22x4+x62dx=[x36x55+x714]01=1615+114=4105 gr cm.

Por lo tanto, el centroide de la placa metálica tendrá coordenadas

x¯=Mxm=Ryδ(x,y)dARδ(x,y)dA=1/241/12=0.5 cm,y¯=Mym=Rxδ(x,y)dARδ(x,y)dA=4/1051/12=16350.4571 cm.

13.10 Aplicaciones estadísticas

13.10.1 Valor medio de una función

En la vimos cómo utilizar la integral definida para calcular el valor medio de una función de una variable en un intervalo unidimensional. Para funciones de dos variables podemos usar la misma idea para calcular el valor medio de una función sobre una región R dividiendo la integral por el área de la región de integración, y para funciones de tres variables dividiendo por el volumen de la región de integración. En general, para calcular el valor medio de una función de varias variables en una región de Rn utilizaremos la integral de la siguiente definición.

Definición 13.6 (Valor medio de una función de varias variables) Dada una función f:RnR integrable sobre una región RRn, el valor medio de f sobre R es

f¯R=Rf(x1,,xn)dARdA

Nota

Obsérvese que la integral del denominador RdA, para una función de dos variables da el área de la región R, mientras que para una función de tres variables da el volumen de la región R.

Ejemplo 13.13 En el vimos que la integral doble de la función f(x,y)=xy sobre la región R={(x,y)R2:0x1,x2yx} era RxydA=124. Como el área de la región R es

RdA=01x2xdydx=01[y]x2xdx=01xx2dx=[x22x33]01=1213=16.

el valor medio de f en R es

f¯R=RxydARdA=1/241/6=14.