Create successful ePaper yourself
Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.
32 Chapitre 2 : Analyse numérique dans un cas simplié<br />
Combining these relations gives<br />
F (y k ) ≤ F (x k ). (2.39)<br />
A similar analysis for the reverse mode also gives F (y k ) ≤ F (x k ).<br />
The components of z k (s) lie on the unit sphere for all choices of s. Consequently,<br />
F ′′<br />
k<br />
(s) is bounded by a nite constant M, uniformly in k and s. Dene the constants<br />
δ = min{ρ, 1/M}<br />
and ¯s k = −δF ′ k(0).<br />
Since ¯s k lies on the interval [0, −ρF k ′ (0)] appearing in (2.6), we have<br />
F (x k+1 ) = F k (s k ) ≤ F k (¯s k ). (2.40)<br />
Expanding in a Taylor series around s = 0, there exists ξ k ∈ [0, s k ] such that<br />
F k (¯s k ) = F k (0) + F k(0)¯s ′ k + 1 2 ¯s2 kF k ′′ (ξ k )<br />
≤ F k (0) + F k(0)¯s ′ k + 1 2 ¯s2 kM<br />
= F k (0) + δF ′ k(0) 2 ( 1 2 δM − 1) ≤ F k(0) − δ 2 F ′ k(0) 2<br />
= F (y k ) − δ 2 F ′ k(0) 2 ≤ F (x k ) − δ 2 F ′ k(0) 2 ,<br />
where the last inequality is (2.39). Combining this with (2.40) gives<br />
F (x k+1 ) ≤ F (x k ) −<br />
Summing this inequality over k yields<br />
F (x k ) ≤ F (x 0 ) −<br />
( δ<br />
2)<br />
F ′ k(0) 2 .<br />
( δ ∑k−1<br />
F i<br />
2) ′ (0) 2 .<br />
i=0<br />
Since the feasible points for (2.1) lie on the unit sphere, the objective function value<br />
is bounded from below. Hence, we have<br />
lim F k(0) ′ = 0. (2.41)<br />
k→∞<br />
By the denition of z k and the fact that s k ∈ [0, −ρF k ′(0)] where F k ′ (0) approaches<br />
0, we also conclude that<br />
‖x k+1 − y k ‖ ≤ c|F k(0)| ′ (2.42)<br />
where c is a constant which is independent of k.<br />
Step 2. The change in the multiplier µ.<br />
The rst-order optimality conditions associated with the subproblems (2.2) can<br />
be stated in the following way : There exist scalars λ k1 , µ k1 , λ k2 , and µ k2 such that<br />
forward [ H 1 0<br />
reverse<br />
] [ ] [ ]<br />
yk1 λk1 y<br />
= k1 + µ k1 x k2<br />
,<br />
0 H 2 y k2 µ k2 y k1 + λ k2 y k2<br />
[ ] [ ]<br />
H1 0 yk1<br />
=<br />
0 H 2 y k2<br />
[ ]<br />
λk1 y k1 + µ k1 y k2<br />
.<br />
µ k2 x k1 + λ k2 y k2<br />
(2.43)