18.11.2014 Views

ϵ2 - Le Cermics - ENPC

ϵ2 - Le Cermics - ENPC

ϵ2 - Le Cermics - ENPC

SHOW MORE
SHOW LESS

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)

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!