Phase shift correction between 2 signals using cross-correlation

Question

Sebsatien 2015-2-18

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/179098-phase-shift-correction-between-2-signals-using-cross-correlation

评论： dpb 2022-5-17

采纳的回答： dpb

在 MATLAB Online 中打开

Hello

I'm in need of some help if you please.

I have 2 sinsoidal signals with a phase shift of 2Pi/3 between them. Basically, I need to eliminate the phase shift and obtain 2 signals with the same phase. I'm looking forward to use the cross-correlation, that represents the differnce between both signals (the phase shift, in our case), to do so.

My problem is to know how to use the cross-correlation to "correct" the shifted signal by eliminating the phase shift indiquated by the cross-correlation. Which operation to use to do so?

Here is a portion of the code for further understanding :

t = 1 : 1000;
phaseshift= 2*pi/3;
s1 = sin (2*pi*t/500);
s2 = sin (2*pi*t/500+ phaseshift);                        
figure(1); clf;
plot (t,s1, t,s2);
Axis_x = [-length(s1)+1 : 1 : length(s1)-1];               
cross_corr = xcorr(s1,s2,'coeff');                       
figure(2); clf;
plot( Axe_x, coss_corr,'r');

Thanks in advance

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

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

Answer 1

dpb 2015-2-18

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/179098-phase-shift-correction-between-2-signals-using-cross-correlation#answer_168445

编辑：dpb 2015-2-18

在 MATLAB Online 中打开

[c,lags]=xcorr(s1,s2);               % compute cross correlation; keep lags vector
[~,iLag]=max(c(find(lags==0):end));  % find the max in one-sided
s3=circshift(s2,[0 iLag]);           % correct for the shift

For your case I get

>> [~,iLag]=max(c(1000:end))
iLag =
 162
>> 2*pi*iLag/500
ans =
  2.0358
>> phaseshift
phaseshift =
  2.0944

This isn't perfect; illustrates even with pure sine

8 个评论
显示 6更早的评论隐藏 6更早的评论

dpb 2015-10-30

在 MATLAB Online 中打开

Again would have to see the actual dataset to be able to tell for certain.

To understand what you get from the xcorr function for that dataset, try the following--

plot(x1), hold on,plot(x2,'r')  % show the two signals...
[c,l]=xcorr(x1,x2);             % compute correl, save lags
figure,plot(l,c)                % plot correlations vs lag

Now align the two figures vertically one under the other so both are visible and note where the peak is on the correlation at lag==4. Now note that if you count from 1 and move the LH signal over it aligns with the other at position 5 which is FOUR steps (lags) over from the original position. Hence owing to the indices in the signal arrays counting from one while the lags count from zero, there's an "off by one" error built in.

Now, this is for a deterministic signal; if there is noise and/or multiple subsets within a pair of signals that do show linear correlation the estimation will show that and the actual peak isn't so definitive as in that example.

ZIYI LIU 2022-5-17

Hi dpb,

What if I don't know the functions of these two lines, but only two set of data like these two lines. Then I can have the lags, how do I know the phase shift? Thanks!

dpb 2022-5-17

You know the number of lags and the sampling time, don't you?

请先登录，再进行评论。

Answer 2

Sebsatien 2015-2-19

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/179098-phase-shift-correction-between-2-signals-using-cross-correlation#answer_168537

Thanks a lot dpb. I tried to keep the lags vector before and I haven't succeded doing so. Now Thanks to you, it's clear. I shall therefore test the code.