Calculus·Unidad 4

Calc-04 — Integration essentials

La integración importa sobre todo por probabilidad: la probabilidad de un intervalo es el área bajo la densidad (integral de la PDF), y la expectation (media) de una variable continua es una integral. También aparece en normalizar distribuciones (que integren a 1). No necesitas ser un virtuoso integrando; necesitas entender qué significa una integral.

Why this matters for ML

La integración importa sobre todo por probabilidad: la probabilidad de un intervalo es el área bajo la densidad (integral de la PDF), y la expectation (media) de una variable continua es una integral. También aparece en normalizar distribuciones (que integren a 1). No necesitas ser un virtuoso integrando; necesitas entender qué significa una integral.

Concepts covered

  • Antiderivada (integral indefinida)
  • Integral definida como área/acumulación
  • Fundamental Theorem of Calculus (conecta derivada e integral)
  • Reglas básicas y sustitución (idea)
  • Uso en probabilidad: densidades y expectation

Intuition first

🎬 3Blue1Brown Essence of Calculus cap. 8 ("Integration and the fundamental theorem") y cap. 9–10. Idea: integrar es sumar infinitas rebanaditas → área acumulada; es la operación inversa de derivar.

Theory & key results

Antiderivada: FF es antiderivada de ff si F=fF'=f. La integral indefinida f(x)dx=F(x)+C\int f(x),dx = F(x)+C.

Integral definida: abf(x)dx\int_a^b f(x),dx = área con signo bajo ff entre aa y bb. Se define como límite de sumas de Riemann (rebanadas f(x)Δxf(x)\Delta x).

Fundamental Theorem of Calculus (FTC): conecta ambos mundos: abf(x)dx=F(b)F(a),donde F=f.\int_a^b f(x),dx = F(b)-F(a),\quad\text{donde } F'=f. Derivar e integrar son operaciones inversas.

Reglas básicas:

  • xndx=xn+1n+1+C\int x^n dx = \frac{x^{n+1}}{n+1}+C (n1n\neq-1)
  • exdx=ex+C\int e^x dx = e^x + C, 1xdx=lnx+C\int \frac1x dx = \ln|x| + C
  • Linealidad: (af+bg)=af+bg\int (af+bg) = a\int f + b\int g.
  • Sustitución (u-substitution): el "chain rule al revés".

Conexión con probabilidad (lo que de verdad usarás):

  • Una densidad (PDF) p(x)p(x) cumple p(x)dx=1\int_{-\infty}^{\infty}p(x),dx = 1 (toda la probabilidad suma 1).
  • P(aXb)=abp(x)dxP(a\le X\le b) = \int_a^b p(x),dx (área bajo la curva).
  • Expectation: E[X]=xp(x)dx\mathbb{E}[X]=\int_{-\infty}^{\infty}x,p(x),dx; y E[g(X)]=g(x)p(x)dx\mathbb{E}[g(X)]=\int g(x)p(x)dx.

En la práctica de ML muchas de estas integrales o se conocen en forma cerrada (gaussiana) o se estiman por muestreo/Monte Carlo — pero la idea de "integral = valor esperado / área de probabilidad" es imprescindible.

Worked example

02(3x2)dx\int_0^2 (3x^2),dx. Antiderivada: x3x^3. FTC: 2303=82^3 - 0^3 = 8.

Verificación de densidad: la uniforme en [0,2][0,2] es p(x)=1/2p(x)=1/2. 0212dx=122=1\int_0^2 \frac12 dx = \frac12\cdot2 = 1. ✓ Su media: 02x12dx=12x2202=122=1\int_0^2 x\cdot\frac12 dx = \frac12\cdot\frac{x^2}{2}\big|_0^2 = \frac12\cdot2 = 1, el centro del intervalo, como esperabas.

Notebook exercises (by hand)

  1. (4x32x+1)dx\int (4x^3 - 2x + 1),dx.
  2. 1e1xdx\int_1^e \frac1x,dx (resultado limpio).
  3. 01e2xdx\int_0^1 e^{2x},dx (sustitución u=2xu=2x).
  4. Verifica que p(x)=2xp(x)=2x en [0,1][0,1] es una densidad válida (01p=1\int_0^1 p = 1).
  5. Calcula la expectation de esa densidad: 01x2xdx\int_0^1 x\cdot 2x,dx.
  6. Explica en palabras por qué E[X]=xp(x)dx\mathbb{E}[X]=\int x,p(x),dx es "el promedio ponderado por probabilidad".

