UTEC · FICA

Matemática IV

Características de una Integral Doble

Índice

Teoría matemática

Definición sobre un rectángulo

Sea f:R⊂R2→Rf:R\subset\mathbb{R}^2\to\mathbb{R} continua en R=[a,b]×[c,d]R=[a,b]\times[c,d]. La integral doble es el límite de sumas de Riemann

∬Rf(x,y) dA=lim⁡m,n→∞∑i=1m∑j=1nf(xi∗,yj∗) ΔA,ΔA=Δx Δy.\iint_R f(x,y)\,dA =\lim_{m,n\to\infty}\sum_{i=1}^{m}\sum_{j=1}^{n}f(x_i^*,y_j^*)\,\Delta A, \qquad \Delta A=\Delta x\,\Delta y.

Si f≥0f\ge 0, la integral representa el volumen bajo z=f(x,y)z=f(x,y) y sobre RR. En aplicaciones, también puede acumular densidad de carga, potencia, masa o energía.

Teorema de Fubini

Si ff es continua sobre el rectángulo RR, la integral doble puede calcularse mediante cualquiera de los dos órdenes:

∬Rf(x,y) dA=∫ab ⁣(∫cdf(x,y) dy)dx=∫cd ⁣(∫abf(x,y) dx)dy.\iint_R f(x,y)\,dA =\int_a^b\!\left(\int_c^d f(x,y)\,dy\right)dx =\int_c^d\!\left(\int_a^b f(x,y)\,dx\right)dy.

La integral interior acumula sobre una franja; la exterior suma las contribuciones de todas las franjas.

Regiones tipo I y tipo II

  • Tipo I. Región descrita por a≤x≤ba\le x\le b y g1(x)≤y≤g2(x)g_1(x)\le y\le g_2(x). Integral iterada ∫ab ⁣∫g1(x)g2(x)f(x,y) dy dx\displaystyle\int_a^b\!\int_{g_1(x)}^{g_2(x)}f(x,y)\,dy\,dx.
  • Tipo II. Región descrita por c≤y≤dc\le y\le d y h1(y)≤x≤h2(y)h_1(y)\le x\le h_2(y). Integral iterada ∫cd ⁣∫h1(y)h2(y)f(x,y) dx dy\displaystyle\int_c^d\!\int_{h_1(y)}^{h_2(y)}f(x,y)\,dx\,dy.

En una región no rectangular, invertir el orden exige redescribir geométricamente los límites. No basta con intercambiar dxdx y dydy.

Propiedades básicas

  • Linealidad. ∬R(αf+βg) dA=α∬Rf dA+β∬Rg dA\iint_R(\alpha f+\beta g)\,dA=\alpha\iint_Rf\,dA+\beta\iint_Rg\,dA.
  • Aditividad. Si R=R1∪R2R=R_1\cup R_2 sin solape interior, ∬Rf=∬R1f+∬R2f\iint_Rf=\iint_{R_1}f+\iint_{R_2}f.
  • Monotonía. Si f≤gf\le g en RR, entonces ∬Rf dA≤∬Rg dA\iint_Rf\,dA\le\iint_Rg\,dA.
  • Constante. ∬Rc dA=c area⁡(R)\iint_Rc\,dA=c\,\operatorname{area}(R).

Idea-fuerza

Antes de integrar, dibuja la región y decide cómo la recorre la integral interior. En rectángulos, Fubini permite elegir el orden más cómodo; en regiones generales, los límites codifican la geometría.

Visualización

  • Figura 1 — Integral iterada doble.

    Sumas de Riemann bajo una superficie sobre un rectángulo

    Figura 1. La integral doble se interpreta como el límite de sumas de Riemann con partición cada vez más fina de la región rectangular.

  • Figura 2 — Orden de integración equivalente.

    Dos órdenes de integración sobre el mismo rectángulo

    Figura 2. Sobre un rectángulo, intercambiar dx dydx\,dy por dy dxdy\,dx recorre la misma región con el mismo valor final cuando se cumplen las hipótesis de Fubini.

  • Figura 3 — Región rectangular y sólido asociado.

    Región rectangular <span class="katex"><span class="katex-mathml"><math xmlns="http://www.w3.org/1998/Math/MathML"><semantics><mrow><mi>R</mi></mrow><annotation encoding="application/x-tex">R</annotation></semantics></math></span><span class="katex-html" aria-hidden="true"><span class="katex-base"><span class="katex-strut" style="height:0.6833em;"></span><span class="mord mathnormal" style="margin-right:0.0077em;">R</span></span></span></span> en el plano <span class="katex"><span class="katex-mathml"><math xmlns="http://www.w3.org/1998/Math/MathML"><semantics><mrow><mi>x</mi><mi>y</mi></mrow><annotation encoding="application/x-tex">xy</annotation></semantics></math></span><span class="katex-html" aria-hidden="true"><span class="katex-base"><span class="katex-strut" style="height:0.625em;vertical-align:-0.1944em;"></span><span class="mord mathnormal">x</span><span class="mord mathnormal" style="margin-right:0.0359em;">y</span></span></span></span> y el sólido bajo la superficie

    Figura 3. La región base RR del plano xyxy y el sólido cuya altura local es f(x,y)f(x,y) sobre cada punto.

Ejemplo Forma 1: Integral sobre un rectángulo

Problema

Una placa cuadrada de capacitor tiene lado L=0,04 mL=0{,}04\ \text{m} y densidad superficial

