06.01.2015 Views

The Steepest Descent Algorithm for Unconstrained Optimization and ...

The Steepest Descent Algorithm for Unconstrained Optimization and ...

The Steepest Descent Algorithm for Unconstrained Optimization and ...

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

It is elementary to show that<br />

Proposition 6.1 h ′ (0) < 0.<br />

q.e.d.<br />

′<br />

h (α) = ∇f (¯ x + αd) ¯ T d. ¯<br />

Because h(α) is a convex function of α, wealsohave:<br />

Proposition 6.2 h ′ (α) is a monotone increasing function α. of<br />

q.e.d.<br />

Figure 5 shows an example of convex function of two variables to be<br />

optimized. Figure 6 shows the function h(α) obtained by restricting the<br />

function of Figure 5 to the line shown in that figure. Note from Figure 6 that<br />

h(α) is convex. <strong>The</strong>re<strong>for</strong>e its first derivative h ′ (α) will be a monotonically<br />

increasing function. This is shown in Figure 7.<br />

Because h ′ (α) is a monotonically increasing function, we can approxiα,<br />

¯ the point that satisfies h ′ (¯ α) = 0, by a suitable bisection<br />

mately compute<br />

method. Suppose that we know a value α ˆ that h ′ (ˆ α) > 0. Since h ′ (0) < 0<br />

<strong>and</strong> h ′ (ˆ α) > 0, the mid-value α ˜ = 0+ 2α ˆ is a suitable test-point. Note the<br />

following:<br />

• If h ′ (˜α) = 0, we are done.<br />

• If h ′ (˜ α) > 0, we can now bracket α ¯ in the interval (0, α). ˜<br />

• If h ′ (˜ α) < 0, we can now bracket α ¯ in the interval ( α, ˜ α). ˆ<br />

This leads to the following bisection algorithm <strong>for</strong> minimizing h(α) = f (¯x +<br />

¯ αd) by solving the equation h ′ (α) ≈ 0.<br />

Step 0. Set k = 0. Set α l := 0 <strong>and</strong> α u := ˆα.<br />

α = αu+α l<br />

Step k. Set ˜ <strong>and</strong> compute h ′ 2<br />

(˜ α).<br />

18

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

Saved successfully!

Ooh no, something went wrong!