Ising模型Python中的自相关函数

Question

我正在尝试按照此 this 并且图表应类似于 此具有较小的标准偏差。这是我的代码：

def computeAutocorrelation(Q_list,taunum,knum):
    A_list_tau = []
    A_std_tau = []
    
    for tau in range(taunum):
      
        A_list_k = [(Q_list[k]*Q_list[k+tau]- np.mean(Q_list[:knum-tau])**2)/variance(Q_list[:knum-tau]) for k in range(knum-tau)]

        A_list_tau.append(np.mean(A_list_k))
        A_std_tau.append(np.std(A_list_k))
        
        if A_list_tau[-1] <=0.01:
            A_list_tau.extend([np.nan]*(taunum-len(A_list_tau)))
            A_std1_tau.extend([np.nan]*(taunum-len(A_std_tau)))
            break

    return A_list_tau, A_std1_tau

但是，当我绘制这个函数时，我的图表看起来像 this 具有大标准偏差。

我在计算自相关函数或标准差时是否犯了任何错误？