Posgrado → Análisis numérico → Integración y diferenciación numérica ↓
Métodos de cuadratura
Los métodos de cuadratura, comúnmente conocidos como integración numérica, son un conjunto de algoritmos para calcular aproximadamente la integral definida de una función. Las integrales definidas son fundamentales en matemáticas, particularmente en cálculo, y son esenciales para encontrar áreas, volúmenes, puntos centrales, entre muchas otras cosas. Los métodos numéricos son necesarios cuando una función es complicada o no es integrable analíticamente.
El objetivo de los métodos de cuadratura es estimar la integral definida de una función en un intervalo. Para una función dada f(x)
, la integral definida entre los límites a
y b
se expresa como:
∫ a b f(x) dx
Una ilustración visual simple de este concepto se puede entender considerando el gráfico de la función f(x)
:
El área sombreada bajo la curva f(x)
entre a
y b
representa la integral definida. En muchos casos, es difícil o imposible encontrar el área analíticamente, por lo que los métodos numéricos entran en juego.
1. Regla del punto medio
La regla del punto medio estima la integral tomando el valor de la función en el punto medio del intervalo y multiplicándolo por el ancho del intervalo. Se representa como:
I ≈ f((a + b) / 2) × (b – a)
Supongamos que queremos aproximar la integral de f(x) = x^2
en el intervalo [1, 3]. Usando la regla del punto medio:
Punto medio: ((1 + 3) / 2) = 2
f(2) = 2^2 = 4
Por lo tanto, la integral aproximada es:
≈ 4 × (3 – 1) = 8
2. Regla del trapecio
La regla del trapecio aproxima el área bajo la curva como un trapecio y calcula su área. La fórmula es:
I ≈ (B – A) × (F(A) + F(B)) / 2
Supongamos que queremos encontrar la integral de la función f(x) = x^3
en [1, 2]:
f(1) = 1^3 = 1
f(2) = 2^3 = 8
Por lo tanto, la integral se estima de la siguiente manera:
I ≈ (2 – 1) × (1 + 8) / 2 = 4.5
3. Regla de Simpson
La regla de Simpson es más precisa y utiliza una parábola para aproximar la curva. La regla de Simpson asume que la función es cuadrática. Se expresa como:
I ≈ (B – A) / 6 × [f(A) + 4f((A + B) / 2) + f(B)]
Probémoslo para f(x) = x^2
en [0, 2]:
f(0) = 0^2 = 0
Punto medio: ((0 + 2) / 2) = 1; f(1) = 1^2 = 1
f(2) = 2^2 = 4
La integral aproximada usando la regla de Simpson es:
I ≈ (2 – 0) / 6 × [0 + 4 × 1 + 4] = 2.67
4. Cuadratura de Gauss
La cuadratura de Gauss mejora la precisión al elegir puntos de muestra óptimos y sus pesos dentro del intervalo. Estos puntos y pesos se obtienen de las raíces de polinomios ortogonales, como los polinomios de Legendre. La fórmula general con n
puntos es:
I ≈ ∑ wi × f( xi )
donde w i
son pesos y x i
son raíces (nodos).
Por ejemplo, usando una cuadratura de Gauss de 2 puntos para ∫ -1 1 x^2 dx
:
Raíces: x 1 = -1/√3, x 2 = 1/√3
Pesos: w 1 = w 2 = 1
Entonces, la integral se puede aproximar como:
I ≈ 1 × ((-1/√3)^2) + 1 × ((1/√3)^2) = 0.6667
Ejemplo de texto
Aquí hay ejemplos simples para mostrar la aplicación de cada método:
Ejemplo 1: Regla del punto medio
Sea f(x) = e^x
en [1, 2]. Aproximar usando la regla del punto medio.
Punto medio: ((1 + 2) / 2) = 1.5
f(1.5) = e^1.5 ≈ 4.4817
Aproximación: 4.4817 × (2 - 1) = 4.4817
Ejemplo 2: Regla del trapecio
Encuentre ∫ 0 π sin(x) dx
usando la regla del trapecio.
f(0) = sin(0) = 0
f(π) = sin(π) = 0
Aproximación: π × (0 + 0) / 2 = 0
Ejemplo 3: Regla de Simpson
Aplica la regla de Simpson para estimar ∫ 0 2 ln(x + 1) dx
.
f(0) = ln(1) = 0
Punto medio f(1) = ln(2) ≈ 0.6931
f(2) = ln(3) ≈ 1.0986
Aproximación: (2/6) × (0 + 4×0.6931 + 1.0986) ≈ 1.7627
Ejemplo 4: Cuadratura de Gauss
Estime ∫ -1 1 x^4 dx
usando cuadratura de Gauss de 2 puntos.
Raíces: x 1 = -1/√3, x 2 = 1/√3
Pesos: w 1 = w 2 = 1
f(x 1 ) = ((-1/√3)^4), f(x 2 ) = ((1/√3)^4)
Aproximación: 1 × f(x 1 ) + 1 × f(x 2 ) = 1/9 + 1/9 = 0.2222
Conclusión
Los métodos de cuadratura son una piedra angular esencial del análisis numérico. Se utilizan ampliamente en una variedad de campos, desde la física hasta las finanzas, donde los cálculos integrales son necesarios pero no es posible realizarlos analíticamente. Estos métodos permiten abordar problemas que involucran funciones complejas y limitan el análisis a resultados numéricos que aún proporcionan información importante sobre los fenómenos subyacentes. Cada método tiene sus propias ventajas y limitaciones, a menudo dependiendo de la naturaleza de la función que se está integrando. Seleccionar el método apropiado puede conducir a cálculos más eficientes y precisos.