Formule di Newton-Cotes

Da Wikipedia, l'enciclopedia libera.
Vai alla navigazione Vai alla ricerca

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.

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

[modifica | modifica wikitesto]

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

[modifica | modifica wikitesto]

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

Nota bene: la notazione è un'abbreviazione di e .

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 uguale a zero).

Formule di Newton-Cotes aperte

[modifica | modifica wikitesto]

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 Regola di Masina
3 Nessun nome

Formule composte

[modifica | modifica wikitesto]

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 applica di volta in volta le formule di Newton-Cotes, e sommando poi i risultati. Questo procedimento è chiamato formula composta.

Si ipotizza di dividere l'intervallo in sottointervalli equidistanti mediante i punti , .

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

Voci correlate

[modifica | modifica wikitesto]

Collegamenti esterni

[modifica | modifica wikitesto]
  • (EN) Formule di Newton-Cotes, su www.math-linux.com. URL consultato il 30 aprile 2008 (archiviato dall'url originale il 12 aprile 2008).
  • (EN) Formule di Newton-Cotes [collegamento interrotto], su hades.ph.tn.tudelft.nl.
  • (EN) Formule di Newton-Cotes su MathWorld
  • (EN) Modulo per l'integrazione secondo Newton-Cotes, su math.fullerton.edu. URL consultato il 30 aprile 2008 (archiviato dall'url originale il 9 giugno 2007).
  Portale Matematica: accedi alle voci di Wikipedia che trattano di matematica