Regola di Cavalieri-Simpson

Per regola di Cavalieri-Simpson o regola di Cavalieri o regola di Simpson si intende un metodo per il calcolo numerico approssimato di integrali definiti della forma:

I:=\int _{a}^{b}f(x)\,dx.

Come tutti i procedimenti per il calcolo approssimato di integrali definiti e per altri calcoli approssimati a partire da funzioni di variabile reale, tale metodo si utilizza per funzioni $\,f\left(x\right)$ delle quali non si conosce la funzione primitiva, oppure della cui primitiva si conoscono solo caratteristiche dalle quali non si riesce a ricavare un'espressione tramite funzioni elementari che possa essere ragionevolmente utilizzata per i calcoli richiesti. Questi metodi approssimati si utilizzano inoltre nei casi in cui non è nota un'espressione analitica della funzione da integrare, ma si conoscono soltanto alcuni suoi valori (ottenuti sperimentalmente o ricavati da altre fonti), oppure quando è noto soltanto il suo diagramma (tracciato con l'ausilio di appositi strumenti o ricavato dalla letteratura).

La formula di quadratura o metodo delle parabole[modifica | modifica wikitesto]

La regola di Cavalieri-Simpson prevede la suddivisione dell'intervallo di integrazione in sottointervalli e la sostituzione in questi sottointervalli della funzione integranda mediante archi di parabola, cioè mediante polinomi quadratici.

Consideriamo dunque $\,\int _{a}^{b}f(x)\,dx.$ Per semplicità supponiamo $f(x)\geq 0$ in tutto l'intervallo di integrazione $\left[a,b\right].$

Suddividiamo $\left[a,b\right]$ in un numero pari $n=2m$ di sottointervalli ciascuno di ampiezza $h:={\frac {b-a}{n}}.$ Poniamo:

x_{i}:=a+i\cdot h,\qquad y_{i}:=f\left(x_{i}\right),\qquad {\text{per }}i=0,1,\ldots ,n,

per gli estremi dei successivi sottointervalli e per i valori che la funzione assume in loro corrispondenza.

Consideriamo anche l'intervallo parziale formato da due sottointervalli consecutivi avente come estremi $x_{0}$ e $x_{2}.$ Oltre a questo consideriamo anche i successivi $m-1$ intervalli parziali aventi come estremi rispettivamente $x_{2}$ e $x_{4},\ldots ,x_{n-2}$ e $x_{n}.$

In ciascuno di questi intervalli parziali ci proponiamo di sostituire $f\left(x\right)$ con una funzione razionale intera di secondo grado. Cominciamo dal primo intervallo parziale e scegliamo un polinomio della forma

y\left(x\right):=A\left(x-x_{1}\right)^{2}+B\left(x-x_{1}\right)+C

in modo che il suo integrale tra $x_{0}$ e $x_{2}$ differisca da quello della funzione originale di una quantità che possa risultare trascurabile.

L'espressione polinomiale sostitutiva rappresenta una generica parabola con asse di simmetria verticale. Per determinare il valore delle costanti $A,B$ e $C$ si impone il passaggio della parabola per i punti di coordinate:

\left(x_{0},\,y_{0}:=f(x_{0})\right),\qquad (x_{1},\,y_{1}:=f(x_{1})),\qquad (x_{2},\,y_{2}:=f(x_{2})).

In questo modo la parabola è univocamente determinata dal seguente sistema di equazioni lineari (nelle incognite $A$ , $B$ e $C$ ):

{\begin{cases}y_{0}=A(x_{0}-x_{1})^{2}+B(x_{0}-x_{1})+C\\y_{1}=C\\y_{2}=A(x_{2}-x_{1})^{2}+B(x_{2}-x_{1})+C,\end{cases}}

da cui risulta:

A={\frac {y_{0}+y_{2}-2y_{1}}{2h^{2}}},\quad B={\frac {y_{2}-y_{0}}{2h}},\quad C=y_{1}.

Detto $J_{1}$ il valore dell'integrale di $y(x)$ nell'intervallo $\left[x_{0},x_{2}\right]$ , si ha:

J_{1}=\int _{x_{0}}^{x_{2}}y(x)\,dx=\int _{x_{0}}^{x_{2}}\left[A\left(x-x_{1}\right)^{2}+B\left(x-x_{1}\right)+C\right]\,dx

=

\left[{\frac {A\left(x-x_{1}\right)^{3}}{3}}+{\frac {B\left(x^{2}-2x_{1}x+x_{1}^{2}\right)}{2}}+Cx\right]_{x_{0}}^{x_{2}}={\frac {2Ah^{3}}{3}}+2Ch.

Sostituendo i valori di $A,B$ e $C$ ricavati dal sistema, si ottiene

\int _{x_{0}}^{x_{2}}f(x)\,dx\simeq J_{1}=\int _{x_{0}}^{x_{2}}y(x)\,dx={\frac {h}{3}}\left(y_{0}+4y_{1}+y_{2}\right).

Operando in modo analogo per il calcolo degli integrali dei polinomi nei successivi $m-1$ intervalli parziali e sommando i valori ottenuti sugli $m$ intervalli parziali e per l'intero intervallo d'integrazione, si ottiene

J=\sum _{k=1}^{m}J_{k}={\frac {h}{3}}\left[\left(y_{0}+4y_{1}+y_{2}\right)+\left(y_{2}+4y_{3}+y_{4}\right)+\ldots +\left(y_{n-2}+4y_{n-1}+y_{n}\right)\right]

e quindi

\int _{a}^{b}f(x)\,dx\simeq J={\frac {h}{3}}\left(y_{0}+4y_{1}+2y_{2}+4y_{3}+2y_{4}+\ldots +2y_{n-2}+4y_{n-1}+y_{n}\right).

Tale formula prende il nome di formula di Cavalieri-Simpson o formula di quadratura delle parabole.

Gli errori[modifica | modifica wikitesto]

Il metodo di quadratura di Cavalieri-Simpson, come ogni metodo di approssimazione numerica, è suscettibile di errori. Oltre all'errore dovuto alla sostituzione della funzione integranda con una successione di funzioni approssimanti, i polinomi quadratici, intrinseco al metodo utilizzato, si riscontrano anche errori dovuti all'arrotondamento dei valori $y_{i}$ che vengono concretamente calcolati con strumenti che inevitabilmente operano con precisione limitata.

Per ridurre al minimo questi ultimi, è consigliabile:

scegliere un passo di integrazione $h$ con un numero finito di cifre decimali;
eseguire i calcoli con un numero di cifre decimali almeno doppio di quello delle cifre che si desiderano esatte nel risultato.

Indicando l'errore intrinseco al metodo $e:=J-I$ , si può dimostrare che:

e\simeq kh^{4},

dove $k$ è una costante che dipende dalla funzione integranda e dall'intervallo d'integrazione. La regola di Cavalieri-Simpson è dunque un metodo del quarto ordine.

Di questo errore può essere molto utile conoscere una maggiorazione. La valutazione accurata di tale maggiorazione non è semplice, poiché richiede di calcolare la derivata quarta della funzione integranda. Per molte funzioni integrande date analiticamente il calcolo della derivata quarta risulta molto oneroso; per funzioni note empiricamente la stessa valutazione della derivata quarta costituisce di per sé un problema di calcolo numerico tendenzialmente oneroso. Di conseguenza in genere per la valutazione dell'errore si preferisce ricorrere a metodi empirici: il più noto e utilizzato è il metodo del dimezzamento del passo.

Da quanto osservato precedentemente segue che, applicando il metodo di Cavalieri-Simpson con un passo di integrazione $h$ , si ottiene l'approssimazione che ora indichiamo con $J_{(1)}$ affetta da un errore che scriviamo $e_{(1)}=J_{(1)}-I\simeq kh^{4}.$

Utilizzando il passo di integrazione ${\frac {h}{2}}$ , si otterrà il valore approssimato $J_{(2)}$ con errore: $e_{(2)}=J_{(2)}-I\simeq k{\frac {h^{4}}{16}}.$

Da tali relazioni segue

e_{(1)}-e_{(2)}=J_{(1)}-J_{(2)}\simeq k\left(h^{4}-{\frac {h^{4}}{16}}\right),

da cui risulta:

k\simeq {\frac {16\left(J_{(1)}-J_{(2)}\right)}{15h^{4}}}.

Poiché, trascurando le approssimazioni da arrotondamento, l'approssimazione migliore è data da $J_{(2)},$ sostituendo il valore di $k$ in $e_{(2)}\simeq k{\frac {h^{4}}{16}},$ si ottiene:

e_{(2)}\simeq {\frac {16\left(J_{(1)}-J_{(2)}\right)}{15h^{4}}}\cdot {\frac {h^{4}}{16}}

da cui

e_{(2)}\simeq {\frac {J_{(1)}-J_{(2)}}{15}}.

Si può quindi assumere il valore assoluto di $e_{(2)}$

{\frac {\left|J_{(1)}-J_{(2)}\right|}{15}}

come maggiorazione dell'errore assoluto

e_{a}=\left|J_{(2)}-I\right|.

È interessante osservare che, se le approssimazioni $J_{(1)}$ e $J_{(2)}$ coincidono per le prime $r$ cifre decimali, risulta:

\left|J_{(1)}-J_{(2)}\right|<10^{-r},

da cui segue

{\frac {\left|J_{(1)}-J_{(2)}\right|}{15}}<{\frac {10^{-r}}{15}}=0,666\ldots \cdot 10^{-(r+1)},

il che equivale a dire che le prime $r$ cifre decimali non sono affette da errore.

Si può quindi concludere che se due approssimazioni di un integrale, di cui la seconda ottenuta dimezzando il passo di integrazione utilizzato per calcolare la prima, coincidono per le prime $r$ cifre decimali, tali cifre si possono ritenere esatte.

Più in generale se si vuole conoscere un'approssimazione di un integrale con la garanzia dell'esattezza per un determinato numero $s$ di cifre decimali, si deve calcolare un certo numero di approssimazioni successive, dimezzando di volta in volta il passo, fino a ottenerne due che coincidono per $s$ cifre.

Osserviamo che può accadere di arrivare alla prima coppia di valori approssimati soddisfacenti che presentano più di s cifre coincidenti. Osserviamo anche che procedendo con la riduzione dell'ampiezza dei sottointervalli, oltre al maggior tempo di calcolo richiesto, si possono avere errori di arrotondamento tutt'altro che trascurabili a causa dell'aumento del numero di operazioni richieste; si può anche arrivare a situazioni che vedono aumentare l'errore complessivo con il ridursi del passo di integrazione. Per tale motivo sono stati concepiti dei metodi di integrazione numerica "adattivi" (o adattativi) che aumentano il numero dei sottointervalli solo nelle zone indicate da un apposito test di errore.

Nota storica[modifica | modifica wikitesto]

La regola è uno dei migliori esempi della legge dell'eponimia di Stigler: pare che questa regola fosse nota già a Torricelli, mentre per certo Cavalieri aveva già dimostrato nel 1635 la sua formulazione geometrica. Anche Keplero l'aveva trattata: infatti molti testi tedeschi la chiamano Keplersche Fassregel 100 anni prima di Simpson. Ciononostante a livello internazionale oggi è riportata secondo la designazione anglosassone come regola di Simpson: anche in questo ambiente era già stata utilizzata per esempio da Gregory. D'altra parte, al contrario dei predecessori, Thomas Simpson oggi è conosciuto essenzialmente per questa regola nonostante ebbe solo il merito di formalizzare un metodo già ampiamente noto.