想要用 sympy 做多变量最小化 [英] Want to do multi-variation minimize with sympy

问题描述

我想通过使用符号化字符的 scipy.optimize 使用最小化

I want to use minimization with scipy.optimize using Symbolized charactors

from scipy.optimize import minimize
from sympy.utilities.lambdify import lambdify
import sympy as sp


x1, x2, x3, x4 = sp.symbols('x1 x2 x3 x4')

FormulaMain = sp.symbols('-2*x1**2*x3+6*x1**2*x4+13*x1**2-3*x1*x2**2+x1*x2+3*x1*x3**2-3*x4+103')
HandleMain  = lambdify((x1,x2,x3,x4),FormulaMain,'numpy')
bnds = ((-1, 1), (-1, 1), (-1, 1), (-1, 1))

PrintParams  = minimize(HandleMain,[1,1,1,1],method='SLSQP',bounds=bnds)

print PrintParams

当我运行代码时，我得到

When I run the code, I get

<lambda>() takes exactly 4 arguments (1 given)

我想我用 [1,1,1,1] 输入了 4 个参数我需要更改代码吗?

I think I have input 4 argument with [1,1,1,1] Are there anything I have to change with the code?

推荐答案

首先:欢迎来到 SO！

First of all: Welcome to SO!

据我所知，lambdify() 不能处理向量.此外，在使用 Sympy 时，确定雅可比很容易.你可以试试:

As far as I know, lambdify() cannot deal with vectors. Furthermore, when using Sympy, determining the jacobian is easy. You could try:

import numpy as np
from scipy.optimize import minimize
from sympy.utilities.lambdify import lambdify
import sympy as sy

sy.init_printing()  # LaTeX like pretty printing for IPython


x1, x2, x3, x4 = sy.symbols('x1 x2 x3 x4')
xx = (x1, x2, x3, x4)
f = -2*x1**2*x3+6*x1**2*x4+13*x1**2-3*x1*x2**2+x1*x2+3*x1*x3**2-3*x4+103
f_n = lambdify(xx, f, modules='numpy')

# Build Jacobian:
jac_f = [f.diff(x) for x in xx]
jac_fn = [lambdify(xx, jf, modules='numpy') for jf in jac_f]


def f_v(zz):
    """ Helper for receiving vector parameters """
    return f_n(zz[0], zz[1], zz[2], zz[3])


def jac_v(zz):
    """ Jacobian Helper for receiving vector parameters """
    return np.array([jfn(zz[0], zz[1], zz[2], zz[3]) for jfn in jac_fn])


bnds = ((-1, 1), (-1, 1), (-1, 1), (-1, 1))
zz0 = np.array([1, 1, 1, 1])

rslts = minimize(f_v, zz0, method='SLSQP', jac=jac_v, bounds=bnds)
print(rslts)

这篇关于想要用 sympy 做多变量最小化的文章就介绍到这了，希望我们推荐的答案对大家有所帮助，也希望大家多多支持IT屋！