3
$\begingroup$

I have to compute the cross-correlation between two transient signals, which could be decomposed into a trend (with non-zero mean value) + noise.

In particular I am interested in the max of the cross-correlation.

It seems that in MATLAB the cross-correlation is more accurate if it is calculated using the function xcov() rather then xcorr(), probably as a consequence of the fact that the function xcov() removes the mean value of the signal.

An example of two signals between which I would have to measure similarity can be this one:

enter image description here

Do you know if this makes sense and why?

$\endgroup$
4
  • 2
    $\begingroup$ covariance and correlation are simply two different things – if you want one, then you don't want the other and vice versa. What is it that you need this for? That defines which you want. $\endgroup$ Nov 19, 2020 at 23:11
  • $\begingroup$ I'd like to compute the max of the cross-correlation as a measure of similarity between the two signals $\endgroup$
    – EmThorns
    Nov 20, 2020 at 0:17
  • $\begingroup$ then, the correlation is indeed what you want, otherwise high-variance signals have a high chance of showing an exceedingly large metric $\endgroup$ Nov 20, 2020 at 9:19
  • $\begingroup$ I agree. But do I have to remove the mean value to optimize the MATLAB routine? If I don't do that, I have a case in which the max of the cross-correlation of two slightly different signals (but not identical) computed with xcorr is higher then the max value of the auto-correlation (so in this case the two signals are identical), still computed using xcorr,. This doesn't make sense to me. If I remove the mean value, this doesn't happen. Does this depend on the MATLAB algorithm? $\endgroup$
    – EmThorns
    Nov 20, 2020 at 13:14

1 Answer 1

1
$\begingroup$

When trying to measure similarity between signals we're basically building a metric.
When doing so we need to ask what we want to be sensitive about.

For instance, if you don't remove the DC Component (The Mean) and use something like an integral to measure similarity (Correlation / Convolution based) then you are sensitive the added DC component.

For instance, take 2 random white noises and add to each a DC component of 1e6.
Measure their similarity by correlation and compare it to 2 Sine Waves with amplitude of 1e-3 and a relative shift phase of 1e-6 [Rad].
In your eyes, which should be more similar? Which is indeed?

I'd say that in most cases in signal processing, similarity should be measured ignoring the DC component.

$\endgroup$
23
  • $\begingroup$ In my case, I don't have an actual DC components, because my signals are transients, so they represent the quick evolution of a variable. Do you think that the mean value should also be removed in this case? In this case, it is even hard to define a mean value, actually there is a trend + superimposed noise. In addition, I suspect that a mean value interferes with the MATLAB algorithm, which uses the FFT to compute the cross-correlation rapidly. Probably this algorithm is optimized for signals with a DC component (more or less constant in time) + zero-mean components/noise. What do you think? $\endgroup$
    – EmThorns
    Nov 21, 2020 at 11:46
  • $\begingroup$ I modified the question, to be more specific. $\endgroup$
    – EmThorns
    Nov 21, 2020 at 11:52
  • $\begingroup$ Is it like rise? Can you model it in a parametric model? $\endgroup$
    – Royi
    Nov 21, 2020 at 15:40
  • $\begingroup$ I think that the shape does not matter so much. As an example, you can imagine a ramp, a successive flat part with a pulse inside and then a descent ... with superimposed noise of course $\endgroup$
    – EmThorns
    Nov 21, 2020 at 15:54
  • 1
    $\begingroup$ DC Component in DSP is basically the Mean of a signal. Regarding to the white noises I wanted to show the case the Mean of 2 signals hides the fast they are basically so different. Hence it makes sense in most cases not to measure it. $\endgroup$
    – Royi
    Nov 30, 2020 at 5:43

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service and acknowledge that you have read and understand our privacy policy and code of conduct.

Not the answer you're looking for? Browse other questions tagged or ask your own question.