问 利用Matlab中的FFT计算自相关

EN

就像你说的，取快速傅立叶变换，然后逐点乘以它的复数共轭，然后使用逆快速傅立叶变换(或者在两个信号互相关的情况下： Corr(x,y) <=> FFT(x)FFT(y)* )

x = rand(100,1);
len = length(x);
%# autocorrelation
nfft = 2^nextpow2(2*len-1);
r = ifft( fft(x,nfft) .* conj(fft(x,nfft)) );
%# rearrange and keep values corresponding to lags: -(len-1):+(len-1)
r = [r(end-len+2:end) ; r(1:len)];