0

对于我的硕士论文，我正在尝试实现一个代码来评估时间和温度的测量数据，然后用最小二乘法计算变量。问题是 Python 似乎没有改变积分边界中的变量。该模型基于热响应测试的移动线源，如果有人熟悉的话。代码如下所示：

def model(time_s, coeffs):
  integral = []
    for t in time_s:

        def f(u):
            return np.exp(-((x0 ** 2 / ((coeffs[0] / cv) + alpha_l * vth))
                            + (y0 ** 2 / ((coeffs[0] / cv) + alpha_t * vth)))
                          * (vth ** 2 / (16 * ((coeffs[0] / cv) + alpha_l) * vth * u)) - u) * (1 / u)

        i, err = quad(f, 0, ((vth ** 2) * t / (4 * ((coeffs[0] / cv) + alpha_l * vth))))
        integral.append(i)

    integral = np.asarray(integral)
    return (heatflow / (4 * m.pi * cv * np.sqrt(((coeffs[0] / cv) + alpha_l * vth) *
                                                ((coeffs[0] / cv) + alpha_t * vth)))) \
           * np.exp((vth * x0) / (2 * ((coeffs[0] / cv) + alpha_t * vth))) * integral\
           + Tb + coeffs[1] * heatflow


def residuals(coeffs, y, time_s):
    return y - model(time_s, coeffs)

guess = [lam_guess, Rb_guess]

x, flag = leastsq(residuals, guess, args=(time_s, temperature_mean))

所以我需要为coeffs[0]and获得一个很好的价值coeffs[1]。问题是 Python 只更改了coeffs[1]. 是否有更好的方法来实现模型女巫改变每个时间步长的积分边界？我想，我对四边形和/或最小二乘的工作方式有误。

4

0 回答 0