wave_transform.py 文件源码

def MAD(x,n=3,fil=1):
        if fil == 1:
            meda = med.median_filter(x,size = (n,n))
        else:
            meda = np.median(x)
        medfil = np.abs(x-meda)
        sh = np.shape(x)
        sigma = 1.48*np.median((medfil))
        return sigma