Cuadratura de Gauss
En análisis numérico un método de cuadratura es una aproximación de una integral definida de una función. Una cuadratura de Gauss n, es una cuadratura construida para obtener el resultado exacto al integrar polinomios de grado 2n-1 o menos. Para esto selecciona los puntos de evaluación xi y los pesos wi de forma conveniente. La regla suele expresarse para una integral en el intervalo [−1, 1], y viene dada por la siguiente expresión:
En el caso particular en que es un polinomio de grado 2n-1 o menos, la cuadratura de Gauss da el valor exacto de la integral. En el caso general, tal cuadratura dará buenas aproximaciones si puede ser bien aproximada por un polinomio de grado 2n-1 o menos, en el intervalo [−1, 1].
Fórmula para calcular
También conocido como método de Gauss-Legendre, los coeficientes están dados por
donde es el polinomios de Legendre de grado n en el intervalo [−1, 1], y los xi son las raíces de dicho polinomio. La siguiente tabla muestra los valores de los xi y los pesos asociados wi, para distintos valores de n.
Número de puntos, n | Puntos, xi | Pesos, wi |
---|---|---|
1 | 0 | 2 |
2 | 1 | |
3 | 0 | 8⁄9 |
5⁄9 | ||
4 | ||
5 | 0 | 128⁄225 |
Cambio de intervalos
Aplicando un cambio de variable lineal, se puede convertir una integral en el intervalo [a, b], en una integral en [−1, 1]:
Luego se puede aplicar la Cuadratura de Gauss para aproximar la integral en [−1, 1]:
Ejemplo
Aproxime la integral de 1 a 5 cuando n = 2 mediante el método de cuadratura de Gauss y después comparelo con el resultado exacto.
Con podemos resolver la integral con exactitud para todos los polinomios de grado igual o menor a 3 para f(x)
Referencias
- https://web.archive.org/web/20100327115656/http://www.tonahtiu.com/notas/metodos/cuadratura_gauss.htm
- Burden, Richard L.; Faires, J. Douglas (1998). «§4.7: Cuadratura gaussiana». Análisis numérico. ISBN 968-7529-46-6.
Enlaces externos
- Esta obra contiene una traducción derivada de «Gaussian quadrature» de Wikipedia en inglés, publicada por sus editores bajo la Licencia de documentación libre de GNU y la Licencia Creative Commons Atribución-CompartirIgual 4.0 Internacional.
Datos: Q767680