2022-04-29 python自学笔记，用python写了一个最优化方法中的黄金分割法

【python学习笔记】写了一个最优化方法中的黄金分割法

一、黄金分割法的原理
二、代码展示

# 利用黄金分割法计算f(x) = x^2 - x + 2在x属于[-1,3]上的解
# x1 = x_down + 0.382 * (x_up - x_down)
# x2 = x_down + 0,618 * (x_up - x_down)

x_down = float(input('请输入f(x)的下限：'))
x_up = float(input('请输入f(x)的上限：'))
x_e = float(input('精度要求：'))

x1 = x_down + 0.382 * (x_up - x_down)
x2 = x_down + 0.618 * (x_up - x_down)

f1 = x1**2 - x1 + 2
f2 = x2**2 - x2 + 2

# 精度要求
while abs(f1 - f2) > x_e:
    if f1 > f2:
        x_down = x1
        x_up = x_up
    if f1 < f2:
        x_down = x_down
        x_up = x2
    if f1 == f2:
        x_down = x1
        x_up = x2
    x1 = x_down + 0.382 * (x_up - x_down)
    x2 = x_down + 0.618 * (x_up - x_down)
    f1 = x1 ** 2 - x1 + 2
    f2 = x2 ** 2 - x2 + 2
    print('重新定义的上限：', round(x_down, 6), '重新定义的上限：', round(x_up, 6), end='\t')
    print('重新定义的x1：', round(x1, 6), '重新定义的x2：', round(x2, 6), end='\t')
    print('f1的值为：', round(f1, 6), 'f2的值为：', round(f2, 6), end='\t')
    print()

x = (x1 + x2) / 2
print('最优解：', round(x, 6), 'f = x^2 - x + 2：', round(x1**2 - x1 + 2, 6))

三、结果显示