Metodo delle tangenti

In matematica, e in particolare in analisi numerica, il metodo delle tangenti, chiamato anche metodo di Newton o metodo di Newton-Raphson, è uno dei metodi per il calcolo approssimato di una soluzione di un'equazione della forma ${\displaystyle f(x)=0}$. Esso si applica dopo avere determinato un intervallo ${\displaystyle [a,b]}$ che contiene una sola radice.

Il metodo consiste nel sostituire alla curva ${\displaystyle y=f(x)}$ la tangente alla curva stessa, partendo da un qualsiasi punto; per semplicità si può iniziare da uno dei due punti che hanno come ascissa gli estremi dell'intervallo ${\displaystyle [a,b]}$ e assumere, come valore approssimato della radice, l'ascissa ${\displaystyle x_{t}}$ del punto in cui la tangente interseca l'asse delle ${\displaystyle x}$ internamente all'intervallo ${\displaystyle [a,b]}$.

Procedendo in modo iterativo si dimostra che la relazione di ricorrenza del metodo è

${\displaystyle x_{n+1}=x_{n}-{\frac {f(x_{n})}{f'(x_{n})}},}$

che permette di determinare successive approssimazioni della radice dell'equazione ${\displaystyle y=f(x)=0}$. Con le ipotesi poste, si dimostra che la successione delle ${\displaystyle x_{n}}$ converge alla radice piuttosto rapidamente.

Più in dettaglio, si dimostra che se ${\displaystyle f\in C^{2}(I)}$ dove ${\displaystyle I}$ è un opportuno intorno dello zero ${\displaystyle \alpha }$ con ${\displaystyle f'(\alpha )\neq 0}$ e se ${\displaystyle x_{0}\in I,}$ allora

${\displaystyle \lim _{n\to +\infty }{\frac {\alpha -x_{n+1}}{(\alpha -x_{n})^{2}}}=-{\frac {f''(\alpha )}{2f'(\alpha )}},}$

cioè la convergenza è quadratica (il numero di cifre significative approssimativamente raddoppia ad ogni iterazione; mentre col metodo di bisezione cresce linearmente), benché locale (cioè non vale per ogni ${\displaystyle I}$). Se invece la radice è multipla, cioè ${\displaystyle f'(\alpha )=0}$ allora la convergenza è lineare (più lenta). Nella pratica, fissata la tolleranza di approssimazione consentita ${\displaystyle \tau }$, il procedimento iterativo si fa terminare quando ${\displaystyle \left|x_{n+1}-x_{n}\right|<\tau \cdot |x_{n+1}|.}$

Il problema di questo metodo è che la convergenza non è garantita, in particolare quando ${\displaystyle f'(x)}$ varia notevolmente in prossimità dello zero. Inoltre, il metodo assume che ${\displaystyle f'(x)}$ sia disponibile direttamente per un dato ${\displaystyle x}$. Nei casi in cui questo non si verifichi e risultasse necessario calcolare la derivata attraverso una differenza finita, è consigliabile usare il metodo della secante.

Il matematico francese François Viète presentò nel 1600^[1] un metodo, già noto nel 1427 da al-Kashi, per la ricerca degli zeri di un polinomio attraverso una perturbazione di una sua soluzione approssimata. Quattro anni dopo Newton venne a conoscenza del metodo di Viète e nel 1669 scoprí autonomamente un metodo per la ricerca degli zeri di un polinomio.

Come esempio mostra la seguente equazione ${\displaystyle f(x)=x^{3}-2x-5=0}$ una cui soluzione ha parte intera ${\displaystyle x_{0}=2}$. Applicando la sostituzione ${\displaystyle x=2+p}$ si ricava il polinomio ${\displaystyle p^{3}+6p^{2}+10p-1=0}$e trascurando i monomi di grado superiore al primo, ossia linearizzando il polinomio, si ottiene ${\displaystyle p=0{,}1}$. Per cui si applica la sostituzione ${\displaystyle p=0{,}1+q}$ e si arriva a ${\displaystyle q^{3}+6{,}3q^{2}+11{,}23q+0{,}061=0}$ e per linearizzazione ${\displaystyle q=-0{,}0054}$. Sostituendo ${\displaystyle q=0{,}0054+r}$ e facendo lo stesso ragionamento si ricava ${\displaystyle r=0{,}00004853}$. Da cui ${\displaystyle x_{3}=x_{0}+p+q+r=2{,}09455147.}$

Si possono fare due osservazioni relative al metodo proposto:

1. ${\displaystyle p=x_{1}-x_{0}=-{\frac {f(x_{0})}{f'(x_{0})}}}$ e ${\displaystyle q=x_{2}-x_{1}=-{\frac {f(x_{1})}{f'(x_{1})}}}$ per cui il metodo trovato da Newton corrisponde al moderno metodo delle tangenti;
2. osservando i valori di ${\displaystyle p}$, ${\displaystyle q}$ e ${\displaystyle r}$ si può notare che il numero di zeri dopo la virgola raddoppia ad ogni passo, allora nell'esempio si ha convergenza quadratica.

Nel 1687, nel Philosophiae Naturalis Principia Mathematica Newton applica per la prima volta il metodo ad un'equazione non polinomiale. È il caso dell'equazione ${\displaystyle x-e\sin(x)=M}$ dove ${\displaystyle M}$ indica l'anomalia media e ${\displaystyle x}$ l'anomalia eccentrica. In questo caso approssimando il seno come somma troncata del suo sviluppo in serie di Taylor Newton ricavava un polinomio e quindi poteva applicare il metodo da lui trovato.

Nel 1690 il matematico Joseph Raphson riuscì a ricavare un metodo iterativo per aggiornare la soluzione approssimata ${\displaystyle x_{k}}$ senza dover calcolare la potenza del monomio completa e nel 1740 Thomas Simpson, nel libro Essays on Several Curious and Useful Subjects in Speculative and Mix's Mathematicks, Illustrated by a variey of Examples ricavò il moderno metodo delle tangenti riconoscendo il ruolo delle derivate prime nell'aggiornamento della soluzione.

Consideriamo una funzione undimensionale ${\displaystyle f\in C^{2}[a,b]}$, e quindi per il teorema di Weierstrass la funzione ammette un minimo ${\displaystyle x}$ da determinare.

Per cui, preso un punto ${\displaystyle x_{0}}$ nell'intervallo, sfruttando la serie di Taylor di ${\displaystyle f}$, si trova che

${\displaystyle f(x)=f(x_{0})+(x-x_{0})f'(x_{0})+{\frac {(x-x_{0})^{2}}{2}}f''(r(x))}$

con ${\displaystyle r(x)}$ compreso tra ${\displaystyle x}$ e ${\displaystyle x_{0}.}$

Per cui, se ${\displaystyle |x-x_{0}|}$ è sufficientemente piccolo, ossia ${\displaystyle (x-x_{0})^{2}\approx 0}$, ponendo ${\displaystyle f(x)=0}$ per trovare l'intersezione della retta tangente nel punto ${\displaystyle (x_{0},f(x_{0}))}$ con l'asse delle ${\displaystyle x}$, si ricava ${\displaystyle x\approx x_{0}-{\frac {f(x_{0})}{f'(x_{0})}}=x_{1}}$. Osservare che l'ultima relazione ha senso solo se ${\displaystyle f'(x_{0})}$ non è nullo. Una volta trovato ${\displaystyle x_{1}}$ si reitera il procedimento. Si è trovato così il seguente algoritmo:

Metodo di Newton Unidimensionale
* Passo 0: Si sceglie un punto ${\displaystyle x_{0}}$ nell'intervallo ${\displaystyle [a,b]}$. Si pone ${\displaystyle k=0.}$
 Per ${\displaystyle k=0,1,\ldots }$ 
* Passo 1: Si determina ${\displaystyle x_{k+1}=x_{k}-{\frac {f(x_{k})}{f'(x_{k})}}}$. 
* Passo 2: Poni ${\displaystyle k=k+1}$ e torna al Passo 1.

Come si è visto, condizione necessaria affinché il metodo sia applicabile è che esista un intervallo ${\displaystyle I=[a,b]}$ in cui ${\displaystyle x_{0}\in I}$ è tale che ${\displaystyle f'(x_{0})\neq 0}$ e ${\displaystyle f''(x_{0})\neq 0}$. Per il teorema degli zeri tale intervallo esiste se e solo se ${\displaystyle f'(x_{0})\neq 0}$ e ${\displaystyle f''(x_{0})\neq 0}$.

Consideriamo una funzione ${\displaystyle f\in C^{2}(\mathbb {R} ^{n})}$ e sia ${\displaystyle x\in \mathbb {R} ^{n}}$ lo zero da determinare. Sfruttando lo sviluppo in serie di Taylor si ha che, preso un generico vettore ${\displaystyle h\in \mathbb {R} ^{n}}$:

${\displaystyle f(x+h)=f(x)+J(x)h+O(\|h\|^{2})}$

dove ${\displaystyle J(x)}$ indica la matrice jacobiana di ${\displaystyle f}$ calcolata nel punto ${\displaystyle x}$.

Per cui, se ${\displaystyle J(x)}$ non è singolare si ottiene un nuovo punto ${\displaystyle y=x+h_{1}}$ dove ${\displaystyle h_{1}}$ è la soluzione del sistema lineare ${\displaystyle J(x)h=-f(x)}$.

Si è trovato così il seguente algoritmo:

Metodo di Newton multidimensionale
* Passo 0: Si sceglie un punto ${\displaystyle x_{0}}$ . Si pone ${\displaystyle k=0}$.
 Per ${\displaystyle k=0,1,...}$ 
* Passo 1: Si risolve ${\displaystyle J(x_{k})h=-f(x_{k})}$ ottenendo il vettore ${\displaystyle h}$. 
* Passo 2: Si determina ${\displaystyle x_{k+1}=x_{k}+h}$.
* Passo 3: Si pone ${\displaystyle k=k+1}$ e si torna al Passo 1.

${\displaystyle f\in C^{2}(I)}$ dove ${\displaystyle I}$ è un opportuno intorno della radice ${\displaystyle \alpha }$ con ${\displaystyle f'(\alpha )\neq 0}$ e se ${\displaystyle x_{0}\in I,}$

• Calcolo di uno zero di una funzione
• Metodo di bisezione
• Metodo delle secanti
• Confronto tra metodo delle secanti e metodo delle tangenti
• Algoritmo di Gauss-Newton

• Wikiversità contiene risorse sul metodo delle tangenti
• Wikimedia Commons contiene immagini o altri file sul metodo delle tangenti

• (EN) Newton’s iterative method, su Enciclopedia Britannica, Encyclopædia Britannica, Inc.
• (EN) Eric W. Weisstein, Newton's Method, su MathWorld, Wolfram Research.
• Presentazione: usare il metodo delle tangenti con Turbo Pascal, da matsoftware.it
• Applet Geogebra per applicare il metodo

Portale Matematica: accedi alle voci di Teknopedia che trattano di matematica