Взаимная корреляция изолированного фазой Sine wave

В этом примере показано, как использовать последовательность взаимной корреляции, чтобы оценить, что фаза отстает между двумя синусоидами. Теоретическая последовательность взаимной корреляции двух синусоид на той же частоте также колеблется на той частоте. Поскольку демонстрационная последовательность взаимной корреляции использует меньше и меньше выборок в больших задержках, демонстрационная последовательность взаимной корреляции также колеблется на той же частоте, но амплитудные затухания как задержка увеличиваются.

Создайте две синусоиды с частотами $$2\pi /10$$ рад/отсчет. Стартовая фаза одной синусоиды 0, в то время как стартовая фаза другой синусоиды $$-\pi$$ радианы. Добавление $$N(0,0.25^2)$$ белый шум к синусоиде с задержкой фазы $$\pi$$ радианы. Установите генератор случайных чисел на настройки по умолчанию для восстанавливаемых результатов.

rng default

t = 0:99;
x = cos(2*pi*1/10*t);
y = cos(2*pi*1/10*t-pi)+0.25*randn(size(t));

Получите демонстрационную последовательность взаимной корреляции в течение двух периодов синусоиды (10 выборок). Постройте последовательность взаимной корреляции и отметьте известную задержку между этими двумя синусоидами (5 выборок).

[xc,lags] = xcorr(y,x,20,'coeff');

stem(lags(21:end),xc(21:end),'filled')

hold on
plot([5 5],[-1 1])

ax = gca;
ax.XTick = 0:5:20;

Вы видите, что peaks последовательности взаимной корреляции в задержке 5 как ожидалось и колеблется с периодом 10 выборок.

Смотрите также

xcorr