Formule di Newton-Coves

Formula dei trapezi[modifica | modifica wikitesto]

Considero una funzione positiva crescente, considero la retta passante per e e calcolo l'area del trapezio sotteso, ponendo . Il trapezio ha area

Osservo che il polinomio interpolante nei nodi e è:
quindi
si approssima con
Supponiamo che , allora
(deriva dall'espressione dell'errore analitico per l'interpolazione di polinomi)
Quindi si ha
Nel caso della formula dei trapezi, ha segno negativo sull'intervallo . Posso quindi applicare il teorema del valor medio.

Teorema 6.1

Supponiamo di avere continua, integrabile, che non cambia segno in . Allora esiste tale che

 

Allora, riconsiderando la formula dell'errore, siccome la funzione peso non cambia segno si applica il teorema del valor medio, esisterà indipendente da x tale che

Formula dei trapezi composita[modifica | modifica wikitesto]

Se l'intervallo è grande, l'errore potrebbe essere elevato, quindi si pensa di passare alla formula composita. Suddivido quindi in sottointervalli, considero la partizione

Pongo (ipotesi semplificativa: intervalli di uguale ampiezza). Si calcola l'area del trapezio su ogni intervallino, e sommando i vari contributi ottengo:
e siccome i nodi interni compaiono due volte posso riscriverlo come
è la distanza tra due nodi successivi, e più la partizione è fine, più l'approssimazione è buona.
Sommando i vari contributi dell'errore analitico, si ha:
(come dimostrato prima è l'integrale del polinomio nodale) e moltiplicando e dividendo per m:
e per il teorema del valor medio discreto, assumendo che sia continuo,
e siccome ottengo il risultato cercato.

Formula di Simpson[modifica | modifica wikitesto]

Considero la funzione , sia il punto medio dell'intervallo. Considerando il polinomio di grado 2 interpolante nei tre nodi equispaziati, si ha

con (passo di equispaziatura dei nodi).
Considerando i polinomi di Lagrange di grado 2 ottengo i pesi e .

Modalità di scrittura dei pesi[modifica | modifica wikitesto]

Osservazione 6.1

Considero un esempio di formula di quadratura con nodi equispaziati. I pesi della formula di quadratura si possono sempre scrivere come dove gli non dipendono dall'intervallo .

 

Posso scrivere i nodi come

Il polinomio generico di grado si scrive come
dove
Considero un cambio di variabili.
tale che .
L'espressione del polinomio di Lagrange diventa
Tenendo conto che , il peso si può quindi riscrivere come
e definisco
in modo che con nel caso di nodi equispaziati.
Inoltre , cioè i pesi sono simmetrici. Questo deriva dal fatto che anche i nodi equispaziati sull'intervallo sono simmetrici.

Metodo dei coefficienti indeterminati[modifica | modifica wikitesto]

Fisso nodi distinti, e impongo che la formula di quadratura abbia grado di precisione .

Inoltre pongo , allora per la formula di quadratura deve valere:

che si può riscrivere come
Ottengo il sistema lineare
dove il termine noto ha componenti
La matrice a primo membro è la matrice di Vandermonde, e se i nodi sono distinti quindi esiste unica la soluzione.

Esempio 6.1 (metodo di Simpson)

Considero un intervallo con nodi . La distanza fra i nodi è .

Il sistema corrispondente è:

Risolvo il sistema:
Siccome , si ha
e si ritrova la formula di Simpson.

 

Grado di precisione[modifica | modifica wikitesto]

Teorema 6.2

Supponiamo , allora

 

Questa formula ha grado di precisione 3 (tutti i polinomi di grado 3 hanno derivata quarta nulla, e quindi errore nullo)

Tutte le formule interpolatorie su nodi equispaziati quando il polinomio interpolante ha grado pari hanno un grado di precisione in più di quanto ci si aspetta, Questo deriva dalla simmetria nei nodi.

Vale in generale che dato come polinomio interpolante, allora il grado di precisione è .

Formula di Simpson composita[modifica | modifica wikitesto]

Definisco , e , . Segue che

(i punti medi vengono contati una volta, gli altri punti due volte)
Per , l'errore decresce.

La formula di quadratura di Simpson e dei trapezi sono dette formule chiuse perché usano gli estremi dell'intervallo . Altre formule vengono dette aperte perché tutti i nodi di quadratura sono interni all'intervallo.

Formula del punto medio[modifica | modifica wikitesto]

Considero la funzione , considero il punto medio e approssimo la funzione con la retta per quel punto. Ottengo quindi:

  1. formula semplice: , dove .
  2. Errore nella formula semplice: , con gradi di precisione .
  3. Formula composita:
  4. Errore nella formula composita: , dove , e .

Considero generico.

dove è il polinomio interpolante di grado , quindi
e .

Osservazione 6.2

Per grado del polinomio interpolante fino a 7, i coefficienti sono tutti positivi, mentre per i pesi possono anche essere negativi (si può quindi avere la possibilità dell'errore di cancellazione).

 

Vale il seguente teorema:

Teorema 6.3

La formula di quadratura con grado di precisione almeno , coincide con la formula di quadratura interpolatoria che usa come nodi di interpolazione .

 
Dimostrazione

è la formula interpolatoria, con polinomio interpolante di grado , segue che il grado di precisione è almeno . Sia il polinomio interpolante di grado di . Allora esiste unica la soluzione del sistema con matrice di Vandermonde, e quindi le due formule di quadratura coincidono.

 

Stime dell'errore a priori[modifica | modifica wikitesto]

Definizione 6.3

L'ordine di infinitesimo di una formula di quadratura rispetto a è il massimo tale che

 


Teorema 6.4

Sia la formula di quadratura interpolatoria di Newton-Coves (nodi equispaziati e polinomio interpolante di grado ), il passo di equispaziatura. Se il grado del polinomio interpolante è pari e allora si dimostra che

(grado di precisione e ordine di infinitesimo ) con
con polinomio nodale, e questa quantità è negativa.

Se invece è dispari e si ha

(grado di precisione e ordine di infinitesimo ) con

 

Ad esempio, nel caso della formula di Simpson (tre nodi, quindi ) dove il polinomio interpolante ha grado 2, il grado di precisione è , e l'ordine di

infinitesimo è .
Teorema 6.5

Per le formule composite, segue che, se è pari,

e il grado di precisione è e l'ordine di infinitesimo è , mentre per dispari,
e il grado di precisione è e l'ordine di infinitesimo è .

In ogni caso

quindi non ci sono problemi di convergenza.

 

Raddoppio del numero dei sottointervalli[modifica | modifica wikitesto]

Consideriamo cosa succede se si raddoppia il numero dei sottointervalli. Nel caso della formula dei trapezi:

  1. Al primo passo, , i nodi sono e , e pongo .
  2. Per , considero il punto medio e chiamo .
  3. Al passo successivo, (raddoppio il numero dei sottointervalli) e divido i due sottointervalli nuovamente a metà, in modo che
    dove sono punti medi dei due sosttointervallini.
  4. Ad ogni livello vengono introdotti nuovi nodi.

Si ha quindi

dove gli sono i nuovi nodi.

Esercizio 6.1

Nel caso della formula di Simpson, raddoppiano i sottointervalli

Infatti,

  1. al passo 1,
  2. Al passo 2, pongo :
    dove sono i punti medi di e .
    cioè
  3. al passo , avendo nodi, si ha:
 
 PrecedenteSuccessivo