Question

我想使用scipy.optimize优化W（（1×9）matrix）中的9个变量。

from scipy.optimize import minimize

def func(W):
    W = W.reshape(1,9) #(1,9)
    Y = df0.values.reshape(49,1) #(49,1)
    X = df1.values.reshape(49,1) #(49,9)
    Z = np.dot(X, W.T) #(49, 1)
    Z = np.abs(Z - Y) #(49, 1)
    Cost = np.sum(Z ,axis=0, keepdims=True)
    return Cost[0][0]  #float

W = np.array([2,3,4,5,6,7,8,9,10])
cons = ({'type': 'ineq', 'fun': W[1]-W[0]})

result = minimize(func, x0=W0, constraints=cons, method="SLSQP")

但是我收到这样的TypeError

'numpy.int32' object is not callable

我将“缺点”更改为

def cons(W):
    return W

cons = (
    {'type': 'ineq', 'fun': cons}
)

然后一切正常，我得到了结果

     fun: 125.4977648197736
     jac: array([26.3666687 , 39.73333454, 46.56666756, 32.76666737, 38.23333454,
       25.20000076,  9.        ,  5.        ,  5.76666737])
 message: 'Optimization terminated successfully.'
    nfev: 332
     nit: 24
    njev: 24
  status: 0
 success: True
       x: array([7.36486798e-03, 8.29918593e-03, 9.61602518e-02, 9.17950729e-03,
       2.98795999e-12, 3.73831662e-12, 3.59100171e-12, 4.73656828e-01,
       1.77345002e+00])

我无法寻求一个好的解决方案。

Answer 1

由于事实，scipy的最小化需要实际功能（而不是其输出）才能正确最小化。参见this question。您的上述代码有效，因为您将函数引用传递给了可调用的函数（cons，而不是cons（W））。您可以尝试为其创建lambda函数，例如：

cons = ({'type': 'ineq', 'fun': lambda *args: W[1]-W[0]})

最小化Cost函数时带有scipy.optimize的TypeError

1 个答案: