Sujet Et Corrigé Optimisation Sans Contraintes
Sujet Et Corrigé Optimisation Sans Contraintes
Sujet Et Corrigé Optimisation Sans Contraintes
Département de mathématiques
L3 mathématiques 2021/2022
Exercice 1 (6 points):
On définit la fonction J : Rn → R par
1
J(x) := ⟨Ax, x⟩ − ⟨b, x⟩,
2
où A ∈ Mn (R) est une matrice symétrique définie positive et b ∈ Rn . On
considère la méthode du gradient à pas optimal pour la minimisation de J.
- Donner l’expression de xk+1 obtenu par cette méthode.
- Montrer que le pas θk s’écrit dans ce cas:
||dk ||2
θk = .
⟨Adk , dk ⟩
- Rappelons que le problème de minimisation de J admet une unique solution
x̄ ∈ Rn caractérisée comme étant l’unique solution de l’équation d’Euler
∇J(x) = 0. Montrer que
⟨A−1 dk , dk ⟩ = 2(J(xk ) − J(x̄)).
Solution:
- L’algorithme de la méthode du gradient à pas optimal s’écrit:
xk+1 = xk + θk dk
dk = −Axk + b
θk est tel que J(xk + θk dk ) = minθ∈R J(xk + θdk )
(1 point).
- Dans l’algorithme du gradient à pas optimal, deux directions successives
sont orthogonales:
⟨dk+1 , dk ⟩ = 0.
En utilisant la formule de xk+1 et de dk on obtient:
0 = ⟨−Axk+1 + b, dk ⟩
= ⟨−A(xk + θk dk ) + b, dk ⟩
= ⟨−Axk + b − θk Adk , dk ⟩
= ⟨dk − θk Adk , dk ⟩
= ||dk ||2 − θk ⟨Adk , dk ⟩
on suppose que dk ̸= 0 et donc ⟨Adk , dk ⟩ =
̸ 0 (car A est définie positive).
Ainsi, θk s’exprime
||dk ||2
θk = .
⟨Adk , dk ⟩
(2 points)
- La solution du problème de minimisation de J est x̄ = A−1 b. Un simple
calcul montre que
∇f (x) = 0 ⇔ Ax = 0.
2
converge pour toute valeur de α ∈]0, 23 [. (1.5 point)
- Le pas optimal de cette méthode est donné par:
2 1
ᾱ = = (1.5 point)
λ1 + λ2 2
Exercice 3 (9 points):
Déterminer et préciser la nature des points critiques des fonctions suivantes:
f (x, y) = x3 + y 3 − 3axy, a ∈ R
g(x, y) = x2 − cos(y)
h(x, y) = y 2 + xy ln(x).
Solution:
1) le gradient le hessien de f sont donnés par:
! !
3x2 − 3ay 2 6x −3a
∇f (x, y) := ∇ f (x, y) =
3y 2 − 3ax −3a 6y
donc pk est un minimum si k est pair et c’est un point selle si k est impair
(3 points).
3) la fonction h est définie pour tout x > 0 et elle est de classe C 2 sur son
domaine de définition. le gradient le hessien de h sont donnés par:
y
! !
y ln(x) + y 2 ln(x) + 1
∇h(x, y) := ∇ h(x, y) = x
2y + x ln(x) ln(x) + 1 2