σ(x,y)=10−6(1+0,5x) C/m2,0≤x≤L,0≤y≤L.\sigma(x,y)=10^{-6}(1+0{,}5x)\ \text{C/m}^2, \qquad 0\le x\le L,\quad 0\le y\le L.

Calcular la carga total QQ.

Resolución

Analítica

La región es rectangular y la densidad no depende de yy:

Q=∫0L ⁣∫0L10−6(1+0,5x) dy dx=10−6L∫0L(1+0,5x) dx=10−6L[x+x24]0L=10−6L(L+L24).\begin{aligned} Q &=\int_0^L\!\int_0^L10^{-6}(1+0{,}5x)\,dy\,dx \\[0.45em] &=10^{-6}L\int_0^L(1+0{,}5x)\,dx \\[0.45em] &=10^{-6}L\left[x+\frac{x^2}{4}\right]_0^L \\[0.45em] &=10^{-6}L\left(L+\frac{L^2}{4}\right). \end{aligned}

Para L=0,04L=0{,}04:

Q=10−6(0,04)(0,0404)=1,616×10−9 C.Q=10^{-6}(0{,}04)(0{,}0404)=1{,}616\times10^{-9}\ \text{C}.

Resultado. La placa almacena Q=1,616 nCQ=1{,}616\ \text{nC} según el modelo.

Diferenciador (Ejemplo Forma 1 · Analítica)
  • Dato. La región tiene límites constantes y forma rectangular.
  • Herramienta. El teorema de Fubini autoriza cualquiera de los dos órdenes.
  • Bifurcación. En el paso de planteamiento se elige integrar primero en yy porque la densidad no depende de esa variable.
Verificación (Ejemplo Forma 1 · Analítica)

Al invertir el orden, la integral en xx vale 10−6(L+L2/4)10^{-6}(L+L^2/4) y la integral posterior en yy multiplica por LL. Se obtiene exactamente el mismo resultado. Además,

10−6L2<Q<1,02×10−6L2,10^{-6}L^2<Q<1{,}02\times10^{-6}L^2,

coherente con una densidad que crece ligeramente desde x=0x=0 hasta x=Lx=L.

Errores comunes (Ejemplo Forma 1 · Analítica)
  • Omitir el factor 10−610^{-6} al sustituir la densidad.
  • Integrar solo una vez y olvidar que la placa tiene dos dimensiones.
  • Mezclar centímetros y metros dentro del mismo cálculo.

Ejemplo Forma 2: Integral sobre una región general

Problema

La zona activa de un sensor es el triángulo con vértices (0,0)(0,0), (0,03,0)(0{,}03,0) y (0,03,0,02)(0{,}03,0{,}02), medidos en metros. La irradiancia es

E(x,y)=500(1+0,3y) W/m2.E(x,y)=500(1+0{,}3y)\ \text{W/m}^2.

Calcular la potencia óptica total recibida.

Resolución

Analítica

La recta inclinada es y=2x/3y=2x/3. Como región tipo I,

R={(x,y):0≤x≤0,03, 0≤y≤2x3}.R=\left\{(x,y):0\le x\le0{,}03,\ 0\le y\le\frac{2x}{3}\right\}.

Por tanto,

P=∫00,03 ⁣∫02x/3500(1+0,3y) dy dx=500∫00,03(2x3+x215)dx=500[x23+x345]00,03=0,1503 W.\begin{aligned} P &=\int_0^{0{,}03}\!\int_0^{2x/3}500(1+0{,}3y)\,dy\,dx \\[0.45em] &=500\int_0^{0{,}03}\left(\frac{2x}{3}+\frac{x^2}{15}\right)dx \\[0.45em] &=500\left[\frac{x^2}{3}+\frac{x^3}{45}\right]_0^{0{,}03} \\[0.45em] &=0{,}1503\ \text{W}. \end{aligned}

Resultado. El sensor recibe aproximadamente P=150,3 mWP=150{,}3\ \text{mW}.

Diferenciador (Ejemplo Forma 2 · Analítica)
  • Dato. La frontera superior depende de xx y la región es triangular.
  • Herramienta. Se describe la región como tipo I antes de plantear la integral.
  • Bifurcación. En el primer paso se obtiene y=2x/3y=2x/3; los límites interiores dejan de ser constantes.
Verificación (Ejemplo Forma 2 · Analítica)

La región tipo II equivalente es

R={(x,y):0≤y≤0,02, 3y2≤x≤0,03}.R=\left\{(x,y):0\le y\le0{,}02,\ \frac{3y}{2}\le x\le0{,}03\right\}.

La integral con orden inverso produce el mismo valor. Como el área triangular es 3×10−4 m23\times10^{-4}\ \text{m}^2, una irradiancia constante de 500 W/m2500\ \text{W/m}^2 daría 0,15 W0{,}15\ \text{W}; el resultado debe ser ligeramente mayor porque EE crece con yy.

Errores comunes (Ejemplo Forma 2 · Analítica)
  • Integrar sobre el rectángulo completo en lugar del triángulo.
  • Intercambiar dxdx y dydy sin redescribir los límites.
  • Usar y=3x/2y=3x/2 en vez de la pendiente correcta 2/32/3.

Errores comunes transversales

  1. Confundir ∬Rf dA\iint_R f\,dA con el área de RR; solo coinciden cuando f=1f=1.
  2. Ignorar las unidades del elemento de área: al pasar de cm a m, el factor es cuadrático.
  3. Suponer que invertir el orden conserva literalmente los límites en una región curva o triangular.
  4. Plantear la integral antes de dibujar la región y comprobar sus intersecciones.