【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))
三、结果显示