最优化第五讲练习题

宇来风满楼

已于 2024-06-21 13:30:28 修改

阅读量834

点赞数 13

分类专栏：优化文章标签：算法

于 2024-06-06 22:36:22 首次发布

本文链接：https://blog.csdn.net/qq_46221910/article/details/139504545

版权

优化专栏收录该内容

2 篇文章 0 订阅

订阅专栏

在这里插入图片描述

def f(x):
    return x*x-4*x+5
a0,b0=1,3

1、均匀搜索
令 $\delta=(b_0-a_0)/N,a_i=a_0+i\delta,i=1,2,3$

while b0-a0>0.1:
    a=np.linspace(a0,b0,5)
    for i in range(1,4):
        if f(a[i-1])>f(a[i]) and f(a[i+1])>f(a[i]):
            a0,b0=a[i-1],a[i+1]
            break
print((a0+b0)/2)

极小点为2
2、黄金区间法
令 $\tau=(\sqrt{5}-1)/2,\lambda=a_0+( 1-\tau)(b_0-a_0),\mu=a_0+\tau(b_0-a_0)$

tao=(pow(5,0.5)-1)/2
lamda=a0+(1-tao)*(b0-a0)
miu=a0+tao*(b0-a0)
while miu-lamda>0.1:
    if f(lamda)<f(miu):
        b0=miu
    else:
        a0=lamda
    lamda=a0+(1-tao)*(b0-a0)
    miu=a0+tao*(b0-a0)
print((lamda+miu)/2)

极小点约等于2.09
3、导数二分法
记区间中点 $\lambda=(a_0+b_0)/2$

def first_order(x):
    return 2*x-4

while True:
    lamda=(a0+b0)/2
    value=first_order(lamda)
    if value==0:
        print(lamda)
        break
    elif value<0:
        a0=lamda
    else:
        b0=lamda