Newton–Cotes-formula

A Wikipédiából, a szabad enciklopédiából
Ugrás a navigációhoz Ugrás a kereséshez

A numerikus analízisben a Newton–Cotes-kvadratúraformulák (amiket Newton–Cotes-szabályoknak is neveznek) olyan képletek csoportja, amelyek numerikus integrálásra (más néven kvadratúrákra) szolgálnak, amik alapjául az integrálni kívánt intervallumot n+1 egyenlő távolságra (lépésközre) osztunk fel. Ezeket a módszereket Isaac Newtonról és Roger Cotesról nevezték el.

A Newton–Cotes-kvadratúraformulák nagyon hatékonyak, ha meghatározottak a függvényértékek az adott ekvidisztáns (egymástól egyenlő távolságra lévő) pontokban. Ha lehetséges más pontokban is meghatározni a függvény értékét, akkor vannak hatékonyabb, célszerűbb módszerek is az integrál kiszámítására, pl.: Gauss-kvadratúra és a Clenshaw–Curtis-kvadratúra.

Leírás[szerkesztés]

Feltesszük, hogy az f függvény értékeit ismerjük az egymástól egyenlő távolságra lévő (ekvidisztáns) xi pontokban, ahol i = 0, …, n. A Newton–Cotes-kvadratúraformuláknak két típusát különböztetjük meg: a "nyitott" típus csak a belső alappontokkal ( x1, x2,…,xn-1 ) számol, míg a "zárt" típus az intervallum kezdő- és végpontját is beleveszi (tehát az alappontok: x0, x1, x2,…, xn-1, xn ). Az n-ed fokú Newton–Cotes-formula általános alakja:

ahol xi = h i + x0, h (más néven lépésköz) egyenlő (xnx0)/n-nel. wi a súlyokat jelöli.

Ahogy az alábbiakban látható, a súlyozást a Lagrange-alappolinomokból kaptuk. Ez azt jelenti, hogy a súlyok csak az xi -től függnek, és nem az adott f függvény értékétől. Legyen L(x) a Lagrange interpolációs polinom megadott alappontjai (x0, f(x0) ), …, (xn, f(xn) ), akkor

A nyitott n-ed fokú Newton–Cotes-formula:

A súlyokat hasonlóan kereshetjük meg, mint a zárt képletnél.

Instabilitás magas alappontszám esetén[szerkesztés]

A Newton–Cotes-kvadratúraformulákat minden n alappontra el lehet készíteni. Igaz ugyan, hogy a Newton-Cotes szabály alkalmazása magas alappontszám esetén néha a Runge-jelenséget vonhatja maga után, ahol a hibatag exponenciálisan növekszik. Más módszerek, mint a Gauss-kvadratúra és a Clenshaw–Curtis-kvadratúra ahol nem egyenlő a lépésköz (és a pontok az intervallum végpontjaiban torlódnak), sokkal stabilabbak, pontosabbak és elfogadottabbak, mint a Newton–Cotes-módszer. Ha ezeket a módszereket nem tudjuk használni, mert a függvényértékek csak ekvidisztáns (egyenlő lépésközű) pontokban ismertek, akkor a Runge-jelenséget az alábbiakban ismertetett összetett formulával lehet elkerülni.

Zárt Newton–Cotes-formulák[szerkesztés]

Ez a táblázat néhány zárt Newton–Cotes-formulát tartalmaz. Az az . jelölést rövidíti.

Zárt Newton–Cotes-formulák
Fok/alappontok száma Elnevezés Képlet Hibatag
1 Trapézszabály
2 Simpson-szabály
3 Simpson 3/8 szabály
4 Boole-szabály, vagy
Bode-szabály (sic)

Boole szabályát egy régebbi sajtóhiba miatt tévesen nevezik Bode szabályának, ami Abramowitz és Stegun egy korai kézikönyvében volt található.[1].

A lépésköz exponenciális növelése a hibatagban, a közelítés hibáját csökkenti. Az f függvény deriváltja a hibatagban megmutatja, hogy melyik polinomot lehet pontosan integrálni (ahol a hibatag egyenlő nullával). Minden második lépésnél a hibatag pontossága kettővel nő. A hibatag (ξ) a és b közé esik.


Nyitott Newton–Cotes-formulák[szerkesztés]

Az alábbi táblázat felsorol néhány nyitott Newton-Cotes kvadratúra formulát.

Nyitott Newton–Cotes-formulák
Fok/alappontok száma Elnevezés Képlet Hibatag
0 Téglalapszabály
1
2
3

Összetett szabályok[szerkesztés]

Ahhoz, hogy a Newton–Cotes-szabály pontos értéket adjon, a h lépésköznek kicsinek kell lennie, ami azt jelenti, hogy az integrációs intervallumnak is kicsinek kell lennie, ami az esetek nagy részében nem fordul elő. Éppen ezért, a numerikus integrálást az adott intervallumon felosztják kisebb részintervallumokra, majd ezeken alkalmazzák a Newton–Cotes-szabályt, végül ezeken a részintervallumokon kapott eredményeket összegzik és kapják a végső eredményt. Ezt nevezzük összetett szabálynak, lásd Numerikus integrálás.

Hivatkozások[szerkesztés]

  • M. Abramowitz és I. A. Stegun, eds. Handbook of Mathematical Functions with Formulas, Graphs, and Mathematical Tables. New York: Dover, 1972. (lásd 25.4.)
  • George E. Forsythe, Michael A. Malcolm, and Cleve B. Moler. Computer Methods for Mathematical Computations. Englewood Cliffs, NJ: Prentice Hall, 1977. (lásd 5.1.)
  • William H. Press, Brian P. Flannery, Saul A. Teukolsky, William T. Vetterling. Numerical Recipes in C. Cambridge, UK: Cambridge University Press, 1988. (lásd 4.1.)
  • Josef Stoer and Roland Bulirsch. Introduction to Numerical Analysis. New York: Springer-Verlag, 1980. (lásd 3.1.)

Külső hivatkozások[szerkesztés]