Méthode de la sécante

En analyse numérique, la méthode de la sécante est un algorithme de recherche d'un zéro d'une fonction f.

La méthode

modifier

La méthode de la sécante est une méthode comparable à celle de Newton, où l'on remplace par On obtient la relation de récurrence :

L'initialisation nécessite deux points x0 et x1, proches, si possible, de la solution recherchée. Il n'est pas nécessaire que x0 et x1 encadrent une racine de f. La méthode de la sécante peut aussi être vue comme une généralisation de la méthode de la fausse position, où les calculs sont itérés.

Démonstration

modifier
La courbe rouge représente la fonction f et le segment en bleu, la sécante.
Illustration des deux premières itérations, pour une autre courbe (ici, la méthode va diverger car x0 et x1 sont choisis trop loin de la solution).

Étant donnés a et b, on construit la droite passant par (a, f(a)) et (b, f(b)). Son équation est :

On choisit c égal à l'abscisse du point d'ordonnée y = 0 de cette droite :

Si l'on extrait c de cette équation, on retrouve la relation de récurrence citée plus haut :

avec

Convergence

modifier

La méthode peut diverger. La suite peut même ne pas être bien définie.

Cependant, si la fonction est dérivable de dérivée continue et que la suite converge, sa limite est bien une racine de f [1].

Sous l'hypothèse que la fonction f soit deux fois continûment différentiable et la solution soit une racine simple de f, on peut démontrer qu'il existe un intervalle autour de la racine tel que, si les valeurs initiales x0 et x1 sont prises dans cet intervalle, la suite est bien définie et converge vers la racine[2]. La méthode aura un ordre de convergence de

qui est le nombre d'or[3].


Aucune de ces deux conditions n'est cependant nécessaire, ni pour appliquer la méthode, ni pour en assurer la convergence. Il n'est pas nécessaire que f soit dérivable : la méthode peut s'appliquer à une fonction continue nulle part dérivable telle que la fonction de Weierstrass[réf. nécessaire].

Voir aussi

modifier

Notes et références

modifier
  1. Jean-Baptiste APOUNG KAMGA, « Résolution approchée d’équations ordinaires (EO): f (x) = 0 », sur université paris Saclay, p. 21
  2. « Équations non linéaires », sur web.univ-ubs.fr (consulté le ) diapos 10-68 et 12-79
  3. Démonstration dans Nikolaï Bakhvalov, Méthodes numériques, Moscou, Éditions Mir, , p. 402-403.

Bibliographie

modifier