为一维优化问题指定边界SciPy - python

我有一个需要优化的一维函数。我的初始值为20,界限为(0,50)

x0=[20]
bounds=(0,50)
sol1=minimize(f,x0,method="SLSQP",bounds=bounds)

但是,这会产生IndexError。

SLSQP Error: the length of bounds is not compatible with that of x0.

我在这里犯什么错误?

python参考方案

正如评论中指出的那样,您可以将边界放入列表中。一个最小的示例如下所示:

from scipy.optimize import minimize


def f(x):
    return (x - 3) ** 2

x0 = [10]
bounds = [(0, 50)]

res = minimize(f, x0, method='SLSQP', bounds=bounds)

然后res.x将给出预期的array([3.])

正如@sascha在评论中指出的那样,对于此类问题,您也可以使用minimize_scalar

只需将上面的导入更改为

from scipy.optimize import minimize_scalar

和使用

res2 = minimize_scalar(f, method='bounded', bounds=bounds)

然后res2.x返回3.0

只需将上面的f替换为您的实际功能即可。

Python sqlite3数据库已锁定 - python

我在Windows上使用Python 3和sqlite3。我正在开发一个使用数据库存储联系人的小型应用程序。我注意到,如果应用程序被强制关闭(通过错误或通过任务管理器结束),则会收到sqlite3错误(sqlite3.OperationalError:数据库已锁定)。我想这是因为在应用程序关闭之前,我没有正确关闭数据库连接。我已经试过了: connectio…

scipy interp1d中的错误 - python

我不了解interp1d报告的结果。我收到NAN,我应该收到号码。In [131]: bb Out[131]: array([ 0. , 1.80286595, 1.87443683, 2.70410611, 3.02764722, 3.11305985, 3.11534355, 3.18695351, 3.20693444]) In [132]: alph…

Python pytz时区函数返回的时区为9分钟 - python

由于某些原因,我无法从以下代码中找出原因:>>> from pytz import timezone >>> timezone('America/Chicago') 我得到:<DstTzInfo 'America/Chicago' LMT-1 day, 18:09:00 STD…

用大写字母拆分字符串,但忽略AAA Python Regex - python

我的正则表达式:vendor = "MyNameIsJoe. I'mWorkerInAAAinc." ven = re.split(r'(?<=[a-z])[A-Z]|[A-Z](?=[a-z])', vendor) 以大写字母分割字符串,例如:'我的名字是乔。 I'mWorkerInAAAinc”变成…

如何打印浮点数的全精度[Python] - python

我编写了以下函数,其中传递了x,y的值:def check(x, y): print(type(x)) print(type(y)) print(x) print(y) if x == y: print "Yes" 现在当我打电话check(1.00000000000000001, 1.0000000000000002)它正在打印:<…