Méthode de Halley

De testwiki
Aller à la navigation Aller à la recherche

Modèle:Redirect

En analyse numérique, la méthode de Halley est un algorithme de recherche d'un zéro d'une fonction utilisé pour les fonctions d'une variable réelle dérivables deux fois et à dérivée seconde continue (i.e. C2). La méthode, présentée par l'astronome Edmond Halley en 1694[1], est une généralisation de la méthode de Newton, à convergence cubique.

Énoncé

Soit f une fonction de classe C² et a un zéro de f. La méthode de Halley consiste à itérer

xn+1=xn2f(xn)f(xn)2(f(xn))2f(xn)f(xn)

à partir d'une valeur x0 proche de a.

Au voisinage de a, la suite vérifie :

|xn+1a|<K|xna|3,

avec K>0 ; ce qui signifie que la convergence est donc (au pire) cubique.

La formulation suivante montre que la méthode de Halley est un raffinement de celle de Newton, et très proche de celle-ci si f(a)=0.

xn+1=xnf(xn)f(xn)f(xn)f(xn)f(xn)2=xnf(xn)f(xn)(1f(xn)f(xn)f(xn)2f(xn))1.

On remarque que dans cette formulation l'expression f(xn)/f(xn) n'est calculée qu'une fois.

Obtention à partir de la méthode de Newton

La formule se déduit de la méthode de Newton appliquée à la fonction g=f/f ; en effet, celle-ci s'écrit  :

xn+1=xng(xn)g(xn),

avec

g(x)=2[f(x)]2f(x)f(x)2f(x)f(x),

d'où le résultat. Si f(a)=0, ceci ne s'applique que si g peut être prolongée par continuité en a.

Application au calcul d'une fonction réciproque

Appliquée à la fonction F(x)=f(x)y, cette méthode conduit à une suite récurrente convergeant vers f1(y).

Ceci peut être utilisé pour le calcul d'une racine n-ième, avec F(x)=xnA ; la suite récurrente définie par xk+1=xk((n1)(xk)n+(n+1)A)(n+1)(xk)n+(n1)A converge vers An.

La méthode de Halley est aussi utilisée pour calculer un logarithme népérien, en inversant la fonction exponentielle.

Notes et références

Modèle:Références

Voir aussi

Liens internes

Liens externes

Modèle:Palette

Modèle:Portail