Question

我需要通过非线性最小二乘回归得到参数（kf，beta1，beta2，gamma）。错误消息是：＆＃34; ValueError：操作数无法与形状（4,7）（0，）和＃34一起广播;

我在收集的下一个数据中做了4次实验：

flujo_ms（x轴数据）：具有7个位置的数组
＆＃34; fri＆＃34;：每个＆＃34; fri＆＃34;是一个由7个职位组成的数组。
brfv：取四个常数值，每个实验一个常数。

我想一起使用四个实验来最小化参数估计的误差。

import numpy as np
from scipy.optimize import leastsq

flux = [flujo_ms, flujo_ms, flujo_ms, flujo_ms]
brfv = [[0.00694] * 7, [0.00972] * 7, [0.0139] * 7, [0.0208]*7]
fr = [fr1, fr2, fr3, fr4]


def foulingRate(parameters, flux, brfv, mlts=8.22):
    kf, beta1, beta2, gamma = parameters
    FR = kf * np.exp(flux * (beta1 * brfv + beta2 * mlts + gamma))
    return FR


def objective(pars, yData, xData, brfv):
    # it will minimize this function
    err = yData - foulingRate(pars, xData, brfv)
    return err

x0 = [5.6 * 10 ** -4, -2.48 * 10 ** 8, 5.1 * 10 ** 4,
      2.81 * 10 ** 6]  # initial values for the parameters

plsq = leastsq(objective, x0, args=(fr, flux, brfv))
print("Fitted parameters = {0}".format(plsq[0]))

Answer 1

我找到了解决方案。我使用模块lmfit：cars9.uchicago.edu/software/python/lmfit

多参数估计。广播问题

1 个答案: