OpenBetaMatemática

← Temas

La integral de Riemann-Stieltjes

La derivada mide cómo cambia algo; la integral mide cómo se acumula. El clímax del capítulo es el Teorema Fundamental, que revela que son operaciones inversas —asombroso, porque sus definiciones (cocientes vs. sumas de áreas) no se parecen en nada. Y la versión de Stieltjes conecta directo con la esperanza: E[f(X)]=fdFE[f(X)] = \int f\, dF.

La integral como límite de sumas

Para medir el área bajo una curva irregular la cortamos en franjas y aproximamos cada una por un rectángulo. La pregunta fina: ¿qué altura le damos si ff varía en la franja? Dos respuestas. Con una partición PP de [a,b][a,b], sean Mi=supfM_i = \sup f y mi=inffm_i = \inf f en el trozo ii:

L(P,f)=imiΔxi    aˊrea    iMiΔxi=U(P,f).L(P, f) = \sum_i m_i\, \Delta x_i \;\le\; \text{área} \;\le\; \sum_i M_i\, \Delta x_i = U(P, f).

La suma inferior subestima, la superior sobreestima, y el área real queda atrapada entre ambas. Al refinar la partición, las inferiores suben y las superiores bajan:

L(P1)L(P2)aˊreaU(P2)U(P1).L(P_1) \le L(P_2) \le \cdots \le \text{área} \le \cdots \le U(P_2) \le U(P_1).

ff es integrable si supPL(P,f)=infPU(P,f)\displaystyle\sup_P L(P,f) = \inf_P U(P,f), y ese valor común es abfdx\int_a^b f\,dx. Equivalentemente: para todo ε>0\varepsilon > 0 existe una partición con U(P,f)L(P,f)<εU(P,f) - L(P,f) < \varepsilon.

Es el espíritu ε\varepsilon del Cap. 4 aplicado al área: si la brecha ULU - L se hace tan chica como quieras, la integral existe. Mirala cerrarse (la banda roja es ULU-L):

Sumas superior e inferior — interactivof(x) = 1 + sin x
n = 8 franjasla banda roja es U − L
4.7231
inferior L
5.5085
superior U
0.7854
brecha U−L
5.1416
∫ (= 2+π)

La suma inferior usa la altura mínima de cada franja (subestima); la superior, la máxima (sobreestima). El área real queda atrapada entre las dos. f es integrable porque al afinar las franjas la brecha U − L tiende a 0: ambas se encuentran en un único valor, la integral.

Continuidad ⇒ integrabilidad. Si ff es continua en [a,b][a,b], es integrable. La razón encadena todo lo anterior: sobre un compacto una función continua es uniformemente continua, así que no varía demasiado en cada franja chica; la brecha MimiM_i - m_i se controla uniformemente y las sumas se aprietan. Funciones con saltos finitos (o discontinuas en un conjunto de medida cero) también pueden ser integrables — el puente exacto a Lebesgue.

Monótona ⇒ integrable, aun con saltos. En una partición uniforme las oscilaciones Mimi=f(xi)f(xi1)M_i - m_i = f(x_i) - f(x_{i-1}) telescopean a f(b)f(a)f(b) - f(a); con la malla chica, esa brecha total multiplicada por el ancho se va a cero. La monotonía limita la oscilación acumulada aunque haya (a lo sumo contables) saltos.

Pero no toda función acotada es integrable. La función de Dirichletf=1f = 1 en los racionales de [0,1][0,1], 00 en los irracionales— tiene en cada subintervalo Mi=1M_i = 1 (hay racionales) y mi=0m_i = 0 (hay irracionales), así que U=1U = 1 y L=0L = 0 para toda partición: las tenazas nunca se cierran. No es Riemann-integrable —y es, exactamente, la función que la integral de Lebesgue sí sabrá integrar. El fracaso de Riemann que motiva toda la teoría de la medida.

La generalización de Stieltjes

