GidaIntegrazioaaurreratua

Gaussen koadraturak

Gaussen koadraturek polinomio ortogonalen bidez nodo eta pisu optimoak nola aukeratzen dituzten: Legendre, Chebyshev, Laguerre eta Hermite.

Ideia orokorra

Newton-Cotesek nodo ekidistanteak finkatzen ditu. Gaussen metodoak alderantziz egiten du: nodoak mugitzen uzten du eta pisuak kalkulatzen ditu erregela ahalik eta gradu handieneko polinomioetarako zehatza izan dadin.

∫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-Legendrek w(x)=1w(x)=1 erabiltzen du [−1,1][-1,1] tartean. Erregela [a,b][a,b]-n aplikatzeko, integrazio-tartea eraldatzen da eta jakobiarrarekin biderkatzen da.

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)
nNodoakPisuak
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

Gaussen beste familiak

Familia pisuari eta eremuari begiratuz aukeratzen da. Integralak jada singularitate bat edo esponentzial bat badu, ez da haren aurka egin behar: pisuan sartzen da.

FamiliaPisua eta eremuaPolinomioak
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)

Bi puntuko Gauss-Legendreren dedukzioa

FrogapenaFrogapena: bi puntuko Gauss-LegendreIkusi orri propio gisa →

  1. [−1,1][-1,1] tartean −r-r eta rr nodoak eta cc pisu berdinak dituen erregela simetrikoa bilatzen dugu, tartea eta w(x)=1w(x)=1 pisua simetrikoak direlako.

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

  2. f(x)=1f(x)=1-erako zehaztasuna eskatzen dugu:

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

  3. xx eta x3x^3 funtzio bakoitiak automatikoki betetzen dira simetriagatik. f(x)=x2f(x)=x^2-rako zehaztasuna eskatzen dugu:

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

  4. r askatuz bi puntuko erregela agertzen da:

    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)