GuíaIntegraciónavanzado

Cuadraturas de Gauss

Cómo las cuadraturas de Gauss eligen nodos y pesos óptimos mediante polinomios ortogonales: Legendre, Chebyshev, Laguerre y Hermite.

Idea general

Newton-Cotes fija nodos equiespaciados. Gauss hace lo contrario: permite mover los nodos y calcula pesos para que la regla sea exacta para polinomios del mayor grado posible.

∫abw(x)f(x) dx≈∑i=1ncif(xi),w(x)>0\int_a^b w(x)f(x)\,dx\approx\sum_{i=1}^{n} c_i f(x_i),\qquad w(x)>0
E=f(2n)(ξ)(2n)!∫abpn2(x)w(x) dxE=\frac{f^{(2n)}(\xi)}{(2n)!}\int_a^b p_n^2(x)w(x)\,dx

Gauss-Legendre

Gauss-Legendre usa w(x)=1w(x)=1 en [−1,1][-1,1]. Para aplicar la regla en [a,b][a,b] se transforma el intervalo de integración y se multiplica por el jacobiano.

x=b−a2y+b+a2,dx=b−a2 dyx=\frac{b-a}{2}y+\frac{b+a}{2},\qquad dx=\frac{b-a}{2}\,dy
∫abf(x) dx≈b−a2∑i=1ncif ⁣(b−a2xi+b+a2)\int_a^b f(x)\,dx\approx\frac{b-a}{2}\sum_{i=1}^{n}c_i f\!\left(\frac{b-a}{2}x_i+\frac{b+a}{2}\right)
nNodosPesos
2-0.577350, 0.5773501.000000, 1.000000
30.000000, -0.774597, 0.7745970.888889, 0.555556, 0.555556
4-0.339981, -0.861136, 0.339981, 0.8611360.652145, 0.347855, 0.652145, 0.347855
50.000000, -0.538469, -0.906180, 0.538469, 0.9061800.568889, 0.478629, 0.236927, 0.478629, 0.236927

Otras familias de Gauss

La familia se elige mirando el peso y el dominio. Si la integral ya contiene una singularidad o una exponencial, conviene no pelearse con ella: se incorpora al peso.

FamiliaPeso y dominioPolinomios
Legendrew(x)=1, [-1,1]p0=1, p1=x, recurrence by Legendre
Chebyshevw(x)=1/sqrt(1-x^2), [-1,1]T0=1, T1=x, Tk=2xT(k-1)-T(k-2)
Laguerrew(x)=e^(-x), [0,+infty)L0=1, L1=1-x, recurrence by Laguerre
Hermitew(x)=e^(-x^2), (-infty,+infty)H0=1, H1=2x, H(k+2)=2xH(k+1)-2(k+1)Hk
∫−11f(x)1−x2 dx≈πn∑i=1nf(xi),xi=cos⁡ ⁣(π2i−12n)\int_{-1}^{1}\frac{f(x)}{\sqrt{1-x^2}}\,dx\approx\frac{\pi}{n}\sum_{i=1}^{n}f(x_i),\qquad x_i=\cos\!\left(\pi\frac{2i-1}{2n}\right)

Deducción de Gauss-Legendre con dos puntos

DeducciónDeducción: Gauss-Legendre de dos puntosVer como página propia →

  1. Buscamos una regla simétrica en [−1,1][-1,1] con nodos −r-r y rr y pesos iguales cc, porque el intervalo y el peso w(x)=1w(x)=1 son simétricos.

    ∫−11f(x)dx≈cf(−r)+cf(r)\int_{-1}^{1}f(x)dx\approx c f(-r)+c f(r)

  2. Exigimos exactitud para f(x)=1f(x)=1:

    2=2c⟹c=12=2c\quad\Longrightarrow\quad c=1

  3. Las funciones impares xx y x3x^3 se satisfacen automáticamente por simetría. Exigimos exactitud para f(x)=x2f(x)=x^2:

    ∫−11x2dx=23=r2+r2=2r2\int_{-1}^{1}x^2dx=\frac{2}{3}=r^2+r^2=2r^2

  4. Despejando r aparece la regla de dos puntos:

    r=13,∫−11f(x)dx≈f ⁣(−13)+f ⁣(13)r=\frac{1}{\sqrt{3}},\qquad \int_{-1}^{1}f(x)dx\approx f\!\left(-\frac{1}{\sqrt{3}}\right)+f\!\left(\frac{1}{\sqrt{3}}\right)