Formula di differenziazione all'indietro

La formula di differenziazione all'indietro (BDF, dall'inglese backward differentiation formula) è una famiglia di metodi impliciti per l'integrazione numerica di equazioni differenziali ordinarie. Sono metodi lineari multipasso che, dati una funzione e un istante temporale, forniscono un valore approssimato della derivata della funzione utilizzando i risultati dei calcoli per istanti temporali precedenti e aumentando così l'accuratezza dell'approssimazione. Questi metodi sono usati specialmente per la soluzione di equazioni differenziali rigide.

Formula generale[modifica | modifica wikitesto]

Le BDF sono usate per risolvere il seguente problema ai valori iniziali:

$y'=f(y,t),\quad y(t_{0})=y_{0}$

Una generica BDF può essere scritta come:

$\sum _{k=0}^{s}a_{k}y_{n+k}=h\beta f(y_{n+s},t_{n+s})$

dove $h$ indica il passo di integrazione e $t_{n}=t_{0}+nh$ . I coefficienti $a_{k}$ e $\beta$ sono scelti in modo tale che il metodo raggiunga l'ordine $s$ , che è il massimo possibile.

I metodi BDF sono impliciti, perciò richiedono la risoluzione di equazioni non lineari ad ogni passo di integrazione. Di solito, per questo scopo si usa un metodo di Newton modificato.

Formule specifiche[modifica | modifica wikitesto]

Le BDF per s da 1 a 6 sono:

BDF1: $y_{n+1}-y_{n}=hf(y_{n+1},t_{n+1})$ (che equivale al metodo di Eulero all'indietro);
BDF2: $y_{n+2}-{\frac {4}{3}}y_{n+1}+{\frac {1}{3}}y_{n}={\frac {2}{3}}hf(y_{n+2},t_{n+2})$ ;
BDF3: $y_{n+3}-{\frac {18}{11}}y_{n+2}+{\frac {9}{11}}y_{n+1}-{\frac {2}{11}}y_{n}={\frac {6}{11}}hf(y_{n+3},t_{n+3})$ ;
BDF4: $y_{n+4}-{\frac {48}{25}}y_{n+3}+{\frac {36}{25}}y_{n+2}-{\frac {16}{25}}y_{n+1}+{\frac {3}{25}}y_{n}={\frac {12}{25}}hf(y_{n+4},t_{n+4})$ ;
BDF5: $y_{n+5}-{\frac {300}{137}}y_{n+4}+{\frac {300}{137}}y_{n+3}-{\frac {200}{137}}y_{n+2}+{\frac {75}{137}}y_{n+1}-{\frac {12}{137}}y_{n}={\frac {60}{137}}hf(y_{n+5},t_{n+5})$ ;
BDF6: $y_{n+6}-{\frac {360}{147}}y_{n+5}+{\frac {450}{147}}y_{n+4}-{\frac {400}{147}}y_{n+3}+{\frac {225}{147}}y_{n+2}-{\frac {72}{147}}y_{n+1}+{\frac {10}{147}}y_{n}={\frac {60}{147}}hf(y_{n+6},t_{n+6})$ .

Metodi con s ≥ 7 non sono utilizzati perché non sono zero-stabili.

Stabilità[modifica | modifica wikitesto]

La regione di stabilità del metodo BDF1 comprende l'intera parte sinistra del piano complesso, quindi questo metodo è A-stabile e può essere utilizzato per la risoluzione di equazioni rigide. Le regioni di stabilità delle BDF di ordine superiore invece contengono la quasi totalità del semipiano sinistro, inclusa la parte negativa dell'asse reale.