En Riemann cada franja pesa por su ancho Δxi\Delta x_i. En Stieltjes, una segunda función α\alpha actúa como peso, y cada franja pesa por Δαi=α(xi)α(xi1)\Delta\alpha_i = \alpha(x_i) - \alpha(x_{i-1}):

abfdα=limif(ti)Δαi.\int_a^b f\, d\alpha = \lim \sum_i f(t_i)\, \Delta\alpha_i.

Por qué importa: cuando α\alpha es una CDF, fdα\int f\,d\alpha es exactamente la esperanza E[f(X)]E[f(X)]. Y unifica lo continuo con lo discreto en una sola notación — subí el salto y mirá aparecer la masa puntual:

Integral de Stieltjes ∫f dα — interactivosalto J en x₀ = 2
f(x) — el integrandof(x₀)α(x) — el peso acumuladosalto J
J = 1.00subí el salto: la masa puntual crece
1.335
parte continua c·∫f dx
1.578
parte salto J·f(x₀)
2.913
∫ f dα (total)

En Stieltjes cada franja pesa por Δα, no por Δx. Donde α sube con pendiente constante, recuperás la integral de siempre (parte continua); donde α salta, ese punto aporta f(x₀)·(tamaño del salto), una masa puntual. Si α es una CDF, esto es exactamente E[f(X)] = ∫ f dF: densidad continua y átomos discretos en una sola fórmula.

  • α\alpha suave (variable continua): dα=α(x)dxd\alpha = \alpha'(x)\,dx recupera f(x)p(x)dx\int f(x)\,p(x)\,dx.
  • α\alpha con saltos (variable discreta): cada salto aporta f(xk)(taman˜o del salto)=f(xk)P(X=xk)f(x_k)\cdot(\text{tamaño del salto}) = f(x_k)\,P(X = x_k).

Una sola fórmula, E[f(X)]=fdFE[f(X)] = \int f\,dF, cubre densidades y masas puntuales. Por eso en estadística teórica la esperanza se escribe como integral de Stieltjes contra la CDF.

El Teorema Fundamental del Cálculo

Derivar e integrar se definieron por separado, con maquinarias distintas. El TFC dice que son inversas, y eso da un atajo gigante: calcular una integral desde la definición es horrible; basta hallar una antiderivada y evaluarla.

Parte 1 (derivar deshace integrar). Si F(x)=axf(t)dtF(x) = \int_a^x f(t)\,dt, entonces F(x)=f(x)F'(x) = f(x).

Parte 2 (la versión práctica). Si F=fF' = f, entonces abf(x)dx=F(b)F(a)\displaystyle\int_a^b f(x)\,dx = F(b) - F(a).

¿Por qué la Parte 1? F(x)F(x) es el área acumulada hasta xx. Si avanzo un pasito hh, gano una franjita de ancho hh y altura f(x)\approx f(x):

F(x+h)F(x)f(x)h    F(x+h)F(x)hf(x),F(x+h) - F(x) \approx f(x)\cdot h \;\Longrightarrow\; \frac{F(x+h) - F(x)}{h} \approx f(x),

y con h0h \to 0 la izquierda es F(x)F'(x). La tasa a la que se acumula el área es la altura de la curva: donde ff es alta el área crece rápido, donde es baja crece lento, donde f=0f=0 el área tiene un extremo. Movés xx y lo ves:

El Teorema Fundamental — interactivof = sin · F = 1 − cos
f(t) — la curvaF(x) = ∫₀ˣ f — el área acumulada
x = 2.20moví x: el área de arriba ES la altura de abajo
El área acumulada hasta x vale F(x) = 1.589. La velocidad a la que crece esa área es la altura de la curva ahí: F '(x) = f(x) = 0.808. Donde f > 0 el área sube; donde f < 0, baja; donde f = 0, F tiene un máximo o mínimo. Derivar el área devuelve la curva: eso es el Teorema Fundamental.

