Estrapolazione di Richardson

Da testwiki.
Vai alla navigazione Vai alla ricerca

LTemplate:'estrapolazione di Richardson è un metodo che permette la costruzione di approssimazioni numeriche di ordine superiore a partire da formule di ordine inferiore. Prende il nome da Lewis Fry Richardson, che introdusse il metodo agli inizi del XX secolo,[1][2] ed è una tecnica importante in analisi numerica, definita da alcuni autori come un metodo per "trasformare il piombo in oro"[3] e "la cui importanza difficilmente può essere sopravvalutata".[4]

Applicazioni dell'estrapolazione di Richardson includono il metodo di Romberg, che applica l'estrapolazione di Richardson alla regola del trapezio, e l'algoritmo di Bulirsch-Stoer per la soluzione di equazioni differenziali ordinarie.

Formulazione

Sia data una funzione A* e una sua approssimazione numerica A(h) con passo h e errore di troncamento di ordine O(hk0)

A*=A(h)+e

il cui residuo e=O(hk0) ha forma

e=C0hk0+C1hk1+C2hk2+

dove Ci sono costanti ignote e ki sono costanti note e tali che ki+1>kii.

Applicando l'approssimazione A con passo ht, dove t>0 è una costante arbitraria, si ottiene

A*=A(ht)+C0(ht)k0+O(hk1).

Moltiplicando entrambi i membri per tk0 e sottraendo la formula originale

A*=A(h)+C0hk0+O(hk1)

il termine di errore di ordine k0 si cancella

tk0A*A*=tk0A(ht)A(h)+tk0C0(ht)k0C0hk0+O(hk1)

ottenendo una nuova formula di ordine k1>k0

A*=tk0A(ht)A(h)tk01+O(hk1)

Il metodo può essere applicato ricorsivamente, costruendo una successione di formule

Ai+1(h)=tk0Ai(ht)Ai(h)tk01

con A0(h)=A(h), tale che A*=Ai(h)+O(hki).

Note

Bibliografia

  • Extrapolation Methods. Theory and Practice by C. Brezinski and M. Redivo Zaglia, North-Holland, 1991.
  • Ivan Dimov, Zahari Zlatev, Istvan Farago, Agnes Havasi: Richardson Extrapolation: Practical Aspects and Applications, Walter de Gruyter GmbH & Co KG, Template:ISBN (2017).

Collegamenti esterni

Template:Portale