黄金分割法求极值点

我们先来了解什么是黄金分割算法：

    黄金分割法也称0.618算法，属于区间收缩法，首先找出包含极小点的初始搜索区间，然后黄金分割点通过对函数值的比较不断缩小搜索区间（当然要保证极小点在搜素区间），当定义域的长度缩小的一定长度时候，就可以用当前区间的端点值的平均近似代替极小值点。

注：适用范围是单谷函数（就是只有一个极大值（转化成求极小值问题）或者极小值点）

通俗点讲就是讲先给定搜索区间比如[a b]，然后另x₁ = a + 0.382(b - a),x₂ = a + 0.618(b - a)，然后把x₁和x₂代入到函数f(x)中比较f(x₁)和f(x₂)的大小，如果f(x₁)>f(x₂)，则让a = x₁，否则b = x₂，然后在新的搜索区间[a b]内，重新找到x₁和x₂重复以上过程，直到b - a<ξ(这个是给出的最小精度)，然后取a,b的平均值近似代替f(x)_min。

代码实现如下：

 1 import numpy as np