Calculate autocorrelation using FFT in Matlab

by

Just like you stated, take the fft and multiply pointwise by its complex conjugate, then use the inverse fft (or in the case of cross-correlation of two signals: 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)];

%# compare with MATLAB's XCORR output
all( (xcorr(x)-r) < 1e-10 )

In fact, if you look at the code of xcorr.m, that’s exactly what it’s doing (only it has to deal with all the cases of padding, normalizing, vector/matrix input, etc…)

More Related Contents:

  1. What is a spectrogram and how do I set its parameters?
  2. Getting FFT peaks from data
  3. Understanding Matlab FFT example
  4. resizing 3D matrix (image) in MATLAB
  5. Sweep / chirp signal ends at incorrect frequency
  6. Proper way to add noise to signal
  7. Matlab Template Matching Using FFT
  8. Analytical Fourier transform vs FFT of functions in Matlab
  9. Remove noise from wav file, MATLAB
  10. How to apply a low-pass or high-pass filter to an array in Matlab?
  11. Is colon notation (:) equal to array(vector) [closed]
  12. To remove blood vessel [closed]
  13. Saving time and memory using parfor?
  14. How to training/testing my own dataset in caffe?
  15. [caffe]: check fails: Check failed: hdf_blobs_[i]->shape(0) == num (200 vs. 6000)
  16. Return Unique Element with a Tolerance
  17. Recursive Anonymous Function Matlab
  18. support vector machines in matlab
  19. How to find all permutations (with repetition) in MATLAB?
  20. How do I retrieve the names of function parameters in matlab?
  21. How to plot a circle in Matlab?
  22. Real time plot in MATLAB
  23. Matlab imshow omit NaN
  24. General method to find submatrix in matlab matrix
  25. find time shift between two similar waveforms
  26. How to abort a running program in MATLAB?
  27. MATLAB debugging: smarter way to stop the code with an specific condition?
  28. Generating all combinations with repetition using MATLAB
  29. Rounding to n significant digits
  30. How do I double the size of a vector in MATLAB with interpolation?

Leave a Comment