2

我有一个平滑函数 f(x) = sin(x / 5) * exp(x / 10) + 5 * exp(-x / 2) 任务是找到一个非平滑函数 h(x) 的最小值= int(f(x)) 在 1 到 30 的区间上。换句话说，f(x) 的每个值都转换为 int 类型，并且该函数只接受整数值。

我正在使用 2 种方法从 scipy.optimize 中找到最小值：最小化和差分进化。最小化给了我-5的结果，而difficult_evolution给了Index error: tuple index out of range 问题是为什么和什么是错的？

这是代码：

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
4

1 回答 1

2

这是因为 的bounds参数differential_evolution需要一个 (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

我希望这有帮助。

于 2017-06-03T19:08:41.910 回答