La zero-stabilità è una condizione di stabilità per i metodi multistep applicati a equazioni differenziali ordinarie. Controlla cosa accade agli errori introdotti nei valori iniziali o nei passi precedenti quando il passo di integrazione tende a zero. È una proprietà strutturale del metodo numerico: senza zero-stabilità, anche un metodo formalmente consistente può fallire nella convergenza.
Il termine “zero” si riferisce al limite h\to 0, cioè al comportamento del metodo quando il passo temporale diventa piccolo e l’equazione differenziale viene vista nella sua parte omogenea dominante.
Metodi multistep
Un metodo lineare a k passi per un problema di Cauchy può essere scritto nella forma
dove h è il passo, y_n approssima la soluzione nel nodo t_n, e i coefficienti \alpha_j,\beta_j definiscono il metodo.
La zero-stabilità riguarda il polinomio caratteristico associato alla parte sinistra:
detto primo polinomio caratteristico.
Condizione delle radici
Il metodo è zero-stabile se tutte le radici \xi_i di \rho soddisfano
e ogni radice con modulo uguale a 1 è semplice. In forma sintetica:
Questa è detta condizione delle radici. Se esiste una radice con modulo maggiore di 1, un piccolo errore può crescere geometricamente. Se una radice sul cerchio unitario è multipla, l’errore può crescere in modo polinomiale e diventare comunque incompatibile con la convergenza uniforme.
Perché è necessaria
Nei metodi a un passo, come il metodo di Eulero, ogni nuovo valore dipende essenzialmente dal valore precedente. Nei metodi multistep, invece, la formula ricicla più valori passati. Questo rende il metodo più efficiente, ma introduce modi numerici parassiti: componenti dell’errore che non rappresentano la soluzione dell’equazione differenziale, ma la dinamica interna dello schema.
La zero-stabilità garantisce che questi modi non vengano amplificati in modo incontrollato quando si raffina la griglia.
Consistenza, stabilità e convergenza
Per i metodi multistep lineari vale il principio di equivalenza di Dahlquist: sotto ipotesi standard, un metodo consistente è convergente se e solo se è zero-stabile. In forma operativa:
La consistenza assicura che l’errore locale tenda a zero; la zero-stabilità assicura che gli errori locali e iniziali non vengano amplificati troppo nel processo di propagazione.
Esempio di lettura della condizione
Se il polinomio caratteristico ha una radice \xi=1, questa è attesa nei metodi consistenti: rappresenta la propagazione della soluzione costante. Il problema non è la presenza di \xi=1, ma la presenza di radici fuori dal disco unitario o di radici multiple sul bordo.
Per esempio, una radice
produce contributi proporzionali a 1{,}2^n, che crescono rapidamente al procedere dei passi. Una radice doppia in \xi=1 può produrre contributi proporzionali a n, incompatibili con una buona propagazione degli errori.
Errori comuni
Un errore frequente è confondere zero-stabilità e stabilità assoluta. La zero-stabilità riguarda il limite h\to 0 e il polinomio \rho; la stabilità assoluta riguarda invece il comportamento del metodo applicato a problemi test, spesso del tipo y'=\lambda y, per valori finiti di h\lambda.
Un secondo errore è pensare che un ordine elevato garantisca convergenza. L’ordine misura la consistenza locale; senza zero-stabilità, l’errore globale può comunque esplodere.
Infine, la condizione va verificata sul metodo, non sul singolo problema differenziale. È una proprietà dello schema multistep e dei suoi coefficienti. Per una visione d’insieme, vedi il formulario di EDO e analisi numerica e la voce sui metodi numerici per EDO.