Nonlinear Opt.: Gold section search 3: Unterschied zwischen den Versionen

Version vom 26. Juni 2013, 13:20 Uhr

Index of contents

1. Theory

2. Mathematical formulation

3. Example for maximum

4. Sources

Theory

Golden section search is a part of unconstrained nonlinear optimization for strictly unimodal functions. This method is used to find a minimum or maximum by narrowing values in a specific interval. The extrema exists in the range. Based on the unimodal function, we assume that the local extrema is simultaneously the global extrema. In this method we do not use derivation, but rather a "step-by-step" approximation to the extrema.

Mathematical formulation (for maximum)

The method of the golden section search has the function Fehler beim Parsen (http://mathoid.testme.wmflabs.org Serverantwort ist ungültiges JSON.): f:\mathbb{R}^{n} \supseteq X \rightarrow \mathbb{R}^{n}

based on an interval (a₀,b₀), where “a” is smaller than ”b”. A termination tolerance of ε is given, with the objective ε to compare the distance between “a” and “b” in each iteration. Through each iteration, the interval will be diminished by the factor λ, which is given with the value 0,618. λ will be cut from the bigger part of the interval. k = n represents the iteration steps.

Step 0: Initialization

Set a₀ := a, b₀:= b, k := 0

       Calculate λk := ak + (1 – δ) * (bk – ak)
       and µk := ak + δ * (bk – ak)
       as well as f(λk) and f(µk)

Step1: Compare the distance with the termination tolerance. If the distance is smaller than the termination tolerance, stop the iteration, because the maximum has been reached.

If b_k – a_k < ε, then STOP.

Step 2: Compare the distance with the termination tolerance. If the distance is greater than the termination tolerance, continue with the following two cases:

Case a)

      If f(λk) < f(µk), 
      set ak+1 := λk,
      bk+1 := bk, 
      λk+1 := µk, 
      µk+1 := ak+1 + δ * (bk+1 – ak+1), 
      f(λk+1) := f(µk), 
      and compute f(µk+1).

Case b)

      If f(λk) ≥ f(µk),
      set ak+1 := ak,
      bk+1 := µk,
      µk+1 := λk,
      λk+1 := ak+1 + (1 – δ) * (bk+1 – ak+1),
      f(µk+1) := f(λk),
      compute f(λk+1)
      set k = k+1, go to Step 1.

The result is:

figure 1: interval demonstration of the two possibilities