Взаимная корреляция с многоканальным входом

Открыть сценарий в реальном времени

Создайте три экспоненциальные последовательности из 11 образцов, задаваемые $^{значениями 0.4n}$ , $^{0.7n}$ и $^{0.999n}$ , с помощью $n≥0$ . Использовать stem3 для построения графика последовательностей бок о бок.

N = 11;
n = (0:N-1)';

a = 0.4;
b = 0.7;
c = 0.999;

xabc = [a.^n b.^n c.^n];

stem3(n,1:3,xabc','filled')
ax = gca;
ax.YTick = 1:3;
view(37.5,30)

Figure contains an axes. The axes contains an object of type stem.

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

[cr,lgs] = xcorr(xabc,'coeff');

for row = 1:3
    for col = 1:3
        nm = 3*(row-1)+col;
        subplot(3,3,nm)
        stem(lgs,cr(:,nm),'.')
        title(sprintf('c_{%d%d}',row,col))
        ylim([0 1])
    end
end

$Figure contains 9 axes. Axes 1 with title c_{11} contains an object of type stem. Axes 2 with title c_{12} contains an object of type stem. Axes 3 with title c_{13} contains an object of type stem. Axes 4 with title c_{21} contains an object of type stem. Axes 5 with title c_{22} contains an object of type stem. Axes 6 with title c_{23} contains an object of type stem. Axes 7 with title c_{31} contains an object of type stem. Axes 8 with title c_{32} contains an object of type stem. Axes 9 with title c_{33} contains an object of type stem.$

Ограничьте расчет лагами от $- 5$ до $5$ .

[cr,lgs] = xcorr(xabc,5,'coeff');

for row = 1:3
    for col = 1:3
        nm = 3*(row-1)+col;
        subplot(3,3,nm)
        stem(lgs,cr(:,nm),'.')
        title(sprintf('c_{%d%d}',row,col))
        ylim([0 1])
    end
end

Вычислить несмещенные оценки автокорреляций и взаимных перекрестных корреляций. По умолчанию задержки выполняются между $- (N-1$ ) и N-1.

cu = xcorr(xabc,'unbiased');

for row = 1:3
    for col = 1:3
        nm = 3*(row-1)+col;
        subplot(3,3,nm)
        stem(-(N-1):(N-1),cu(:,nm),'.')
        title(sprintf('c_{%d%d}',row,col))
    end
end

См. также

Функции

xcorr

Документация по инструментам обработки сигналов

Поддержка

Памятка переводчика

1. Если смысл перевода понятен, то лучше оставьте как есть и не придирайтесь к словам, синонимам и тому подобному. О вкусах не спорим.

2. Не дополняйте перевод комментариями “от себя”. В исправлении не должно появляться дополнительных смыслов и комментариев, отсутствующих в оригинале. Такие правки не получится интегрировать в алгоритме автоматического перевода.

3. Сохраняйте структуру оригинального текста - например, не разбивайте одно предложение на два.

4. Не имеет смысла однотипное исправление перевода какого-то термина во всех предложениях. Исправляйте только в одном месте. Когда Вашу правку одобрят, это исправление будет алгоритмически распространено и на другие части документации.

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