blr_v12.py 文件源码

def integrand(x, lam, sigma, H):
    C2 = lambda H: np.pi / (H*gamma(2*H)*np.sin(H*np.pi))
    A = (lam*sigma)**2 / C2(H)
    return A * (2*(1-np.cos(x))*np.abs(x)**(-2*H -1)) / (lam**2 - 2*lam*(1-np.cos(x)) + 2*(1-np.cos(x)))

##################################################################################