Metodi Numerici per Equazioni Non Lineari

Indice dei contenuti

    I metodi numerici per equazioni non lineari permettono di trovare approssimazioni delle radici di f(x)=0f(x) = 0 quando la soluzione esatta non è disponibile in forma chiusa.

    Metodo delle Secanti

    Generalizzazione di Newton che non richiede il calcolo della derivata. Dati due punti xn1x_{n-1} e xnx_n: xn+1=xnf(xn)xnxn1f(xn)f(xn1)x_{n+1} = x_n - f(x_n)\,\frac{x_n - x_{n-1}}{f(x_n) - f(x_{n-1})}

    Converge con ordine 1.618\approx 1.618 (sezione aurea). Meno rapido di Newton, ma più robusto se ff' è difficile da calcolare.

    Metodo Regula Falsi (Falsa Posizione)

    Come il metodo delle secanti, ma si mantiene sempre un intervallo [a,b][a,b] con f(a)f(b)<0f(a)f(b) < 0: xn+1=anf(bn)bnf(an)f(bn)f(an)x_{n+1} = \frac{a_n f(b_n) - b_n f(a_n)}{f(b_n) - f(a_n)}

    Garantisce la convergenza (come la bisezione) ma converge più lentamente.

    Metodo del Punto Fisso

    Si riscrive f(x)=0f(x) = 0 come x=g(x)x = g(x) e si itera xn+1=g(xn)x_{n+1} = g(x_n).

    Teorema di Banach-Caccioppoli (contrazione): se g:[a,b][a,b]g: [a,b] \to [a,b] è una contrazione (g(x)L<1|g'(x)| \leq L < 1 per ogni xx), allora l’iterazione converge all’unico punto fisso.

    Ordine di Convergenza

    L’errore en=xnxe_n = x_n - x^* soddisfa asintoticamente: en+1Cenp|e_{n+1}| \approx C\,|e_n|^p

    Il valore pp è l’ordine di convergenza:

    • Bisezione: p=1p = 1 (convergenza lineare, lenta ma robusta)
    • Secanti: p1.618p \approx 1.618
    • Newton: p=2p = 2 (convergenza quadratica, se f(x)0f'(x^*) \neq 0)
    • Newton con radici multiple di molteplicità mm: p=1p = 1 (degrada a lineare)

    Ultimo aggiornamento: