Calcolo Numerico
Laboratorio 5:
Metodo di Newton Modificato
Varianti del metodo di Newton
Ángeles Martínez Calomardo
[Link]
[Link]@[Link]
Laurea Triennale Ing. Industriale
Meccanica (matricole dispari) e Aerospaziale
A.A. 2015–2016
Calcolo Numerico (Ingegneria Industriale) (Laurea
Introduzione
Triennale
a Matlab/Octave
Ing. Industriale Meccanica
Ángeles
(matricole
Martínez
dispari)
Calomardo
e Aerospaziale
1 /A.A.
8 2015–2
Convergenza del Metodo di Newton-Raphson
La convergenza del metodo di Newton è innanzittutto locale cioè è garantita
solo se x0 è sufficientemente vicino alla soluzione.
Inoltre, nel caso in cui ξ sia uno zero semplice (cioè f 0 (ξ ) 6= 0) è quadratica
ovvero:
|xk+1 − ξ | ' M |xk − ξ |2
Questa proprietà rende il metodo di Newton assai veloce e, a dispetto della
sua età (343 anni), il più usato per equazioni scalari.
Se invece la molteplicità dello zero è maggiore di uno, il metodo di
Newton converge linearmente.
Detta r la molteplicità di ξ , la convergenza quadratica può essere ripristinata
utilizzando il cosiddetto metodo di Newton modificato:
x0 fissato
f (xk )
xk+1 = xk − r 0 , k≥0
f (xk )
Calcolo Numerico (Ingegneria Industriale) (Laurea
Introduzione
Triennale
a Matlab/Octave
Ing. Industriale Meccanica
Ángeles
(matricole
Martínez
dispari)
Calomardo
e Aerospaziale
2 /A.A.
8 2015–2
Metodo di Newton-Raphson per radici multiple
Esercizio
Si scriva la function newtonmod.m ottenuta a partire dalla function newton.m,
aggiungendo come ultimo parametro di ingresso, oltre a quelli della function
newton, la molteplicità della radice, r, e facendo in modo che implementi il
metodo di Newton modificato per la ricerca di radici multiple.
Calcolo Numerico (Ingegneria Industriale) (Laurea
Introduzione
Triennale
a Matlab/Octave
Ing. Industriale Meccanica
Ángeles
(matricole
Martínez
dispari)
Calomardo
e Aerospaziale
3 /A.A.
8 2015–2
Metodo di Newton-Raphson per radici multiple
Esercizio
Si scriva lo script chiamato scriptmod.m che:
1 usi la function newtonmod.m per approssimare lo zero della funzione
(x − 1)2 log(x) nell’intervallo [1, 3], per i valori r = 1 e r = 3, con tolleranza
tol= 10−9 .
2 confronti la convergenza del Metodo di Newton nei due casi mediante un
grafico semilogaritmico che rappresenti gli scarti ottenuti in entrambi i casi in
funzione del numero di iterazioni. (Si noti che il grafico deve essere unico
con due curve, una per ogni valore di r).
Calcolo Numerico (Ingegneria Industriale) (Laurea
Introduzione
Triennale
a Matlab/Octave
Ing. Industriale Meccanica
Ángeles
(matricole
Martínez
dispari)
Calomardo
e Aerospaziale
4 /A.A.
8 2015–2
Metodo di Newton-Raphson per radici multiple
Un esempio del grafico richiesto è il seguente:
Profilo di convergenza Newton Modificato versus Newton
100
Newton Modificato
Newton
10-2
10-4
10-6
abs(scarto)
10-8
10-10
10-12
10-14
0 10 20 30 40 50
iter
Calcolo Numerico (Ingegneria Industriale) (Laurea
Introduzione
Triennale
a Matlab/Octave
Ing. Industriale Meccanica
Ángeles
(matricole
Martínez
dispari)
Calomardo
e Aerospaziale
5 /A.A.
8 2015–2
Metodo della tangente fissa
Si copi il file contenente la function newton.m con il nome tfissa.m e si
modifichi opportunamente in modo che implementi il metodo della tangente fissa:
f (xk )
xk+1 = xk − 0
, k = 0, 1, . . . , x0 fissato
f (x0 )
Tale function deve avere come dati di ingresso:
la funzione f la sua derivata prima d f , x0 , tol e itmax
e restituire come parametri di output:
il numero di iterazioni impiegate iter, il vettore con le iterate
xk , k = 0, . . . , iter e il vettore degli scarti sk = xk − xk−1 , k = 1, . . . , iter.
Come test di arresto si usi quello basato sullo scarto.
Calcolo Numerico (Ingegneria Industriale) (Laurea
Introduzione
Triennale
a Matlab/Octave
Ing. Industriale Meccanica
Ángeles
(matricole
Martínez
dispari)
Calomardo
e Aerospaziale
6 /A.A.
8 2015–2
Metodo della tangente fissa
Esercizio
Si scriva uno script chiamato scripttfissa.m che chiamando le due function
(newton.m e tfissa.m) risolva con i due metodi l’equazione f (x) = 0 dove
f (x) = e−x + cos(x) − 3,
e con:
x0 = −1
itmax = 100
tol = 10−9 .
Creare un (UNICO) grafico semilogaritmico con il profilo di convergenza dei due
metodi.
Il grafico può essere salvato tramite il comando print nel file [Link].
p r i n t −d p d f g r a f i c o t f . p d f
Calcolo Numerico (Ingegneria Industriale) (Laurea
Introduzione
Triennale
a Matlab/Octave
Ing. Industriale Meccanica
Ángeles
(matricole
Martínez
dispari)
Calomardo
e Aerospaziale
7 /A.A.
8 2015–2
Confronto tra Newton-Raphson e tangente fissa
Un esempio del grafico richiesto è il seguente:
Profilo di convergenza dei due metodi (f(x)= exp(-x) + cos(x) - 3)
10 0
Newton
-1 Tfissa
10
10 -2
10 -3
10 -4
scarto
10 -5
10 -6
10 -7
10 -8
10 -9
10 -10
1 2 3 4 5 6 7 8 9 10 11
iterazioni
Si osserva che il metodo della tangente fissa pur convergendo linearmente è molto
veloce: quanto vale la costante asintotica?
Calcolo Numerico (Ingegneria Industriale) (Laurea
Introduzione
Triennale
a Matlab/Octave
Ing. Industriale Meccanica
Ángeles
(matricole
Martínez
dispari)
Calomardo
e Aerospaziale
8 /A.A.
8 2015–2