- html - 出于某种原因,IE8 对我的 Sass 文件中继承的 html5 CSS 不友好?
- JMeter 在响应断言中使用 span 标签的问题
- html - 在 :hover and :active? 上具有不同效果的 CSS 动画
- html - 相对于居中的 html 内容固定的 CSS 重复背景?
我使用python2.7,需要找到多元标量函数的最大值。
换句话说我有这个功能:
def myFun(a,b,c,d,e,f):
# complex calculation that takes about 30 seconds
return res # res is a float
该函数不是凸函数。
我指定每个参数 a、b、c、d、e 和 f 的最小和最大可能值。我需要找到什么参数组合大约会导致 myFun
的最大值。我将为它提供一个不错的起点。
我尝试进行强力网格搜索,但考虑到我的函数需要多长时间来计算,它是不可行的。
我已经研究了 scipy 包。我特别看到了scipy.optimize.fmin_slsqp
功能。这适合我的问题吗?或者也许scipy.optimize.fmin()
?还有其他适合此的功能/模块吗?
最佳答案
您可能想尝试 CVXPY (http://www.cvxpy.org/en/latest),奇怪的是,它是 CVXOPT(凸求解器)的非凸扩展。然后,您可以使用 CVXOPT 进行凸优化,或者使用 CVXPY 进行非凸优化,只要适合您的问题即可。
Python 中有很多非凸求解器,其中许多都列在这里: https://scicomp.stackexchange.com/questions/83/is-there-a-high-quality-nonlinear-programming-solver-for-python ...但看来您实际上是在询问连续求解器,它可能是局部的或全局的,并且可以处理昂贵的函数。
就个人而言,我建议mystic
( https://pypi.python.org/pypi/mystic )。确实,我是作者,但它已经获得了大约十年的充足资金,并且它可以解决其他软件包几乎无法解决的高度受限的非凸问题。它还可以处理具有非线性约束的基本凸优化。此外,mystic
专为大规模并行计算而构建,因此您可以在多个级别的优化中轻松利用并行计算。如果您有足够的资源,mystic
可以进行集成优化,您可以将其想象为能够进行网格搜索(您可以在其中选择点的初始分布),而不是使用固定点在网格上,mystic
使用并行启动的快速线性求解器。
以下是 mystic
附带的近 100 个示例之一:
'''
Maximize: f = 2*x[0]*x[1] + 2*x[0] - x[0]**2 - 2*x[1]**2
Subject to: x[0]**3 - x[1] == 0
x[1] >= 1
'''
提出了两种解决方案(一种用于线性求解器,一种用于全局求解器):
def objective(x):
return 2*x[0]*x[1] + 2*x[0] - x[0]**2 - 2*x[1]**2
equations = """
x0**3 - x1 == 0.0
"""
bounds = [(None, None),(1.0, None)]
# with penalty='penalty' applied, solution is:
xs = [1,1]; ys = -1.0
from mystic.symbolic import generate_conditions, generate_penalty
pf = generate_penalty(generate_conditions(equations), k=1e4)
from mystic.symbolic import generate_constraint, generate_solvers, solve
cf = generate_constraint(generate_solvers(solve(equations)))
# inverted objective, used in solving for the maximum
_objective = lambda x: -objective(x)
if __name__ == '__main__':
from mystic.solvers import diffev2, fmin_powell
from mystic.math import almostEqual
result = diffev2(_objective, x0=bounds, bounds=bounds, constraint=cf, penalty=pf, npop=40, ftol=1e-8, gtol=100, disp=False, full_output=True)
assert almostEqual(result[0], xs, rel=2e-2)
assert almostEqual(result[1], ys, rel=2e-2)
result = fmin_powell(_objective, x0=[-1.0,1.0], bounds=bounds, constraint=cf, penalty=pf, disp=False, full_output=True)
assert almostEqual(result[0], xs, rel=2e-2)
assert almostEqual(result[1], ys, rel=2e-2)
这是另一个:
"""
Fit linear and quadratic polynomial to noisy data:
y(x) ~ a + b * x --or-- y(x) ~ a + b * x + c * x**2
where:
0 >= x >= 4
y(x) = y0(x) + yn
y0(x) = 1.5 * exp(-0.2 * x) + 0.3
yn = 0.1 * Normal(0,1)
"""
解决方案:
from numpy import polyfit, poly1d, linspace, exp
from numpy.random import normal
from mystic.math import polyeval
from mystic import reduced
# Create clean data.
x = linspace(0, 4.0, 100)
y0 = 1.5 * exp(-0.2 * x) + 0.3
# Add a bit of noise.
noise = 0.1 * normal(size=100)
y = y0 + noise
@reduced(lambda x,y: abs(x)+abs(y))
def objective(coeffs, x, y):
return polyeval(coeffs, x) - y
bounds = [(None, None), (None, None), (None, None)]
args = (x, y)
# 'solution' is:
xs = polyfit(x, y, 2)
ys = objective(xs, x, y)
if __name__ == '__main__':
from mystic.solvers import diffev2, fmin_powell
from mystic.math import almostEqual
result = diffev2(objective, args=args, x0=bounds, bounds=bounds, npop=40, ftol=1e-8, gtol=100, disp=False, full_output=True)
assert almostEqual(result[0], xs, tol=1e-1)
assert almostEqual(result[1], ys, rel=1e-1)
result = fmin_powell(objective, args=args, x0=[0.0,0.0,0.0], bounds=bounds, disp=False, full_output=True)
assert almostEqual(result[0], xs, tol=1e-1)
assert almostEqual(result[1], ys, rel=1e-1)
对于并行计算,mystic
可以利用 pathos
和 pyina
(参见:https://github.com/uqfoundation),您只需传递层次结构配置要用于并行运行的映射函数。这很容易。对于库存问题,它可能不是最快的,但(在我看来)对于无法解决的问题(由于大小或复杂性),它是最佳选择。
关于python-2.7 - 非凸优化器,我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/35767277/
我之前让 dll 注入(inject)器变得简单,但我有 Windows 7,我用 C# 和 C++ 做了它,它工作得很好!但是现在当我在 Windows 8 中尝试相同的代码时,它似乎没有以正确的方
我正在尝试制作一个名为 core-splitter 的元素,该元素在 1.0 中已弃用,因为它在我们的项目中起着关键作用。 如果您不知道 core-splitter 的作用,我可以提供一个简短的描述。
我有几个不同的蜘蛛,想一次运行所有它们。基于 this和 this ,我可以在同一个进程中运行多个蜘蛛。但是,我不知道如何设计一个信号系统来在所有蜘蛛都完成后停止 react 器。 我试过了: cra
有没有办法在达到特定条件时停止扭曲 react 器。例如,如果一个变量被设置为某个值,那么 react 器应该停止吗? 最佳答案 理想情况下,您不会将变量设置为一个值并停止 react 器,而是调用
https://code.angularjs.org/1.0.0rc9/angular-1.0.0rc9.js 上面的链接定义了外部js文件,我不知道Angular-1.0.0rc9.js的注入(in
我正在尝试运行一个函数并将服务注入(inject)其中。我认为这可以使用 $injector 轻松完成.所以我尝试了以下(简化示例): angular.injector().invoke( [ "$q
在 google Guice 中,我可以使用函数 createInjector 创建基于多个模块的注入(inject)器。 因为我使用 GWT.create 在 GoogleGin 中实例化注入(in
我在 ASP.NET Core 1.1 解决方案中使用配置绑定(bind)。基本上,我在“ConfigureServices Startup”部分中有一些用于绑定(bind)的简单代码,如下所示: s
我在 Spring MVC 中设置 initBinder 时遇到一些问题。我有一个 ModelAttribute,它有一个有时会显示的字段。 public class Model { privat
我正在尝试通过jquery post发布knockoutjs View 模型 var $form = $('#barcodeTemplate form'); var data = ko.toJS(vm
如何为包含多态对象集合的复杂模型编写自定义模型绑定(bind)程序? 我有下一个模型结构: public class CustomAttributeValueViewModel { publi
您好,我正在尝试实现我在 this article 中找到的扩展方法对于简单的注入(inject)器,因为它不支持开箱即用的特定构造函数的注册。 根据这篇文章,我需要用一个假的委托(delegate)
你好,我想自动注册我的依赖项。 我现在拥有的是: public interface IRepository where T : class public interface IFolderReposi
我正在使用 Jasmine 测试一些 Angular.js 代码。为此,我需要一个 Angular 注入(inject)器: var injector = angular.injector(['ng'
我正在使用 Matlab 代码生成器。不可能包含代码风格指南。这就是为什么我正在寻找一个工具来“ reshape ”、重命名和重新格式化生成的代码,根据我的: 功能横幅约定 文件横幅约定 命名约定 等
这个问题在这里已经有了答案: Where and why do I have to put the "template" and "typename" keywords? (8 个答案) 关闭 8
我开发了一种工具,可以更改某些程序的外观。为此,我需要在某些进程中注入(inject)一个 dll。 现在我基本上使用这个 approach .问题通常是人们无法注入(inject) dll,因为他们
我想使用 swing、spring 和 hibernate 编写一个 java 应用程序。 我想使用数据绑定(bind)器用 bean 的值填充 gui,并且我还希望它反射(reflect) gui
我有这段代码,当两个蜘蛛完成后,程序仍在运行。 #!C:\Python27\python.exe from twisted.internet import reactor from scrapy.cr
要点是 Spring Batch (v2) 测试框架具有带有 @Autowired 注释的 JobLauncherTestUtils.setJob。我们的测试套件有多个 Job 类提供者。因为这个类不
我是一名优秀的程序员,十分优秀!