因此,我尝试编写一个scipy.optimize.minimize的最小工作示例,其中不止一个示例。
基本上,我的示例适用于一个变量的lambda函数,但是一旦添加另一个变量,它就会崩溃。
lamX = lambda x: (x-2)**2
q0X = np.ones(1)
solX = optimize.minimize(lamX, x0=q0X)
lamXY = lambda x,y: (x-2)**2 + y**2
q0XY = np.ones(2)
solXY = optimize.minimize(lamXY, x0=q0XY)
前三行执行没有错误并给出了正确的结果,但后三行给出了以下错误
solXY = optimize.minimize(lamXY, x0=q0XY)
File "/usr/lib/python3/dist-packages/scipy/optimize/_minimize.py", line 444, in minimize
return _minimize_bfgs(fun, x0, args, jac, callback, **options)
File "/usr/lib/python3/dist-packages/scipy/optimize/optimize.py", line 913, in _minimize_bfgs
gfk = myfprime(x0)
File "/usr/lib/python3/dist-packages/scipy/optimize/optimize.py", line 292, in function_wrapper
return function(*(wrapper_args + args))
File "/usr/lib/python3/dist-packages/scipy/optimize/optimize.py", line 688, in approx_fprime
return _approx_fprime_helper(xk, f, epsilon, args=args)
File "/usr/lib/python3/dist-packages/scipy/optimize/optimize.py", line 622, in _approx_fprime_helper
f0 = f(*((xk,) + args))
File "/usr/lib/python3/dist-packages/scipy/optimize/optimize.py", line 292, in function_wrapper
return function(*(wrapper_args + args))
TypeError: <lambda>() missing 1 required positional argument: 'y'
任何人都可以给我提示我做错了什么事
答案 0 :(得分:1)
您的lambda需要为x
使用类似数组的对象。我使用它来工作:
>>> lamXY = lambda x: (x[0]-2)**2 + x[1]**2
>>> q0XY = np.ones(2)
>>> solXY = optimize.minimize(lamXY, x0=q0XY)
>>> solXY
fun: 3.865407235741147e-16
hess_inv: array([[0.75, 0.25],
[0.25, 0.75]])
jac: array([-9.04871520e-09, -1.62848344e-08])
message: 'Optimization terminated successfully.'
nfev: 12
nit: 2
njev: 3
status: 0
success: True
x: array([ 1.99999999e+00, -1.55929978e-08])
虽然对于python中的lambda通常不是这样,但scipy似乎希望lambda具有单个输入变量。