用二分法求函数极小值点

二分法求函数极小值点的代码如何实现（最好matlab，python也可以）；多谢啦~

def fun(t):
    return t ** 4 / 4 + t ** 2 * 5 - 20 * t + 3

a, b = 1, 2
e = 10e-6
x0 = fun(a)
x1 = fun(b)

while a <= b:
    m = (a + b) / 2
    x = fun(m)
    if abs(x - x1) <= e:
        result = m
        break
    
    if abs(x - x0) > abs(x - x1):
        x0 = x
        a = m
    else:
        x1 = x
        b = m

print(m, fun(m))
# result 1.59423828125 -14.561859965720927

对比图如下：

先计算1和2时表达式的值，然后计算1.5时表达式的值，看一下哪一段更小，则在这一段继续二分

• 你可以看下这个问题的回答https://ask.csdn.net/questions/675730
• 你也可以参考下这篇文章：最小错误率贝叶斯分类器（包含matlab程序）
• 除此之外, 这篇博客: 详细介绍遗传算法的原理以及最值问题代码实现（MATLAB/Python）中的 （1）单变量最值问题 部分也许能够解决你的问题, 你可以仔细阅读以下内容或跳转源博客中阅读:
  

• 求解函数y=sin(x)+pi*cos(x)-exp(x.^2+3)，x∈[-1,1] 的最大值