Python lab

import numpy as np, sympy as sp

x = sp.symbols("x")
print(sp.integrate(3*x**2, (x, 0, 2)))            # 8
print(sp.integrate(1/x, (x, 1, sp.E)))            # 1
print(sp.integrate(2*x, (x, 0, 1)))               # 1 (densidad válida)
print(sp.integrate(x*2*x, (x, 0, 1)))             # 2/3 (expectation)

# integral numérica (Monte Carlo) de ∫_0^1 2x dx = 1
n = 1_000_000
u = np.random.rand(n)
print("MC estimate:", (2*u).mean())               # ~1.0

Examen final 📝

Intenta cada nivel antes de abrir las soluciones.

🟡 Medio

  1. 03(2x+1)dx\int_0^3 (2x+1),dx.
  2. 121xdx\int_1^2 \frac{1}{x},dx.

🟠 Medio-difícil

  1. Verifica que f(x)=32x2f(x)=\frac{3}{2}x^2 en [0,c][0,c] es densidad; halla el cc que la normaliza a integrar 1.
  2. Para la densidad exponencial f(x)=λeλxf(x)=\lambda e^{-\lambda x} (x0x\ge0), muestra que 0f(x)dx=1\int_0^\infty f(x),dx=1 y calcula E[X]=0xf(x)dx\mathbb{E}[X]=\int_0^\infty x,f(x),dx (integración por partes).

🔴 Difícil

  1. Demuestra que para la gaussiana estándar, 12πex2/2dx=1\int_{-\infty}^\infty \frac{1}{\sqrt{2\pi}}e^{-x^2/2},dx=1 (usa el truco del cuadrado y coordenadas polares, o al menos esboza el argumento).
  2. La constante de normalización de una distribución (el "partition function" ZZ) es una integral. Explica por qué en ML a veces ZZ es intratable y cómo se conecta esto con la idea de estimar integrales por Monte Carlo (Prob-08).
✅ Soluciones
  1. [x2+x]03=9+3=12[x^2+x]_0^3=9+3=12.
  2. [lnx]12=ln20.693[\ln x]_1^2=\ln 2\approx0.693.
  3. 0c32x2dx=12c3=1c3=2c=21/31.26\int_0^c \frac32 x^2,dx=\frac12 c^3=1\Rightarrow c^3=2\Rightarrow c=2^{1/3}\approx1.26.
  4. 0λeλxdx=[eλx]0=1\int_0^\infty \lambda e^{-\lambda x}dx=[-e^{-\lambda x}]_0^\infty=1. E[X]=0xλeλxdx\mathbb{E}[X]=\int_0^\infty x\lambda e^{-\lambda x}dx; por partes =[xeλx]0+0eλxdx=0+1λ=1λ=[-xe^{-\lambda x}]_0^\infty+\int_0^\infty e^{-\lambda x}dx=0+\frac1\lambda=\frac1\lambda.
  5. Sea I=ex2/2dxI=\int e^{-x^2/2}dx. Entonces I2=e(x2+y2)/2dxdyI^2=\iint e^{-(x^2+y^2)/2}dx,dy; en polares =02π ⁣0er2/2rdrdθ=2π1=2π=\int_0^{2\pi}!\int_0^\infty e^{-r^2/2}r,dr,d\theta=2\pi\cdot1=2\pi. Luego I=2πI=\sqrt{2\pi}, y con el factor 12π\frac{1}{\sqrt{2\pi}} integra a 1.
  6. Z=p(x)dxZ=\int p̃(x),dx(x)dx sobre todo el espacio; en modelos de alta dimensión (p. ej. energy-based models) esta integral no tiene forma cerrada y el espacio es enorme ⇒ intratable. Se estima con muestreo (Monte Carlo / MCMC): promediar muchas muestras aproxima la integral por la LLN.