Anar al contingut

Quadratura de Gauss

De L'Enciclopèdia, la wikipedia en valencià

En anàlisis numèric un método de quadratura és una aproximació d'una integral definida d'una funció. Una quadratura de Gauss n, és una quadratura construïda per a obtindre el resultat exacte en integrar polinomis de grau 2n-1 o menys. Per a açò selecciona els punts d'evaluació xi i els pesos wi de forma convenient. La regla sol expressar-se per a una integral en l'interval [−1, 1], i ve donada per la següent expressió:

∫−11f(x)dx≈∑i=1nwif(xi)

En el cas particular en que f(x) és un polinomi de grau 2n-1 o menys, la quadratura de Gauss dona el valor exacte de l'integral. En el cas general, tal quadratura donarà bones aproximacions si f(x) pot ser ben aproximada per un polinomi de grau 2n-1 o menys, en l'interval [−1, 1].

Fòrmula per a calcular wi

També conegut com a método de Gauss-Legendre, els coeficients estan donats per

wi=2(1−xi2)[P'n(xi)]2;

a on Pn és el polinomis de Legendre de grau n en l'interval [−1, 1], i els xi són les raïls de dit polinomi. La següent taula mostra els valors dels xi i els pesos associats wi, per a distints valors de n.

Número de punts, n Punts, xi Pesos, wi
1 0 2
2 ±1/3 1
3 0 8⁄9
±3/5 5⁄9
4 ±(3−26/5)/7 18+3036
±(3+26/5)/7 18−3036
5 0 128⁄225
±135−210/7 322+1370900
±135+210/7 322−1370900

Canvi d'intervals

Aplicant un canvi de variable llineal, es pot convertir una integral en l'interval [a, b], en una integral en [−1, 1]:

∫abf(x)dx=b−a2∫−11f(b−a2u+a+b2)du

Després es pot aplicar la Quadratura de Gauss per a aproximar l'integral en [−1, 1]:

∫abf(x)dx≈b−a2∑i=1nwif(b−a2ui+a+b2)

Referències