Как можно реализовать функцию с помощью numpy.correlate()?

Есть реализация функции оценки выборочного коэффициента автокорреляции:

def comp_acf(N, l):
    res = []
    for i in range(l):
        j = len(N) - i
        op_1 = np.sum(N[:len(N)-i] * N[i:])
        op_2 = np.sum(N[:len(N)-i])
        op_3 = np.sum(N[i:])
        op_4 = np.sum( N[:len(N)-i]**2)
        op_5 = np.sum( N[i:]**2)
        g = ((j * op_1) - (op_2 * op_3))/((np.sqrt((j * op_4) - pow(op_2, 2))) * (np.sqrt((j * op_5) - pow(op_3, 2))))
        res.append(g)
    return res

Как можно её реализовать с помощью функции numpy.correlate()? Формула оценки представлена на следующем скриншоте: