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ù)，對應(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冻河，感謝??
