斐波那契數(shù)列搜索，參考Edwin《最優(yōu)化導(dǎo)論》第四版7.3章節(jié)，算法采用go語(yǔ)言實(shí)現(xiàn)。

/****************************************** FileName : fibonacci_search.go* Author : fredric* Date : 2017.09.01* Note : 斐波那契數(shù)列搜索算法* History : *****************************************/ package search import("fmt" )func _get_fibonacci(i int) int{if i == 1 {return 1}else if i == 2{return 2}else{return _get_fibonacci(i - 1) + _get_fibonacci(i - 2)} }func _test_func(x float64) float64 {return x*x*x*x - 14*x*x*x + 60*x*x - 70*x }func _test_func_01(x float64) float64 {return (x - 1) * (x - 1) }/* * 基于黃金分割的思路對(duì)分割的比例系數(shù)p進(jìn)行優(yōu)化 * p采用斐波那契數(shù)列，即 * p1 = 1 - FN/FN+1 * p2 = 1 - FN-1/FN * ... * PN = 1 - F1/F2 * 總的壓縮比：p1*p2*..pN = 1/FN+1 * 因此 F N+1 需要能夠滿足壓縮比 */ func DoFibonnaciSearch(){fmt.Println("DoFibonnaciSearch")//最小區(qū)間為0.2//此時(shí)需要斐波那契的壓縮比 1 + 2e/F N + 1 <= 最小區(qū)間長(zhǎng)度／初始長(zhǎng)度//取e是一個(gè)很小的整數(shù)，如0.05//則N等于第五次迭代可以滿足要求a0 := 0.0b0 := 2.0p := 0.0delta := 0.05 //最后一次增加一個(gè)小整數(shù)做偏移for i := 5; i >=1; i-- {if i != 1 {fmt.Printf("a0 = %f, b0 = %f \n", a0, b0)//獲取斐波那契數(shù)列p1 := _get_fibonacci(i)p0 := _get_fibonacci(i - 1)p = 1 - float64(p0)/float64(p1)a1 := a0 + p * (b0 - a0)b1 := a0 + (1 - p) * (b0 - a0)f_a1 := _test_func_01(a1)f_b1 := _test_func_01(b1)if f_b1 > f_a1 {b0 = b1}else{a0 = a1}fmt.Printf("a1 = %f, b1 = %f f_a1 = %f, f_b1 = %f p = %f\n", a1, b1, f_a1, f_b1, p)}else{a1 := a0 + (1/2 - delta) * (b0 - a0)b1 := a0 + (1/2 - delta) * (b0 - a0)f_a1 := (a1 - 1) * (a1 - 1)f_b1 := (b1 - 1) * (b1 - 1)if f_b1 > f_a1 {b0 = b1}else{a0 = a1}}}//for i := 5; i >=1; i-- {fmt.Printf("a0 = %f, b0 = %f", a0, b0) }

?

轉(zhuǎn)載于:https://www.cnblogs.com/Fredric-2013/p/7531603.html

總結(jié)

以上是生活随笔為你收集整理的最优化算法-斐波那契数列搜索的全部?jī)?nèi)容，希望文章能夠幫你解決所遇到的問(wèn)題。

如果覺(jué)得生活随笔網(wǎng)站內(nèi)容還不錯(cuò)，歡迎將生活随笔推薦給好友。