最優化算法-黃金分割搜索

/*****************************************
 * FileName  : golden_section_search.go
 * Author    : fredric
 * Date      : 2017.09.01
 * Note      : 黃金分割算法
 * History   :
*****************************************/
package search 

import(
    "fmt"
)

const(
    GOLDEN_SECTION_VALUE = 0.382
)

/*
* 黃金分割法的用途在于求救一元函數的極值，前提要求該函數在區間內有唯一極值，即單峰
* 取作用域[a0,b0]
* 每次迭代取[a1, b1]，其中a1 - a0 = b0 - b1 = p*(b0 - a0)
* 通過判斷如果
* : f(b1) > f(a1) ---> 則極值在區間a0 ~ b1之間
* : f(b1) < f(a1) ---> 則極值在區間a1 ~ b0之間
* 迭代繼續，直至滿足循環條件
* 通過計算可得p = 0.382
*/
func DoGoldenSectionSearch(){

    //假設方程為法f(x) = (x -1)*(x-1)，區間為0 ～ 2
    a0 := 0.0
    b0 := 2.0

    delta := 0.02 
    count := 0

    a1 := 0.0
    b1 := 2.0

    for b1 - a1 > delta {

        count++

        a1 = a0 + GOLDEN_SECTION_VALUE * (b0 - a0)
        b1 = a0 + (1 - GOLDEN_SECTION_VALUE) * (b0 - a0)

        f_a1 := (a1 - 1) * (a1 - 1)
        f_b1 := (b1 - 1) * (b1 - 1)

        if f_b1 > f_a1 {

            b0 = b1

        }else{

            a0 = a1
        }
    }

    fmt.Printf("a1 = %f, b1 = %f count=%d", a1, b1, count)

}