def xcorr(a,b,mode ='coef'):
    fa = np.fft.fft(a,n=len(a)*2-1)
    fb = np.fft.fft(b,n=len(a)*2-1)
    xx = fa*np.conj(fb)
    xcc = np.fft.fftshift(np.fft.ifft(xx))
    aa = fa*np.conj(fa)
    xaa = np.fft.fftshift(np.fft.ifft(aa))
    bb = fb*np.conj(fb)
    xbb = np.fft.fftshift(np.fft.ifft(bb))
    if mode == 'coef':
        xcc = xcc/np.sqrt(xaa[len(a)-1]*xbb[len(b)-1])
    return np.real(xcc)

以上為政宏老哥寫的xcorr函數(shù)，對(duì)應(yīng)numpy或者scipy的為

cc = np.correlate(lope1, lope, mode='full')/np.sqrt(np.correlate(lope,lope,mode='valid')[0]
                                                                *np.correlate(lope1,lope1,mode='valid')[0])

速度提升巨大，lgnb，感謝??