Strange results with xcorr

Question

Guilherme T Silva 2016-10-13

1
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/307199-strange-results-with-xcorr

评论： Michael King 2019-3-1

Hi,

I'm running the following script to calculate the cross-correlation between two signals:

 s1 = sig.signals.values(:,1,5);
 s2 = ref.signals.values(:,1,5);
 C = xcorr(s1,s2);
 figure; subplot(2,1,1); plot([s1 s2]); subplot(2,1,2); plot(C);

Unfortunately, instead of the expected result I get maximum cross-correlation on the frame size (1024 samples), like shown in the figure:

Any clue on what's going on?

Thanks, Guilherme

1 个评论
显示 -1更早的评论隐藏 -1更早的评论

Michael King 2019-3-1

Hi Guilherme T Silva,

Did you ever find a suitable solition to your question? I have found that detrending with the detrend function did not help with the issue you presented above that I was also having.

Michael

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Erick Oberstar 2018-7-25

3
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/307199-strange-results-with-xcorr#answer_330336

XCORR wants zero mean inputs - remove the DC offset and I suspect your results will be more inline with what you expect. You may also try passing in the 'unbiased' flag

3 个评论
显示 1更早的评论隐藏 1更早的评论

ES_Thorny 2018-12-6

I have a similar problem.

I have to compute the cross-correlation between two transient signals with non-zero mean. As they represent transient phenomena, it doesn't make sense to me to subtract the mean value The objective is to compute the time delay between the maximum of the two signals. The signals are not exactly correlated-similar, but I guess this is always the case ...

If i try to compute the time delay using XCORR, I get close results to what I expect (i.e the time delay checked visually by checking where the maximum of the two signals are) only using the UNBIASED options.

Why is that? Does the unbiased routine subtract the mean values from both my signals?

Please let me know what you think.

Regards,

E

Michael King 2019-3-1

Hi Emanuele Spinosa,

Did you find a suitable solution to your question? I have the same issue that Guilherme T Silva posted on 14 Oct 2016. And I have tried to use the detrend function without success, which I think was suggested by Erick Oberstar on 25 Jul 2018 and Pistachio on 14 Nov 2018

请先登录，再进行评论。

Answer 2

Image Analyst 2016-10-13

1
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/307199-strange-results-with-xcorr#answer_238888

I'm not sure what you were expecting. It looks reasonable to me. What were you expecting? As the signals slide past, they first start to overlaps at -1/2 the frame width, then the overlap gets more and more until they completely overlap, at a shift of the frame width, then the start to separate and the sum of the products decreases until they finally completely separate at 1.5 times the frame width. Then, since MATLAB can't have negative indexes, the shift of -1/2 the frame width will show up at index 1, and the correlation at a shift of 1.5 times the frame width will happen at twice the frame width. Were you expecting something different? Like you wanted to see ripples on the triangle or didn't expect to even see a triangle at all? Please explain what is strange to you.

2 个评论
显示无隐藏无

Guilherme T Silva 2016-10-14

I was expecting to see a peak around 35, which is the number of samples one signal appears to be lagged from the other. My main goal is to detect this number.

Image Analyst 2016-10-14

Try subtracting off the mean before you correlate. The problem is that your signal is on huge offsets, so essentially you're practically convolving two very tall step (rect) functions. And as you know, the convolution (correlation) of two rect functions is a triangle. That's like the first thing they teach when they teach convolution to you.

请先登录，再进行评论。