Update Algorithme_de_newton.md

This commit is contained in:
Gergaud 2020-11-06 09:14:22 +01:00 committed by GitHub
parent f62dcf7ec3
commit 52362a78ec
No known key found for this signature in database
GPG key ID: 4AEE18F83AFDEB23

View file

@ -8,12 +8,12 @@ La fonction ``f`` étant ``C^{2}`` , on peut remplacer ``f`` au voisinage de l
On choisit alors comme point ``x_{k+1}`` le minimum de la quadratique q lorsquil existe et
est unique, ce qui nest le cas que si ``\nabla^{2} f (x)`` est définie positive. Or le minimum de q est
réalisé par ``x_{k+1}`` solution de : ``\nabla^{2} f (x_{k+1}) = 0`` , soit :
``\nabla f\left(x_{k}\right)+\nabla^{2} f\left(x_{k}\right)\left(x_{k+1}-x_{k}\right)=0``
réalisé par ``x_{k+1}`` solution de : ``\nabla q (x_{k+1}) = 0`` , soit :
``\nabla f\left(x_{k}\right)+\nabla^{2} f\left(x_{k}\right)\left(x_{k+1}-x_{k}\right)=0,``
ou encore, en supposant que ``\nabla^{2} f (x_{k})`` est définie positive :
``x_{k+1}=x_{k}-\nabla^{2} f\left(x_{k}\right)^{-1} \nabla f\left(x_{k}\right)``
``x_{k+1}=x_{k}-\nabla^{2} f\left(x_{k}\right)^{-1} \nabla f\left(x_{k}\right).``
La méthode ne doit cependant jamais être appliquée en utilisant une inversion de la
matrice Hessienne (qui peut être de très grande taille et mal conditionnée), mais plutôt en utilisant :
@ -34,14 +34,14 @@ est définie positive (par continuité de ``\nabla^{2} f``).
#### Données:
f , ``x_{0}`` première approximation de la solution cherchée, ``\epsilon > 0`` précision demandée.
``f , x_{0}`` première approximation de la solution cherchée, ``\epsilon > 0`` précision demandée.
#### Sorties
une approximation de la solution du problème ``\min _{x \in \mathbb{R}^{m}} f(x)`` .
une approximation de la solution du problème ``\min _{x \in \mathbb{R}^{n}} f(x)`` .
#### 1.Tant que le test de convergence est non satisfait
a. Calculer d k solution du système : ``\nabla^{2} f (x_{k}) d_{k} = - \nabla f (x_{k})``
a. Calculer ``d_k`` solution du système : ``\nabla^{2} f (x_{k}) d_{k} = - \nabla f (x_{k})``
b. Mise à jour : ``x_{k+1} = x_{k}+ d_{k} , k = k + 1``
#### 2.Retourner : ``x_{k}``.