我有一个平滑函数 f(x) = sin(x/5) * exp(x/10) + 5 * exp(-x/2) 任务是找到非平滑函数 h 的最小值(x) = int(f(x)) 在 1 到 30 的区间内。换句话说,f (x) 的每个值都转换为 int 类型,并且该函数仅接受整数值。
我使用两种方法从 scipy.optimize 中找到最小值:minimum 和 Differential_evolution。最小化给出了-5的结果,而difference_evolution给出了索引错误:元组索引超出范围 问题是为什么以及出了什么问题?
这是代码:
import math
import numpy as np
from scipy.optimize import minimize
from scipy.optimize import differential_evolution
from scipy.linalg import *
import matplotlib.pyplot as plt
def f(x):
return np.sin(x / 5.0) * np.exp(x / 10.0) + 5 * np.exp((-x / 2.0))
def h(x):
return f(x).astype(int)
x = np.arange(1, 30)
y = h(x)
plt.plot(x, y)
plt.show()
#res = minimize(h, 30, method='BFGS')
#print res
res = differential_evolution(h, [1, 30])
print res
最佳答案
这是因为 difference_evolution
的 bounds
参数需要一个 (min, max) 对序列。鉴于您只有一对 x 的最小值和最大值,您可以执行以下操作:
res = differential_evolution(h, [(1, 30)])
print res
这将导致:
fun: -11.0
message: 'Optimization terminated successfully.'
nfev: 92
nit: 5
success: True
x: array([ 25.76747524])
更多信息请查看官方文档:https://docs.scipy.org/doc/scipy-0.19.0/reference/generated/scipy.optimize.differential_evolution.html
我希望这会有所帮助。
关于python - Scipy.optimize Differential_evolution 索引错误 : tuple index out of range,我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/44347216/