12.07.2015 Views

Esercizi svolti di esame ed esonero - Dipartimento di Matematica

Esercizi svolti di esame ed esonero - Dipartimento di Matematica

Esercizi svolti di esame ed esonero - Dipartimento di Matematica

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

Nella fattorizzazione <strong>di</strong> Cholesky, invece,da cuiH −1 =( )1 05 10, (H t ) −1 = (H −1 ) t ,K ∞ (H) = 15 , K ∞ (H t ) = 15.Il numero globale <strong>di</strong> con<strong>di</strong>zionamento associato alla soluzione dei due sistemi linearié quin<strong>di</strong> <strong>di</strong> 3780 per la fattorizzazione LU e <strong>di</strong> 225 (ben piú basso, quin<strong>di</strong>) per lafattorizzazione <strong>di</strong> Cholesky. Per confronto, il numero <strong>di</strong> con<strong>di</strong>zionamento della matriceA vale K ∞ (A) = 169.<strong>Esercizi</strong>o 3.c) Anche se non é strettamente necessario per la convergenza del metodo, conviene porsiin un intervallo in cui la funzione cambi <strong>di</strong> segno. La scelta piú semplice é quelladell’intervallo [1, 2]. Calcoliamo ora la costante <strong>di</strong> contrazione del metodo. Poichéa = 1 e b = 2, il metodo si scrive x k+1 = g(x k ) cone derivando g si hag(x) = x −b − af(b) − f(a) f(x) = x − 1 3 (x2 − 2),g ′ (x) = 1 − 2x 3che essendo una funzione lineare, assume il suo massimo modulo in un estremo. Quin<strong>di</strong>L = sup |g ′ (x)| = max(|g ′ (1)|, |g ′ (2)|) = 1[1,2]3(osserviamo che questo giustifica a posteriori la scelta dell’intervallo [1, 2]). Infine,poiché la ra<strong>di</strong>ce é interna all’intervallo [a, b], supponendo che anche x 0 lo sia, si ottiene|x k − √ 2| ≤ L k |x 0 − √ 2| ≤ L k (b − a) = L k ,e la precisione richiesta si raggiunge quando 3 k ≥ 1000, ovvero alla settima iterazione.59

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!