Modelo mental: si ff es el caudal (litros/seg), FF es el agua acumulada; el balde se llena a la velocidad a la que entra el agua, F=fF' = f. Un caso central: la CDF F(x)=xp(t)dtF(x) = \int_{-\infty}^x p(t)\,dt es el área acumulada de la densidad, y F(x)=p(x)F'(x) = p(x): densidad y CDF son antiderivada/derivada una de otra. Cuando no hay antiderivada cerrada, volvés a las sumas de rectángulos refinadas —la cuadratura numérica (trapecio, Simpson, Gauss). El TFC para lo analítico; las sumas para lo numérico.

La Parte 2, con rigor. Tomá una partición; el valor medio da en cada trozo un tit_i con F(xi)F(xi1)=f(ti)ΔxiF(x_i) - F(x_{i-1}) = f(t_i)\,\Delta x_i. Sumando, la izquierda telescopea a F(b)F(a)F(b) - F(a), y la derecha es una suma de Riemann atrapada entre L(P)L(P) y U(P)U(P); como ff es integrable, ambas convergen a abf\int_a^b f, forzando F(b)F(a)=abfF(b) - F(a) = \int_a^b f.   \;\blacksquare El TFC es, además, el caso 1-D del teorema de Stokes (Cap. 10): lo que se acumula dentro de [a,b][a,b] se lee desde su frontera {a,b}\{a, b\}.

Propiedades

ab(cf+g)dx=cabfdx+abgdx(linealidad)\int_a^b (cf + g)\,dx = c\int_a^b f\,dx + \int_a^b g\,dx \qquad \text{(linealidad)} acfdx=abfdx+bcfdx(aditividad, a<b<c)\int_a^c f\,dx = \int_a^b f\,dx + \int_b^c f\,dx \qquad \text{(aditividad, } a<b<c\text{)} fg    fg(monotonıˊa),abfdxabfdx(acotacioˊn).f \le g \;\Rightarrow\; \int f \le \int g \qquad \text{(monotonía)}, \qquad \left| \int_a^b f\,dx \right| \le \int_a^b |f|\,dx \qquad \text{(acotación)}.

No son tecnicismos: la linealidad es E[aX+Y]=aE[X]+E[Y]E[aX + Y] = aE[X] + E[Y]; la monotonía es XYE[X]E[Y]X \le Y \Rightarrow E[X] \le E[Y]; la acotación es la desigualdad triangular integral, base de cotas como Jensen.

Ejercicios resueltos

La receta del TFC es siempre la misma: (1) hallá una antiderivada FF (subiendo el exponente con xndx=xn+1n+1\int x^n\,dx = \tfrac{x^{n+1}}{n+1}); (2) evaluá F(b)F(a)F(b)-F(a).

1. 03x2dx=[x33]03=2730=9.\displaystyle \int_0^3 x^2\,dx = \left[\frac{x^3}{3}\right]_0^3 = \frac{27}{3} - 0 = 9.

2. Con linealidad, 12(x3+2x)dx=[x44+x2]12=(4+4)(14+1)=274=6.75.\displaystyle \int_1^2 (x^3 + 2x)\,dx = \left[\frac{x^4}{4} + x^2\right]_1^2 = (4 + 4) - \left(\tfrac14 + 1\right) = \frac{27}{4} = 6.75.

3. Polinomio completo, 02(3x24x+1)dx=[x32x2+x]02=(88+2)0=2.\displaystyle \int_0^2 (3x^2 - 4x + 1)\,dx = \left[x^3 - 2x^2 + x\right]_0^2 = (8 - 8 + 2) - 0 = 2.

4. Aditividad sin calcular nada: si 01f=3\int_0^1 f = 3 y 12f=5\int_1^2 f = 5, entonces 02f=01f+12f=8.\displaystyle \int_0^2 f = \int_0^1 f + \int_1^2 f = 8.

Las reglas de derivación del Cap. 5, leídas al revés, se vuelven la maquinaria para integrar. No hacen falta dos juegos de herramientas: una operación deshace a la otra.