python目标函数_使用limfit.minimize在Python中最小化目标函数

我对lmfit.minimize最小化程序包有问题.实际上,我无法为我的问题创建正确的目标函数.

问题定义

>我的函数：yn = a_11 * x1 ** 2 a_12 * x2 ** 2 … a_m * xn ** 2,其中xn-未知数,a_m-

??????????????系数. n = 1..N,m = 1..M

>就我而言,对于x1,..,x5,N = 5；对于y1,y2,y3,M = 3.

我需要找到最佳值：x1,x2,…,x5,以便可以满足y

我的问题：

>错误：ValueError：操作数不能与形状(3,)(3,5)一起广播.

>我是否在Python中正确创建了问题的目标函数？

我的代码：

import numpy as np

from lmfit import Parameters, minimize

def func(x,a):

return np.dot(a, x**2)

def residual(pars, a, y):

vals = pars.valuesdict()

x = vals['x']

model = func(x,a)

return y - model

def main():

# simple one: a(M,N) = a(3,5)

a = np.array([ [ 0, 0, 1, 1, 1 ],

[ 1, 0, 1, 0, 1 ],

[ 0, 1, 0, 1, 0 ] ])

# true values of x

x_true = np.array([10, 13, 5, 8, 40])

# data without noise

y = func(x_true,a)

#************************************

# Apriori x0

x0 = np.array([2, 3, 1, 4, 20])

fit_params = Parameters()

fit_params.add('x', value=x0)

out = minimize(residual, fit_params, args=(a, y))

print out

if __name__ == '__main__':

main()