Formule di Newton-Cotes: differenze tra le versioni

Da Wikipedia, l'enciclopedia libera.
Vai alla navigazione Vai alla ricerca
Contenuto cancellato Contenuto aggiunto
Matti90 (discussione | contributi)
Riga 32: Riga 32:
''Nota bene: la notazione <math>f_i</math> è un'abbreviazione di <math>f(x_i)</math>.
''Nota bene: la notazione <math>f_i</math> è un'abbreviazione di <math>f(x_i)</math>.
{| {{prettytable|text-align=center}}
{| class="wikitable" style="margin:1em auto 1em auto; background:white"
|-
|+ '''Formule di Newton-Cotes chiuse'''
! Grado !! Nome comune !! Formula !! Termine d'errore
! Grado !! Nome comune !! Formula !! Termine d'errore
|-
|- align="center"
| 1 || [[Regola del trapezio]] || <math> \frac{h}{2} (f_0 + f_1) </math> || <math>-\frac{h^3}{12}\,f^{(2)}(\xi)</math>
| 1 || [[Regola del trapezio]] || <math> \frac{h}{2} (f_0 + f_1) </math> || <math>-\frac{h^3}{12}\,f^{(2)}(\xi)</math>
|-
|- align="center"
| 2 || [[Regola di Cavalieri-Simpson]] || <math> \frac{h}{3} (f_0 + 4 f_1 + f_2) </math> || <math>-\frac{h^5}{90}\,f^{(4)}(\xi)</math>
| 2 || [[Regola di Cavalieri-Simpson]] || <math> \frac{h}{3} (f_0 + 4 f_1 + f_2) </math> || <math>-\frac{h^5}{90}\,f^{(4)}(\xi)</math>
|-
|- align="center"
| 3 || Regola di Cavalieri-Simpson <br /> (con fattore 3/8) || <math> \frac{h}{8} (3 f_0 + 9 f_1 + 9 f_2 + 3 f_3) </math> || <math>-\frac{3h^5}{80}\,f^{(4)}(\xi)</math>
| 3 || Regola di Cavalieri-Simpson<br />(con fattore 3/8) || <math> \frac{h}{8} (3 f_0 + 9 f_1 + 9 f_2 + 3 f_3) </math> || <math>-\frac{3h^5}{80}\,f^{(4)}(\xi)</math>
|-
|- align="center"
| 4 || Regola di Boole || <math> \frac{h}{90} (7 f_0 + 32 f_1 + 12 f_2 + 32 f_3 + 7 f_4) </math> || <math>-\frac{8h^7}{945}\,f^{(6)}(\xi)</math>
| 4 || [[Regola di Boole]] || <math> \frac{h}{90} (7 f_0 + 32 f_1 + 12 f_2 + 32 f_3 + 7 f_4) </math> || <math>-\frac{8h^7}{945}\,f^{(6)}(\xi)</math>
|}
|}



Versione delle 16:24, 16 apr 2012

In analisi numerica, le formule di Newton-Cotes sono un gruppo di formule adoperate nell'integrazione numerica (detta anche quadratura) che si basano sulla valutazione dell'integrando in n+1 punti equidistanti. Le formule sono chiamate così in onore di Isaac Newton e Roger Cotes.

Le formule di Newton-Cotes possono essere utili se il valore dell'integrando nei punti equidistanti è noto. Se è possibile modificare i punti dove è valutato l'integrando, allora è preferibile usare altri metodi come la quadratura di Gauss.

Descrizione

Si assume che il valore di una funzione f è noto nei punti equidistanti xi, per i = 0, …, n. Esistono due tipi di formule di Newton-Cotes: la forma "chiusa", che valuta il valore della funzione in tutti i punti, e la forma "aperta", che non considera i valori della funzione nei suoi estremi. La formula di Newton-Cotes chiusa di grado n è definita come:

dove xi = h i + x0, con h (chiamato passo) uguale a (xnx0)/n. Le wi sono chiamate pesi.

Come si può notare nella seguente derivazione, i pesi provengono dai polinomi di Lagrange. Ciò significa che dipendono solo dalle xi e non dalla funzione f. Sia L(x) il polinomio di Lagrange di interpolazione per i punti noti (x0, f(x0) ), …, (xn, f(xn) ); allora:

La formula di Newton-Cotes aperta di grado n è definita invece come:

I pesi si ricavano in maniera simile alla formula chiusa.

Instabilità per gradi molto grandi

Le formule di Newton-Cotes si possono ricavare per qualsiasi grado n. Tuttavia, per n molto grandi, la formula può risentire del fenomeno di Runge, dove l'errore cresce in maniera esponenziale se n è elevato. In tal caso si adottano solitamente metodi molto stabili, come ad esempio la quadratura di Gauss con punti di integrazione non equidistanti. Se poi non possono essere usati neanche questi, perché l'integrando è definito solo su punti equidistanti, allora il fenomeno di Runge può essere evitato usando una formula composta, come descritto in seguito.

Formule di Newton-Cotes chiuse

Nella tabella sono elencate alcune formule di Newton-Cotes di tipo chiuso.

Nota bene: la notazione è un'abbreviazione di .

Grado Nome comune Formula Termine d'errore
1 Regola del trapezio
2 Regola di Cavalieri-Simpson
3 Regola di Cavalieri-Simpson
(con fattore 3/8)
4 Regola di Boole

L'esponente del passo h nel termine d'errore mostra il tasso di decremento dell'errore di approssimazione. La derivata di f nel termine d'errore mostra quali polinomi possono essere integrati esattamente (cioè con errore pari a zero).

Formule di Newton-Cotes aperte

Nella tabella sono elencate alcune formule di Newton-Cotes di tipo aperto.

Formule di Newton-Cotes aperte
Grado Nome comune Formula Termine d'errore
0 Regola del rettangolo
1 Nessun nome
2 Nessun nome
3 Nessun nome

Formule composte

Per ottenere una certa accuratezza dalle formule di Newton-Cotes, il passo h deve essere piccolo; ciò significa che l'intervallo di integrazione dovrà essere anch'esso piccolo, il che non è sempre vero. Per questo motivo, di solito si sceglie di calcolare l'integrale dividendo l'intervallo in tanti piccoli sottointervalli, ai quali si applicano di volta in volta le formule di Newton-Cotes, e sommando poi i risultati. Questo procedimento è chiamato formula composta.

Voci correlate

Collegamenti esterni

  Portale Matematica: accedi alle voci di Wikipedia che trattano di matematica